from DYL

b2mn.dat

- 1. physics modules

*label          (lblmn: character*60)
  'b2mn.dat template for SOLPS-ITER' # 添加标签，会在包fstate中体现
*b2cmpa  basic parameters
*b2cmpb  boundary conditions
*b2cmpt  transport coefficients
*cflim    cflme    cflmi    cflmv     ...
  '-1'      0.3     10.0      0.5       0        0        0        0        0 # 对于2M cflmi选择 0.3，ITER选择为10
*endphy

- 2. RUN modules

*RUN
'b2mndr_ntim'     '5000' # 运行步数
'b2mndr_dtim'     '1e-5' # 步长

'b2mndt_nstg0'    '1' # 外层迭代次数
'b2mndt_nstg1'    '1'
'b2mndt_nstg2'    '10' # 最内层迭代次数

'b2stbc_boundary_namelist' '1' #指定从b2.boundary.parameter中读取边界条件，而不是b2ah.DNcase，强制从b2.boundary.parameter中读取边界
'b2stbr_neutrals_namelist' '1' # 指定从b2.neutrals.parameter 中读取流体中性信息，中性信息会被input.dat覆盖
'b2tqna_transport_namelist' '1' # 指定从b2.transport.parameter中读取输运系数，而不是b2ah
'b2tqna_inputfile' '1' # 指定从b2.transport.inputfile中读取输运系数

'b2mndr_na_min'   '1e4' # 每个网格所有粒子密度的下限，不需保证小于na_new，也要小于b2ai.dat的所有初始密度。（eps）
'b2mndr_na_new'   '1e14' # 每个网格，所有新粒子的初始密度（假如没有被初始价态文件覆盖）

'b2mndr_savecpu'  '3600.0' # 隔多长cpu时间，保存一个plasmastate.***,与b2co.exe结合，从crash中恢复。

'b2mndr_eirene'                   '1' # 非零与eirene耦合，必须用-DB25_EIRENE选项编译


'b2mndr_rescale_neutrals'         '1' # 第一个时间步，中性密度的乘子
'b2mndr_rescale_neutrals_sources' '1e-10' # 中性源乘子（eirene计算所得或者对应的比率系数）

• 1. nstg 指定内部迭代的次数

     
     image.png

---

- 2. OUTPUT module

*OUTPUT

*'tallies_netcdf'  '1' # 不等于零，产生‘b2tallies.nc’,包含区域标签信息；如果b2mndr_stim 小于零，在原有包tallies.nc上追加(>>)，否咋覆盖（>）
*'b2stbr_b2wall_netcdf'  '0' # 不等于零，产生b2wall.nc.包含wall标签信息。b2mndr_stim 小于零追加，否则强制重定向
*'balance_netcdf'  '0' # 非零产生balance.nc,包含所有平衡后处理的行信息
*'balance_average' '0' # 非零，对balance.nc的所有行取平均（总步数）

'b2mndr_b2time'   '1' # 指定写包time.nc 与时间依赖文件的时间步，大于零写最后的时间步（没明白），产生b2time.nc
'b2mndr_tally'    '10' #  指定tallies 的时间步长，大于零在最后的时间步写（明白了，就是隔多少部写一次信息）

-3. GEOMETRY module

*GEOMETRY

'b2mwti_jxi'  '61' # 内中平面的index，默认为（24+（72-24）/4=36）；DN（两个inner cut中间）；线性装置，在nx/4处
'b2mwti_jxa'  '38' # 外中平面index，SN（72-（72-24）/4=60）；DN(rightcut(1)+rightcut(2))/2;线性装置3nx/4;

-4. 数值，物理，原子 module

*NUMERICS


*PHYSICS

*ATOMIC

-5. BODONG module

*BUDONG

'b2mndr_min_areshe'     ' 1.00E-03'
'b2mndr_min_areshi'     ' 1.00E-03'
'b2mndr_min_aresco'     ' 1.00E-03' # 运行停止残差的最小值，所有标准同时满足，才结束一步运行。
'b2stbc_sna0ep'   '1.0e-10' # 指定网格中源项的小的等级
'b2trcl_lluciani' '1' # 指定lluciani模型的相关参数。1，采用标准连接长度计算公式。2，反向兼容SOLPS4.0。3，spb的新的lluciani系数
'b2news_potit'    '100' # 电势方程迭代的最大次数；potitmin指定迭代的最小次数。要保证最小值小于最大值。

'b2news_recalculate_contributions' '1'  # 非零时，当源项的错误流标签被激活时，从新计算

'eirene_repeat_first_call'      '1' # 大于零，在eirene_mc中重复第一次调用。在feedback喷气中很有用（到达eirene的流被忽略，采用eirene input file的值）
'eirene_extrap' '0' # 置1，最外层网格数值通过周围的真实网格计算；置零不改变。（疑问，靶板index的选择）

-5. DRIFT module

*DRIFT
*'b2news_facdrift_dec' '1.0'
*'b2news_facdrift_inc' '1.0'
*'b2news_facdrift_start' '1.0'
*'b2news_facdrift_target' '1.0'
# facdrift的梯度参数，乘到抗磁漂移项和惯性电流项。代码从facdrift=facdrift_start的时间步开始运行。如果facdrift_target不等于facdrift_start，每一个时间步长，乘以facdrift_inc;如果该时间步不收敛，facdrift通过facdrift_dec下降；中性——粒子摩擦流也需要抗磁漂移和E$\times$B其中有一个打开。
'b2news_facExB_dec'  '1.0'
'b2news_facExB_inc'  '1.0'
'b2news_facExB_start' '1.0'
'b2news_facExB_target' '1.0'
# 乘到E$\times$B dirft 的梯度参数，过程与抗磁漂移相同。 


*'b2news_facvis_dec'  '1.0'
*'b2news_facvis_inc'  '1.0'
*'b2news_facvis_start' '1.0'
*'b2news_facvis_target' '1.0'
# 漂移粘滞项项乘子（viscosity），指定facvis的梯度信息，处理方式同上。
*CURRENT
*'b2tfhe_vis_per'  '1.0' # 垂直粘滞电流项乘子。非零，电势方程通过b2nnp7子例程（七点模板）
*'b2tfhe_vis_par'  '1.0' # 平行粘滞驱动电流乘子。置零，粘滞驱动电流开关关闭，否则打开。
*'b2tfhe_neutral'  '1.0' # 粒子——中性电流项乘子。置零，粒子——中性电流项开关关闭。否则开启；粒子——中性电流项也需要开启抗磁漂移或者E*B漂移
*'b2tfhe_vis_q'    '1.0' # 热粘滞效应驱动电流乘子

-6. VISCOSITY TERM

*VISCOSITY TERM
'b2siav_addvis'   '1.0' # 动量方程中，粘滞张量对热流散度贡献的散度乘子。
'b2siav_addvis1'  '1.0' # 非零，加上粘滞张量磁场x-变化的散度贡献
'b2mndr_hz'       '1.0' # hz已经被引进到新的平行动量方程；置零用旧的方程；非零用新方程。
'b2sigp_style'    '2' # 置一，使用SPB格式的压强梯度项（动量平衡方程右侧）

-7. FLUX LIMIX

*FLUX LIMIT
'b2trno_flux_limit_to_dpa'  '1' # 置一，流限制对中性对dpa0的贡献被采用（b2tlc0.F中采用的输运系数），推荐置一；b2tlc0.F包含流限制参数，通过b2tlc0z_alpha/gamma 和 b2fnb_alpha/gamma指定，不能同时和零不同(不能同时不等于零）
'b2tlc0_alpha'              '1.0' 
'b2tlc0_gamma'              '1.0'
# 指定关于dpa0的流限制参数（压强驱动中性系数）；alpha是经典流限制的乘子；gamma(?)是流限制公式使用的指数项。若alpha置零，不采用流限制。
'b2tfhe_lim_flux'           '0' # 置零，流限制不直接作用于电子热流，而是通过输运系数传递。‘b2trcl_donductive_limit’置一，则推荐该开关置零。
'b2tfhi_lim_flux'           '0' # 离子热流。同上。
'b2trcl_conductive_limit'   '1' # 置一，平行电子和离子热流的流限制被输运系数采用。推荐置一，若置零，以上两个开关必须置零。

-8. *NEW manner OF CALCULATING DRIFT

*NEW manner OF CALCULATING DRIFT
'b2tfnb_drift_style'   '1' #  置零，漂移速度在网格中心计算。置一，在网格界面计算。推荐置一

-9. NEW FORM of SOME TERMS

*NEW FORM of SOME TERMS
'b2tfnb_mdf'    '1' # 置一，采用新的SPB格式，计算粒子流。推荐在开漂移时置一。
'b2tfhe_mdf'    '1' #新的SPB格式计算电子热流，漂移置一
'b2tfhi_mdf'    '1' # 离子热流，同上

-9. ABOUT SOLVER

*ABOUT SOLVER
'b2ux9p_style'          '2' # 选择九点模板矩阵计算器。置二，其他0,1返回错误。
'b2ux9p_nltrsol'        '0' # 矩阵计算器的输出标签。更大的数值意味着增加输出等级。

'b2ux7p_style'          '2' # 七点矩阵计算器，同九点要求。
'b2ux7p_nltrsol'        '0' # 同上上
'b2tfhe_vis_per'        '1' # 垂直粘滞电流的乘子；非零，电势方程通过b2npp7使用7点模板计算；推荐置一，在开漂移的时候，此时出垂直粘滞和对应的粘滞电流考虑在内。

'b2ux5p_style'          '2' # 5点计算器
'b2ux5p_nltrsol'        '0' # 同上

-10. SOME RECOMMENDED KEYS ABOUT DRIFT

*SOME RECOMMENDED KEYS ABOUT DRIFT
'b2tqce_fke_Zhdanov'    '1' # 置一，电子热传导采用Zhdanov 表达式，推荐置一。只有b2qce_model 和b2fhe_pTe 同时置一，该开关才激活。
'b2tfhe_fch_pTe'        '1' # 推荐置一，采用新的电子热流计算。
'b2tfhe_alfTeEh'        '1' # 推荐置一，采用旧的电子热流计算。
'b2sicf_phm0'           '1'  # 推荐置一，离心力项乘子
'b2trcl_lluciani'       '3' # 非零，网格上的输运系数属于闭合磁场线通过lucini模型。置一，连接长度。置二，SOLPS4.0。置三，spb's新的Lciani系数。

-11. SOME SWITCH FORM

*SOME SWITCH FROM BENCHMARK
'b2stbc_istyle_fchi'    '1' # 置一，从边界条件明确采用粒子流表达式（bx*cs*na）而不是fna。使用边界条件 BCPOT=11.
'b2tqna_model_sig'      '1' # 置一，使用芯部常数密度计算反常传导sig0=dfsig*qe*ne(nmdpl,-1)(nmdpl-中平面index)，置零，sig0=dfsig*qe*ne(x,y)(dfsig用于计算扩散对应的反常电传导)

'b2mndt_av'                          '1'             # 大于零，计算平均
'b2mndt_av_continue'                 '1'             # 大于零continue from averaged solution