It is used to 'correct' the CX data 0 => do not fix 1 => only fix H if CX data is < 1e-40 [default] 2 => fix if CX data is < 1e-40 3 => fix H 4 => fix all At the moment the only species with CX data in ADAS is C. Be careful with options that use the fit formula other than for H. See the comments in ratstr.F for the origin of the fit formula used to 'fix' the CX data.
When set to 1, disables the correction of the hydrogen atomic rate data with the Weisheit data using ratwei. Must be set to 1 when using b2ardr_amjhydhel. *** Use with caution! ***
When not set to 0, reads reaction cards from the AMJUEL/HYDHEL database files to overwrite rates for neutral hydrogenic species. The list of reaction cards to be used should be given in the b2ar.dat file. Cannot be used in conjunction with no_weisheit.eq.0.
When set to 1, disables the tail smoothing for of rate coefficients to represent the effect of a power-law tail on the distribution function. It has been added for compatibility with SOLPS4.3
Switch for choosing extrapolation method of ADAS data: If adas_extrap.eq.0, a simple 2-D linear extrapolation in log-log space is used (may lead to spurious results). If adas_extrap.eq.1, use the extrapolated result, but growing no faster than the parameter variation to the 3/2 power (default). If adas_extrap.eq.2, use the nearest domain edge value.
When fix_recomb is changed from 0, and if the ADAS option is in use, then uses the PRB file to add the contribution arising from is-->is-1 processes. If the AMJUEL/HYDHEL rates are in use and fix_recomb.eq.0, then the recombination energy is substracted to the electron recombination energy costs. May lead to negative rates! This term includes the bremsstrahlung. For ADAS, the '0' option only contains the bremsstrahlung for the is-->is-1 process. The correction is only applied for isonuclear sequences with nuclear charge up to zmax_recomb, inclusive. Use of this option is still experimental!
Number of AMJUEL/HYDHEL reactions you explicitly specify in b2ar.dat, when b2ardr_amjhydhel = 1. A set of default reactions is taken but can be overwritten by the reactions you specify here.
The number of intervals in the discretisation of the density (rtnn) and temperature (rtnt) ranges for the atomic physics tables.
Modify the atomic rates by including a hot electron population of temperature te_hot (in eV) and a population faction ne_hot_frac. Still experimental.
Lower temperature bound (in eV) used for rates imported from AMJUEL/HYDHEL database for neutral hydrogenic species. Rates below this bound are given the same value as they have at t_min. Only operational if within the bounds provided by tlohi in the b2ar.dat file header.
Upper temperature bound (in eV) used for rates imported from AMJUEL/HYDHEL database for neutral hydrogenic species. Rates above this bound are given the same value as they have at t_max. Only operational if within the bounds provided by tlohi in the b2ar.dat file header.
Lower density bound (in m^{-3}) used for rates imported from AMJUEL/HYDHEL database for neutral hydrogenic species. Rates below this bound are given the same value as they have at n_min. Only operational if within the bounds provided by nlohi in the b2ar.dat file header.
Upper density bound (in m^{-3}) used for rates imported from AMJUEL/HYDHEL database for neutral hydrogenic species. Rates above this bound are given the same value as they have at n_max. Only operational if within the bounds provided by nlohi in the b2ar.dat file header.
contains the file name of the geometry file to be read. The file will be looked for in the run directory, in the ../baserun directory, in $SOLPSTOP/data.local/meshes, and in $SOLPSTOP/modules/Carre/meshes/$DEVICE. To be used in b2ag.dat.
Cell index, on the basis mesh, of the outer midplane. Default value depends on geometry: Lower Single-null : jxa=rightcut(1)-(rightcut(1)-leftcut(1))/4, i.e. three quarters of the way between the two cuts. Upper Single-null : jxa=leftcut(1)+(rightcut(1)-leftcut(1))/4, i.e. one quarter of the way between the two cuts. Stellarator island : jxa=leftcut(1)+(rightcut(1)-leftcut(1))/2, i.e. midway between the two cuts. Double-null : jxa=(rightcut(1)+rightcut(2))/2, i.e. halfway between the two outer cuts. Limiter geometry : jxa=nx/2 Straight geometry : jxa=nx/2 Superseded by the RZOMP definition from b2.user.parameters.
Cell index, on the basis mesh, of the inner midplane. Default value depends on geometry: Lower Single-null : jxi=leftcut(1)+(rightcut(1)-leftcut(1))/4, i.e. one quarter of the way between the two cuts. Upper Single-null : jxi=rightcut(1)-(rightcut(1)-leftcut(1))/4, i.e. three quarters of the way between the two cuts. Stellarator island : jxi=leftcut(1)+(rightcut(1)-leftcut(1))/4, i.e. one quarter of the way between the two cuts. Double-null : jxi=(leftcut(1)+leftcut(2))/2, i.e. halfway between the two inner cuts. Limiter geometry : jxi=nx/4 Straight geometry : jxi=nx/4 Superseded by the RZIMP definition from b2.user.parameters.
Flux surface index, on the basis mesh, of the active separatrix. Can only be provided for slab geometries, otherwise deduced from geometry.
If pbs_from_basis_mesh.eq.0, the pbs array is computed from values on the actual working mesh (old treatment). If pbs_from_basis_mesh.eq.1, the pbs array is computed from values on the (possibly) finer basis mesh from the geometry file and the contributions from each basis cell are added to obtain the value on the working mesh. In both cases, mind the value of 'b2news_area_fix', which should be also declared in b2ag.dat if the non-default behaviour is wanted.
When b2agdr_redef_pbs.eq.1, geometrical quantities are adjusted so as to ensure that the poloidal flux between two flux surfaces remains constant.
Obsolete. Should use b2agdr_redef_pbs instead.
periodic_bc specifies if this is either an island or limiter geometry. If periodic_bc.eq.1 then island/limiter treatment is turned on. We differentiate between the two case through nncut: nncut.eq.0 = limiter case nncut.ge.1 = island divertor case (there should be nncut islands then) The case periodic_bc.eq.-1 is used to remove tranverse physics in 1-D cases. If periodic_bc.eq.0 then we have the usual treatment. No other values are currently allowed. To be used in b2ag.dat.
isymm specifies the type of symmetry of the geometry: isymm.eq.0 implies a slab geometry, isymm.eq.1.or.isymm.eq.2 imply rotational symmetry about the crx=0 axis, and isymm.eq.3.or.isymm.eq.4 indicate rotational symmetry about the cry=0 axis. Other values are not allowed. isymm.eq.1.or.isymm.eq.3 mean a toroidal geometry case, i.e. with a toroidal (out-of-plane) magnetic field component. isymm.eq.2.or.isymm.eq.4 mean a cylindrical geometry case, i.e. no toroidal magnetic field component. To be used in b2ag.dat.
coreregno is the boundary index of the core boundary in the input files b2ah.dat and b2mn.dat. For standard single-null and double-null cases, the default value of coreregno is 1. For a straight geometry or limiter case, the default value of coreregno is 0. coreregn2 is the boundary index of the second core boundary in case of a double-null geometry, and is not used otherwise.
pfrregno1 is the boundary index of the first half of the private flux region boundary on which feedback is being applied (through b2stbc_ndes or b2stbc_private_flux_puff) in the input files b2ah.dat and b2mn.dat. pfrregno2 is the second half of the same private flux region. For a standard single null, the pair should be (0,2). For a standard double-null, the lower private flux region is labelled with the pair (0,5), and the upper private flux region with the pair (2,3). One can use only one private flux region by setting pfrregno1.eq.pfrregno2.
solregno is the boundary index of the SOL North boundary in the input files b2ah.dat and b2mn.dat. For a standard single-null case, solregno is 3. For a standard double-null case, solregno for the outer SOL North boundary would be 7.
For 1-D cases, gives the width of the domain (in m) in the third (toroidal) dimension.
When the first flight model is turned on, this number must be greater than or equal to the num ber of boundary cells on the mesh. Only in cases with special geometries (i.e. limiters, islands, ...) need it be increased beyond the default value above. This variable is only used within the first flight module.
Distance used as the matching criterion when reading the geometry file.
If the mesh is offset (See b2agfs_xoffset, b2agfs_yoffset below) and Bt_adjust.eq.1, then the magnetic field is recomputed for toroidal geometries, keeping the pitch constant but changing the toroidal field magnitude according to Bt = Bt0/R. To be used in b2ag.dat.
The magnetic field will be multiplied by Bt_rescale. All components of the field are scaled together. To be used in b2ag.dat.
The magnetic field line pitch will be multiplied by pit_rescale.
This means that the poloidal field component is multiplied by pit_rescale, while the toroidal field component is left unchanged. To reverse the plasma current direction, set pit_rescale to -1.0.
The sign convention used is that a positive poloidal field points in the direction of increasing
If Bt_reversal.eq.1, then the sign of the toroidal field component only is reversed. The sign convention used is that a positive toroidal field leads to a Bx(grad(B)) direction pointing down, so B is out of the page when looking at the usual view of the poloidal plane with the center line on the left. To be used in b2ag.dat.
Minimum allowed value for the pitch angle (in degrees) at the plates.
xoffset and yoffset are offsets of the basis mesh in the x- and y-direction respectively. The mesh will be translated by (xoffset,yoffset) if any of the two is non-zero from its position given in the local_sonnet (See 'b2agfs_geometry' above) file. To be used in b2ag.dat.
Rescaling factors of the x- and y- coordinates of the basis mesh. To be used in b2ag.dat.
Number of topological cuts in geometry. Needed only to force a different value from the one computed automatically by b2ag.
Parameters describing the position and extent of the first topological cut, i.e. corresponding to the first X-point. Needed only to force different values from the ones computed automatically by b2ag. leftcut is the poloidal index of cells directly to the RIGHT of the left branch of the cut. rightcut is the poloidal index of cells directly to the RIGHT of the right branch of the cut. bottomcut is the radial index of cells directly below the cut. topcut is the radial index of cells directly above the cut.
Parameters describing the position and extent of the second topological cut, i.e. corresponding to the second X-point. Needed only to force different values from the ones computed automatically by b2ag. leftcut2 is the poloidal index of cells directly to the RIGHT of the left branch of the second cut. rightcut2 is the poloidal index of cells directly to the RIGHT of the right branch of the second cut. bottomcut is the radial index of cells directly below the second cut. topcut is the radial index of cells directly above the second cut.
Range of an optional isolated region to be included in the geometry. The isolated region will extend over the cell range spanned by (nxiso1:nxiso2,nyiso1:nyiso2). If only some of the four coordinates are given, the region is made to extent from the given coordinates to the end of the grid. Neighbourhood arrays and region indices are automatically adjusted.
Switch for conversion with b2us; can be specified in the b2mn.dat file read by b2us during case conversion to unstructured format. When set to 1, unstructured b2fgmtry_us file is written by b2us. When set to 0, it is not.
Switch for conversion with b2us; can be specified in the b2mn.dat file read by b2us during case conversion to unstructured format. When set to 1, unstructured b2fstati_us file is written by b2us. When set to 0, it is not.
Switch for conversion with b2us; can be specified in the b2mn.dat file read by b2us during case conversion to unstructured format. When converting a structured case with b2us and conv_triangles_old_co = 1, the coordinates of the nodes of the triangles from the original fort.33 file will be used in the unstructured format. Otherwise, the triangular mesh will be exactly matched to the B2.5 grid. Setting the switch to 1 may be required when converting triangle grids that include extreme slender/skewed triangles (typically manifested through FOLNEUT errors at runtime after case conversion).
Switch for conversion with b2us; can be specified in the b2mn.dat file read by b2us during case conversion to unstructured format. When set to 1, b2us will automatically try to reconstruct outer and inner midplane from the geometry (magnetic field and cell coordinates), and write the corresponding coordinates in b2.user.parameters. When set to 0, the IMP and OMP definition will be based on the values of b2mwti_jxi, b2mwti_jxa, b2mwti_jsep.
If abs(cos(alpha)*bx).ge.b2us_prep_Qalfmin, the beta angle is assigned to zero during writing of b2fgmtry. If abs(cos(alpha)*bx).lt.b2us_prep_Qalfmin, boundary faces are treated like field-aligned faces in BCMOM=13, BCCON=14, BCENE=15 (style=1), BCENI=15 (style=0) and BCPOT=11. Leakage coefficients (if they are needed) are set equal to b2us_prep_Qalfmax.
If abs(cos(alpha)*bx).ge.b2us_prep_Qalfmin .and. abs(cos(alpha)*bx).lt.b2us_prep_Qalfmax, fcPbs is set equal to b2us_prep_Qalfmax*fcS for boundary faces during writing of b2fgmtry.
Obsolete. Removed from code.
Obsolete. Removed from code.
If ramp_slow.eq.0 (default), the ExB and diamagnetic drifts multipliers are ramped on each call of b2news, otherwise they are only changed on the first call of the nstg loop on the timestep.
Master switch to turn on the 9-point stencil numerical treatment. If set to 0, reverts to the old 5-point stencil. Only usable with the SOLPS5.2 physics model. The perpendicular viscosity current terms are not yet available within this treatment.
Minimum residuals for stopping a run. All criteria must be met simultaneously for the run to stop. The continuity equation residual applies to the "ismain" species (See "Run" section).
Threshold of the maximum residual among the equations solved below which the run is stopped. Only meaningful when compiled with FIXED_POINT and with adjoint AD code.
Minimum residuals for an internal solution loop to stop. Used in conjunction with (not in replacement of) the nstg(0:2) numbers. All nonzero criteria must be met simultaneously. The continuity equation residual applies to the 'ismain' species (See 'Run' section).
Allows for solving only for the species range [nsmin:nsmax-1].
Upper and lower bounds used when writing out data for Eirene.
If the Mach number based on the total fluid velocity of a background species exceeds eirene_M_max, all velocity components passed to Eirene are restricted to make the Mach number equal to eirene_M_max. It is an additional restriction on top of eirene_ua_min and eirene_ua_max.
Print min and max values of te, ti, na & ua
If eirene_extrap.eq.1, then guard cell plasma parameter values are calculated from the neighbouring real cell values. If eirene_extrap.eq.0, the guard cell values are used unchanged.
If ank_mods.ne.0, then uses an additional scheme to ensure particle balance as the B2.5 solution evolves, due to the internal iteration scheme, away from the plasma background on which the Eirene sources were originally computed at the beginning of the time step. The user is referred to the text in $SOLPSTOP/doc/Source_Scaling_in_B2.pdf for a full description of the method used. Cannot be used in conjunction with 'eirene_ionising_core'.
If dpc_fix.eq.1, then uses the true particle source from the stratum. If dpc_fix.eq.2, sets this particle source to zero. The equivalent of the old behaviour is dpc_fix.eq.0 and is wrong!
If b2stbr_eir_src_nhist.gt.1, then Eirene sources are accumulated and a moving average is computed. Note that only one of b2stbr_eir_src_nhist, eirene_neutr_avg, or eirene_underrelax may be nonzero.
If eirene_neutr_avg.gt.0, then Eirene sources are accumulated and averaged until eirene_neutr_avg+1 steps, at which point the averaging process is reset. Note that only one of b2stbr_eir_src_nhist, eirene_neutr_avg, or eirene_underrelax may be nonzero.
If eirene_underrelax.gt.0, then an underrelaxation scheme is used for the Eirene sources, with an underrelaxation ratio of 1/eirene_underrelax. Note that only one of b2stbr_eir_src_nhist, eirene_neutr_avg, or eirene_underrelax may be nonzero.
If eirene_uub_style.eq.0, compute poloidal velocities passed to Eirene from ua, vadia, wadia, vaecrb arrays. If eirene_uub_style.eq.1, compute poloidal velocities passed to Eirene from fna and fna_eir fluxes.
Obsolete. Removed from code. New implementation supersedes treatment of area_fix=3.
If ion_vlct_restrict.ne.0, the parallel speed of ions and fluid neutrals is limited to be no more than M times the global plasma sound speed. In addition, if ion_vlct_restrict.eq.2, the ExB ion speeds are limited to be no more than M the local ion sound speed. M is given by b2npmo_ion_vlct_restrict_M.
Determines M for ion_vlct_restrict.ne.0.
Puts a hard limit on the absolute value of the parallel velocity, but applies only to fluid neutrals and in case the spatially hybrid fluid-kinetic model is active (spatial_hybrid.ne.0). This switch was introduced when b2npmo_ion_vlct_restrict only applied to ions. Now that b2npmo_ion_vlct_restrict also applies to fluid neutrals, the removal of b2mndr_ua_max can be considered. For now it is kept for backward compatibility with published results.
If b2nppo_restr_po.ne.0.0, the electric potential is restricted such that its ratio (in absolute value) with the electron temperature (in eV) does not exceed b2nppo_restr_po. If b2nppo_restr_po is negative, the limit only applies to regions of negative electric potential. Only applies on open field lines.
If b2upht_rte_min.gt.0.0, the electron to ion temperature ratio is forced to not be lower than rte_min (i.e. Te/Ti >= rte_min).
If b2sihs_shivis.eq.0, then shivc and shiva (ion heat source due to viscosity) will be zeroed out in boundary-volume cells.
If b2sihs_shi.eq.0, then all b2sihs_shi* and b2sihs_shivis keys will be assigned to zero.
If b2sihs_shidu.eq.0, then shidu (ion heat source due to divergence of velocity) will be zeroed out in boundary-volume cells.
If b2sihs_shidd.eq.0, then shidd (ion heat source due to divergence of drift velocity) will be zeroed out in boundary-volume cells.
If b2sihs_shedu.eq.0, then shedu (electron heat source due to divergence of velocity) will be zeroed out in boundary-volume cells.
If b2sihs_shedd.eq.0, then shedd (electron heat source due to divergence of drift velocity) will be zeroed out in boundary-volume cells.
vis_per is a multiplier to the perpendicular viscosity current. If vis_per.ne.0.0, the electric potential equation is solved by the subroutine b2npp7 using a 7-point stencil. The value '1.0' is recommended for runs with drifts when perpendicular viscosity and corresponding perpendicular viscosity current is taken into account. Cannot be used in conjunction with the 9-point stencil numerical treatment. If used with the 5-point stencil treatment and with fluid neutrals, must be equal to b2tfnb_vis_per for numerical stability reasons.
If sheath_drift_fix.eq.0, then the drift velocity used in the sheath boundary conditions is the sum of diamagnetic and ExB contributions (old, but likely wrong, treatment). If sheath_drift_fix.eq.1, then the drift velocity used in the sheath boundary conditions is the ExB velocity only (recommended).
This switch concerns a 9-point correction for the Neumann type boundary conditions (BCCON/BCMOM/BCENE/BCENI/BCPOT/BCENK/BCENZ = 2) for non-orthogonal cells. If b2stbc_BC2_cor9.eq.0.0, gradients tangential to the boundary are ignored. The actual normal gradient then deviates from the one requested by the user. If b2stbc_BC2_cor9.eq.1.0, an additional term based on the vertex values of the boundary face is added to correct for the tangential gradient, to get the correct (up to interpolation error) normal gradient.
If fix_fch_in_fhe_sheath.eq.0, then the wrong (old) implementation of the FCH contribution to FHE used; If fix_fch_in_fhe_sheath.eq.1, then the wrong (second) implementation of the FCH contribution to FHE used; If fix_fch_in_fhe_sheath.eq.2, then the correct implementation of the FCH contribution to FHE used;
If potential_at_guard_cell.eq.1, the electric potential used for the computation of the incident energy for sputtering yields is taken as the value in the guard cell. If potential_at_guard_cell.eq.0, the value from the neighbouring real cell is used instead.
If temperature_at_guard_cell.eq.1, the temperatures used for recycling and reflection are taken as the value in the guard cell. If temperature_at_guard_cell.eq.0, the temperatures are interpolated between the guard cell and the neighbouring real cell.
When set to '1', flux limit of the parallel electron and ion heat fluxes is applied to transport coefficients. It is recommended '1'. If 'b2trcl_conductive_limit' '0' then the keys 'b2tfhe_lim_flux' and 'b2tfhi_lim_flux' must be '0'.
If core_cond_limit.eq.0, then the heat flux limits due to chvemx and chvimx are not applied in the core. Switch only active if 'b2trcl_conductive_limit' .ne. 0.
If '0', use the SOLPS5.0 scheme for neutral heat conductivity flux limits. If '1', SPb form to calculate parallel and perpendicular neutral heat conductivity flux limit (does not preserve symmetry, kept for backward compatibility reasons). If '2', modifies the SPb treatment for the flux limits to be applied on the transport coefficients directly. It is recommended '2'.
Switch to choose between various interpolation schemes for transport coefficients (Does not apply to pinch velocities vla or to the temperature-driven electric conductivity alf). mode.eq.-1: harmonic averaging mode.eq.0 : geometric averaging mode.eq.1 : arithmetic averaging (SOLPS5.0 formulation) mode.eq.2 : arithmetic averaging (SOLPS4.0 formulation)
Switches to toggle between the standard 5.0 hybrid scheme (no_hybr.eq.0) and the old 4.0 upwind scheme (no_hybr.eq.1). These apply separately to the electron heat, ion heat and particle conservation equations. If hybr2.ne.0, then a mixed scheme is used where only the conductivity is computed according to the 4.0 scheme. Obsolete. Not available for WG. Replaced by the 'discr_meth' switches.
If upwind.ne.0, then the SOLPS4 upwind scheme in b2tfhe_ and b2tfhi_ to obtain the heat fluxes. Obsolete. Not available for WG. Replaced by the 'discr_meth' switches.
If pflux_cor.eq.1 and fna_mdf is used, enforces that the integral particle flux across the domain boundaries computed from fna is equal to the same integral computed from the fna_mdf fluxes. Not yet available for WG.
Artificial coefficient providing faster convergence to the neoclassical electric field but giving a distortion of the flows in the SOL as a side-effect. Can be applied <>1 during the convergence and turned off for the final stage of calculations. Use with caution.
Number of expected nonzero matrix elements per matrix row.
Multipliers to the number of nonzero elements in the solution matrix for workspace arrays in the matrix solver.
Choose the type of 5-point stencil matrix solver. Style.eq.0 = iluter, Style.eq.1 = slv5pt, Style.eq.2 = MA28. NOTE: Only style.eq.2 will give good results. Other values are NOT recommended! style.eq.0 and style.eq.1 are only applicable to linear geometries with no cuts and no isolated regions.
Choose the type of 7-point stencil matrix solver. Style.eq.0 = iluter, Style.eq.1 = slv5pt, Style.eq.2 = MA28copy3. Style.eq.3 = SDRV from YSMP. NOTE: Only style.eq.2 will give good results. Applying style.eq.3 should be corrected and is no longer recommended. It might give slow convergence or even divergence of the potential equation. Other values are NOT recommended! style.eq.0 .or. style.eq.1 will return an error.
Choose the type of 9-point stencil matrix solver. Style.eq.0 = iluter, Style.eq.1 = slv5pt, Style.eq.2 = MA28copy. NOTE: Only style.eq.2 will give good results. Other values are NOT recommended! style.eq.0 and style.eq.1 will return an error.
Paremeter needed for iluter matrix solver.
Multiplier to the neoclassical current and diamagnetic heat flux convective boundary conditions, also multiplied by facdrift. Adds a poloidal variation consistent with neoclassics and diamagnetic contributions to the heat flux boundary conditions. Should only be used with the 5.0 model and heat flux boundaries. Allows use b2stbc_integral_current if fchy_dia.eq.0, forbids it otherwise. Not yet available for WG.
If fchy_dia_coreonly.eq.1, the neoclassical current and diamagnetic heat flux convective boundary conditions are only applied on the core boundary. If fchy_dia_coreonly.eq.0, those are applied to all boundaries, including PFR, North SOL edge and diverted island center boundary. Not yet available for WG.
Real parameter which establishes boundary conditions for radial component of the current for North and South boundaries when these lie on closed flux surfaces. If b2stbc_neoclassical is 0 then the radial component of the current is zero. If b2stbc_neoclassical is not 0 the radial component of the current is given by the neoclassical value. b2stbc_neoclassical is superseded if facdrift.ne.0. This cannot be used in conjunction with b2stbc_integral_current below.
Multiplier to the ExB velocity for sheath boundary conditions in b2stbc_spb and BCMOM=13 case of b2stbc_phys. Not yet available for WG.
If integral_current.gt.0, corrects current sources so that there be no surface current at the North and South edges of the edge plasma. The value of integral_current multiplies the correction term added to the current source. This setting is incompatible with the normal use of diamagnetic drift terms (facdrift.gt.0.0 .and. b2stbc_fchy_dia.ne.0.0) or neoclassical boundary conditions (b2stbc_neoclassical.gt.0.0). Not yet available for WG.
Step length (in metres) for computing the first flight model chords.
Number of chords started from each start point in the first flight model.
Workspace size given to the first flight table.
Small level sources introduced in all cells.
Frequency (in number of calls to b2stbc_phys) at which the constant value of the potential at the boundaries where BCPOT=16 is applied will be recomputed. Not yet available for WG.
Stabilizing term to the source associated with boundary condition BCENE=15, when bcene_15_style.eq.1 (default). To help suppress electron temperature oscillations, this parameter can be set to a positive value. The range 1-100 is recommended.
Stabilizing term to the source associated with boundary condition BCENI=15. To help suppress ion temperature oscillations, this parameter can be set to a positive value. The range 1-100 is recommended.
Stabilizing term to the heat sources associated with ion-electron energy exchange.
Minimal density maintained in all cells for all species (in m^{-3}). It must be true that na_min is smaller than na_new, as well as any of the initial densities provided in b2ai.dat. This switch replaces 'b2mndr_na0eps' from SOLPS5.x. It is also used to avoid division by 0 when calculating kinetic neutral velocities from the kinetic flux densities in b2tinnt.F.
Initial density (in m^{-3}) put in all cells for all new species if not overwritten by initial state file. This switch replaces 'b2mndr_na0eps' from SOLPS5.x.
Sets the minimum electron temperature value (in eV).
Sets the maximum electron temperature value (in eV).
Sets the minimum ion temperature value (in eV).
Sets the maximum ion temperature value (in eV).
Sets the minimum fluid neutral temperature value (in eV).
Sets the maximum fluid neutral temperature value (in eV).
Sets the minimum turbulent ExB kinetic energy value (in eV).
Sets the maximum turbulent ExB kinetic energy value (in eV).
Sets the minimum turbulent ExB enstrophy value (in s^{-2}).
Sets the maximum turbulent ExB enstrophy value (in s^{-2}).
Prohibited switch. As of v3.2.0, there are no faces between guard cells anymore.
Parameters for spatial variation of the diamagnetic, ExB and viscosity drift term multipliers. These are allowed to have tanh profiles centered at tanh_a and of width tanh_b, where both quantities are in metres and measured along the outer midplane. The maximum value of facdrift and fac_ExB is controlled by their respective _start and _target switches, see Run section for details. NOTE: none of them have been converted to WG! For consistency with the naming convention of related variables: b2news_facExB_tanh_a is an alias for b2news_fac_ExB_tanh_a, b2news_facExB_tanh_b is an alias for b2news_fac_ExB_tanh_b, b2news_facvis_tanh_a is an alias for b2news_fac_vis_tanh_a, and b2news_facvis_tanh_b is an alias for b2news_fac_vis_tanh_b. Additionally, for all deep core cell rings up to and including the iy_nocoreExB-th, fac_ExB(ix,iy) is set to zero. This means that ExB contributions and flows will be zeroed out in those cells. For consistency with the boundary condition treatment, the value iy_nocoreExB.eq.1 is NOT recommended!
When set to '2', SPb form of adding currents on the South core boundary is included by using BCPOT=12, and on the South PFR and North boundaries by using BCPOT=13. It is recommended '2', in conjunction with the BCPOT=12 and BCPOT=13 settings, respectively. When set to '1', SPb form of adding currents on South and North boundaries is done explicitly in addition to other existing boundary conditions for the potential equation. The old 5.0 calculation is recovered by using the value '0'.
If '1', explicitly use the expression (bx*cs*na) of particle flux (fna) from boundary condition instead of fna, for boundary condition type BCPOT=11. Not yet available for WG.
If no_b2sral_call.eq.0, add an additional call to recompute the sources in b2news_ to reproduce the behaviour from SOLPS4.
If do_2nd_b2npco_call.eq.1, perform a second call to the density equation solve to improve particle balance, as done in SOLPS4.
If re_eval_prtls_fluxes.eq.1, the particle fluxes are recomputed at the end of b2news_. This is necessary for rescaling of the Eirene sources during coupled runs, so this switch is superceded by use_eirene, and also needed to reproduce SOLPS4 runs.
If 0, refers to the obsolete SOLPS5.0 physics model, no longer supported. If 1, the SOLPS5.2 physics model is used, and b2news_ is called to compute the new plasma state. (default) If 2, the SOLPS5.2 physics model is used, and b2news_m is called instead to run the code in fully self-consistent time-dependent mode. B2news_m uses a different order of solving equations: first the density equations, then the parallel velocity, electrical potential and Te and Ti equations (and optionally also the Tn and k-epsilon equations, if requested).
When set to a value of 'K', residuals at each K-th step will be written in b2ftrace. It helps to avoid a huge b2ftrace files during long time run.
Feedback relaxation parameters. See Physics section for individual feedback quantities.
Additional feedback strength adjustment parameters for BCCON=13 boundary condition, density adjusted to match a specific flux. The density is adjusted by a factor of: (1.0_R8+CONPAR(IS,IB,2)*(DFS-FS)/ (abs(DFS)+abs(FS)+abs(type13_norm)+abs(NAS*type13_fac))) where FS is the current particle flux, DFS is the desired particle flux, and NAS the current volume-averaged density. The desired flux DFS is computed as CONPAR(IS,IB,1)+CONPAR(IS,IB,3) where the first number is provided directly by the user in b2.boundary.parameters and the second is internally set to correspond to the re-entry of ionised neutrals that have crossed into the core for cases run with Eirene where eirene_ionising_core is activated. See also the description of 'eirene_ionising_core'. Not yet available for WG.
Parameters related to solving the temperature equation for the target plate elements in depth. If linlog.eq.1, a linear subdividing of the plate element is used. If linlog.eq.2, a logarithmic subdividing is used, with the surface layer being ratio times thinner than the last layer in the bulk. ratio must be larger than 1.
If '1', flux limit to neutrals contribution to dpa0 - the diffusion coefficient is applied in b2tlc0.F. It is recommended '1'. b2tlc0.F has the flux limit parameters α and γ which are given by 'b2tlc0_alpha' and 'b2tlc0_gamma'. 'b2tfnb_alpha' and 'b2tlc0_alpha' cannot be different from zero simultaneously. 'b2tfnb_alpha' gives another form of flux limit which is applied to the whole particle flux.
If '1', flux limit to neutrals contribution to vsa0 - the viscosity is applied in b2tlv0.F. It is recommended '1'. b2tlv0.F has the flux limit parameters α and γ which are given by 'b2tlv0_alpha' and 'b2tlv0_gamma'.
Experimental: allows for density feedback of the number of electrons accompanying species 'isfb' at location (ixfb,iyfb). Note that ixfb and iyfb are real numbers!
Switches off the stabilization of the source coefficients.
Main under-relaxation parameter. Should be set to 1.0 (default) if running in time-dependent mode (b2mndt_style.eq.2). Otherwise, default value is 0.5. Setting this switch to zero prevents updating of the principal code quantities.
Under-relaxation parameters for the parallel momentum, continuity, potential and energy equations respectively. Setting these switches to zero prevents updating of the associated code quantities.
Specifiy multipliers to Rhie-Chow corrections from particle equations Not used (and default set to zero) if running in time-dependent mode (b2mndt_style.eq.2). b2npco_pcm0 is a multiplier for the update on the parallel velocity following the pressure-correction equation. b2npco_pcm1 is a multiplier for the Rhie-Chow terms in the particle fluxes.
rxg specifies a special under-relaxation parameter. A term of the form "abs(residual)/rxg" is added to the diagonal of the matrix of the correction equation, with the effect of limiting the computed correction. rxg is dimensionless and of order unity; smaller values of rxg imply stronger damping. Not used if running in time-dependent mode (b2mndt_style.eq.2).
Specifies multipliers to Rhie-Chow corrections from heat equations. Not used (and default set to zero) if running in time-dependent mode (b2mndt_style.eq.2). b2npht_pcm0 is a multiplier for the update on the parallel velocity following the update of Ti and Te. b2npht_pcm1 is a multiplier for the Rhie-Chow terms in the flow fields of ion, electron, and neutral heat.
rxg specifies a special under-relaxation parameter. A term of the form "abs(residual)/rxg" is added to the diagonal of the matrix of the correction equation, with the effect of limiting the computed correction. rxg is dimensionless and of order unity; smaller values of rxg imply stronger damping. Not used if running in time-dependent mode (b2mndt_style.eq.2).
When set to '1', SPb form of the program b2sihs_ is called. b2npht_style.eq.0 no longer available for WG.
Normalisation factor for the parallel momentum equation. Not used when running in time_dependent mode (b2mndt_style.eq.2).
If poteq.eq.0, the potential equation is jumped over and not solved. If poteq.eq.2, the potential is set to 3.1*Te/qe as per SOLPS4.0. If poteq.eq.1, the potential equation is solved according to the no_solve switch settings. If poteq.ne.1, then 'b2tfhe_no_current' must be set to '1'.
When set to '1', the total friction force cancel is not calculated at the guard boundary cells. It is recommended '1'.
style specifies the precise form of interpolation used in the computation of flcb and cvcb. When set to '1', SPb form of the transport terms in the momentum correction equation is used. It is recommended '1'.
When set to '1', SPb form of an expression that occurs in the electron-atom thermal force is used. Only used if b2sigp_style.ne.2. It is recommended '1'.
Specify the type of linearisation used in the thermal force term.
Switch to choose between various discretization schemes when computing the poloidal electric field. See code for details.
Switch to choose between various discretization schemes when computing the radial electric field. See code for details.
When set to '1', SPb form of the pressure gradient term on the right hand of the momentum balance equation is used. When set to '2', the parallel current term contains a correction due to impurities, when Z_eff is not equal to the average plasma ion charge. Only usable if the potential equation is solved simultaneously (b2news_poteq.eq.1). The default value of '2' is recommended for cases with one main hydrogenic species (or at least where the main species is the lightest one) and any amount of impurities. It is not guaranteed to yield correct results for hydrogen isotopic mixtures or non-hydrogenic plasmas.
When set to '1', a CFL-like restriction is applied to the pressure gradient term for minority ions and for neutrals species.
If 1 then zero passed sources in dead regions, If 2 zero passed sources in dead regions and core boundary cells If 0 then SKIP
rf0-4 are numerical parameters, see the code. All must be non-negative and of order unity; larger values are meant to provide more stabilisation. Some or all may be set to zero and removed eventually. I have had some experience to indicate that rf0, rf2, rf4 may have a useful function; however, later changes in the code may have made all these parameters superfluous.
Multipliers to the density, parallel momentum, heat, potential, and electron particle time-derivative source terms, respectively. 'phm4' is normally zero since the potential equation contains no source terms.
(rf0/1/2/3 for numerical stabilisation)
(rg0 for numerical stabilisation; needs experiments.)
Specifies the type of linearisation used in the charge exchange momentum source term.
Specifies the type of linearisation used in the sources derived from the Monte Carlo neutrals. If linearisation.eq.0, all sources are fed to the constant term, if linearisation.eq.1, positive terms go in the constant term, and negative terms in the proportional term. 'eirene_mc_linearization' is an alias for this switch.
Same as above, but applied to the sources coming for an externally coupled code providing plasma sources, for instance a kinetic treatment of trace impurity ions. 'b2stbm_linearization' is an alias for this switch.
Specifies the frequency (in units of full B2.5 timesteps) at which an externally coupled code providing additional source (for instance a kinetic treatment of trace impurity ions such as IMPGYRO) should be called.
Damping coefficient applied to the atomic physics sources.
Specifies the type of linearisation used in the charge exchange momentum source term.
Multiplier to the pressure force term.
If '1', SPb new form of calculating electron heat flux is used. It is recommended '1' for runs with drifts.
If no_current.eq.1, all currents are set to zero. The setting no_current.eq.1 must be used if the potential equation is not explicitly solved for, i.e. for b2news_poteq.ne.1.
Multiplier to the parallel current.
If '1', SPb new form of calculating ion heat flux is used. It is recommended '1' for runs with drifts.
When set to '0', drift velocities are calculated in cell centers (NO LONGER SUPPORTED IN WG CODE). When set to '1', drift velocities are calculated in cell faces. It is recommended '1'.
Multiplier to the poloidal electrical drift velocity.
Multiplier to the radial electrical drift velocity.
Multiplier to the poloidal diamagnetic drift velocity.
Multiplier to the radial diamagnetic drift velocity.
If '1', SPb new form of calculating particle flux is used. It is recommended '1' for runs with drifts.
(stylec is a numerical switch, needs experiments)
Linearisation constant.
If enable.ne.0, then the sources for particle, parallel momentum, electron and ion heat are rescaled, for non-boundary cells, if the local timestep exceeds the physics timescale implied by those sources. enable.ne.0 not yet available for WG.
Integer switch to choose the numerical treatment of the ion pressure gradient calculation for the ion-neutral current. If b2tinnt_style.eq.0, then the chain rule is first applied to the product of density and temperature, and (grad n)/n is rewritten as grad(log(n)). If b2tinnt_style.eq.1 (default), then the gradient of the product of density and temperature is calculated directly.
Determines the linearisation of the boundary condition heat source due to neutral recycling/reflection when erecyc in used. Not used for advanced fluid neutral (AFN) model. lin_shi0.eq.0: put source in shi0(iCv,0). lin_shi0.eq.1: put source/ti in shi0(iCv,1).
Determines the discretization method for the calculation of the fluxes. discr_meth.eq.0: central scheme discr_meth.eq.1: upwind scheme discr_meth.eq.2: SOLPS5 hybrid scheme discr_meth.eq.3: SOLPS4 continuity scheme. Not yet available for WG.
Determines the interpolation weights to get the parallel velocities at the cell faces. style_int_vel.eq.0: interpolation based on cell volumes. style_int_vel.eq.1: interpolation based on cell connector lengths. Only the combination b2tfnb_style_int_vel=0 and b2sihs_style_int_vel=0 works robustly for now. Other combinations gave rise to oscillations. More research needed.
Determines the interpolation weights to get the parallel velocities at the cell faces for the calculation of divergence ua. style_int_vel.eq.0: interpolation based on cell connector lengths. style_int_vel.eq.1: interpolation based on cell volumes. Only the combination b2tfnb_style_int_vel=0 and b2sihs_style_int_vel=0 works robustly for now. Other combinations gave rise to oscillations. More research needed.
If larger than zero, specifies the number of optimization iterations between writes of the intermediate state file b2fstate_optim.XXXX.
Percentage value at which drifts are reset at each function/gradient evaluation to avoid simualtion crash.
When resetting drifts, the drift increase multiplier (e.g. fac_exb_inc) will be computed so that drifts are back at 100% in ntim/b2optim_reset_iter iterations.
Specifies the number of timesteps between writes of the b2time.nc time-dependent file. If b2time.gt.0, always writes out on the last timestep. If b2mndr_stim.lt.0, data from the current run is appended to the existing b2time.nc file, otherwise the file is overwritten.
Specifies the number of timesteps between writes of tallies. If tally.gt.0, always writes out on the last timestep.
Controls detailed monitoring output. moitlv.eq.-1 means no output, moitlv.eq.0 means output once per timestep, moitlv.eq.1 means output once per nstg0 iteration, moitlv.eq.2 means output once per nstg1 iteration and moitlv.eq.3 means output once per nstg2 iteration. See 'b2mndr_nstg?' description in Run section for more details.
Controls quick trace output. moqtlv.eq.-1 means no output, moqtlv.eq.0 means output once per timestep, moqtlv.eq.1 means output once per nstg0 iteration, moqtlv.eq.2 means output once per nstg1 iteration and moqtlv.eq.3 means output once per nstg2 iteration. See 'b2mndr_nstg?' description in Run section for more details.
Specifies the maximum number of instances at which movie data will be output.
Specifies the number of timesteps between movie output. If greater than 1 and the last timestep does not correspond to a planned frame, an extra endstate frame is added.
Specifies the maximum number of instances at which extra writes of b2fplasmf.xxxx will occur.
Specifies the number of timesteps between b2fplasmf.xxxx writes.
If non-zero, the movie file b2movies.nc is generated, containing several 2d quantities. Gives the real-time interval between writes of movie frames in b2movies.nc. If non-zero, must be at least as large as the simulation timestep.
Controls the quantities outputted to b2movies.nc. If cdfmovie_fields.ge.1 then only the main state variables (na, ne, Te, Ti, etc.) are written with each write to b2movies.nc. If cdfmovie_fields.ge.2 then additional quantities (e.g. fluxes, particle and energy sources) are also included.
Another option for plasma state file output. If greater than 0, gives the number of B2.5 full interations between successive write-ups of plasmastate files, in addition to the ones written at regular CPU intervals (see b2mndr_savecpu and b2mndr_plasmatim). Should be used, at the exclusion of other plasmastate write-up frequency settings, in conjunction with the Leuven Monte Carlo averaging scheme.
If b2mndt_av.gt.0, turns on computation of running averages.
If b2mndt_av_continue.gt.0, continuously perform running averages.
If ntim_batch.gt.0, number of iterations used to compute batch averages, written out in 'b2batch.nc'.
If ntim_run.gt.0, number of iterations used for writing running averages.
If b2mndt_av_batch_all.gt.0, produces standard output for batch averages.
Another option for plasma state file output. If greater than 0, gives the real-time interval between successive write-ups of plasmastate files, in addition to the ones written at regular CPU intervals (see b2mndr_savecpu and b2mndr_ntim_save). Must then be at least as large as the simulation timestep.
Frequency, in timesteps, at which an IDS is saved containing the plasma state. Can be used in conjunction with b2mndr_ids_time and b2mndr_ids_av. If used, an IDS time slice will be written at the end of run. If negative, a single IDS will be written at the end of the run. Will not work if running a "0 timesteps" simulation, in which case the use of b2_ual_write is recommended.
Frequency, in simulation seconds, at which an IDS is saved containing the plasma state. Can be used in conjunction with b2mndr_ids_save and b2mndr_ids_av. If non-zero, must be at least as large as the simulation timestep. If used, an IDS time slice will be written at the end of run. Will not work if running a "0 timesteps" simulation, in which case the use of b2_ual_write is recommended.
Frequency, in number of batch averages, at which an IDS is saved containing the averaged plasma state. Can be used in conjunction with b2mndr_ids_save and b2mndr_ids_time. The averaged plasma state is saved as a separate occurrence of the edge_profiles and edge_sources IDSs.
If iout.eq.1, a large set of *.dat output files will be produced containing the values of a variety of code quantities, evaluated at the end of the B2.5 iteration. If iout.eq.4, an even larger set of output files will be produced, for the purposes of a full run analysis, evaluated in the individual routines where the quantities are used. For even more detailed debugging analysis, one should instead use "procedure_name"_iout.eq.1. The files and their content are fully described in the Output_description.pdf file in the $SOLPSTOP/doc directory and Appendix G of the SOLPS-ITER manual.
If get_residuals.eq.1, then the residuals of unsolved equations will be outputted into b2ftrace. Equations are unsolved if 'b2news_no_solve' is nonzero or any elements of the solveco, solvemo, solveee, solveei, or solvepo are set to .false. . By default, get_residuals.eq.0 and the residuals for unsolved equations are set to zero. The residual of the total momentum equation is nonzero when all momentum equations for parallel velocities are solved together. The residual of the total energy equation is nonzero when Te and Ti equations are solved together.
If append.ne.0, the output files produced by the b2wdat_iout.eq.4 switch are appended upon every write, instead of being rewritten every time.
Specifies the number of significant digits with which the *.dat files are written out. Defaults to 6 in normal mode and 15 in debug mode. Must be positive.
If old_style.gt.0, old-fashioned (SOLPS4 style) output is added at the end of the b2mn.prt file.
If b2mndr_av_read.gt.0, the initial plasma is the averaged solution read from the b2faveri file. Check $SOLPSTOP/doc/SOLPS-ITER_Eirene_averaging.pdf for details.
The five switches above are safeguards numbers for when printing changes after a time step. The change is computed as: deltaX = abs((X(t)-X(t-1))/(X(t)+X_eps))
Produces a numbered 'tran' file every trantim real-time seconds. An endstate file is written if it falls between scheduled write-up times. Only available within the -DJET environment. If non-zero, must be at least as large as the simulation timestep.
The diagnostic values from b2time.nc use guard cell values if target_offset.eq.0, and values from the neighbouring real cell if target_offset.eq.1. Fluxes are not affected.
Controls additional output to b2time.nc. If eq 1 then a few 2d arrays (ne, Te, Ti) are written with each write to b2time.nc. If eq 2 then fluxes, electric potential, kinetic energy, and fluid particle and energy source terms are also included (e.g., rsana, rsahi, rqrad).
Index of the species used to create the 'dp3d?.last10' diagnostic files. Defaults to the neutral species associated with ismain. If no such species is declared, the default will be set to ismain. Is expected to either correspond to a neutral species or to ismain. Can also be used in b2ar.dat to specify the index of the main neutral CX species when resorting to the AMNS database (starting from version 1.3.2). If ismain is also defaulted, then will be 0.
Specifies the amount of data that is written out to b2ftrace. See the manual (Section on b2yq) for full details.
If tallies_netcdf.ne.0, the file 'b2tallies.nc' is created, which contains the regional tallies in CDF format. If b2mndr_stim.lt.0, data from the current run is appended to the existing b2tallies.nc file, otherwise the file is overwritten. If b2mndt_av_ntim_batch.gt.0, the file 'b2batch.nc' is created, which contains the batch averages of the tallies. If b2wall_netcdf.ne.0, the file 'b2wall.nc' is created, which contains the wall tallies in CDF format [written every b2wall_netcdf main calls]. If b2mndr_stim.lt.0, data from the current run is appended to the existing b2wall.nc file, otherwise the file is overwritten. If balance_netcdf.ne.0, the file 'balance.nc' is created, which contains all of the arrays required by the balance post-processing routines, in CDF format. If balance_average.ne.0, the balance arrays are averaged over the number of completed timesteps (itim).
If ank_tracing.ge.1, additional tracing output from Andrei Kukushkin is produced, in files to be found in the tracing/ directory inside the run directory. The traces will be written every ank_tracing iteration. If b2mndr_stim.lt.0, data from the current run is appended to the existing files, otherwise the files are overwritten.
Controls level of output in b2stbc and subservient routines. Level 1 (diagno.ge.1) output includes wrong_flow warnings, transport coefficients at the (nx/2,ny) basis mesh position, the solution found for the given heat flux at constant temperature bc, and all feedback diagnostics. Level 2 output (diagno.ge.2) includes the integral current correction, the boundary sources on the first call, excessive sound speeds, and the iterative progress to the given heat flux at constant temperature bc. Level 3 output (diagno.ge.3) includes the boundary sources on each call, and the boundary.prt file listing all boundary cells.
Output flag for the first_flight model. If output.ge.1, the files start, end, and fort.88 are created, which contain information about the chords used in the model.
Output flag to control diagnostics related to use of ion_vlct_restrict switch. If vlct_diagno.eq.0, the velocity restrictions are applied silently. If vlct_diagno.eq.1 (default), the user only gets a count of how many times the velocity restriction has been applied for each species, if any. If vlct_diagno.eq.2, the user gets the full details of where and how large the applied velocity restriction was, if any.
For coupled runs, the code will create two files named fort.30 and fort.31, containing the plasma geometry and plasma state, respectively, for use by Eirene. This happens irrespective of the setting of the two switches above. If however, the user wishes to obtain these files for a non-coupled run, then these two switches can be set to 1 and the files will then be produced at the end of the run. If savef31.lt.0, then a fort.31.###### file is produced after every B2.5 time step. Additionally, if savef31.ne.0, the fort.31 file will be updated before every Eirene call for coupled runs.
If inverse_ua.eq.1, the code will produce a 'b2fstati_with-ua' file that contains the same plasma information, but with the opposite sign convention for the parallel velocity. The run will then proceed with the new velocities with their sign inverted.
New number of species desired when reading a new atomic rates file using b2yr.exe. Must be specified within b2yr.dat. To be used for checking atomic rates after bundling and before any plasma state files have been created or converted, i.e. when the number of species in the atomic rates file b2frates does not match the number of species declared in the run parameters file b2fpardf.
Controls level of output in b2srsm. If diagno.ge.2, output will be given on every call. If diagno.eq.1, output will be given on main calls only.
If cpu.gt.0, prints out the time spent in the matrix solver.
Output flag for the matrix solvers. Larger numbers mean increasing output level.
If nonzero, prints sources computed by Eirene. If greater than 1, prints them on every use of the recycling sources, not just when they are computed.
If nwrite.gt.0, prints the content of the sparse matrix in the b2_matrix file, for the first nwrite calls.
If the iteration number is equal to ncallout, then several output files 'b2ne_npmo', 'b2ne_xppb', 'b2ne_npco', 'b2ne+nppo' which detail the convergence behaviour and residuals.
If diagno.gt.0, outputs details of the calculations of neutral transport coefficients when using the new_df0 model. See switch b2tqna_new_df0 for more details.
idout0 and idout1 specify the desired selection of output segments; idout0 for output before the calculation and idout1 for output at the conclusion of the calculation. The generic description of either of these 'idout' variables follows. idout specifies the desired selection of output segments for printed and graphical output. Associated with each segment is a four-character word, which is identified on the printed output. (The same word may be associated with several segments.) A particular segment will be produced only if its associated word has a match, as defined by the routine strmas, in idout.
This switch is to be provided within b2yt.dat to indicate to the b2yt converter the format of the original Eirene input.dat file to be modified. The accepted values are (case-insensitive): 'old' for input files from SOLPS4.0 and SOLPS5.0 runs using 'old' Eirene_96 'new' for input files from SOLPS4.0 and SOLPS5.0 runs using 'new' Eirene_99 'facelift' for input files from SOLPS5.1 runs 'juelich' for input files from Juelich Eirene versions (2008 and younger) 'iter' for input files from SOLPS4.2/4.3 and SOLPS-ITER runs
This switch is to be provided within b2yt.dat to indicate to the b2yt converter the format of the original fort.44 file to be modified. The accepted values are (case-insensitive): '4.3' for a fort.44 file produced from SOLPS4.3 runs '5.0' for a fort.44 file produced from SOLPS5.0 runs '5.1' for a fort.44 file produced from SOLPS5.1 runs '5.2' for a fort.44 file produced from SOLPS5.2 runs 'iter' for input files from SOLPS-ITER runs (no conversion necessary)
This switch is to be used when running the b2_ual_write utility program to obtain an IDS from the SOLPS run results. If set to 1, the IDS will indicate that the results come from a converted SOLPS4.3 run with no additional calculation. By default, the IDS will indicate a SOLPS-ITER run.
If AFN_out.ne.0, boundary flux components of the advanced fluid neutral models are written out. See code for details.
Multiplier to heat flux contribution to divergence of viscosity tensor in the momentum equation.
When not equal to '0.0', adds contribution to divergence of viscosity tensor coming from x-variations in B.
If style_qip.eq.1, adds a classical ion heat conductivity term to the heat flux used to compute the heat viscosity current (see manual for full details).
If b2sigp_style is set to '2', this switch has no effect. When set to '1', the new correct form of the friction force is used, applicable for non-hydrogenic plasmas or hydrogenic mixtures. The value '0' corresponds to the old SOLPS5.0 treatment. Not recommended unless wanting to recover older 5.0 results.
When set to '1', there is no radial contribution in Joule heating because there are no physical reasons to take it into account. The value '0' corresponds to the old SOLPS5.0 treatment.
Multiplier to the parallel momentum source term associated with the anomalous current. It is recommended '1.0'. Only active if both the ExB and diamagnetic drifts are turned on. The value '0.0' corresponds to the old SOLPS5.0 treatment.
Multiplier of the centrifugal force term. It is recommended '1.0'. The value '0.0' corresponds to the old SOLPS5.0 treatment.
Multiplier of the centrifugal force correction term due to linearization. It is recommended '1.0'. The value '0.0' corresponds to the old SOLPS5.0 treatment.
Real parameter which determines anomalous current. Both b2tfhe_anomalous and b2tanml_anomalous point to the same switch (backward compatibility with SOLPS5.2): b2tfhe_anomalous is checked first and only if it is not found is b2tanml_anomalous looked for. If b2tfhe_anomalous is 0 then anomalous current is switched off otherwise anomalous current is switched on, and the numerical value of b2tfhe_anomalous acts as a multiplier to the anomalous electrical conductivity and thermo-electric coefficients provided.
Real parameter which multiplies ExB flows. If b2news_ExB.eq.0 and b2news_facExB_start.eq.0 then ExB flows are switched off. If b2news_ExB is nonzero, then ExB flows are multiplied by that constant throughout the run. See also Run section on switches b2news_facExB_... for more details. A spatial fac_ExB profile is also possible, see Numerics section for details.
Real parameter which multiplies the viscous drift flows. If b2news_vis.eq.0 and b2news_facvis_start.eq.0 then the viscous drift flows are switched off. If b2news_vis is nonzero, then the viscous drift flows are multiplied by that constant throughout the run. See also Run section on switches b2news_facvis_... for more details. A spatial fac_vis profile is also possible, see Numerics section for details.
Real parameter which multiplies the ion inertial current.
Real parameter which multiplies the diamagnetic current.
Real parameter which multiplies the convective heat flux due to grad B-drift of guiding centers in non-modified heat fluxes of electrons and ions.
Real parameter which multiplies the ion-neutral current. If b2tfhe_neutral is 0 then the ion-neutral current is switched off otherwise the ion-neutral current is switched on. The ion-neutral current also requires either the diamagnetic or ExB drifts to be turned on as well.
Integer switch to turn off or on the ion-neutral current in the core region (applies to coupled runs only). This is recommended in cases where the neutral densities are very low in the core region and the ion-neutral current is likely to vary widely from one iteration to the next as a result of Monte Carlo noise. If fchin_in_core.eq.0 (default), then fchin is set to zero in the core. If fchin_in_core.eq.1, then fchin is unchanged. fchin_in_core.eq.1 not yet available for WG.
Real parameter which multiplies the Pfirsch-Schlueter electron heat flux and conductivity.
Real parameter which multiplies the current driven by parallel viscosity. If b2tfhe_vis_par is 0 then the viscosity-driven current is switched off otherwise the viscosity-driven current is switched on.
Real parameter which multiplies the current driven by heat viscosity effects.
Real parameter which turns on stochastic current. If b2tfhe_stochastic is 0 then stochastic current is switched off otherwise stochastic current is switched on. Not yet available for WG.
Width of the stochastic current layer (in metres), measured from the separatrix inward, along the poloidal index ixref (given by b2tqna_ixref). If b2tfhe_stochastic.ne.0, then b2tstch_delta must be greater than zero. Not yet available for WG.
Multiplier to the magnetic field line stochastic diffusion coefficient, describing the stochastic conductivity. Not yet available for WG.
Multiplier to the stochastic conductivity. Not yet available for WG.
If lluciani.ne.0, then transport coefficients on cells belonging to closed field lines are modified according to the Luciani model. If lluciani.eq.1, the standard connection length formulation is used. If lluciani.eq.2, the old SOLPS4.0 formulation is used for backward compatibility. If lluciani.eq.3, Spb's new form Luciani's coefficient.
If lthf21.ne.0, then the heat transfer model is modified according to the 21-moment approach. This should only be used for multi-fluid runs. Not yet available for WG.
If lvis21.ne.0, then the viscosities are modified according to the 21-moment approach. This should only be used for multi-fluid runs. This switch should not be used along with 'b2tfhe_vis_par' to avoid double-counting of classical viscosity effects. Not yet available for WG.
If art_rad.ne.0, then an artificial radiation loss term is added to the electron cooling rate. Art_rad represents the percent fraction of impurities in the plasma.
The 20 switches above are all feedback switches and require that b2stbc_feedback.ne.0 (See Run section for details) or that lfeedback=.true. in b2.boundary.parameters. Feedback schemes are now all applied in the same way in a dedicated module b2us_feedback (see FEEDBACK_CONTROL namelist). Switches in b2mn.dat are still available and will override whatever is specified in b2.feedback_control.parameters, but it is recommended to not use switches anymore. This will allow more flexibility in applying such schemes. The feedback is done as a boundary condition and is under-relaxed using the b2stbc_...._alpha switches (See Numerics section for details). fheycore is the radial electron heat flow entering the core boundary. Will apply fb_type=11, fb_rescale=1, fb_actuator=4 on region defined by coreregno. fhiycore is the radial ion heat flow entering the core boundary. Will apply fb_type=12, fb_rescale=1, fb_actuator=5 on region defined by coreregno. fchycore is the radial current entering the core boundary. Will apply fb_type=13, fb_rescale=1, fb_actuator=4 on region defined by coreregno. fnaycore is the radial flux of species "isfeedback" entering the core boundary. Will apply fb_type=10, fb_rescale=1, fb_actuator=3 on region defined by coreregno. If fhiycore_kinetic_energy.eq.1, the parallel kinetic energy flux is counted as part of the ion energy flux. For double-null cases, the f..ycore variables contain the total flow through both core boundaries taken together (to be checked for WG code!). The core boundaries are identified by the coreregno and coreregn2 switches (see Geometry section). These switches need to be manually adjusted for WG cases as now they indicate the boundary where the feedback is applied. For switches like b2stbc_xxxycore or private_gas\_puff the user has to specify in b2.feedback_control.parameters} where the feedback is calculated (see b2us.regions.parameters produced at conversion), while coreregno still indicates where it is applied. pfrregno and solregno lose meaning in kinetic cases if the actuator is the gas puff. Coreregn2 and pfrregno2 are not yet converted/used. The other seven switches are density feedback switches: nesepm is the outer midplane separatrix electron density. (See rzomp switch in b2.user.parameters for setting of midplane position). This is controlled through particle input of main plasma species isfeedback through the core boundary (See Geometry section for the core boundary indices coreregno and coreregn2). nesepm_sol is the outer midplane separatrix electron density, but this time controlled through a gas puff outside the SOL region. (See b2mwti_jxa switch in Geometry section for setting of midplane position). This is done through particle input of neutrals species isfeedback through the SOL North boundary (See Geometry section for the SOL North boundary index solregno). If nesepm_overshoot.gt.1.0, the gas puff is turned off when the midplane separatrix density reaches nesepm*nesepm_overshoot. Minimum and maximum values of the feedback gas puff can be set using b2stbc_nesepm_minpuff and b2stbc_nesepm_maxpuff. It will apply fb_type=3, fb_rescale=1, fb_actuator=1 on region defined by solregno. nepedm_sol is similar to nesepm_sol but the feedback is exerted on the electron density at the top of the pedestal, along the outer midplane, at the icped-th ring (where icped=iyped+2) counting from the inner core boundary. It will apply fb_type=16, fb_rescale=1, fb_actuator=1 on region defined by solregno. ndes_sol is the total particle content from the isonuclear sequence of species 'isfeedback' over the entire simulation domain. It is governed, like nesepm_sol, by nesepm_overshoot and nesepm_alpha and corresponds to a SOL boundary feedback. Will apply fb_type=14, fb_rescale=1, fb_actuator=1 on region defined by solregno. WARNING: in the structured code this switch was asking for BCCON=12, so density BC for fluid neutrals feedback but then it switches to gas puff for kinetic neutrals. Now it is ONLY gaspuff. volrec_sol is the total particle source due to volume recombination coming from Eirene. Its feedback is governed by volrec_overshoot and volrec_alpha and volrec_beta parameters (See Numerics section for the latter two). Will apply fb_type=15, fb_rescale=1, fb_actuator=1 on region defined by pfrregno1. WARNING: not available yet for fluid neutrals. nesepm_pfr is the outer midplane separatrix electron density, but this time controlled through a gas puff outside the PFR region. (See b2mwti_jxa switch in Geometry section for setting of midplane position). This is done through particle input of neutrals species isfeedback through the private flux boundaries (See Geometry section for the PFR South boundary indices pfrregno1 and pfrregno2). will apply fb_type=3, fb_rescale=1, fb_actuator=1 on region defined pfrregno1 ndes is the total particle (neutrals + ions) content of the species whose neutral species has index isfeedback. This is controlled through particle input of neutrals species isfeedback through the private flux boundaries (See Geometry section for the private flux region boundary indices pfrregno1 and pfrregno2). will apply fb_type=14, fb_rescale=1, fb_actuator=1 on region defined by pfrregno1. WARNING: in the structured code this switch was asking for BCCON=12, so density BC for fluid neutrals feedback but then it switches to gas puff for kinetic neutrals. Now it is ONLY gas puff. private_flux_puff NOT YET CONVERTED TO WG CODE! is the particle flow of species isfeedback through the private flux region boundary. This is controlled through particle input of neutrals species isfeedback through the private flux boundaries (See Geometry section for the private flux region boundary indices pfrregno1 and pfrregno2). All but the first of these feedback schemes ("nesepm") can also be run though an Eirene neutral gas puff. In that case, one should set an appropriate gas puffing stratum in the Eirene input file input.dat, as well as in the b2.neutrals.parameters file. The code will assume, unless told otherwise through the use of "eirene_nesepm_istra", that the first stratum of type 'C' found in b2.neutrals.parameters is the stratum used for the feedback. If "eirene_nesepm_istra" is set to a positive number, the code will consider the feedback stratum to be the stratum in the "eirene_nesepm_istra"-th position in b2.neutrals.parameters. b2stbc_iyped is the index in the OMP CV list (defined with RZOMP in b2.user.parameters) of the pedestal cell (defaults to icsepomp/2).
These switches are associated with the boundary condition of type 13, 16, 17, 20 or 21 for particle flux, or electron and ion energy fluxes, i.e. constant density or temperature giving a prescribed flux through a radial surface, located type13_ref or type16_ref radial steps away from the boundary of the computational domain, respectively, and applied type20_ref steps away. See code and description in b2cdcn for details. Type 21 also applies to the electric potential boundary condition. It scales the feedback strength according to the value of the variable on the ring type21_ref steps inside. When type16_ref.gt.1, then the position where the core flux tallies are computed moves along with the boundary condition. If type16_kinetic_energy.eq.1, then the parallel kinetic energy flux is included in the ion heat flux component of type 16 and 17 boundary conditions.
If bcene_15_style.eq.0, the electron sheath boundary condition computes the electron current as the difference between the ion flow and the poloidal current. If bcene_15_style.eq.1 (recommended, default), the electron current at the sheath is computed directly from the sheath potential drop. This method has proven numerically more stable, especially in cases with drifts.
If bceni_15_style.eq.0 (default), use the ion sheath transmission factor specified by ENIPAR(IB,1). If bceni_15_style.eq.1, use the self-consistent ion sheath transmission factor implied by the truncated drifting Maxwellian ion distribution determined by the Eirene sheath model (ENIPAR(IB,1) is not used). See Eirene manual Section 1.5.1 for the detailed expression.
If bcene_16_style.eq.1, the boundary condition is enforced using the modified fluxes fhe_mdf. Recommended when running cases with drifts. bcene_16_style.eq.1 not yet available for WG.
If bceni_16_style.eq.1, the boundary condition is enforced using the modified fluxes fhi_mdf. Recommended when running cases with drifts. bceni_16_style.eq.1 not yet available for WG.
If secmodel.eq.1, then the secondary electron emission coefficient at the plates is computed locally according to a kinetic model, otherwise the default values of cbsch(7,ireg) and/or gammae are used. Not yet available for WG.
If boundary_sources.ne.0, allow for the specification of boundary values for all equations, using BC type 7, provided in the bv_na.dat, bv_ua.dat, bv_te.dat, bv_ti.dat, and bv_po.dat files. All files must be present.
If equation_sources.ne.0, allow for the specification of source term values for all equations, per unit volume, provided in the art_sna.dat, art_smo.dat, art_she.dat, art_shi.dat, and art_sch.dat files. All files must be present.
If use_mms.ne.0, specifies a run using the Method of Manufactured Solutions. The exact solutions to match are provided by the b2mndr_set_XX_numerical switches. An exact solution is available for comparison if the corresponding set_XX_numerical switch is nonzero. The exact solutions are provided, respectively, in the exact_na.dat and exact_ne.dat, exact_ua.dat, exact_te.dat, exact_ti.dat, and exact_po.dat files. The exact_na.dat and exact_ne.dat files must be provided simultaneously.
Sputtering model switches. Not yet available for WG. See the code in b2stbr and the b2mod_sputter module for specific implementation details and references. The default values represent the case of Graphite plates. Sput_src is the atomic number of the plasma species which causes chemical sputtering. Sput_dst (when .ge.0) is the species index of the destination species produced by physical sputtering, chemical sputtering and/or RES. If either of the latter two processes are included, the chemical sputtering and RES production rate calculations provided in the code assume that sput_dst points to a Carbon species. Sput_dst2 and sput_dst3 (when .ge.0) represent other species produced by wall interactions for mixed materials scenarios. Sput_frac_flag is the switch to turn on mixed materials scenarios (when sput_frac_flag.eq.1). Plate_model.eq.0 means the 0-D time-independent plate heating model while plate_model.eq.1 indicates the 1-D time-dependent plate heating treatment. Plate_model.eq.2 gives acces to a 2-D time-dependent model. Plate_option chooses the initialisation of the plate temperature profile. If plate_option.eq.1, the profile is set to the constant given in plate_temp. If plate_option.eq.2, the profile is computed to be the 0-D equilibrium profile. If plate_option.eq.3, the profile is read from results of the previous run. Sput_phys_model is a switch for choosing between the TRIM tables (model 1, default) or an empirical formula (model 0). When TRIM data is not available, the empirical formula is automatically used. Extrapolations of low and high energy ranges beyond the TRIM table data is done using the same physical dependencies as the empirical formula. Sput_chem_model is a switch for the chemical sputtering model used. If sput_chem_model.eq.0 (default), the empirical formula is used. If sput_chem_model.eq.1, a constant with a low energy cut-off is used. The cutoff occurs at approximately sput_chem_cutoff_alpha and the width is determined by sput_chem_cutoff_beta (the larger the value, the narrower the width over which the transition from 0 to 1 occurs). For model 0, sput_frc is a multiplier to the empirical formula, while for model 1, sput_frc is the constant chemical sputtering yield. Sput_frc is superseded by the chem_sput array from b2.neutrals.namelist if the latter is used. For neutrals species, we add a factor of α*na*vbar to the particle flux to the plate (used to compute chemical and RES sputtering), where vbar is the average neutral particle speed. Sput_phys turns on physical sputtering when sput_phys.gt.0.0 and contains a (real) multiplier for the physical sputtering rate looked up in the TRIM tables. It is superseded by the phys_sput array in b2.neutrals.parameters if the latter is used. sput_phys_col indicates which column to use in the TRIM table; the default (3) corresponds to 30 degrees incidence. The columns are for 0 15 30 45 55 65 75 80 85 degrees. This angle is also used in the empirical formula. Plate_temp is the temperature (in Kelvin) of the COOLED end of the divertor plates. The surface temperature is computed self-consistently (using a 1-D description) according to the incident heat fluxes from the plasma. The plate is assumed to have a thickness of plate_thick, measured in metres. When plate_thick.eq.0 (default), the plate surface temperature is assumed to be the same as plate_temp. Sput_res turns on radiation enhanced sublimation (RES) when .gt.0.0 and is a multiplier to the RES rate. The switch sputter_energy_on turns on the energy contribution from the sputtered particles (physical and chemical). If sputter_energy_on.eq.0, the particles are sputtered cold, with no energy contribution to the ion heat equation. If sputter_energy_on.ne.0 (default), the chemically sputtered particles are re-introduced into the plasma with the thermal energy corresponding to the surface temperature of the material, while physically sputtered particles are re-introduced with the TRIM-calculated energy sputtering yield (if the TRIM model is used). Therm_evap turns on thermal evaporation when .gt.0.0 and is a multiplier to the thermal evaporation rate. Redep_alpha is the multiplier to the reduction of the sputtering yield due to prompt redeposition. The promptly redeposited species is counted as eroded-then-deposited species.
Reflection model switches. If the recycling coefficients (particle and energy) are set to zero, then the code attempts to use, for surfaces in contact with walls, the TRIM reflection data for particle and energy reflection coefficients, when refl_model.eq.1. If refl_model.eq.0, an empirical formula is used instead. The model switch is subservient to reflection_on.ne.0. If reflection_on.eq.0, then no reflection is used. If reflection_on.eq.1 (default), then the incident species is reflected unchanged. Otherwise, the reflected flux is of the associated neutral species to the incident particles. Not converted for WG. It is recommended to use the advanced fluid neutral (AFN) models instead.
Turns on the coronal model if coronal.ne.0. I.e., the reaction rates of all reactions are made density-independent, and the value computed at the lowest density (given in b2ar.dat) is used.
hz has been introduced into the new form of the parallel momentum balance equation. If fac_hz = 0.0 then hz = 1 and old form of equations is used. If fac_hz = 1.0 then new form of equations is used.
When set to '1.0', the old form of the electron heat flux calculation is used and b2tfhe_fch_pTe should be set to '0.0'. It is recommended to use 0.0. Cannot be set to 1.0 when using the SOLPS5.2 physics model.
When set to '1.0', the new form of the electron heat flux calculation is used and b2tfhe_alfTeEh should be set to '0.0'. It is recommended to use 1.0. Cannot be set to 0.0 when using the SOLPS5.2 physics model.
Xcur is a multiplier to the parallel viscosity, perpendicular viscosity, heat viscosity, ion inertial and anomalous current contributions to the ion poloidal flows (particle and energy). If set to 0 (default), these currents are carried by electrons. If set to 1, these poloidal currents are carried by the main hydrogenic ions. The latter option can only be used if the plasma has a single majority hydrogenic species. Otherwise, an error will be returned.
Ycur is a multiplier to the parallel viscosity, perpendicular viscosity, heat viscosity, ion inertial and anomalous current contributions to the ion radial flows (particle and energy). If set to 0, these currents are carried by electrons. If set to 1 (default), these radial currents are carried by the main hydrogenic ions. The latter option can only be used if the plasma has a single majority hydrogenic species. Otherwise, an error will be returned.
vis_per is a multiplier to the perpendicular viscosity current contributions to the ion poloidal flows (particle and energy). Subservient to b2tfnb_xcur and b2tfnb_ycur. Cannot be used in conjunction with the 9-point stencil numerical treatment. If used with the 5-point stencil treatment and with fluid neutrals, must be equal to b2tfhe_vis_per for numerical stability reasons.
vis_q is a multiplier to the heat viscosity current contributions to the ion poloidal flows (particle and energy). Subservient to b2tfnb_xcur and b2tfnb_ycur.
When set to '1', the Zhdanov expression is used in the electron thermal conductivity. It is recommended to use 'b2tqce_fke_Zhdanov' '1'. This switch is only active if, simultaneously, one has b2tqce_model.eq.1 and b2tfhe_fch_pTe.eq.1.0.
When set to '0', b2tqce calculates all classical contributions to the electron transport coefficients on all cells. When set to '1', the electrical conductivity (sig) and the thermo-electric coefficient (alf) are not calculated in the guard cells that lie on open flux surfaces.
The last switch is subservient to b2tqna_user_transport, and only used if user_transport.ne.0. When this is the case, the electron perpendicular heat diffusivity and anomalous perpendicular particle transport coefficient are scaled so as to maintain the required_te_gradient (in units of eV/m), at the iyref-th radial ring along the outer midplane. By default, this is the ring immediately inside the separatrix. The transport_eta switches are activated when user_transport.eq.3 and allow modifications of the electron heat conductivity according to a constant eta (Grad Te / Grad Ne) model. See routine set_transport_eta code for details. A model for disruption transport coefficients is available with user_transport.eq.7. See routine set_transport_disruption code for details. Not yet available for WG. The value user_transport.eq.8 activates the neoclassical model for computation of the transport coefficients in the core. See routine set_transport_neo and subroutines called for details. Implemented from the NEOART package.
Upper and lower limits of the diffusivity coefficients df0 computed in b2tqna for each neutral species.
When new_df0.eq.1, the neutral diffusivity is computed according to the local charge exchange, ionisation and elastic collision rates, instead of the standard form with constant cross-sections.
Ballooning model switches. If ballooning.ne.0, then all transport coefficients in cell (i) are rescaled by a factor of ballooning_rescale*abs(bb_ref/bb(i))**ballooning . The default value for bb_ref is the arithmetic average of the total magnetic field strength over the entire computational domain. If b2tqna_ballooning_sig.ne.0 (default), the ballooning factor is also applied to the anomalous radial electrical conductivity and the anomalous radial thermo-electric coefficient.
Scaling factor for all ion and electron transport coefficients inside private flux regions.
Scaling factor for all ion and electron transport coefficients inside divertor SOL regions.
Multiplier of the friction term between charged species (only ions in case b2sigp_style.eq.2).
If b2sigp_style.eq.2, multiplier to the friction force term between electrons and ions. Otherwise, multiplier of the ehxp term in the older expression for the thermal force.
Multiplier of the electron thermal gradient term in the thermal force term and parallel current
Multiplier of the ion thermal gradient term in the thermal force term.
Parameters for the computation of the thermal force term.
Parameters for the computation of the thermal force term.
If lambda is positive, the Coulomb logarithm is set to lambda. If lambda is negative, then abs(lambda) is the lower bound to be used when computing the Coulomb logarithm according to the NRL formulae. The computation of the Coulomb logarithm takes place in b2tlnl. The default formula is to use the Braginskii formulation, which will give the same answer for electron-electron (ee), electron-ion (ei), and ion-ion (ii) collisions. Respectively, non-default values of '1', '2' and '3' for b2tlnl_ii, _ee, and _ei, respectively, make use of the calculation according to Wesson.
Multiplier to the Pfirsch-Schlueter flows.
The poloidal convective heat flux contains a prefactor of (3/2+BoRiS). The default gives us the internal energy equation. The value BoRiS = 1.0 gives us the total energy equation. Only for use to benchmark with codes using the total energy equation. When used, the meaning of fhe and fhi changes from internal energy fluxes to total energy fluxes.
When conduction_only.eq.1, convective heat transfer terms are turned off. Only to be used when benchmarking against pure conduction cases.
Parameters for the flux limit to the convective neutral flow. Alpha is a multiplier to the classical flux limit value. The larger α is, the weaker the flux limit is. γ is the exponent used in the flux-limiting formula. The smaller γ is, the stronger the flux limit is. If alpha.eq.0, no flux limit is applied. flux_limit_min_ti specifies the minimum ti to be used (in eV).
Parameters for the flux limit to dpa0 - pressure driven neutral diffusion. Alpha is a multiplier to the classical flux limit value. γ is the exponent used in the flux-limiting formula. If alpha.eq.0, no flux limit is applied.
Parameters for the flux limit to the heat conductivity of the neutrals. Alpha is a multiplier to the classical flux limit value. The larger α is, the weaker the flux limit is. γ is the exponent used in the flux-limiting formula. The smaller γ is, the stronger the flux limit is. If alpha.eq.0, no flux limit is applied. flux_limit_min_ti specifies the minimum ti to be used (in eV).
Parameters for the flux limit to the viscosity of the neutrals. Alpha is a multiplier to the classical flux limit value. The larger α is, the weaker the flux limit is. γ is the exponent used in the flux-limiting formula. The smaller γ is, the stronger the flux limit is. If alpha.eq.0, no flux limit is applied.
If style = 0 then the original form of the viscosity flux limit is used, otherwise the SPb flux limit is used.
EXPERIMENTAL SWITCH GROUP, USE WITH CARE. If set to 1, enables gradual modification of the electron heat flux limiter towards low density regions. Above b2tlhe_ne_max (default 1.0e18 m^{-3}), the standard value cflme is used. Below b2tlhe_ne_min (default 1.0e17 m^{-3}), the value b2tlhe_cflme_min is used (default 0.2). For densities in between b2tlhe_ne_min and b2tlhe_ne_max, the heat flux limiter is computed based on linear interpolation between the limiting values.
EXPERIMENTAL SWITCH GROUP, USE WITH CARE. If set to 1, enables gradual modification of the ion heat flux limiter towards low density regions. Above b2tlhi_ni_max (default 1.0e18 m^{-3}), the standard value cflmi is used. Below b2tlhi_ni_min (default 1.0e17 m^{-3}), the value b2tlhi_cflmi_min is used (default 1.0e1). For densities in between b2tlhi_ni_min and b2tlhi_ni_max, the heat flux limiter is computed based on linear interpolation between the limiting values.
EXPERIMENTAL SWITCH GROUP, USE WITH CARE. If set to 1, enables gradual modification of the parallel viscosity flux limiter towards low density regions. Above b2tlmv_ni_max (default 1.0e18 m^{-3}), the standard value cflmv is used. Below b2tlmv_ni_min (default 1.0e17 m^{-3}), the value b2tlmv_cflmv_min is used (default 0.5). For densities in between b2tlmv_ni_min and b2tlmv_ni_max, the flux limiter is computed based on linear interpolation between the limiting values.
Multiplier of the contribution to electron heat sources from divergence(ue,ve). This term is superseded by the BoRiS switch if invoked.
Multiplier of the contribution to ion heat sources from divergence(ua,va). This term is superseded by the BoRiS switch if invoked.
Multiplier of the contribution to ion heat sources from viscous heating due to poloidal velocity differences. This term is superseded by the BoRiS switch if invoked.
Multiplier of the contribution to electron heat sources from Joule heating (electron-ion friction).
Multiplier of the contribution to ion heat sources from atom-atom friction. This term is superseded by the BoRiS switch if invoked.
Multiplier of the contribution to electron heat sources from divergence of the electrical drift. That contribution is also multiplied by facExB.
Multiplier of the contribution to ion heat sources from divergence of the electrical drift. That contribution is also multiplied by facExB.
Multiplier of the contribution to heat sources from friction due to diamagnetic velocities. Normally already included in 'phm3' term above.
Multiplier of the contribution to heat sources from viscous heating due to radial velocity differences. This term is only included when using b2nph9_style.eq.1 or b2npht_style.eq.1 (defaults). This term is superseded by the BoRiS switch if invoked.
Multiplier of the contribution to electron and ion heat sources from divergence of the gradTxB heat flux. That contribution is also multiplied by facdrift. This term is no longer needed as it has been transferred from being a source to being included in the flows.
When set to '0' or '2', the code calls the standard b2stbc routine, which uses the particle flux with drift terms included in the expression of the electron particle flux (fne). Option '1' is obsolete. It is recommended '2'.
phm0 : Multiplier to the charge exchange rate coefficient. Subservient to styl0: only used if styl0.eq.0 or is.lt.2. If styl0.ne.0, the charge exchange rate coefficients for species of index 2 and above are set to zero.
Multiplier to the ionisation rate coefficient.
Multiplier to the recombination rate coefficient.
Multiplier to the heat loss rate coefficient.
Multiplier of the sources due to atomic physics.
If '0', flux limit is not applied directly to the electron heat flux but is applied through transport coefficients if 'b2trcl_conductive_limit' is '1'. It is recommended '0'. b2tfhe_lim_flux.eq.1 not yet available for WG.
If '0', flux limit is not applied directly to the ion heat flux but is applied through transport coefficients if 'b2trcl_conductive_limit' is '1'. It is recommended '0'. b2tfhi_lim_flux.eq.1 not yet available for WG.
Multiplier to the temperature equipartition term.
Multiplier for the classical parallel viscosity.
If model.eq.1, use the Balescu formulation from SOLPS5.2 classical parallel ion heat diffusivity. If model.eq.2, use the older Braginskii SOLPS4.0 model. Note: old option model.eq.3 removed, replaced with model.eq.1., but numerical treatment w.r.t. factor 4/3 according to old model.eq.3.
If model.eq.1, use the fitted Balescu formulation from SOLPS5.0 classical parallel electron heat diffusivity. If model.eq.2, use the older Braginskii SOLPS4.0 model. If model.eq.3, use the 21-moment Balescu results.
If '1', use constant density in core region to calculate anomalous conductivity sig0=dfsig*qe*ne(nmdpl,-1), nmdpl - number of midplane cell. If '0', use sig0=dfsig*qe*ne(x,y)
If '1', the anomalous contributions in the parallel direction to the electrical conductivity csig and the thermo-electric coefficient calf are zeroed out.
Integer switch to choose the expression for the ion-neutral conductivity. If '1', the ion-neutral conductivity is computed from hydrogen gas diffusivity (SOLPS5.2 treatment). If '0' (default), it is computed from the local CX rate (SOLPS5.0/5.1 and AFN treatments).
If pol_anom_scale.ne.1, the poloidal anomalous transport coefficients are scaled by the number. To get the same treatment as setting the old switch (radial_only) to 1, set the new switch to 0.0. This multiplication is to only take place for charged species.
Deprecated. Use "DEFINE_LINES" in block 12 of Eirene input file instead.
If > 0 then repeats the first call to Eirene in eirene_mc so many times. Useful if a density feedback gas puff is being used because the flux passed to Eirene on its first call is ignored (instead the value in the Eirene input file is used).
If > 0 use recyceir (non species dependent) to specify the recycling coefficients, else if 0 use recyc (species dependent).
If <> 0 then recycles the neutral flux having crossed the core boundary within Eirene as ions. The recycling is modulated as per the poloidal density distribution of the ions, and neutrals come back as fully-stripped ions. 'eirene_ionizing_core' is an alias for this switch. If the value = 1, then the flux is added by direct modification of the sources in the guard cells --- this will only work if a standard flux boundary condition is applied at that boundary. If the value is < 0, then the absolute value specifies which boundary in b2.boundary.parameters is to be used. This will only work for the BCCON = 13 type boundary condition, for which a negative value for ionizing_core is required. Cannot be used in conjunction with 'eirene_ank_mods'. Use of this switch is not allowed if the core boundary condition is not of a flux type. Moreover, certain boundary conditions expressly include the flux of neutrals already (BCCON = 6, 19, 21, 23, 26, or 27), in which case use of this switch is redundant and should be avoided.
If eirene_background.eq.0, the ion velocities passed to Eirene to be used for the collisions are based on grad-B and ExB drifts (vadia + vaecrb). If eirene_background.eq.1, these velocities contain the full diamagnetic and ExB drifts (wadia + vaecrb). Note: recycling fluxes are always computed based on grad-B and ExB drifts only (and are not affected by this switch), because diamagnetic drift flows largely close within the sheath.
If eirene_sheath_pot.eq.1, the sheath potential drop as computed by B2.5 (i.e. including effects of parallel currents, secondary electron emission, etc.) is passed to EIRENE to compute ion acceleration in the sheath. If eirene_sheath_pot.eq.0, the sheath potential drop is recomputed by EIRENE, usually assuming zero current and secondary electron emission.
If atomic_physics_rescale is not 0, then the code will rescale the rates stored in b2frates according to the multipliers given in the b2.atomic_physics_rescale.parameters inputfile before making use of them.
..set the contribution ic in NEOART 0 --- classical particle flux 1 --- banana plateau contribution 2 --- Pfirsch-Schlueter contribution 3 --- both banana and PS 4 --- all contributions mind that B2.5 already calculates the classical transport ! avoid double transport, 0+4 for cross-checks only !
Main switch controlling the use of the k(-eps)-model. In this model, anomalous transport coefficients are determined self-consistently through the solution of an additional transport equation for the turbulent kinetic energy k (variable "kt" in the code) and possibly the turbulent enstrophy zeta (variable "zt" in the code). Note that both the k and k-eps models are currently under active development/testing, and should be used with care. In particular, the k-eps model (option '2') is not suitable for routine testing yet. If '0': k(-eps)-model not used. If '1': one additional transport equation solved for k. If '2': two additional transport equations solved for k and zeta. Warning: under development - expert use only! Note: solutions to these equations are only used to compute the transport coefficients if the switch b2tqna_transport_keps is set.
Compute transport coefficients based on k(-eps)-model. If '0': transport coefficients not computed based on k(-eps)-model. If '1': transport coefficients computed using k-model (i.e. based on kt). If '2': transport coefficients computed using k-eps model (i.e. based on kt and zt). Warning: under development - expert use only!
If '0': use Larmor radius at OMP separatrix in k(-eps)-model. If '1': use local Larmor radius for each cell in k(-eps)-model.
If '1': print output for k(-eps)-model.
Multiplier for the diffusion coefficient in the k(-eps) model. Recommended range: 0.1 ... 1.0
Multiplier for the shear contribution to the diffusion coefficient in the k(-eps) model. Only affects results if ExB drifts are present.
Multiplier for the electron heat conductivity (w.r.t. diffusion coefficient).
Multiplier for the ion heat conductivity (w.r.t. diffusion coefficient).
Multiplier for the ion viscosity (w.r.t. diffusion coefficient).
Multiplier for the anomalous electrical conductivity (w.r.t. diffusion coefficient x qe x ne(omp)).
Multiplier for the anomalous thermo-electric coefficient (w.r.t. diffusion coefficient x sqrt(qe/te) x ne(omp)).
Multiplier for the k-conductivity (w.r.t. diffusion coefficient).
Multiplier for the diffusivity in z-equation (w.r.t. diffusion coefficient).
fflokt and fconkt are fudge factors to approximate enhanced parallel transport of k due to term related to (fluctuations in parallel current) * (fluctuations in potential) In principle this models the same effect as skt_diss (thus: needs b2sikt_fac_sheath_* = 0.0).
fflozt and fconzt are fudge factors to the transport of zt. See code for details.
Multiplier to the perpendicular current due to Reynolds-stress in k-model.
Mechanism to enable a smooth transition between standard anomalous transport coefficients and k(-epsilon)-model coefficients. The final coefficients (dna0, vsa0, hcib, hce0, sig0, alf0) are linearly interpolated between the standard and the k(-epsilon) values based on keps_fac. keps_fac starts out as b2tqna_keps_init and the value is multiplied by b2tqna_keps_inc at each iteration.
Multiplier for the parallel transport term of kt.
Model to compute the sources in the energy equations due to the k-eps model. Currently only a single model is available.
Currently only a single computation style for the sources in the energy equations due to the k-eps model is available.
If '0': use sound speed, connection length, and Larmor radius at OMP in k(-eps)-model. If '1': use local sound speed, connection length, and Larmor radius.
Numerical stabilization of sources. See code for details.
Multiplier for the dissipation term of kt in the SOL region.
Multiplier for the dissipation term of kt in the core region.
Multiplier for the sheath loss term of kt in the SOL region.
Multiplier for the sheath loss term of kt in the core region.
If '0': use connection length in expression for dissipation term of kt. If '1': use b2sikt_fac_diss_lpar in expression for dissipation term of kt.
If '1': sets minimum source for k to zero.
If '1': add anisothermal contribution to transport coefficients from k-(epsilon)-model.
Multiplier for Reynolds-stress turbulent viscosity in k-model equations. Only active if ExB drifts are present.
Multiplier for transport of k associated with turbulent ExB heat fluxes.
Main switch to turn on the AFN model for hydrogenic neutrals. If this switch is turned on, some default values are changed (see description of the corresponding switches).
Switch to turn on the spatially hybrid fluid-kinetic neutral model. For multispecies simulations, the hybrid approach can only be used for the hydrogenic neutrals. For other species, you can use either a fully kinetic or purely fluid approach.
When set to 1, the AFN recycling boundary conditions are used for hydrogenic neutrals. When b2mn_afn.ne.0 and there is no specification of b2stbr_recycle_afn, the switch is automatically set to 1.
When set to 1, the coarse representation of the integrated TRIM reflection coefficients is used for the AFN boundary conditions, which is expected to be sufficient. In case you want to use the original fine representation, use b2stbr_afn_bcs_use_coarse = 0. However, the fine representation is only present for D on Be, C and W.
When set to 1, the AFN transport coefficients are used for hydrogenic neutrals. When b2mn_afn.ne.0 and there is no specification of b2stbr_transport_afn, the switch is automatically set to 1.
When set to 1, a collision frequency related to neutral-neutral collisions is added to the AFN transport coefficients. This is an ad-hoc model for fluid neutral transport in regions where plasma density is extremely low and CX effectively absent (e.g. below dome, far-SOL,...). When set to 0, the correction is not added.
When set to 1, perpendicular transport in the AFN model due to neutral-neutral collisions is modelled as density diffusion rather than pressure diffusion. When set to 0, it is modelled as pressure diffusion.
Advanced hybrid option. Forbidden to use for multispecies simulations (ns > 1) and only in use for the strata with HYB_TYPE='A'. The local charge-exchange Knudsen number Kn is calculated at the boundary. If (Kn < auto_spatial_hyb_Kn_1) then the recycled atom is treated as a fluid. If (Kn > auto_spatial_hyb_Kn_2) then the recycled atom is treated kinetically. (otherwise) linear combination between fluid and kinetic.
Macroscopic length scale (in m) used to calculate the local Knudsen number.
Should be used in combination with MAXW in b2.neutrals.parameters (see description MAXW).
Should be used in combination with MAXW in b2.neutrals.parameters (see description MAXW).
Should be used in combination with use_auto_spatial_hyb (see description use_auto_spatial_hyb).
Should be used in combination with use_auto_spatial_hyb (see description use_auto_spatial_hyb).
When set to 1, the energy required for the Franck-Condon dissociation of the thermally released particles at the surfaces in the AFN boundary conditions is subtracted from the electron energy equation. It is not recommended to turn this switch on, because it seems to lower the electron temperature to unrealistic values.
Potential drop used for BCENI/E = 15 when the potential equation is not solved (pot_eq.ne.1), in units of the local electron temperature.
Determines the method for applying flux limiters to fluid neutrals. style.eq.0: old treatment where the radial and poloidal conductive fluxes are limited separately. style.eq.1: improved treatment (isotropic flux limiters) where the total conductive flux is limited, preventing artificial rotation of the flow.
Minimum transport coefficients for use in k-epsilon model. It is also possible to use these minima for other transport models by setting b2tqna_limit_coeff.eq.1.
b2tqna_limit_coeff.ne.0: the minimum transport coefficients dna_min, vsa_min, hci_min and hce_min are for whatever transport model is selected. b2tqna_limit_coeff.eq.0: the minimum transport coefficients dna_min, vsa_min, hci_min and hce_min are only used for the k-epsilon model.
Constant c_q1 in calculation of additional viscosity terms calculated in b2siav.
Multiplier for the pressure gradient term. Setting b2sigp_phm0.ne.1.0 should only be done for specific testing purposes.
Multiplier for the convective term due to new form of ion viscosity term in the momentum correction equation. Applies to hydrogenic ions and fluid neutrals.
Multiplier for the convective term due to new form of ion viscosity term in the momentum correction equation. Applies only to hydrogenic fluid neutrals.
Multiplier for radial diamagnetic current in BCPOT=12.
Multiplier for radial inertial current in BCPOT=12.
internal_energy_sources.eq.0: the externally provided energy sources are total energy sources and a conversion to internal energy sources is performed. internal_energy_sources.eq.1: the externally provided energy sources are already internal energy sources.
Multiplier to the anomalous current.
The correct switch is now b2tfhe_anomalous. b2tanml_anomalous is kept for backward compatibility. b2tanml_anomalous will only be read if b2tfhe_anomalous.eq.1.0.
Anomalous velocity in parallel momentum balance equation.
Under-relaxation factors related to feedback schemes for boundary conditions. Not yet available for WG.
These switches serve as identification for the simulation. They can be inherited from the SOLPS-GUI: IDS path : Only available if using IMAS Access Layer 5. The directory path where the IMAS data entry will be saved. Supersedes the pulse (previously shot) and run numbers. May begin with "$HOME", "$IMASDIR", or "$SOLPSTOP", that will be parsed by the code. If it begins with a '/', is understood as an absolute path, otherwise will be prefixed with the $IMASDIR environment variable. Defaults to the last line in the 'shotnumber.history' file if present, and to $IMASDIR/${pulse_number}/${run_number} otherwise or if $IMASDIR is modified from its default value by the switches below. Pulse number : Pulse number identifying the run (supersedes b2mndr_shot_number). Defaults to the last number found in shotnumber.history, or 0 if the file is not found. Must be positive and non-zero if an IMAS data entry is requested. If using IMAS Access Layer 4, cannot exceed 214748. Run number : The number of the run. Must be positive or zero. If using IMAS Access Layer 4, limited to 5 digits (i.e. span from 0 to 99999). Database : IMAS IDS database name (supersedes b2mndr_device). Defaults to 'solps-iter' if $DEVICE if undefined. If $DEVICE is 'iter', the code internally changes it to 'ITER' to follow the IMAS convention. Can be used to rebuild the IMASDIR path. User : The user who ran the simulation. Can be used to rebuild the IMASDIR path.
If read_b2fstate.ne.0, b2aidr will read an existing b2fstate file and append to it a flat state description for additional new species being introduced in the run.
Number of timesteps desired in this run. If ntim.eq.0, the code simply reads the input file and creates the b2fplasma and b2fstate files corresponding to the present plasma state.
Timestep (in seconds).
One of the dynamic timestep control parameters. Only active if both delta_max and delta_min are nonzero. The timestep is increased by a factor dt_change_inc or decreased by a factor dt_change_dec, within the bounds [dt_min, dt_max], according to whether the largest instantaneous change in all equations is smaller than delta_min or larger than delta_max, respectively.
stim specifies the initial time --- default 0. If set to a positive or zero value, this overwrites the time value read from b2fstati. If set to a negative value, the run continues from the time read in b2fstati. In that case, the tracing data is appended to the existing files, otherwise the tracing files are overwritten.
etim specifies the end time. Only active if etim > stim.
Time-dependent mode and iterative mode switches. The basic code timestep proceeds as follows: ..test input arguments ..compute auxiliary quantities ..prepare source computation ..do i0=1,nstg0 ..compute log-log linearised rate coefficients ..do i1=1,nstg1 ..compute source linearisation ..do i2=1,nstg2 ..perform one inner iteration ..re-compute auxiliary quantities ..produce monitoring output ..enddo ..enddo ..enddo This can be completed by the the 'b2mndt_nstg_ares??' switches. See 'Numerics' section for details.
If no_solve.eq.1, then the code is run without actually solving any equations. All secondary and dependent data is computed. The nstg(0:2) array is overwritten to '1's. The simulation time will not be updated. The code will compute fluxes, sources, transport coefficients, etc... 'ntim' times but not update the basic plasma quantities. Additionally, if no_solve is set to a negative value, then only some equations are solved (simulation time and nstg arrays will behave normally). no_solve.eq.-1 will activate the parallel momentum equations only. no_solve.eq.-2 will activate the density equations only. no_solve.eq.-4 will activate the potential equation only. (Active only if poteq.eq.1, otherwise, the behaviour dictated by the poteq setting applies). no_solve.eq.-8 will activate the heat equations only. These can be combined. For example, no_solve.eq.-3 will activate the parallel momentum and particle conservation equations. If the no_solve switch requires an equation not to be solved, the settings in the corresponding solvexx arrays from b2.numerics.parameters are moot. The latter are reserved for fine-tuning numerical diagnostics. If wishing to toggle whether to solve or not to solve each equation separately, one should set, in the b2.numerics.parameters input file, the SOLVEMO, SOLVEMT, SOLVECO, SOLVEE, SOLVEEI, SOLVEET and SOLVEPO arrays, respectively, for each of the parallel momentum equations, the total momentum equation, each density equation, the electron and ion heat equations, and the total energy equation. When running in time-dependent mode, the total energy and total momentum equations are not solved. The table below shows which combination of equations will be solved depending on the value of switch 'b2news_no_solve' *---------------------------------------------------- * | value | co | mo | he | hi | po | *---------------------------------------------------- * | 0 | + | + | + | + | + | *---------------------------------------------------- * | 1 | - | - | - | - | - | *---------------------------------------------------- * | -1 | - | + | - | - | - | *---------------------------------------------------- * | -2 | + | - | - | - | - | *---------------------------------------------------- * | -3 | + | + | - | - | - | *---------------------------------------------------- * | -4 | - | - | - | - | + | *---------------------------------------------------- * | -5 | - | + | - | - | + | *---------------------------------------------------- * | -6 | + | - | - | - | + | *---------------------------------------------------- * | -7 | + | + | - | - | + | *---------------------------------------------------- * | -8 | - | - | + | + | - | *---------------------------------------------------- * | -9 | - | + | + | + | - | *---------------------------------------------------- * | -10 | + | - | + | + | - | *---------------------------------------------------- * | -11 | + | + | + | + | - | *---------------------------------------------------- * | -12 | - | - | + | + | + | *---------------------------------------------------- * | -13 | - | + | + | + | + | *---------------------------------------------------- * | -14 | + | - | + | + | + | *----------------------------------------------------
CPU limit. If cpu.eq.0.0, no CPU limit is enforced. Otherwise, the code stops smoothly after it has been running for at least 'cpu' CPU seconds.
Elapsed time limit. Is elapsed.eq.0.0, no time limit is enforced. Otherwise, the code stops smoothly after is has been running for at least 'elapsed' wall clock seconds.
If savecpu.gt.0.0, CPU time interval after which save files plasmastate.xxxx are written. Other additional options for writing of the checkpointing files are available by means of the b2mndt_ntim_save and/or b2mndr_plasmatim switches. These can be used, in conjunction with the b2co.exe instruction, to recover a crashed run.
ismain identifies the index of the main plasma species. It must hold that ismain is not a neutral species.
Ramping parameters for facdrift, which multiplies the diamagnetic terms as well as the ion inertia current. The code is started on the first time step with facdrift=facdrift_start. If facdrift_target.ne.facdrift_start, then, on each time step, facdrift is multiplied by facdrift_inc. If the code does not converge on the timestep, facdrift is decreased by facdrift_dec. A facdrift profile is also possible, see Numerics section for details. The ion-neutral friction current requires either facdrift or fac_ExB to be turned on as well.
Ramping parameters for facExB, which multiplies the ExB terms. Same treatment as above for facdrift. Can be superseded by b2news_ExB to set a constant value throughout the run.
Ramping parameters for facvis, which multiplies the drift viscosity terms. Same treatment as above for facdrift. Can be superseded by b2news_vis to set a constant value throughout the run.
Namelist file indicators. If xxx_namelist.eq.1, then the code looks for the file b2.xxx.parameters, in which it expects to find the xxx namelist. In the case of boundary_namelist, defaults to 1 for double-null and snowflake geometries, 0 otherwise.
Profiles file indicator. If b2sral_inputfile.eq.1, then the code looks for the file b2.sources.profile in which it expects to find the 'profile' namelist. This namelist will contain profiles of sources measured from the outer midplane separatrix location, which is set with the [RZ]OMP variables in b2.user.parameters.
Profiles file indicator. If b2tqna_inputfile.eq.1, then the code looks for the file b2.transport.inputfile, in which it expects to find the transport namelist. This namelist will contain profiles of transport parameters measured from the outer midplane separatrix location, which is set with the [RZ]OMP variables in b2.user.parameters.
Turns on coupling with the Eirene Monte Carlo neutral code if nonzero. To be used, the code must be compiled with the -DB25_EIRENE option.
Turns on coupling with the ASTRA core transport code if nonzero. To be used, the code must be compiled with the -DASTRA option.
Multiplier to the neutral sources (either computed by Eirene or the corresponding rate coefficients).
Multiplier to the neutral density on the first timestep.
Multiplier of all densities on the first timestep.
Multiplier to the totally ionised species sources at the core boundary.
Feedback on the total heavy particle density. If density_control.ne.0, the sum of all densities is kept constant.
If feedback.eq.1, turns on feedback mode for the boundary conditions. See Physics section for more details. Equivalent to LFEEDBACK=.true. in the b2.boundary.parameters namelist.
Timestep index after which the feedback in b2stbc is activated. NOT YET CONVERTED IN WG CODE
If first_flight.ne.0, turns on the first flight model. See Physics section for additional details.
New number of species desired when creating a new initial state file using b2yt.exe. Must be specified within b2yt.dat.
Radial index of separatrix in the converted grid. Must be specified within b2yt.dat. Only applies when converting linear geometries.
New resolution of target depth desired when creating a new initial state file using b2yt.exe. Must be specified within b2yt.dat.
Rescaling of neutral densities by rescale_neutrals when creating a new initial state file using b2yt.exe. Must be specified within b2yt.dat.
When set to 1, b2yt will attempt to interpolate a plasma solution and create new input files for a new grid which is not commensurate to the original grid, but still has the same topology. Must be specified within b2yt.dat.