py4sci

Table Of Contents

Previous topic

DCON README

Next topic

PYPEC Source documentation

This Page

Namelist Variables

COIL Namelist Inputs

COIL_CONTROL

ceq_type = “efit”
This should be identical to the eq_type specified in equil.in
machine = “d3d”
As of IPEC 3.00 supported machines include ‘nstx’, ‘d3d’, ‘kstar’, and ‘iter’.
ip_direction = “positive”
Set as ‘positive’ (default) or ‘negative’ for CCW or CW from a top down view respectively.
bt_direction = “negative”
Set as “positive” (default) or “negative” for CCW or CW from a top down view respectively.
coil_num = 3
Total number of coil sets to be activated.

cmpsi = 64

cmtheta = 480

cmzeta = 40

coil_name(1) = “iu”
Array values should be specified for each of the coil arrays to be activated. For example “coil_name(1) = c”. The supported coil sets (and number of coils in each) as of IPEC 3.00 include, - for nstx: “rwmef” (6), “tfef” (12), “pf5ef” (2), “ppu” (12), “ppl” (12), “psu” (12), “psl” (12), “vpu” (12), “vpl” (12), “vsu” (12), “vsl” (12), “hhfw” (24) - for d3d: “iu” (6), “il” (6), “c” (6), “tbm_solenoid” (1), “tbm_racetrack” (1) - for kstar: “fecu” (4), “fecm” (4), “fecl” (4) - for iter: “efcu” (6), “efcm” (6), “efcl” (6), “bl2u” (9), “bl2m” (9), “bl2l” (9), “avvu” (9), “avvm” (9), “avvl” (9)
coil_cur(1,1) = 2178.0
Array of Ampere coil current for the nc-th coil in the set corresponding to coil_name(nc) (default = 0).

coil_cur(1,2) = 1577.0

coil_cur(1,3) = -601.0

coil_cur(1,4) = -2178.0

coil_cur(1,5) = -1577.0

coil_cur(1,6) = 601.0

coil_name(2) = “il”

coil_cur(2,2) = -2178.0

coil_cur(2,3) = -1577.0

coil_cur(2,4) = 601.0

coil_cur(2,5) = 2178.0

coil_cur(2,6) = 1577.0

coil_cur(2,1) = -601.0

coil_name(3) = “c”

coil_cur(3,1) = -1551.0

coil_cur(3,2) = 268.0

coil_cur(3,3) = 1744.0

coil_cur(3,4) = 1551.0

coil_cur(3,5) = -268.0

coil_cur(3,6) = -1744.0

COIL_OUTPUT

ipec_interface = .TRUE.

IPEC Namelist Inputs

IPEC_INPUT

jac_in = “pest”
Coordinate system of 3D field input. Choose from hamada, pest, boozer, equal_arc or other
jsurf_in = 1
Set to 1 for surface weighted spectrum Common for surfmn inputs
tmag_in = 0
0 for ordinary toroidal angle 1 for magnetic coordinate angle
data_flag = .TRUE.
Set true to read a 3D field spectrum from the file specified by the infile variable and apply it to the plasma boundary.
data_type = “surfmn”
Type of vacuum 3D error field input file. Choose from ‘vac3d’ or ‘surfmn’.
infile = “/p/gpec/users/jpark/data/3d_field/d3d/surfmn.out.cossin.144182.03750_Intrins”
Path to vacuum 3D field input file of type specified by data_type. This file must specify the 3D vacuum field on the plasma boundary as determined in DCON.
nmin = 0
The minimum toroidal mode number included in the infile data (typically 0).
nmax = 64
The maximum toroidal mode number included in the infile data.
mmin = -64
The minimum poloidal mode number included in the infile data (typically -64).
mmax = 64
The maximum poloidal mode number included in the infile data (typically 64).
harmonic_flag = .FALSE.
If true Tesla amplitudes must be assigned to variables cosmn(int) and sinmn(int), where int is the applied 3D field m harmonic on the plasma surface. The toroidal mode number of the applied field is set by the DCON mode being analyzed.

cosmn(3) = 0.0001

displacement_flag = .FALSE.
Take all boundary inputs as normal displacements rather than magnetic fields
fixed_boundary_flag = .FALSE.
Sets the plasma boundary as defined in the equilibrium input file fixed in all IPEC calculations.

coil_flag = .TRUE.

IPEC_CONTROL

resp_index = 0

sing_spot = 0.0005

reg_flag = .TRUE.
Regularize unintegrable pmodb in psi
reg_spot = 0.05
dB ~ (m-nq)^2/[(m-nq)^2+reg_spot^2)] if reg_flag
chebyshev_flag = .TRUE.
Set true for chebyshev polynomials of the 3D perturbations in (psi, artheta).
nche = 20
Number of coefficients used in the chebyshev polynomials (typicaly ~30).

IPEC_OUTPUT

jac_out = “hamada”
Coordinates of IPEC output files. Choose from hamada, pest, boozer, equal_arc or other
jsurf_out = 0
Set to 1 for a surface weighted spectrum b_{mn}=ointleft(deltamathbf{B}cdothat{n}

ight)left( heta,phi ight)e^{-ileft(m heta-nphi ight)}Jcdot|nablapsi|d heta dphi (an invariant flux on rational surfaces). Set to 0 for unweighted spectrum b_{mn}=ointleft(deltamathbf{B}cdothat{n} ight)left( heta,phi ight)e^{-ileft(m heta-nphi ight)}d heta dphi.

tmag_out = 1
0 for ordinary toroidal angle. 1 for magnetic coordinate angle
resp_flag = .FALSE.
Output ipec_response_n#.out file.

singcoup_flag = .FALSE.

singfld_flag = .FALSE.
Output the q, Psi_{N}, resonant mathbf{B}, singular current, island width and Chirikov parameter at rational surfaces inside the plasma, including the plasma response.
vsingfld_flag = .FALSE.
Output the psi_{N}, resonant mathbf{B}, singular current, island width and Chirikov parameter at rational surfaces inside the plasma, using vacuum calculations. Requires coil_flag=True.
pmodb_flag = .FALSE.
Plasma |\delta B| on flux surfaces.
xbnormal_flag = .FALSE.
Output the normal components of the displacement, covarient perturbed magnetic field, and contravarient perturbed field (area weighted -> can compare between coordinates) with plasma response included.
vbnormal_flag = .FALSE.
Output the normal components of the displacement, covarient perturbed magnetic field, and contravarient perturbed field (area weighted -> can compare between coordinates) from vacuum calculations. Requires coil_flag=True.
fun_flag = .FALSE.
Outputs in functional format (on a coordinate grid) by back transforming the harmonic information (WARNING: this will result in larger data output).
flux_flag = .FALSE.
New flux surfaces psi(r+delta r,z+delta z,phi) for visualization of perturbation.
eqbrzphi_flag = .FALSE.
Output the original equilibrium field on a (r,z,phi) grid.
brzphi_flag = .FALSE.
Output the (r,z,phi) components of the perturbed magnetic fields everywhere on the (r,z,phi) grid. - IPEC 1 only outputs perturbed magnetic fields valid inside the plasma on the (r,z,phi) grid. - IPEC 2.0 and later includes the plasma response and vacuum fields when the coil_flag is the vacuum input.
xrzphi_flag = .FALSE.
Output the plasma displacement on a (r,z,phi) grid.
vbrzphi_flag = .FALSE.
Output the field due to the surface current defining the IPEC final solution boundary condition at the plasma surface (original external field boundary condition and plasma response) on an r,z grid. This does not represent a true field.
vvbrzphi_flag = .FALSE.
Outputs the field due to the surface current defining the external 3D field boundary condition on a (r,z,phi) grid. This field does not represent a true field.
bin_flag = .TRUE.
Output binary files for use with the xdraw command when appropriate.
bin_2d_flag = .TRUE.
Output binary files containing 2D contours for use with the xdraw command when appropriate.

vsbrzphi_flag = .FALSE.

ss_flag(7) = .FALSE.

ss_flag(8) = .FALSE.

xbrzphifun_flag = .FALSE.

arzphifun_flag = .FALSE.

ntv_flag = .TRUE.
Output first order quantities to bin files for interfacing with PENT code

IPEC_DIAGNOSE

div_flag = .FALSE.

radvar_flag = .TRUE.

PENT Namelist Inputs

PENT_INPUT

kinetic_file = “/p/gpec/users/nlogan/data/a10/kinetics/kinetic_a10_MARSbenchmark.dat”
1 header line followed by column data for psi_n, ni(m^-3), ne(m^-3), ti(eV), te(eV), omega_E(rad/s).
o1fun_file = “ipec_order1_fun_n1.bin”
IPEC output generated using ntv_flag or ascii B(psi,theta).
o1mn_file = “ipec_order1_n1.bin”
IPEC output generated using ntv_flag, or ascii B(m,n).

idconfile = “euler.bin”

o1native = .TRUE.
o1fun_file and o1mn_file are IPEC output bin files.
imass = 2
Mass of main species in atomic mass units (0 for electrons).
icharge = 1
Charge of main species in elementry charge units.
zimp = 6
Charge of impurity species in elementry charge units.
zmass = 12
Mass of impurity species in atomic mass units.
nl = 6
Number of bounce harmonics calculated (-nl<=l<=nl).
collision = “harmonic”
Collision operator. Choose from zero,small,krook, or harmonic.

PENT_CONTROL

wefac = 1
Multiplier applied to omega_E profile.
wdfac = 1
Multiplier applied to oemga_D.

wpfac = 1

ptfac = 0.001
Fractional cutoff in Lambda space near trapped-passing boundry.
kolim_flag = .FALSE.
Infinite rotation limit.
neorot_flag = .FALSE.
Approximate numerator of resonant operator using neoclassical rotation.
divxprp_flag = .TRUE.
Include arc-length change in perturbed Action.
xlsode_flag = .TRUE.
Use LSODE package to integrate in energy space.
lmdalsode_flag = .TRUE.
Use LSODE package to integrate in pitch angle.
lsode_rtol = 1e-09
Relative tolarance in LSODE.
lsode_atol = 1e-12
absolute tolarance in LSODE.
xmax = 72
Upper limit of energy integration.
ximag = 0
Step contour of integration off of real axis to avoid poles when collision is zero.
nx = 128
Number of energy spline points if xlsode_flag false.
nt = 128
Number of theta points used in bounce averaging.
nlmda = 128
Number of Lambda points used in splines.
maskpsi = 2
Step size taken by PENT on dcon psi grid (determined by idconfile).

PENT_OUTPUT

passing_flag = .FALSE.
Calculate torque and energy contributed by passing particles.
gen_flag = .TRUE.
Perform full calculation in general geometry.
rlar_flag = .TRUE.
Perform Reduced LAR approximate calculation.
clar_flag = .FALSE.
Perform Circular LAR approximate calculation.
kinetic_flag = .FALSE.
Output kinetic profiles.
bounce_flag = .FALSE.
Output sample bounce averaging profiles.
pitch_flag = .FALSE.
Output sample pitch angle profiles.
energy_flag = .FALSE.
Output sample energy profiles.
psiout(1) = 0.1
Include this surface in sampled outputs.
psiout(10) = 1.0
Include this surface in sampled outputs.

PENT_ADMIN

fmnl_flag = .FALSE.
Do not change unless you know what you are doing.
diag_flag = .FALSE.
Do not change unless you know what you are doing.