Cleaned up the input file and deleted superseded input files in preparation to... (f83d6415) · Commits · Kumar, Atul / LinearFokkerPlanck Axisymmetric

InputFiles/xp_IonSlowingDown.in

+33 −15

Original line number	Diff line number	Diff line
		&in_nml
		in%fileDescriptor = '2021_01_29e',
		! Simulation name:
		! ===============
		in%fileDescriptor = '2021_01_29f',

		! Magnetic field input data:
		! =========================
		in%rootDir ="/home/nfc/myRepos/LinearFokkerPlanck_Axisymmetric",
		in%BFieldFileDir = "/BfieldData", ! Location of file
		in%BFieldFileDir = "/BfieldData",
		in%BFieldFile = "/Bfield_b.txt",
		in%nz = 501, ! Number of points in BFieldFile

		! Simulation conditions:
		! =====================
		in%Nparts = 50000, ! Total number of particles
		in%Nsteps = 5000, ! Total number of time steps
		in%dt = 0.5E-7,
		in%threads_request = 10, ! Time step in [s]
		in%particleBC = 2, ! 1: source, 2: periodic
		in%dt = 0.5E-7, ! Time step in [s]
		in%zmax = +5.0,
		in%zmin = -5.0,
		in%threads_request = 10,

		in%iSave = .true.,
		in%iPush = .true., ! Enable RK4 particle pusher
		@@ -19,23 +26,30 @@ in%iHeat = .false., !
		in%iPotential = .false., ! Enable electric potential
		in%iDrag = .false., ! Enable plasma drag

		in%jstart = 1, !
		in%jend = 5000, !
		in%jincr = 10, !
		! Time steps to record:
		! ==============
		in%jstart = 1, ! start frame
		in%jend = 5000, ! end frame
		in%jincr = 10, ! increment

		! Collision operator conditions:
		! ==============================
		in%Te0 = 10., ! eV
		in%Ti0 = 10., ! eV
		in%ne0 = 1.0E+19, ! /m**3
		in%zmax = +5.0,
		in%zmin = -5.0,

		in%Aion = 1.,
		in%ne0 = 1.0E+19, ! m**-3
		in%Aion = 1., ! Ion mass
		in%Zeff = 1., ! Zeff for ion/impurity mixture
		in%Zion = 1., ! Charge number for main ion species
		in%species_a = 2, ! =1 for electrons, =2 for ions
		in%elevel = 10, ! Energy at which slowing down power is calculated E = elevel*Te0
		in%CollOperType = 2, ! 1 for Boozer-Only, 2 for Boozer-Kim

		! Particle boundary conditions:
		! ============================
		in%particleBC = 3, ! 1: Isotropic plasma source, 2: NBI, 3: Periodic

		! Initial conditions:
		!===================
		in%zp_InitType = 2, ! 1: uniform load, 2: gaussian load
		in%zp_init = 0.0, ! Mean particle injection
		in%zp_init_std = 0.3, ! STD of particle injection
		@@ -43,6 +57,8 @@ in%Ep_init = 1000., !
		in%Tp_init = 10.; ! Thermal kinetic energy
		in%xip_init = 0.924, ! Mean pitch angle where xip = cos(theta) = vpar/v

		! RF heating operator conditions:
		! ===============================
		in%f_RF = 28.0E+9 ! RF frequency
		in%zRes1 = 2.4 ! Defines interval where RF is present
		in%zRes2 = 3.0 ! Defines interval where RF is present
		@@ -51,6 +67,8 @@ in%kpar = 5864. !
		in%Ew = 10000 ! EBW E-minus wave electric field magnitude
		in%n_harmonic = 4 ! Cyclotron harmonic

		! Electric potential conditions:
		! =============================
		in%s1 = 0,
		in%s2 = 1.3,
		in%s3 = 4.0,

InputFiles/xp_Xray.in

deleted100644 → 0

+0 −65

Original line number	Diff line number	Diff line
		&in_nml
		in%fileDescriptor = '2021_01_27j',

		in%rootDir ="/home/nfc/myRepos/LinearFokkerPlanck_Axisymmetric",
		in%BFieldFileDir = "/BfieldData", ! Location of file
		in%BFieldFile = "/Bfield_b.txt",
		in%nz = 501, ! Number of points in BFieldFile

		in%Nparts = 50000, ! Total number of particles
		in%Nsteps = 500, ! Total number of time steps
		in%dt = 0.5E-7,
		in%threads_request = 10, ! Time step in [s]
		in%particleBC = 2, ! 1: source, 2: periodic

		in%iSave = .true.,
		in%iPush = .true., ! Enable RK4 particle pusher
		in%iColl = .false., ! Collisions or no collisions
		in%iHeat = .false., ! Enable RF heating operator
		in%iPotential = .false., ! Enable electric potential
		in%iDrag = .false., ! Enable plasma drag

		in%jstart = 1, !
		in%jend = 500, !
		in%jincr = 1, !

		in%Te0 = 10., ! eV
		in%Ti0 = 10., ! eV
		in%ne0 = 1.0E+19, ! /m**3
		in%zmax = +5.0,
		in%zmin = -5.0,

		in%Aion = 1.,
		in%Zeff = 1., ! Zeff for ion/impurity mixture
		in%Zion = 1., ! Charge number for main ion species
		in%species_a = 2, ! =1 for electrons, =2 for ions
		in%elevel = 10, ! Energy at which slowing down power is calculated E = elevel*Te0
		in%CollOperType = 2, ! 1 for Boozer-Only, 2 for Boozer-Kim

		in%zp_InitType = 2, ! 1: uniform load, 2: gaussian load
		in%zp_init = 0.0, ! Mean particle injection
		in%zp_init_std = 0.3, ! STD of particle injection

		in%kep_InitType = 2, ! 1: Maxwellian, 2: Beam at E = kep_init
		in%kep_init = 1000., ! Initial energy of EEDF

		in%xip_InitType = 2, ! 1: Isotropic, 2: Beam at xip = xip_init
		in%xip_init = 0.707, ! Mean value of xip = upar/u
		in%xip_init_std = 0.98, ! STD in xip

		in%f_RF = 28.0E+9 ! RF frequency
		in%zRes1 = 2.4 ! Defines interval where RF is present
		in%zRes2 = 3.0 ! Defines interval where RF is present
		in%kper = 29920. ! Perpendicular wave number of EBW
		in%kpar = 5864. ! Parallel wave number of EBW
		in%Ew = 10000 ! EBW E-minus wave electric field magnitude
		in%n_harmonic = 4 ! Cyclotron harmonic

		in%s1 = 0,
		in%s2 = 1.3,
		in%s3 = 4.0,
		in%phi1 = 0.,
		in%phi2 = 0,
		in%phi3 = 0.,

		/

src/Main.f90

+2 −2

Original line number	Diff line number	Diff line
		@@ -275,7 +275,7 @@ TimeStepping: do j = 1,in%Nsteps

		if (in%particleBC .EQ. 1) then
		call ReinjectParticles(zp(i),kep(i),xip(i),in,ecnt2,pcnt2)
		else if (in%particleBC .EQ. 2) then
		else if (in%particleBC .EQ. 3) then
		zp(i) = in%zmin
		end if

		@@ -283,7 +283,7 @@ TimeStepping: do j = 1,in%Nsteps

		if (in%particleBC .EQ. 1) then
		call ReinjectParticles(zp(i),kep(i),xip(i),in,ecnt1,pcnt1)
		else if (in%particleBC .EQ. 2) then
		else if (in%particleBC .EQ. 3) then
		zp(i) = in%zmax
		end if