所有输入参数

Basic 模块

参数名	取值	描述
nlextr	true(false)	Turn on(off) extrapolation method.
nrun	1000	Maximum number of time steps.
leqtyp	11	Determines which equilibrium package is to be used      for equilibrium flux surfaces computed with the fixed boundary moments code VMEC from S. P. Hirshman.
		for xxx
reread	1000	Time [sec] to reread namelist.
lsmord	2	Order of binomial filter, smoothing can be bypassed by setting korder < 0.
theta	0.5	Degree of implicitness in Crank—Nicholson differencing schemes.
thetap	0.5	Degree of implicitness for iteration in predictor—corrector scheme.
nsedit	0	Short printout frequency.
nplot	0	Plot frequency.
errmax	1.0E-6	Maximum allowed relative extrapolated error in any zone when extrapolation is used.
delmax	1.0E-6	Maximum relative change allowed in any zone over one time step without reducing the time step.
ntrans	1	If NTRANS < 1 - all diffusion coefficients are set to zero.
nbound	1	Type of boundary condition.
rastar	1	Normalization for empirical model above.
nedit	0	Generate major printout every NEDIT time steps.
nskip	0	Skip NSKIP time steps between major printouts.
splot	0	Plot every splot sec.
sedit	0	Long printout every 5.0 sec.
ntty	0	Frequency of output to terminal.
tinit	0.0	Initial time.
tmax	1.0	Maximum time.
dtinit	0.01	Initial timestep.
dtmax	0.01	Maximum timestep.
nzones	55	Number of radial zones.
cpvelc	0.0	Energy per particle carried by convection.
cpvion	0.0	Energy per particle carried by convection.
ntype	1	Local alpha plasma heating model.
rcwals	1.0	Major radius to alpha scrape off.
rdwals	0.0	Minor radius to alpha scrape off ( 1.5 * rminor ).
einoiz	0.0	Energy loss for each hydrogen ioniztion (1804).

一维数组

参数名	index	取值	描述
nlomt2	13	true(false)	(Not) Skip sbrtn cmpres.
	21	true(false)	(Not) Skip sbrtn cmpres.
	26	true(false)	(Not) skip sbrtn pwso_update.
limprd	32	2	For computes CMEAN, C2MEAN, WEIRS, and DWEIRS in the ghost zone just outside the plasma.
		3	xxx
lnumer	31	1/0	Smooth on fast ion densities if lnumer(31)>0.

Impurity 模块

参数名	取值	描述
ntaomc	2	Non-equilibrium impurity radiation model.
extzef	0.0	z effective.
extf	0.0	ne/ni dilution factor.
grecyc	0.0	Recycling off wall.
ngprof	0.0	Recycling off wall.
ngpart	0	Recycling off wall.
recflx	0	Recycling factor for diffusive losses to the wall.
recscr	0	Recycling factor for scrape-off losses.
nwcool	1100	Print local and vol rad power, electron cooling rate.
nwline	1000	Print local and vol rad power, line radiation.
nwioni	0	Ionization losses.
nwreko	0	Recombination radiation.
nwchex	0	Charge exchange radiation.
nwbrem	0	Bremsstrahlung radiation.

一维数组

参数名	index	取值	描述
cfutz	200	1	If cfutz(200)<0, prescribed impurity influxes are considered ionized. Setting cfutz(200)>0, activates the neutral impurity influx model.
	201	0.04	Prescribed incoming kinetic energy for neutrals of first impurity species species.
	202	0.04	Neutral impurity influx can only be prescribed for the first two impurity species.
	203	0.18	Electron energy loss per initial ionization of a neutral atom of first impurity species.
	204	0.18	Same as cfutz(203), but for second impurity species.
	207	1	N+2-cfutz(207)=index of outermost radial zone allowing ionization of neutrals of first impurity neutrals of first impurity species.
	208	1	Same as cfutz(207), but for second impurity species.
	209	0.9	Prescried maximum penetration depth for neutrals of first impurity species.
	210	0.9	Same as cfutz(209), but for second impurity species.
	281	1.2	Switch for ${\chi }_i$.
	295	1	$e_{An,19}\to e_{An,19}+v_{e^{\ast }}/(cfutz(95)+v_{e^{\ast }}),A=D,i,e,v$ for all $n$.
	410	1	Multiplier on electron—ion equilibration rate.
	480	1	xxx
	481	1	xxx
xfutz	2	1.732	Maximum rate of increase of time step.
	6	0.5	≠ 0. → Override relative extrapolation error method of decreasing time step, and decrease time step when necessary by xfutz(6).
cimprd	10	0.1	xxx

Transport 模块

一维数组

参数名	index	取值	描述
cfutz	1	0	x pseudoclassical contribution to $D_{HH},H=1,2$.
	2	0	x pseudoclassical contribution to $D_{II},I=l,h$.
	3	0	x pseudoclassical contribution to ${\chi }_e$.
	4	0	x pseudoclassical contribution to ${\chi }_i$.
	5	0	x Bohm contribution to $D_{HH},H=1,2$.
	6	0	x Bohm contribution to $D_{II},I=l,h$.
	7	0	x Bohm contribution to ${\chi }_e$.
	8	0	x Bohm contribution to ${\chi }_i$.
	9	1	x simplified neoclassical contribution to $D_{HH}$, diagonal terms only; $H=1,2$.
	10	0	x simplified neoclassical contribution to $D_{II}$, diagonal terms only; $I=l,h$.
	11	1	x simplified neoclassical contribution to ${\chi }_e$.
	12	1	x simplified neoclassical contribution to ${\chi }_i$.
	19	1	x $v_{Ware}^{H,e}$ (factor multiplying Ware pinch of hydrogen isotopes and electron energy).
	21	0	$D_{HH}$ → + cfutz(21); $H=1,2$.
	23	0	${\chi }_e$ → + cfutz(23).
	24	0	${\chi }_e$ → + cfutz(24).
	29	0	${\chi }_e$ → + cfutz(29) $(r/r_{scr}{)}^{cfutz(30)}$.
	30	2	Modified by “soft Bohm limit”.
	35	0.0002	$D_{HH}$ → + cfutz(35) $(r/r_{scr}{)}^{cfutz(36)};H=1,2$.
	36	-6	$D_{HH}$ → + cfutz(35) $(r/r_{scr}{)}^{cfutz(36)};H=1,2$.
	37	0	$D_{HH}$ → + cfutz(37)/$n_e;H=1,2$.
	38	0	$D_{II}$ → + cfutz(38) $(r/r_{scr}{)}^{cfutz(39)};H=l,h$.
	39	3	$D_{II}$ → + cfutz(38) $(r/r_{scr}{)}^{cfutz(39)};H=l,h$.
	42	0	$D_{II}$ → + cfutz(42)/$n_e;I=l,h$.
	43	500000	$D_{aa}$, maximum allowed value; $a=1,2,l,h$.
	44	500000	${\chi }_j$, maximum allowed value; $j=i,e$.
	57	0	${\chi }_e$ → + 0, for cfutz(58)=0.
	58	0	${\chi }_e$ → + 0, for cfutz(58)=0.
	70	1	Note that even if cfutz(70) > 0, there will be no He volume source unless ${}^{4}He$ is specified as an impurity species.
	71	0.95	Fraction of helium ions crossing out at $r=r_{wall}$ which is recycled as a volume source.
	74	0.12	Electron energy lost per $H{e}^0$ ionization.
	76	0.04	(Unique) energy with $H{e}^0$ atoms are launched.
	81	4	Flow bias parameter: If cfutz(81)≠0, zone centers for some of the differencing formulas are offset by $(v\mathrm{\Delta }t)/(cfutz(81)\mathrm{\Delta }r)$, where $v$ is the convective velocity.
	110	1	cfutz(110)>0. Activities “nearly exact” neoclassical particle diffusivities between each ionic species and all others, and replaces simplified off-diagonal neoclassical diffusivities.
	111	1	Factor multiplying classical contribution to “nearly exact” neoclassical particle diffusivities.
	112	1	Factor multiplying banana-plateau contribution to “nearly exact” neoclassical particle diffusivities.
	113	1	Factor multiplying Pfirsch—Schluter contribution to “nearly exact” neoclassical particle diffusivities.

Transport model 模块

参数名	取值	描述
fdrint	0	Electron-ion energy interchange coefficient.
fdr values	0,0,0	xxx
fig values	0.8,0.8,0.8	xxx
fec values	0,0,0	xxx
fkb values	1,0.65,0.65	xxx
frb values	1,1,1	xxx
fmh values	0,0,0	xxx

一维数组

参数名	index	取值	描述
cthery	1	0	xxx
	2	1.013	xxx
	6	0	xxx
	7	1	xxx
	8	3.5	xxx
	12	-4	xxx
	13	-4	xxx
	14	-4	xxx
	15	-4	xxx
	17	0	xxx
	50	10	xxx
	51	10	xxx
	52	10	xxx
	53	10	xxx
	54	10	xxx
	78	1	xxx
	79	1	xxx
	81	-4	xxx
	82	0	xxx
	83	-4	xxx
	85	2	xxx
	86	0.15	xxx
	111	0	xxx
	112	0	xxx
	113	0	xxx
	114	0	xxx
	121	1	xxx
	123	1	xxx
	124	0	xxx
lthery	3	4	xxx
	4	1	xxx
	5	2	xxx
	6	1	xxx
	7	27	xxx
	8	21	xxx
	9	2	xxx
	13	3	xxx
	14	1	xxx
	15	0	xxx
	16	0	xxx
	17	1	xxx
	21	1	xxx
	25	0	xxx
	26	400	xxx
	27	1	xxx
	29	5	xxx
	30	1	xxx
	31	1	xxx
	32	-2	xxx
ltrnsp	20	82	xxx
	21	50	xxx

Initialization 模块

参数名	取值	描述
rmajor	814	Major radius in cm.
rminor	280	Minor radius in cm.
curent	21000	Current in kA.
bz	56800	Toroidal field in Gauss.
ti0	3	Initial ti(r) parabola^etifit.
te0	3	Initial te(r) parabola^etefit.
ti1	3	Initial ti(r) parabola^etifit.
te1	3	Initial te(r) parabola^etefit.
etifit	1	Initial ti(r) parabola^etifit.
etefit	1	Initial te(r) parabola^etefit.
ngas	-2,-3	Sets the hydrogen isotopes for species a. If wtgas(a)=0, ngas=1 → H, ngas=-2 → D, ngas=-3 → T.
denga0	2.0354e13,2.0354e13	Initial central density of hydrogen isotope a; ignore if dengas(a,1)>0.
denga1	2.1481e12,2.1481e12	Density at all times in outer "dummy zone" for hydrogen isotope; ignored if dengas(a,1)>0.
denim0	7.3680e11,4.6050e11	Initial central density of impurity species a; ignored if denimp(a,1)>0.
denim1	7.7760e10,4.8600e10	Density at all times in outermost radial zone for impurity species a; ignored if denimp(a,1)>0.
nimp	6,2	Impurity species a.
ehfit	1,1	xxx
eifit	1,1,1	Outer exponents.
eehfit	2,2	Initial density of hydrogen isotope a is $n_a^H=n_{1a}^H+(n_{0a}^H-n_{1a}^H)[1-(r/r_{wall}{)}^{x_a^H}{]}^{y_a^H}$.
eeifit	2,2,2	Initial density of impurity species a is $n_a^I=n_{1a}^I+(n_{0a}^I-n_{1a}^H)[1-(r/r_{wall}{)}^{x_a^I}{]}^{y_a^I}$.

Boundary condition 模块

参数名	取值	描述
tcoldp	1	Energy of gas puffed neutrals [keV].
tcold	1	Temperature of recycled neutral hydrogen.
bdtime	0.0,90.0,94.0,98.0,102.0,106.0,110.0,114.0,118.0,122.0,126.0,150.0,250.0,500.0,900.0	xxx
bdtee	3,3,4.206,4.7,5,5,5,5,5,5,5,5,5,5,5	xxx
bdtie	3,3,4.206,4.7,5,5,5,5,5,5,5,5,5,5,5	xxx
flgas	0,0	Specifies a gas puffing influx rate for neutral hydrogen isotope a, linearly interpolated between times $t_j^H$. This is independent of density monitoring.

一维数组

参数名	index	取值	描述
bdhyde(i,j)	1	2.1481E+12,2.1481E+12,1.0224E+13,1.5375E+13,1.8450E+13,2.1525E+13,2.4600E+13,2.7675E+13,3.0750E+13,3.3825E+13,3.3825E+13,3.3825E+13,3.3825E+13,3.3825E+13,3.3825E+13	xxx
	2	2.1481E+12,2.1481E+12,1.0224E+13,1.5375E+13,1.8450E+13,2.1525E+13,2.4600E+13,2.7675E+13,3.0750E+13,3.3825E+13,3.3825E+13,3.3825E+13,3.3825E+13,3.3825E+13,3.3825E+13	xxx
bdimpe(i,j)	1	7.7760E+10,7.7760E+10,3.5200E+11,5.0000E+11,6.0000E+11,7.0000E+11,8.0000E+11,9.0000E+11,1.0000E+12,1.1000E+12,1.1000E+12,1.1000E+12,1.1000E+12,1.1000E+12,1.1000E+12	xxx
	2	4.8600E+10,4.8600E+10,7.2000E+11,1.8750E+12,2.2500E+12,2.6250E+12,3.0000E+12,3.3750E+12,3.7500E+12,4.1250E+12,4.1250E+12,4.1250E+12,4.1250E+12,4.1250E+12,4.1250E+12	xxx

Gas puffing and pellet injection 模块

参数名	取值	描述
npuff	2	Gas puffing parameters, monitor line average density.
npel2	1	Milora-foster pellet model ornl-tm-5776.
gflmax	7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19,7e19	Maximum gas puffing rate which can be forced by density monitoring linearly interpolated between times $t_j^G$.
denmon	1,2.7e13,2.7e13,4.0e13,5.0e13,6.0e13,7.0e13,8.0e13,9.0e13,10.0e13,11.0e13,11.0e13,11e13,11e13,11e13,11e13	Target density for feedback monitoring of $n_e$ at $t_j^G,j=2,3,...,20$. denmon(1)=1.
gfract	0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5	Fraction of hydrogen isotope i, in gas puffing due to density monitoring, linearly interpolated between times $t_j^G$.
timp	0,90,94,98,102,106,110,114,118,122,126,130,140,150,160,170,180,200,500,999	Times to change specified impurity influx.
ftzeff	1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5,1.5	xxx

一维数组

参数名	index	取值	描述
flimp	1	1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11,1e11	Specified the flux of impurity ion or neutral species a across $r=r_{wall}$, linearly interpolated between times $t_j^I$.
rpa	1	0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0	Minimum major radius along the flight path of pellet launched at time $t_j^P$.
ypa	1	0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5,0.5	xxx
tpela	1	999.0,999.0,999.0,999.0,999.0,999.0	Times to launch pellets. For $j>1$, pellet is ignored if $t_j^P\le 0$.
vpela	1	0.5e5,0.5e5,0.7e5,1.0e5,1.3e5,1.5e5	Velocity of pellet launched at time $t_j^P$.
rpela	1	0.08,0.08,0.08,0.08,0.08,0.08	Radius of pellet launched at time $t_j^P$.
npelga	1	-2,-2,-2,-2,-2,-2,-2,-2,-2,-2	Hydrogen isotope in solid pellet launched at time $t_j^P$. If npelga does not follow the numbering conventions of ngas the pellet does not contribute to the hydrogen density.

NBI heating 模块

参数名	取值	描述
cjbeam	1	Beam driven current.
hnchek	0.1	Max frac chnge in ne befor recomp beam.
htchek	0.3	Max frac chnge in te befor recomp beam.
nipart	20000	Number of Monte Carlo particles used.
niprof	20	Max number of steps between beam comp.

一维数组

参数名	index	取值	描述
hton	1	90,90,90,90,90,90	Time at which beam b is turned on.
htoff	1	999,999,999,999,999,999	Time at which beam b is turned off.
hangle	3	28.5,28.5,28.5,28.5,28.5,28.5	Angles defining orientation of the line around which the distribution of beam neutrals has symmetry.
hrmaj	1	800,800,800,800,800,800	If hrmaj(b)≤0, we set $R_b^{MAJ}$=rmajor(1).
hrmin	1	500,500,500,500,500,500	If $R_b^{MIN}\le 0$, we set $R_b^{MIN}$=rminor(1).
haper	1	50,50,50,50,50,50	Width of aperture normal to the beamline at pivot point.
haperv	1	150,150,150,150,150,150	Height of aperture normal to the beamline at pivot point.
hdiv	1	0,0,0,0,0,0	Beam divergence.
hwidth	1	50,50,50,50,50,50	Width of injector face.
hfutz	1	0.25	“Approximation factors” for beam treatment.
height	1	150,150,150,150,150,150	Height of injector face.
hfocl	1	1e3,1e3,1e3,1e3,1e3,1e3	Horizontal focal length of beam.
hfoclv	1	1.4e3,1.4e3,1.4e3,1.4e3,1.4e3,1.4e3	Vertical focal length of beam.
hlenth	1	900,900,900,900,900,900	Distance along the beamline from injector face to pivot point.
nhaper	1	2,2,2,2,2,2	1 → circular port, 2 → rectangular port.
nhshap	1	2,2,2,2,2,2	1 → circular injector face, 2 → rectangular injector face.
hebeam	1	1000,1000,1000,1000,1000,1000	Energy of full energy component of injected neutrals.
hpowmw	1	16.5,16.5,0,0,0,0	xxx
nhbeam	1	-2,-2,-2,-2,-2,-2	Sets the injected species.

二维数组

参数名	index(i)	index(j)	取值	描述
hfract(i,j)	1	1	0.998,0.001,0.001,0.998,0.001,0.001,0.998,0.001,0.001,0.998,0.001,0.001,0.998,0.001,0.001,0.998,0.001,0.001	Current fractions of neutral atoms at full, half and one-third energy.

AEGIS 模块

Grid type parameters

参数名	取值	描述
eq_type	efit	type of equilibrium file: 'aegis', 'efit', 'gato' and 'toqdcon'.
psi_type	pest	type of psi: 'pest' and 'original'.
theta_type	hamada	type of theta: 'eqarc', 'hamada', 'pest' and 'general'.
a_fixed	1.0	not use.
s_sp	0.0	not use.
nest_sp	405	not use.

Physical effects

参数名	取值	描述
lg_amass	.true.	apparent mass effect.
lg_flr	.false.	FLR effect.
lg_denread	.false.	read density profile in a file named "profiles" .
zm	2	Parameter in FLR effect.
fac_ome_n	1	Parameter in FLR effect.
fac_ome_t	0.7	Parameter in FLR effect.

Grid setting

参数名	取值	描述
nthet2_p	801	number of theta grid in plasma.
npsi_sht_p	1001	number of psi grid in plasma.
chipct_efit	0.9998	cutting-off at edge.
n_eqstep	1	keeping seperation points.
n_surf_regmin	30	regions with resonances.
gd_power	3	compacking surface near the resonance.
psi_min_p	0.0002	cutting-off in the core.
nreg_add	10	additional number of shooting regions (first form).
psi_add(0-10)	0.0~1.0	additional region boundary.
nreg_add_mid	10	additional number of shooting regions (second form).
psi_add_mid(0-10)	0.0~1.0	additional region boundary.

Mode features

参数名	取值	描述
n_tor	1	toroidal mode number.
m_pol_plus	-4	number of "poloidal mode number": m_max=int(n_tor*qmax)+5-m_pol_plus.
m_pol_minus	-4	number of "poloidal mode number": m_min=int(n_tor*qmin)-4-m_pol_minus
comega	(-0.0,1e-3)	Re: uniformly rotation, Im: gamma
omg_rot	0.0	non-uniformly rotation profile used with rot_type.
bc_axis	.true.	initial conditions.

Shooting

参数名	取值	描述
tol_shoot	1.d-6	tol_shoot*factol: tolerance of integrating.
hp_fac	0.05	setting up high precision shooting regions.
fac_tol	1.d-2	tol_shoot*factol: tolerance of integrating.
n_tol		not use.
relabs	“R”	the type of error control.

Equilibrium profile $n$, $T$ and $P$

参数名	取值	描述
ntp_profile_type	'n' or 't' or 'p'	number density, temperature profiles
den_0	1.d20	number density at the magnetic axis
(d/t)coef(1-10)		ntp_profile='n': set n profile, and T=p/n; ntp_profile='t': set T profile, and n=p/T; ntp_profile='p': n=T=$\sqrt{p}$.

rotation profile

参数名	取值	描述
rot_type	'polynomial' or 'pressure'	rotation profile, 'pressure': use pressure profile
rotcoef(i)		used with rot_type='polynomial'

Diagnostics control parameters

参数名	取值	描述
diag_readeq	'true'	check grad shafranov eq.
diag_regrid	'true'	check q_toq and q.
diag_matrix	'false'	save for matrix
diag_shoot	'true'	plot independent solutions
diag_eigen	'false'	not use.

Spectrum

参数名	取值	描述
lg_spectrum	.false.	spectrum analysis.
lg_wrtplot	.true.	output grid point.
============================================

Vacuum wall type

参数名	取值	描述
nwall_shape	1 or 2	1 for conformal wall, 2 efit wall (reconstruct wall).
nwall_type	0 or 1	0: conducting wall, 1: resistive or liquid metal wall.
relabs_shoot_v	'R'	the type of error control.
tol_shoot_v	1.d-7	tolerance of integrating.
bwall	1.2	position of wall
ss_wx	0.01	wall spline for efit case (nwall_shape=2).
ss_wz	0.2	wall spline for efit case (nwall_shape=2).
nthet2_v	401	number of theta grid in vacuum.
npsi_v	40	number of psi grid in vacuum.

Thin liquid metal wall

参数名	取值	描述
vlm_0	0	liquid metal wall matrix.
thet_vel_in	0	liquid metal wall matrix.
thet_vel_out	1	liquid metal wall matrix.
omwall_over_omalf	1.d-4	complex determinant for dispersion.

Thick wall

参数名	取值	描述
nthe2_w	301	number of theta grid in wall.
wthick	0	number of psi grid in wall.
rlambda_w	(0,0)	included wall shoot matrix.
relabs_shoot_w	'R'	the type of error control.
tol_shoot_w	1.d-5	tolerance of integrating.

Diagnostic

参数名	取值	描述
wallreconst_lg	.true.	reconstruct wall.
diag_wallmat	.false.	cylinder check.
diag_wallshoot	.false.	cylinder check.
plot_pest	.true.	plot the real PEST Fourier component.
skip_green	.false.	preparision for green func.calculation.
v_shoot	.true.	shooting.
lg_wrtplot	.true.	output results.

Reread 模块

一维数组

参数名	index	取值	描述
reread	Reread1	90.2	Time [sec] to reread namelist.
	Reread2	245	Time [sec] to reread namelist.
	Reread3	999	Time [sec] to reread namelist.
dtmax	Reread2	0.5	Maximum allowed time step.
	Reread3	0.2	Maximum allowed time step.
	Reread4	0.5	Maximum allowed time step.

二维数组

参数名	index1	index2	取值	描述
cimprd	Reread1	10	1	xxx
	Reread2	10	1	xxx

Baldur control 模块

参数名	取值	描述
qmin	1.01	Minimum initial q value.
mombnd	2	Number of moments specified for equilibrium boundary.
mom	3	Number of moments computed in equilibrium interior.
ntbkmx	9	Number of breakpoint times.
shift	0.05	Initial shift Raxis - Rmid.
mflxs	25	Number of equilibrium flux surfaces.
swton	10.07	Sawtooth time on
swtoff	230	Sawtooth time on

一维数组

参数名	index	取值	描述
lsweq	1	3	xxx
	4	0	xxx
	5	1	xxx
	6	3	xxx
	7	1	xxx
	8	0	xxx
	10	0	xxx
	11	0	xxx
lsawth	1	1	xxx
	2	1	xxx
	3	10	xxx
	5	1	xxx
	6	0	xxx
	7	1	xxx
	10	2	xxx
csawth	1	0.9	xxx
	5	0	xxx
	9	1	xxx
tbk	1	10.0,20.0,30.0,40.0,50.0, 75.0, 100.0, 300.0, 999	xxx
neqdt	1	10,10,10,10,10,25,25,30,30	xxx
rmajbt	1	8.14,8.14,8.14,8.14,8.14,8.14,8.14,8.14,8.14	xxx
rminbt	1	2.8,2.8,2.8,2.8,2.8,2.8,2.8,2.8,2.8	xxx
elngbt	1	1.6,1.6,1.6,1.6,1.6,1.6,1.6,1.6,1.6	xxx
trngbt	1	0.24,0.24,0.24,0.24,0.24,0.24,0.24,0.24,0.24	xxx
btort	1	5.68,5.68,5.68,5.68,5.68,5.68,5.68,5.68,5.68	xxx
rtort	1	8.14,8.14,8.14,8.14,8.14,8.14,8.14,8.14,8.14	xxx
curmat	1	21,21,21,21,21,21,21,21,21	xxx

first part

first line of jobxdat:

F 0 0 0

分别是程序中的变量 nlres, nresum,nledge,nrec

• nlres 是否是继续上一次模拟，目前建议只取 .F.，nresum，nledge，nrec: 在该版本中没有用到。

jobxdat 中第 2-5 行，对该算例进行描述。

run1 namelist

基本参数，按照 commen block 划分

basic system parameters

1. 浮点数型变量

altime:

cptime:

stime:

versno: version number of the code.

reread: time to reread namelist

2. 整数型变量

nresum,nrec,nledge: not used.

ndiary,npunch,nin,nread,nonlin,nout,nprint: channel to read/write messages to teletype (do not change).

nstep: time step number.

*nrun: maximum number of time steps

3. 逻辑型变量

nlend： label for exit of the run.

nlres: label for resuming a previous run. .T.: resume a run. .F.: start a new run.

4. 字符型变量 label1,label2, label3, label4, label5, label6, label7, label8: page headers

control switches

1. 浮点数型变量 dimension(32)

*cnumer

cinout

cstart

ctrnsp

cemprc

cbeams

cdivrt

cauxil

*cimprd

cnugas

cneocl

cfsion

cstabl

2. 整数型变量 dimension(32)

*lnumer

linout

lstart

*ltrnsp

lemprc

lbeams

lstabl

ldivrt

lauxil

*limprd

lnugas

lneocl

lfsion

standard neutral gas variabl

1. 浮点数型变量

gblosi(2):

gfdni(2):

gflowi(2):

gfluxi(2):

gflx0i(2):

gfrac1(2):

gnchek: not used.

gnel1(2):

gneut(2):

*`grecyc`: fraction of ions or neutral leaving outermost zone which are recycled as neutral gas.

`gsputs(4)`:

`gtchek`: not used.

`gtflwi(20)`:

`gtprfi`:

`gvsrci(2)`:

`gwmin`: Minimum weight below which particles are no longer followed in Monte Carlo neutral algorithm.

`gxdome(2)`: Ratio to ion temperature of minimus and maximum energies used in least squares temperature calculation based on neutral outflux spectrum.

`gxmaxe`: maximum energy for neutral outflux spectrum in eV.

`gxmine`: minimum energy for neutral outflux spectrum in eV.

`gxphi(4)`: Defines azimuthal angle bins used in computation of neutral outflux spectrum.

`gxthet(4)`: Defines bins for angle between neutral outflux detector and normal to the wall, used in computation of neutral outflux spectrum.

`rhons(2,55)`:

`tns(2,55)`:

2. 整数型变量

*`ngpart`: number of particles in Monte Carlo computation of each source of neutral gas.

*`ngprof`: neutral profile every `NGPROF` time steps.

`ngsplt`: Number of splitting surfaces for neutral gas. (Section 2.9.1); must be less than or equal to 10.

`ngsprf`: 

`ngxene`: Number of energy bins for neutral outflux spectrum. must be less than or equal to 100. Not used.

`ngzone`: A minimum for the number of spatial zones in neutral gas computatoin.

3. 逻辑型变量

`nlglim`: `.TRUE.` $\rightarrow$ Include reflection off limiter; Ignore `TCOLD` and `NLGMON` for recycleing neutrals, and launch these neutrals with a Maxwellian energy distribution at the ion temperature of the outer dummy zone. Not that the gas puffing neutrals are not affected by `NLGLIM`.

`nlgmon`: `.TRUE.` $\rightarrow$ Monoenergetic, (`.FALSE.` $\rightarrow$ Maxwellian) recycled neutrals with energy(temperature) `TCOLD` and gas puff neutrals with energy (temperature) `TCOLDP`.    

`nlgpin`:     

`nlgref`:`.TRUE.` $\rightarrow$ Neutrals loss energy due to wall reflection; assumes a stainless steel wall.

`nlgspt`: `.TRUE.` $\rightarrow$ Sputtering of iron by neutrals crossing out over $r=r_{wall}$ is computed. [Sputtered iron is added to impurities as an influx analogous to `FLIMP` if `NIMP(1)`=26 or `NIMP(2)`=26]