Replaced all Nparts with NC and Nsteps with NS (1de26b21) · Commits · Kumar, Atul / LinearFokkerPlanck Axisymmetric

InputFiles/xp_IonSlowingDown.in

+53 −53

Original line number	Diff line number	Diff line
		&in_nml
		&params_nml
		! Simulation name:
		! ===============
		in%fileDescriptor = '2021_02_17c',
		params%fileDescriptor = '2021_02_17c',

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

		! Simulation conditions:
		! =====================
		in%Nparts = 10000, ! Total number of particles
		in%Nsteps = 50000, ! Total number of time steps
		in%dt = 1.5E-7, ! Time step in [s]
		in%zmax = +2.0,
		in%zmin = -2.0,
		params%NC = 10000, ! Total number of particles
		params%NS = 50000, ! Total number of time steps
		params%dt = 1.5E-7, ! Time step in [s]
		params%zmax = +2.0,
		params%zmin = -2.0,

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

		! Time steps to record:
		! ==============
		in%jstart = 1, ! start frame
		in%jend = 50000, ! end frame
		in%jincr = 200, ! increment
		params%jstart = 1, ! start frame
		params%jend = 50000, ! end frame
		params%jincr = 200, ! increment

		! Collision operator conditions:
		! ==============================
		in%Te0 = 250., ! eV
		in%Ti0 = 250., ! eV
		in%ne0 = 3.0E+19, ! m**-3
		in%Aion = 1., ! Ion mass
		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
		params%Te0 = 250., ! eV
		params%Ti0 = 250., ! eV
		params%ne0 = 3.0E+19, ! m**-3
		params%Aion = 1., ! Ion mass
		params%Zion = 1., ! Charge number for main ion species
		params%species_a = 2, ! =1 for electrons, =2 for ions
		params%elevel = 10, ! Energy at which slowing down power is calculated E = elevel*Te0
		params%CollOperType = 2, ! 1 for Boozer-Only, 2 for Boozer-Kim

		! Particle boundary conditions:
		! ============================
		in%BC_Type = 2, ! 1: Isotropic plasma source, 2: NBI, 3: Periodic
		in%BC_zp_mean = 0.0, ! Mean particle injection
		in%BC_zp_std = 0.3, ! STD of particle injection
		in%BC_Ep = 20000., ! Drift kinetic energy
		in%BC_Tp = 10, ! Thermal kinetic energy
		in%BC_xip = 0.707, ! Mean pitch angle where xip = cos(theta) = vpar/v
		params%BC_Type = 2, ! 1: Isotropic plasma source, 2: NBI, 3: Periodic
		params%BC_zp_mean = 0.0, ! Mean particle injection
		params%BC_zp_std = 0.3, ! STD of particle injection
		params%BC_Ep = 20000., ! Drift kinetic energy
		params%BC_Tp = 10, ! Thermal kinetic energy
		params%BC_xip = 0.707, ! Mean pitch angle where xip = cos(theta) = vpar/v

		! Initial conditions:
		!===================
		in%IC_Type = 2, ! 1: uniform load, 2: gaussian load
		in%IC_zp_mean = 0.0, ! Mean particle injection
		in%IC_zp_std = 0.3, ! STD of particle injection
		in%IC_Ep = 20000., ! Drift kinetic energy
		in%IC_Tp = 10, ! Thermal kinetic energy
		in%IC_xip = 0.707, ! Mean pitch angle where xip = cos(theta) = vpar/v
		params%IC_Type = 2, ! 1: uniform load, 2: gaussian load
		params%IC_zp_mean = 0.0, ! Mean particle injection
		params%IC_zp_std = 0.3, ! STD of particle injection
		params%IC_Ep = 20000., ! Drift kinetic energy
		params%IC_Tp = 10, ! Thermal kinetic energy
		params%IC_xip = 0.707, ! 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
		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
		params%f_RF = 28.0E+9 ! RF frequency
		params%zRes1 = 2.4 ! Defines interval where RF is present
		params%zRes2 = 3.0 ! Defines interval where RF is present
		params%kper = 29920. ! Perpendicular wave number of EBW
		params%kpar = 5864. ! Parallel wave number of EBW
		params%Ew = 10000 ! EBW E-minus wave electric field magnitude
		params%n_harmonic = 4 ! Cyclotron harmonic

		! Electric potential conditions:
		! =============================
		in%s1 = 0,
		in%s2 = 1.3,
		in%s3 = 4.0,
		in%phi1 = 0.,
		in%phi2 = 0,
		in%phi3 = 0.,
		params%s1 = 0,
		params%s2 = 1.3,
		params%s3 = 4.0,
		params%phi1 = 0.,
		params%phi2 = 0,
		params%phi3 = 0.,

		/

InputFiles/xp_Test.in

+53 −53

Original line number	Diff line number	Diff line
		&in_nml
		&params_nml
		! Simulation name:
		! ===============
		in%fileDescriptor = '2021_02_01a',
		params%fileDescriptor = '2021_02_01a',

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

		! Simulation conditions:
		! =====================
		in%Nparts = 100, ! Total number of particles
		in%Nsteps = 600, ! Total number of time steps
		in%dt = 0.5E-7, ! Time step in [s]
		in%zmax = +5.0,
		in%zmin = -5.0,
		params%NC = 100, ! Total number of particles
		params%NS = 600, ! Total number of time steps
		params%dt = 0.5E-7, ! Time step in [s]
		params%zmax = +5.0,
		params%zmin = -5.0,

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

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

		! Collision operator conditions:
		! ==============================
		in%Te0 = 100., ! eV
		in%Ti0 = 100., ! eV
		in%ne0 = 5.0E+19, ! m**-3
		in%Aion = 1., ! Ion mass
		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
		params%Te0 = 100., ! eV
		params%Ti0 = 100., ! eV
		params%ne0 = 5.0E+19, ! m**-3
		params%Aion = 1., ! Ion mass
		params%Zion = 1., ! Charge number for main ion species
		params%species_a = 2, ! =1 for electrons, =2 for ions
		params%elevel = 10, ! Energy at which slowing down power is calculated E = elevel*Te0
		params%CollOperType = 2, ! 1 for Boozer-Only, 2 for Boozer-Kim

		! Particle boundary conditions:
		! ============================
		in%BC_Type = 2, ! 1: Isotropic plasma source, 2: NBI, 3: Periodic
		in%BC_zp_mean = 0.0, ! Mean particle injection
		in%BC_zp_std = 0.3, ! STD of particle injection
		in%BC_Ep = 20000., ! Drift kinetic energy
		in%BC_Tp = 10, ! Thermal kinetic energy
		in%BC_xip = 0.707, ! Mean pitch angle where xip = cos(theta) = vpar/v
		params%BC_Type = 2, ! 1: Isotropic plasma source, 2: NBI, 3: Periodic
		params%BC_zp_mean = 0.0, ! Mean particle injection
		params%BC_zp_std = 0.3, ! STD of particle injection
		params%BC_Ep = 20000., ! Drift kinetic energy
		params%BC_Tp = 10, ! Thermal kinetic energy
		params%BC_xip = 0.707, ! Mean pitch angle where xip = cos(theta) = vpar/v

		! Initial conditions:
		!===================
		in%IC_Type = 2, ! 1: uniform load, 2: gaussian load
		in%IC_zp_mean = 0.0, ! Mean particle injection
		in%IC_zp_std = 0.3, ! STD of particle injection
		in%IC_Ep = 1000., ! Drift kinetic energy
		in%IC_Tp = 10.; ! Thermal kinetic energy
		in%IC_xip = 0.707, ! Mean pitch angle where xip = cos(theta) = vpar/v
		params%IC_Type = 2, ! 1: uniform load, 2: gaussian load
		params%IC_zp_mean = 0.0, ! Mean particle injection
		params%IC_zp_std = 0.3, ! STD of particle injection
		params%IC_Ep = 1000., ! Drift kinetic energy
		params%IC_Tp = 10.; ! Thermal kinetic energy
		params%IC_xip = 0.707, ! 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
		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
		params%f_RF = 28.0E+9 ! RF frequency
		params%zRes1 = 2.4 ! Defines interval where RF is present
		params%zRes2 = 3.0 ! Defines interval where RF is present
		params%kper = 29920. ! Perpendicular wave number of EBW
		params%kpar = 5864. ! Parallel wave number of EBW
		params%Ew = 10000 ! EBW E-minus wave electric field magnitude
		params%n_harmonic = 4 ! Cyclotron harmonic

		! Electric potential conditions:
		! =============================
		in%s1 = 0,
		in%s2 = 1.3,
		in%s3 = 4.0,
		in%phi1 = 0.,
		in%phi2 = 0,
		in%phi3 = 0.,
		params%s1 = 0,
		params%s2 = 1.3,
		params%s3 = 4.0,
		params%phi1 = 0.,
		params%phi2 = 0,
		params%phi3 = 0.,

		/

InputFiles/xp_TestCollisions.in

+53 −53

File changed.

Preview size limit exceeded, changes collapsed.

InputFiles/xp_TestDrag.in

+2 −2

Original line number	Diff line number	Diff line
		@@ -12,8 +12,8 @@ params%nz = 501,

		! Simulation conditions:
		! =====================
		params%Nparts = 70000, ! Total number of particles
		params%Nsteps = 10000, ! Total number of time steps
		params%NC = 70000, ! Total number of particles
		params%NS = 10000, ! Total number of time steps
		params%dt = 0.5E-7, ! Time step in [s]
		params%zmax = +5.0,
		params%zmin = -5.0,

src/Modules.f90

+13 −14

Original line number	Diff line number	Diff line
		@@ -220,7 +220,7 @@ TYPE paramsTYP
		CHARACTER*150 :: BFieldFile, BFieldFileDir
		REAL(r8) :: Ti0, Te0, ne0, dt
		REAL(r8) :: Aion, Zion
		INTEGER(i4) :: Nparts, Nsteps, nz, species_a, species_b
		INTEGER(i4) :: NC, NS, nz, species_a, species_b
		INTEGER(i4) :: jstart, jend, jincr
		INTEGER(i4) :: threads_given
		LOGICAL:: iDrag, iPotential, iSave, iPush, iHeat, iColl
		@@ -238,7 +238,6 @@ TYPE paramsTYP
		END TYPE paramsTYP

		TYPE plasmaTYP
		REAL(r8) :: NC, Nsteps, dt
		REAL(r8) , DIMENSION(:), ALLOCATABLE :: zp, kep, xip, a
		INTEGER(i4), DIMENSION(:), ALLOCATABLE :: f1 , f2 , f3 , f4
		REAL(r8) , DIMENSION(:), ALLOCATABLE :: dE1, dE2, dE3, dE4, dE3_hat
		@@ -264,8 +263,8 @@ SUBROUTINE InitPlasma(plasma,params)
		! Declare local variables:
		INTEGER(i4) :: NC, NS

		NC = params%Nparts
		NS = params%Nsteps
		NC = params%NC
		NS = params%NS

		! Allocate memory:
		ALLOCATE(plasma%zp(NC) ,plasma%kep(NC),plasma%xip(NC),plasma%a(NC) )