Scaling factors for rtsa: ionisation rates of species (is).
Scaling factors for rtra: recombination rates of species (is).
Scaling factors for rtqa: electron cooling rates of species (is).
Scaling factors for rtcx: charge exchange rates of species (is).
Scaling factors for rtrd: line radiation rates of species (is).
Scaling factors for rtbr: bremsstrahlung radiation rates of species (is).
Number of boundary segments.
Specifies the nature of the boundary segment N = 'North' boundary S = 'South' boundary W = 'West' boundary E = 'East' boundary X = 'X' boundary used for specifying a fixed value on a row of cells Y = 'Y' boundary used for specifying a fixed value on a column of cells
Contains parameters helping to define the boundary conditions for the continuity equation of species (is). See description of BCCON below for details.
Contains parameters helping to define the boundary conditions for the parallel momentum equation of species (is). See description of BCMOM below for details.
Contains parameters helping to define the boundary conditions for the electron energy equation. See description of BCENE below for details.
Contains parameters helping to define the boundary conditions for the ion energy equation. See description of BCENI below for details.
Contains parameters helping to define the boundary conditions for the potential equation. See description of BCPOT below for details.
For North, South or X boundary conditions, it specifies the row index; for West, East and Y boundary conditions, it specifies the column index.
For North, South or X boundary conditions, it specifies the start column index; for West, East and Y boundary conditions, it specifies the start row index.
For North, South X boundary conditions, it specifies the end column index; for West, East and Y boundary conditions, it specifies the end row index.
Contains the size of the list of cells where a boundary condition is applied.
Contains the X-coordinate of the cells where boundaries conditions are applied.
Contains the Y-coordinate of the cells where boundaries conditions are applied.
Specifying the type of density boundary condition for each segment and species (fastest varying index is species); makes use of CONPAR to specify additional information, as indicated:
Specifying the type of parallel momentum or velocity boundary condition for each segment and species (fastest varying index is species); makes use of MOMPAR to specify additional information, as indicated
Specifying the type of electron energy or temperature boundary condition for each segment; makes use of ENEPAR to specify additional information, as indicated
Specifying the type of ion energy or temperature boundary condition for each segment; makes use of ENIPAR to specify additional information, as indicated
Specifying the type of electric potential or current boundary condition for each segment; makes use of POTPAR to specify additional information, as indicated
Adiabatic coefficient multiplying the ion temperature when computing the sound speeds.
Secondary electron emission coefficient.
If .true. will also call b2stbc_bas boundary condition routine after the end of the physics boundary condition routine (governed by the data from b2ah.dat and b2mn.dat).
Indicates whether a feedback scheme is used. Obsolete. Superceded by 'b2stbc_feedback'.
Number of dead (or isolated) regions.
Density of species (is) in (m-3) to impose in isolated regions.
Electron temperature (in eV) to impose in isolated regions.
Ion temperature (in eV) to impose in isolated regions.
Electric potential (in V) to impose in isolated regions.
Column number of bottom left corner of the dead region (iiso).
Column number of top right corner of the dead region (iiso).
Row number of bottom left corner of the dead region (iiso).
Row number of top right corner of the dead region (iiso).
Name of the next file to use for reading a new /BOUNDARY/ namelist.
When the B2.5 run simulation time, in seconds, modulo(BOUNDARY_TIME_MOD), reaches or exceeds BOUNDARY_TIME_SWITCH, reads the new namelist from BOUNDARY_FILENAME. Also switches to the new namelist as the ELM count (here time/BOUNDARY_TIME_MOD) changes. Only active if BOUNDARY_TIME_MOD is greater than 0.
If BOUNDARY_TIME_MOD is greater than 0, time (in seconds) within an ELM cycle at which a new BOUNDARY namelist is read. Otherwise, B2.5 simulation time, in seconds, when to read the next BOUNDARY namelist file from BOUNDARY_FILENAME. Only active if BOUNDARY_TIME_SWITCH is greater than 0.
List of core boundary segments, according to their numbering in the current /BOUNDARY/ namelist. For limiter, SN, and DN topologies, the code will attempt to build the array from the topology information available, i.e., it will contain the list of 'South' boundaries forming a closed contour. For linear cases, in case a "core" type boundary condition is desired, then LCBS must be specified explicitly.
If .true. (default), writes the content of the namelist to stdout after it has been read.
List of boundary profile files to be imposed as boundary conditions for the density.
List of boundary profile files to be imposed as boundary conditions for the parallel velocity.
List of boundary profile files to be imposed as boundary conditions for the ion energy.
List of boundary profile files to be imposed as boundary conditions for the electron energy.
List of boundary profile files to be imposed as boundary conditions for the electric potential.
If > 0, allows for a communication of particle fluxes across vacuum regions. This option only applies to neutrals. The density boundary condition is based on the difference between the average pressure and the local pressure.
Number of communicating vacuum regions.
Option for resorbing the pressure difference.
Radial index of the ring on which the neutral pressure is computed for region IREG.
Poloidal lower bound of the range over which the neutral pressure is computed for region IREG.
Poloidal upper bound of the range over which the neutral pressure is computed for region IREG.
Parameter for setting the pressure correction. See above.
Parameter for setting the pressure correction. See above.
Sets the target value of the quantity associated with species (ISPECIES) for the feedback scheme defined by NA_FEEDBACK_CHOICE(ISPECIES).
Sets the time of reference (in s) for the feedback of species (ISPECIES).
Choice of quantity on which the feedback is computed:
Option for computing the new fedback quantity.
Choice for the actuator used for the feedback.
Factor by which the rescaling is slowed. Rescaling factor is : Option 1: (1 + alpha*target/current) / (1 + alpha) Option 3: 2**(tanh(log(x)/beta)*log(alpha)/log(2)) Options 4 and 5: Corresponds to parameter F
Factor by which the rescaling is slowed. See above. Options 4 and 5: Corresponds to parameter V (ffb_rtvn)
Ratio of current to target feedback value at which linear rescaling transitions to tanh slowed rescaling for feedback options 7 and 8.
Options 4 and 5: Corresponds to parameter C. If negative, C is computed as the initial total particle content of the sequence.
Lower poloidal bound for the region over which the quantity used for the feedback of species (ISPECIES) is computed, if applicable.
Upper poloidal bound for the region over which the quantity used for the feedback of species (ISPECIES) is computed, if applicable.
Lower radial bound for the region over which the quantity used for the feedback of species (ISPECIES) is computed, if applicable.
Upper poloidal bound for the region over which the quantity used for the feedback of species (ISPECIES) is computed, if applicable.
Index of boundary condition through which the feedback is being applied. If negative refers to the EIRENE stratum index instead. Specifically for actuator 4 it refers to the index of special surface group.
Index of divertor target to which feedback is applied.
Minimum value of actuator being applied.
Maximum value of actuator being applied.
Minimum gas puff being applied. Obsolete.
Maximum gas puff being applied. Obsolete.
If the density is larger than target*overshoot, the gas puff is turned off.
If .true. both feedback target and feedback current are inverted. If .false. actuator value increases in order to increase current feedback value. If .true. actuator value decreases in order to increase current feedback value.
Last value used for the core electron energy radial flux feedback.
Last value used for the core ion energy radial flux feedback.
Last value used for the core particle radial flux feedback. Corresponds to 'isfeedback' B2.5 species.
Last value used for the core radial electric current feedback.
Last value used for the radial inner PFR boundary flux feedback. Corresponds to 'isfeedback' B2.5 species.
Last value used for the radial outer PFR boundary flux feedback. Corresponds to 'isfeedback' B2.5 species.
Last value used for the radial particle flux feedback in the SOL. Corresponds to 'isfeedback' B2.5 species.
Last value used for the feedback control actuator for each isonuclear sequence.
Number of neutral sources, or 'strata'. Must not be larger than DEF_NSTRA from the [$SOLPSTOP/modules/B2.5/]src/modules(.local)/b2mod_dimensions.F file. Need not include the time-dependent stratum. If a time-dependent stratum is used, is incremented internally as needed. The incremented value is referred to as NSTRAT below.
Position in the B2.5 grid of the strata. Similar use as BCPOS from /BOUNDARY/ namelist.
Start coordinate on the B2.5 grid the strata. Similar use as BCSTART from /BOUNDARY/ namelist.
End coordinate on the B2.5 grid of the strata. Similar use as BCEND from /BOUNDARY/ namelist.
Contains the size of the recycling boundary lists. Similar use as BC_LIST_SIZE from /BOUNDARY/ namelist.
Contains the X-coordinate of the cells where recycling boundaries are applied. Similar use as BC_LIST_X from /BOUNDARY/ namelist.
Contains the Y-coordinate of the cells where recycling boundaries are applied. Similar use as BC_LIST_Y from /BOUNDARY/ namelist.
Identifies the base material(s) of this stratum wall. The number corresponds to the B2.5 species index produced by sputtering. Check b2cdci for details. If the bulk species is not sputtered, should contain -1.
Indicates whether chemical sputtering is allowed from this wall stratum.
Stores the particle recycling coefficients of species (is) on stratum (istra). Defaults to 1.0 for non-carbon species and 0.0 for carbon species. Particles recycle into the neutral species associated to their isonuclear sequence. Applies to B2.5 neutral fluid species. Also multiplies Eirene recycling fluxes if 'eirene_use_recyceir' is set to 0 (see b2cdci for details).
Stores the parallel momentum recycling coefficients of species (is) on stratum (istra). The parallel momentum is carried back by the recycled neutrals. Applies only to B2.5 neutral fluid species.
Stores the energy recycling coefficients of species (is) on stratum (istra). Defaults to 0.3 for non-carbon species and 0.0 for carbon species. The energy is carried back by the recycled neutrals. Applies only to B2.5 neutral fluid species.
Stores the particle-into-ion recycling coefficients of species (is) on stratum (istra). Particles recycle into the next ionised species associated to their isonuclear sequence. Applies only to B2.5 neutral fluid species.
Multiplier to the Eirene recycling fluxes from stratum (istra). Only used if 'eirene_use_recyceir' is set to 1 (default).
The element (istra,1) contains the strength of constant gas puff strata (type 'C') for Eirene in particles/second. Its values overwrite the values provided by the FLUX variables from the strata description in block 7 of the Eirene input file.
Contains the type of stratum for Eirene. Possible options include: 'N','S','W','E' - topological mesh boundaries (same use as BCCHAR in /BOUNDARY/ namelist) 'V' - volume recombination source (related RCPOS, RCSTART, RCEND and RC_LIST variables are ignored). 'C' - constant or feedback gas puff source (related RCPOS, RCSTART, RCEND and RC_LIST variables are ignored). Gas puffs for B2.5 fluid neutrals are handled via b2.boundary.parameters. Unless specified otherwise by use of 'eirene_nesepm_istra', it is the first 'C' stratum that is used for feedback puff schemes. See b2cdci for details. 'T' - time-dependent source for Eirene (related RCPOS, RCSTART, RCEND and RC_LIST variables are ignored). Long-lived Eirene neutrals are stored in this stratum. See EIRENE_STEP_DT below.
Relaxation parameters multiplying the Eirene particle, parallel momentum, electron energy and ion energy sources respectively before use in B2.
Stores the physical sputtering multiplier for B2.5 species (is) sputtering on the wall of stratum (istra). Must be non-zero for physical sputtering to occur. Only applies if using B2.5 fluid neutral model.
Stores the chemical sputtering multiplier for B2.5 species (is) sputtering on the wall of stratum (istra). Must be non-zero for chemical sputtering to occur. Only applies if using B2.5 fluid neutral model.
Length of CPU time devoted to Eirene calls after the first one (in s). If zero, the value given in the Eirene input file will be used.
Physical time (in s) the Eirene particles are to be followed before being passed to the time-dependent stratum. If zero, the value given by DTIMV in block 13 of the Eirene input file will be used.
Frequency of Eirene calls. Eirene is called every EIRENE_MOD full B2.5 iterations.
Initial volume recombination rate from stratum (istra) to be used for the volume recombination scaling regulation algorithm. Rendered obsolete by 'eirene_dpc_fix'.
Maximum rate of increase of the volume recombination source strength. The volume recombination regulation scaling scheme is activated if VOLRECINC is neither 0 nor 1. Rendered obsolete by 'eirene_dpc_fix'.
Weight of new volume recombination strength relative to that of previous iteration in the volume recombination regulation scaling scheme. Rendered obsolete by 'eirene_dpc_fix'.
Specifies the index of the first B2.5 species involved in stratum (istra) for EIRENE.
Specifies the index of the last B2.5 species involved in stratum (istra) for EIRENE.
Name of the next file to use for reading a new /NEUTRALS/ namelist.
When the B2.5 run simulation time, in seconds, modulo(NEUTRALS_TIME_MOD), reaches or exceeds NEUTRALS_TIME_SWITCH, reads the new namelist from NEUTRALS_FILENAME. Also switches to the new namelist if the ELM count (here time/NEUTRALS_TIME_MOD) changes. Only active if NEUTRALS_TIME_MOD is greater than 0.
If NEUTRALS_TIME_MOD is greater than 0, time (in seconds) within an ELM cycle at which a new NEUTRALS namelist is read. Otherwise, B2.5 simulation time, in seconds, when to read the next NEUTRALS namelist file from NEUTRALS_FILENAME. Only active if NEUTRALS_TIME_SWITCH is greater than 0.
If l_neutrad >= 0, then the radiated power due to the neutrals atoms is taken directly from the Eirene calculation, instead of being recomputed by B2.
If l_neutflux >=0, then correct treatment of the incident fluxes in B2.5 and b2plot; if <0, then old (approximate) treatment
Indicates with which Eirene atomic species to scale the Eirene stratum (istra) (0 means no scaling, default). Atomic species 0 stands for electrons.
Contains the index of the Eirene atomic species corresponding to the B2.5 species (is).
Contains the isonuclear sequence index of the B2.5 species (is).
Contains the index of the B2.5 neutral fluid species corresponding to the Eirene atomic species (iatm).
Contains the index of the first B2.5 fluid for each species.
Contains the index of the B2.5 isonuclear sequence with which the Eirene atomic species (IATM) should be scaled.
Contains the index of the Eirene atomic species with which the Eirene molecular species (IMOL) should be scaled.
Contains the number of atoms from Eirene atomic species IATM within each molecule of Eirene molecular species IMOL.
Contains the index of the Eirene atomic species with which the Eirene test ion species (IION) should be scaled. If not provided, LIONSCL will map to LMOLSCL for molecular ions that match an existing declared molecule, or to the first species listed in the name of the molecular ion as declared in the Eirene input file if no matching molecule is found.
Contains the indices of Eirene non-standard surfaces (block 3A in input.dat) corresponding to core boundaries.
Contains the indices of Eirene non-standard surfaces (block 3A in input.dat) corresponding to target boundaries.
Contains the indices of Eirene surfaces related to the recycling sources.
Contains the number of Eirene surfaces for each Eirene recycling stratum.
Specifies the fraction of Eirene atomic species (iatm) in one particle puffed from stratum (istra).
Debug output control for eirene_mc routine. See code for usage.
Non-zero elements yield additional print-out data for debugging B2-Eirene coupling. See code for specific uses. Many slots still available for user-specific needs. If debug_flags(81).gt.0, will print out target fluxes in files 'targb2.datv', 'targb2n.datv' and 'targb2pl.datv'. If debug_flags(90).gt.0, will print output about Eirene non-standard surfaces [debug_flags(90):debug_flags(91)] and strata [debug_flags(92):debug_flags(93)] (0 means sum over all strata). Additional output can be obtained for strata [debug_flags(94):debug_flags(95)]. If only debug_flags(90) is specified, the output will include all individual Eirene strata, for all non-standard surfaces, starting from debug_flags(90).
Maximum number by which Eirene neutral sources may be scaled in either direction within the ank-mods scheme. See explaining text in $SOLPSTOP/doc/Source_Scaling_in_B2.pdf for details.
Specifies the mixed material species index related to B2.5 species (is).
Indicates whether mixed material species (itrack) participates in chemical sputtering. Defaults to .true. for first tracked species if sput_dst is carbon, and to .false. for all other cases.
Multiplier to the chemical sputtering coefficient of carbon for computing the rate used for redeposited carbon.
Indicates whether the presence of beryllium is to impact on the chemical sputtering yield of carbon.
The chemical sputtering yield of carbon is multiplied by (1.0-C/2*(tanh((frac-A)/B)-tanh((-A)/B))) where frac is the fractional content of Be in the surface layer material.
See above.
See above.
Number of special surfaces groups for diagnostics.
List of surface segments (non-default standard surfaces [NDSS] or additional surfaces in Eirene notation) included in the groups. Negative numbers correspond to NDSS.
Labels of groups of special surfaces.
Index of the first Eirene surface belonging to a special surface group in the L_SPCSRF list.
Index of the last Eirene surface belonging to a special surface group in the L_SPCSRF list.
Particle absorption on the surface. If positive, will supercede the setting from the Eirene input file: RECYCT = 1-SPS_ABSR. Ignored if negative.
Surface transparency in positive direction. If positive, will supercede the setting from the Eirene input file: TRANSP(1,) = SPS_TRNO. Ignored if negative.
Surface transparency in negative direction. If positive, will supercede the setting from the Eirene input file: TRANSP(2,) = SPS_TRNI. If negative, the setting from SPS_TRNO is used.
Surface material in Eirene notation (e.g., 1206 for C). If positive, will supercede the setting from the Eirene input file: ZNML = SPS_MTRI. Ignored if negative.
Surface material in human notation (e.g., 'C').
Surface temperature (in eV). If set larger to 1.e10, will be ignored. Otherwise, will supercede the setting from the Eirene input file: EWALL = SPS_TMPR.
Fudge factor for physical sputtering. If positive, will supercede the setting from the Eirene input file: RECYCS = SPS_SPPH. Ignored if negative.
Fudge factor for chemical sputtering. If positive, will supercede the setting from the Eirene input file: RECYCC = SPS_SPCH. Ignored if negative.
Sputtered particle species flag for chemical sputtering. If positive, will supercede the setting from the Eirene input file: ISRC = SPS_SGRP. Ignored if negative.
If .true. (default), writes the content of the namelist to stdout after it has been read.
Obsolete. Maximum number of points needed to describe a region contour in Eirene.
Indicates whether a stratum is a time-dependent gas puff. Should only be .true. for strata defined as gas puffs.
Indicates whether a stratum is using a time-dependent gas puff function. Should only be .true. for strata defined as gas puffs. For use with TIME_DEP_PUFF_CASE and TIME_DEP_PUFF_PARAM.
Used to select the functional form of the time-dependent gas puff function. For use with TIME_DEP_PUFF_FUNC and TIME_DEP_PUFF_PARAM.
Used to set the parameters for the functional form of the time-dependent gas puff function. For use with TIME_DEP_PUFF_FUNC and TIME_DEP_PUFF_CASE.
Number of points over which the time profile of the strength of gas puff (istra) is given. Defaults to 0. Current maximum value is 100.
For (i,ik,istra) in (1:NGPDATA(istra),1:2,1:NSTRAT), GPDATA(i,1,istra) contains the time point (i) for stratum (istra) (in s). GPDATA(i,2,istra) contains the gas puff strength of stratum (istra) at time point (i) (in particles/s). The gas puff strength before the first time point is given by USERFLUXPARM(istra,1). The gas puff strength after the last time point is given by the GPDATA value of the last time point. Otherwise, the gas puff strength is linearly interpolated from the given data.
Passed to Eirene. Chemical sputter yield of wall surface (ilim). Element 0 corresponds to the default setting for all surfaces.
Nuclear charge of atomic species causing the sputtering. If set to its default value of 0, then the chemical sputtering switches from this block are not used and the settings from the Eirene input file are used.
Passed to Eirene. Eirene atomic species index of the sputtered particle. If igass_chemical > natmi, the data from chemical_sputter_yield is not used and the yield from the Eirene surface blocks is used instead.
Passed to Eirene. Eirene type index of the sputtered particle. Atoms = 1, Ions = 4. Defaults to 0, meaning 1 eV atom chemical sputtering.
Passed to Eirene. Mass*100 + nuclear charge of the sputtered particle. Defaults to 0, internally changed to 1206 = carbon.
Contains the last value of the strength of volume recombination sources from stratum (istra).
Multiplier to the time used in solving the continuity equation of species (is) in region (ireg).
Multiplier to the time used in solving the parallel momentum equation of species (is) in region (ireg).
Multiplier to the time used in solving the electron heat equation in region (ireg).
Multiplier to the time used in solving the ion heat equation in region (ireg).
Multiplier to the time used in solving the continuity equation of species (is) in the far SOL defined by iy >= iy_SOL_dtXX.
Multiplier to the time used in solving the parallel momentum equation of species (is) in the far SOL defined by iy >= iy_SOL_dtXX.
Multiplier to the time used in solving the electron heat equation in the far SOL defined by iy >= iy_SOL_dtXX).
Multiplier to the time used in solving the ion heat equation in the far SOL defined by iy >= iy_SOL_dtXX.
Defines the region of the SOL (iy >= iy_SOL_dtXX) where the time multipliers dt[XX]_SOL should be applied. If iy_SOL_dtXX.eq.-2 (default), no time multiplication is applied for any equation independently of dt[XX]_SOL values. If iy_SOL_dtXX points to a closed flux surface (i.e. if iy_SOL_dtXX.le.jsep), the code produces an error message and stops.
Indicates whether the continuity equation for species (is) is to be solved in region (ireg). Subservient to 'no_solve' switch.
Indicates whether the parallel momentum equation for species (is) is to be solved in region (ireg). Subservient to 'no_solve' switch.
Indicates whether the total parallel momentum equation is to be solved in region (ireg). For the total parallel momentum equation to be solved in a region, all momentum equations in the range defined by 'b2news_nsmin' and 'b2news_nsmax' must have SOLVEMO true in that region. Subservient to 'no_solve' switch. It is important to note that the 'b2news_no_solve' switch allows for solving or not all momentum equations together as a whole, along with the total equation. Defaults to .true. unless using the time-dependent mode, in which case it is set to .false. and should not be used.
Indicates whether the potential energy equation is to be solved in region (ireg). Subservient to 'no_solve' switch.
Indicates whether the electron energy equation is to be solved in region (ireg). Subservient to 'no_solve' switch.
Indicates whether the ion energy equation is to be solved in region (ireg). Subservient to 'no_solve' switch.
Indicates whether the total energy equation is to be solved in region (ireg). For the equation to be solved in region (ireg), both SOLVEEE and SOLVEEI must be true in that region. Subservient to 'no_solve' switch. It is important to note that the 'b2news_no_solve' switch allows to solve or not to solve the electron and ion heat equations together with total energy equation, as a group. Defaults to .true. unless using the time-dependent mode, in which case it is set to .false. and should not be used.
Minimum time scale of evolution allowed for all equations. Only active is 'b2srsm_enable' is set to 1.
De-multiplier to the timestep in the core. Only active if less than 1. Should be larger than 0. Applies fully to the innermost core ring of cells (IY .eq. -1). See CORE_DT_FACTOR for further use.
Multiplier to the timestep in the core. Only active is less than 1. Should be larger than 0. Multiplies each successive core ring of cells (increasing IY) by CORE_DT_FACTOR, until the local time step multiplier is equal to 1.
Pressure correction speed-up parameter α_a, acting on the density contribution from species a. Will apply in the core region, except for the nrings outer surfaces. See nrings_for_no_speedup_averaging for details. See Pressure_correction_speed-up.pdf in $SOLPSTOP/doc for a full description. Should be roughly equal to corr_core_dt below. Does not apply to neutral species.
Pressure correction speed-up parameter α_T, acting on the temperature contributions. Will apply in the core region, except for the nrings outer surfaces. See nrings_for_no_speedup_averaging for details. See Pressure_correction_speed-up.pdf in $SOLPSTOP/doc for a full description. Should be roughly equal to corr_core_dn above.
Multiplier to the neutral source differential between puffing (plus sputtering plus core boundary flux plus any external sources) and pumping, which is further added to the ionization source for code speed-up. One multiplier must be provided per isonuclear sequence. Such settings might be relevant for e.g. a mixture of main ion and trace impurity. See E. Kaveeva et al., Nucl. Fusion 58 (2018) 126018 for details and keep in mind that in that paper the scheme was described in its original species-unresolved form.
Maximum allowed fraction of ionization source for a given isonuclear sequence from its neutral state (coming from Eirene) to be added or subtracted to the real Eirene source in the method of effective sources for code speed-up.
Parameter which controls exclusion/inclusion of confined region from/in the method of effective sources for code speed-up for each isonuclear sequence. If .true., then the flux through the core boundary of that species is taken into account during the source correction computation, and the source correction is applied on closed flux surfaces (default behavior). If .false., then the particle imbalance is computed as puffing plus sputtering plus any external sources minus pumping without account of the flux through the core boundary, and no source correction is applied on closed surfaces (such a setting might be helpful if the core flux is much bigger than the puffing rate and the core is isolated from the SOL by a strong transport barrier).
Name of the next file to use for reading a new /NUMERICS/ namelist.
When the B2.5 run simulation time, in seconds, modulo(NUMERICS_TIME_MOD), reaches or exceeds NUMERICS_TIME_SWITCH, reads the new namelist from NUMERICS_FILENAME. Also switches to the new namelist as the ELM count (here time/NUMERICS_TIME_MOD) changes. Only active if NUMERICS_TIME_MOD is greater than 0.
If NUMERICS_TIME_MOD is greater than 0, time (in seconds) within an ELM cycle at which a new NUMERICS namelist is read. Otherwise, B2.5 simulation time, in seconds, when to read the next NUMERICS namelist file from NUMERICS_FILENAME. Only active if NUMERICS_TIME_SWITCH is greater than 0.
If .true. (default), writes the content of the namelist to stdout after it has been read.
If .true. (default), adds dte(ix,iy)/qe to the potential correction after the internal energy balance equations are solved, where dte(ix,iy) is the electron temperature correction on the time step. Individually set for each region index.
Number of points over which the source profile of (kind_data,kind_source,is) is defined. Should not exceed NY+2. If KIND_DATA=1, the data is expressed as a profile in physical distance from the separatrix (in metres) along the outer midplane. The user can change this default reference location by use of the 'i[xy]ref_profile' switches. If KIND_DATA=2, the data is expressed as a profile in flux distance from the separatrix (not yet available).
For (i,ir,ik,is) in (1:2,1:NY+2,1:NKIND_SOURCE,0:NS), SDATA(1,ir,:,:) contains the radial location of the profile point (ir). SDATA(2,ir,:,:) contains the source profile value at point (ir). KIND_SOURCE=1 means particle source of species (is) (in particles/m3) KIND_SOURCE=2 means parallel momentum source for species (is) (in kg.m/s/m3) KIND_SOURCE=3 means electron heat source (in Watts/m3) KIND_SOURCE=4 means ion heat source (in Watts/m3) KIND_SOURCE=5 means electric charge source (in Coulombs/m3) KIND_SOURCE=6 means non-ambipolar electron particle source (in e/m3)
Number of points over which the axial profile of (kind_data,kind_source,is) is defined. Should not exceed NX+2. Defaults to 0. If KIND_DATA=1, the data is expressed as a profile in physical distance, here connection length, rescaled from 0.0 to 1.0. For closed field lines, the reference location for the zero of distance is set by use of the 'i[xy]ref_profile' switches. If KIND_DATA=2, the data is expressed as a profile in (ix) cell index, again normalized from 0.0 to 1.0 to match the [0:nx-1] interval.
Multiplier to the poloidal source profile in the axial direction. Same convention for KIND_SOURCE as above.
Additional divertor ion heat source (in Watts/m3)
Poloidal location along which the sources radial profile is measured. If not provided, the value from 'set_transport_ixref' is used.
Radial location along which the sources axial profile is measured. If not provided, the value from 'set_transport_iyref' is used.
Name of the next file to use for reading a new /PROFILE/ namelist. Quantities not present in the new file will be inherited from the old one.
When the B2.5 run simulation time, in seconds, modulo(SOURCES_TIME_MOD), reaches or exceeds SOURCES_TIME_SWITCH, reads the new namelist from SOURCES_FILENAME. Also switches to the new namelist if the ELM count (here time/SOURCES_TIME_MOD) changes. Only active if SOURCES_TIME_MOD is greater than 0.
If SOURCES_TIME_MOD is greater than 0, time (in seconds) within an ELM cycle at which a new PROFILE namelist is read. Otherwise, B2.5 simulation time, in seconds, when to read the next PROFILE namelist file from SOURCES_FILENAME. Only active if SOURCES_TIME_SWITCH is greater than 0.
Array dimension NX used in this namelist.
Array dimension NY used in this namelist.
Array dimension NS used in this namelist.
Contains the chemical sputtering yield of species (is) at position (ix,iy) in element (:,:,:,1) and the physical sputtering yield in element (:,:,:,2).
Contains the energy chemical sputtering yield of species (is) at position (ix,iy) in element (:,:,:,1) and the energy physical sputtering yield in element(:,:,:,2).
Number of points over which the source profile of (kind_data,kind_coef,is) is defined. Should not exceed NY+2. If KIND_DATA=1, the data is expressed as a profiles in physical distance from the separatrix (in metres). If KIND_DATA=2, the data is expressed as a profile in flux distance from the separatrix (not yet available).
For (i,ir,ik,is) in (1:3,1:NY+2,1:NKIND_COEFF,0:NS), TDATA(1,ir,:,:) contains the radial location of the profile point (ir). By default, these are measured at the outer midplane, as distance to the separatrix in metres. The user can change this by setting the 'set_transport_i[xy]ref' switches to choose a different location for the distance reference. TDATA(2,ir,:,:) contains the transport profile value at point (ir). TDATA(3,ir,:,:) contains the ELM transport profile value at point (ir). KIND_COEFF=1 means density-driven particle diffusivity for species (is) KIND_COEFF=2 means pressure-driven particle diffusivity for species (is) KIND_COEFF=3 means ion heat diffusivity for species (is) KIND_COEFF=4 means electron heat diffusivity KIND_COEFF=5 means poloidal pinch velocity for species (is) KIND_COEFF=6 means radial pinch velocity for species (is) KIND_COEFF=7 means viscosity for species (is) KIND_COEFF=8 means field-driven radial current conductivity KIND_COEFF=9 means temperature-driven radial current conductivity
If ADDSPEC(is,ikind,spec).ge.0, then the transport coefficient profile of type (ikind) from species (spec) is also used for species with index ADDSPEC(is,ikind,spec).
Name of the next file to use for reading a new /TRANSPORT/ namelist.
When the B2.5 run simulation time, in seconds, modulo(TRANSPORT_TIME_MOD), reaches or exceeds TRANSPORT_TIME_SWITCH, reads the new namelist from TRANSPORT_FILENAME. Also switches to the new namelist if the ELM count (here time/TRANSPORT_TIME_MOD) changes. Only active if TRANSPORT_TIME_MOD is greater than 0.
If TRANSPORT_TIME_MOD is greater than 0, time (in seconds) within an ELM cycle at which a new TRANSPORT namelist is read. Otherwise, B2.5 simulation time, in seconds, when to read the next TRANSPORT namelist file from TRANSPORT_FILENAME. Only active if TRANSPORT_TIME_SWITCH is greater than 0.
If region_flags(ireg,ikind) is .true. (default), then the transport parameters profiles for region (ireg) and kind (ikind) are used.
If .true., transport coefficients profiles are not implemented in the private flux regions.
If .true., transport coefficients profiles are not implemented in the divertor regions.
If .true., then the transport coefficients profiles from the current b2.transport.inputfile are increased by a factor of 1.0+Gaussian where Gaussian is a Gaussian profile in the poloidal direction of amplitude SCALING_STRENGTH extending from SCALING_IX_BEGIN to SCALING_IX_END inclusively. The profile has a decay length of SCALING_WIDTH (in units of the number of poloidal cells).
See above.
See above. Defaults to about 1/3 of the interval over which the scaling is to be done.
See above.
See above.
Time (in seconds) modulo ELM_TIME_PERIOD at which the ELM phase begins and the ELM data must be used.
Time (in seconds) modulo ELM_TIME_PERIOD at which the ELM phase ends and the ELM data is no longer used.
Indicates the real frequency of simulated ELMs. See above for usage. If zero, no ELM profiles are used.
Poloidal position at which the ELM profile starts to be applied.
Poloidal position at which the ELM profile ends being applied.
All flags follow: FLAG=0: Use model from b2ah.dat and b2mn.dat. Default. FLAG=1: Constant, value set to PARM. FLAG=2: 1/N model, multiplied by PARM. FLAG=3: Bohm model, multiplied by PARM. FLAG=4: flux-scaled model, multiplied by PARM.
Transport model flag for density-driven diffusion.
Transport model flag for pressure-driven dffusion.
Transport model flag for pinch velocity.
Transport model flag for viscosity.
Transport model flag for ion heat diffusivity.
Transport model flag for electron heat diffusivity.
Transport model flag for field-driven current radial conductivity.
Transport model flag for temperature-driven current radial conductivity.
Parameter for the density-driven particle diffusion coefficient for species (is).
Parameter for the pressure-driven particle diffusion coefficient for species (is).
Parameter for the anomalous radial pinch velocity for species (is).
Parameter for the viscosity for species (is).
Parameter for the heat diffusivity coefficient for species (is).
Parameter for the electron heat diffusivity coefficient.
Parameter for the anomalous radial field-driven current conductivity.
Parameter for the anomalous radial temperature-driven current conductivity.
Name of the next file to use for reading a new /TRANSPORT/ namelist.
When the B2.5 run simulation time, in seconds, modulo(TRANSPORT_TIME_MOD), reaches or exceeds TRANSPORT_TIME_SWITCH, reads the new namelist from TRANSPORT_FILENAME. Also switches to the new namelist if the ELM count (here time/TRANSPORT_TIME_MOD) changes. Only active if TRANSPORT_TIME_MOD is greater than 0.
If TRANSPORT_TIME_MOD is greater than 0, time (in seconds) within an ELM cycle at which a new TRANSPORT namelist is read. Otherwise, B2.5 simulation time, in seconds, when to read the next TRANSPORT namelist file from TRANSPORT_FILENAME. Only active if TRANSPORT_TIME_SWITCH is greater than 0.
For the electron, ion, and viscous heat flux limit coefficients, if the main multiplier is not specified but CFLM*_CORE and CFLM*_SOL are set to equal non-zero values, their value will be used. If CFLM*_CORE and CFLM*_SOL are different and non-zero, then a radial tanh profile, centered at radial OMP position CFLM*_TANH_A (in metres) and of width CFLM*_TANH_B is applied. The flux limit coefficients in the PFR are set to a constant value equal to that at the separatrix. The OMP is located at the poloidal index provided by 'b2news_fac_ref', or, if not specified, 'b2mwti_jxa'. For all flux limit multipliers, setting them to zero turns off the corresponding flux limiter. Otherwise, the multipliers must be positive.
Main multiplier to the electron heat flux limit.
Main multiplier to the ion heat flux limit.
Main multiplier to the viscous heat flux limit.
Multiplier to the thermo-electric coefficient flux limit.
Multiplier to the friction force flux limit.
Limit to multiplier to the electron heat flux limit in the deep core.
Limit to multiplier to the electron heat flux limit in the far SOL.
Radial location, measured (in metres) as a radial distance along the OMP away from the separatrix, for the center of the tanh profile multiplying the electron heat flux limit.
Radial width, measured (in metres) as a radial distance along the OMP, for the tanh profile multiplying the electron heat flux limit. Must be positive.
Limit to multiplier to the ion heat flux limit in the deep core.
Limit to multiplier to the ion heat flux limit in the far SOL.
Radial location, measured (in metres) as a radial distance along the OMP away from the separatrix, for the center of the tanh profile multiplying the ion heat flux limit.
Radial width, measured (in metres) as a radial distance along the OMP, for the tanh profile multiplying the ion heat flux limit. Must be positive.
Limit to multiplier to the viscous heat flux limit in the deep core.
Limit to multiplier to the viscous heat flux limit in the far SOL.
Radial location, measured (in metres) as a radial distance along the OMP away from the separatrix, for the center of the tanh profile multiplying the viscous heat flux limit.
Radial width, measured (in metres) as a radial distance along the OMP, for the tanh profile multiplying the viscous heat flux limit. Must be positive.
If .true. (default), writes the content of the namelist to stdout after it has been read.
The four *_CNV parameters are only activated if 'b2tqna_user_transport' is set to '6'.
Value of the "blob" convection velocity for species (is) at the outer grid edge (in m/s).
Exponent in radial profile of the "blob" velocity for species (is).
First exponent in poloidal profile of the "blob" velocity for species (is).
Second exponent in a poloidal profile of the "blob" velocity for species (is).
Used by the user specified set_transport_eta transport model. See code for details.
List of surface indices (EIRENE notation) which are used for calculation of helium enrichment.
B2.5 x-cell index corresponding to the first edge of the bypass between the inner and outer divertor (Western divertor, edge closest to target). If non-positive, counted backwards from the X-point location in the lower PFR, from the inner upper target in the upper PFR, from the lower outer target in the outer SOL, and from the inner upper target in the inner SOL.
B2.5 x-cell index corresponding to the first edge of the bypass between the inner and outer divertor (Eastern divertor, edge closest to target). If non-positive, counted backwards from the outer lower target in the lower PFR, from the lower outer target in the outer SOL, from the upper outer target in the upper PFR, and from the inner upper target in the inner SOL.
B2.5 x-cell index corresponding to the second edge of the bypass between the inner and outer divertor (Western divertor, edge furthest to target). If non-positive, counted backwards from the X-point location in the PFR, from the outer target in the outer SOL for a SN, from the top targets in the SOL for a DN.
B2.5 x-cell index corresponding to the second edge of the bypass between the inner and outer divertor (Eastern divertor, edge furthest to target). If non-positive, counted backwards from the X-point location in the PFR, from the outer target in the outer SOL for a SN, from the top targets in the SOL for a DN.
Species index of the helium atoms in Eirene. The code attempts to find a match by default.
Species index of the neon atoms in Eirene. The code attempts to find a match by default.
Species indices of the hydrogen isotopes in Eirene. If using only one hydrogen species, only the first element needs to be provided. Otherwise the 3 elements correspond to H/D/T. The code attempts to find a match by default.
Number of hydrogen isotope nuclei for molecules in Eirene. If using only one hydrogen species, only the first element needs to be provided. Otherwise the 3 elements correspond to H/D/T. The code attempts to find a match by default.
Fusion power occuring in core (including neutrons, in Megawatts).
Ratio of He to DT pumping speeds (typically, 0.8).
Location of the pump. 0 no pump at all (default), 1 - SN of lower DN PFR, 2 - (outer) SOL, 3 - upper DN PFR, 4 - inner SOL for DN
Actual number of surface groups for PFR flows.
List of surface segments for PFR flows.
First positions in this list for each group.
Last positions in this list for each group.
Labels for surface groups for PFR flows.
Actual number of surface groups for neutral data.
List of surface segments for neutral data.
First positions in this list for each group.
Last positions in this list for each group.
Labels for the neutral data surface groups.
Nominal pumping speed.
Electron temperature at the target (in eV) which is used a threshold identifying the transition of the local magnetic flux tube to detachment. The width of the zone with Te less than TE_DET_THRESHOLD (counted from the strike point towards the SOL side) is outputted into user_SPb.trc. This distance monotonically increases with increasing of radiated power fraction or with decreasing of peak target energy loads, therefore it may be used as a measure of the degree of detachment.
Switches to turn on/off the tracing files controlled by the ank_tracing switch. The files are defined in b2mod_diag, in order:
Filename where /USER/ namelist is stored.
Dimension NDEPTH used for the arrays within this namelist. Represents the number of depth layer discretising the wall elements for the wall model. If using the 0-D model or the time-independent 1-D model, will contain 1 (default). Should not exceed the value of the parameter NDEPTH declared in b2mod_wall.F.
Indicates the (ix) position in the grid of wall element (iwall). Default ordering of the wall elements proceeds clockwise in physical space from the lower left corner.
Indicates the (iy) position in the grid of wall element (iwall). Default ordering of the wall elements proceeds clockwise in physical space from the lower left corner.
XYMAP(ix,iy) contains the wall index (iwall) of the wall element located a grid locaion (ix,iy). If there is no wall element at this position, contains 0. Default ordering of the wall elements proceeds clockwise in physical space from the lower left corner.
Contains the filename from which to extract the surface material properties for wall element (iwall). The files are to be found in $SOLPSTOP/data.local/Surface_properties/ or $SOLPSTOP/modules/B2.5/Database/Surface_properties/.
Contains the filename from which to extract the bulk material properties of the eventual coating for wall element (iwall).
Contains the filename from which to extract the bulk material properties for wall element (iwall). The files are to be found in $SOLPSTOP/data.local/Bulk_properties/ or $SOLPSTOP/modules/B2.5/Database/Bulk_properties/.
Contains the filename from which to extract the material properties for alloy (nalloy) which may be present in mixed materials deposited layers. Not yet operational. The files are to be found in $SOLPSTOP/data.local/Bulk_properties/, $SOLPSTOP/data.local/Surface_properties/, $SOLPSTOP/modules/B2.5/Database/Surface_properties/, or $SOLPSTOP/modules/B2.5/Database/Bulk_properties/.
Contains the temperature (in Kelvin) of wall element (iwall) at depth layer (idepth). If plate_option.eq.1, will be set to backplate_temp(iwall). If plate_option.eq.2 and empty, will be set to equilibrium 1-D profile deduced from plasma incoming fluxes. If plate-option.eq.3, must be set.
Indicates whether wall element (iwall) is inertially cooled instead of actively cooled.
Contains the temperature (in Kelvin) maintained by cooling at the back end of wall element (iwall).
Contains the thickness (in metres) of wall element (iwall).
Contains the thickness (in metres) of the eventual coating on wall element (iwall).
Multiplier to the time for the equations for temperature and composition evolution of wall element (iwall).
Contains the amount of deposited material (in atoms) from species (itrack) onto wall element (iwall).
Contains the amount of eroded material (in atoms) of species (itrack) from wall element(iwall).
Contains the chemical sputter yield of 'sput_dst' or 'sput_dst_bulk' species on wall element (iwall) caused by B2.5 species (is).
Contains the physical sputter yield of 'sput_dst' or 'sput_dst_bulk' species on wall element (iwall) caused by B2.5 species (is).
Contains the fraction of returned energy carried by sputtered particles of species 'sput_dst' or 'sput_dst_bulk' species from wall element (iwall) caused by B2.5 species (is).
Contains the RES sputter yield of 'sput_dst' or 'sput_dst_bulk' species on wall element (iwall) caused by B2.5 species (is).
Contains the rate of thermal evaporation of species (is) (in particles/second) from wall element (iwall).
Contains the backscattering fraction for incoming B2.5 species (is) onto wall element (iwall).
Contains the backscattered energy fraction for incoming B2.5 species (is) onto wall element (iwall).
Indicates how much simulation time has elapsed for wall element (iwall) (in seconds).
Indicates the wall area (in square metres) for wall element (iwall). Defaults to the area computed from the B2.5 grid information.
Contains the amount of deposited material (in monolayers) from species (itrack) onto wall element (iwall).
Contains the amount of eroded material (in monolayers) of species (itrack) from wall element (iwall).
Indicates the number of elemental constituents within the surface layer of wall element (NWALL).
Contains the atomic numbers Z of the various elements present within the surface layer of wall element (iwall). Defaults to 6 for the first value, 0 otherwise.
Contains the relative atomic abundances of the various elements present within the surface layer of wall element (iwall). Defaults to 1.0 for the first value, 0.0 otherwise. Replaces LAYER_NRELCONSTITUENTS.
Name of the experiment being modelled.
Shot number identifying the run. Defaults to the last number found in shotnumber.history, or 0 if the file is not found.
Time point of the experimental shot being simulated.
Label for the run.
If .true. (default), saves data from b2time.nc.
If .true. (default), saves data from b2fplasma.
If .true. (default), saves data from b2tallies.nc.
If .true. (default), saves data from b2movies.nc.
Indicate the shot number to overwrite (to be used only when updating an already saved run with 'resave_mds' script). Defaults to the last number found in shotnumber.history, or 0 if the file is not found.
Number of coils in the device. For each coils the switches R1, Z1, R2, Z2, CURRENT, windings and nlayers have to be provided.
Inner radius of the coil in units of m.
Axial start position of the coil in units of m.
Outer radius of the coil in units of m.
Axial end position of the coil in units of m.
Current in the coil in units of A.
Number of windigs of each coil.
Number of layers to divide each coil in the radial direction.
Initial spacial step size used in the calculation of the field in units of m.
Toroidal magnetic field in units of T.
Two column file with the profile of axial cell size in the numerical space. First line of the file determine the number of rows to read.
Two column file with the profile of radial cell size in the numerical space. First line of the file determine the number of rows to read.
Radial position in m of the northeastmost point of the plasma region.
Axial position in m of the northeastmost point of the plasma region.
Radial position in m of the southeastmost point of the plasma region.
Axial position in m of the southeastmost point of the plasma region.
Radial position in m of the southwestmost point of the plasma region.
Axial position in m of the southwestmost point of the plasma region.
Radial position in m of the northwestmost point of the plasma region.
Axial position in m of the northwestmost point of the plasma region.
Starting radial point in m for the distribution of cells provided by fn_r.
End radial point in m for the distribution of cells provided by fn_r.
Starting axial point in m for the distribution of cells provided by fn_z.
End axial point in m for the distribution of cells provided by fn_z.
Guard cell size for the North boundary in m.
Guard cell size for the South boundary in m.
Guard cell size for the West boundary in m.
Guard cell size for the East boundary in m.
Number of external regions to be used in the grid.
Number of nodes per region. Minimum of 3 is required per region.
Radial positions in m of the nodes defininig the external region
Axial positions in m of the nodes defininig the external region
Radial starting index of the node on the B2.5 boundary for the region.
Axial starting index of the node on the B2.5 boundary for the region.
Radial final index of the node on the B2.5 boundary for the region.
Axial final index of the node on the B2.5 boundary for the region.
If .true., writes the content of the namelist to stdout after it has been read.