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.
In the WG code, BCHAR is always equal to 'A' (automatic) and a boundary list is defined by the mesh boundary faces that have fcLbl.ge.BCSTART and fcLbl.le.BCEND.
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. Obsolete for WG.
Starting label index (fcLbl) of the cells belonging to this boundary condition.
Ending label index (fcLbl) of the cells belonging to this boundary condition.
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. Obsolete for WG.
Contains the Y-coordinate of the cells where boundaries conditions are applied. Obsolete for WG.
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
Specifying the type of boundary condition for the turbulent kinetic energy (kt) for each segment; makes use of ENKPAR to specify additional information, as indicated
Specifying the type of boundary condition for the turbulent enstrophy (zt) for each segment; makes use of ENZPAR to specify additional information, as indicated
Adiabatic coefficient multiplying the ion temperature when computing the sound speeds.
Secondary electron emission coefficient.
Obsolete.
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 cases with closed field lines, 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.
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.
Number of feedback schemes applied.
Choice of quantity on which the feedback is computed:
Option for computing the new fedback quantity.
Choice for the actuator used for the feedback.
Sets the target value of the quantity associated with feedback (IFB) for the feedback scheme defined by FB_TYPE(IFB).
Specifies the species used to calculate the controlled variable and to which he actuator is applied. It follows the B2.5 plasma species enumeration. In general, controlled quantities are calculated on the whole isonuclear sequence to which fb_species(ifb) belongs to. Then if fb_actuator = 1, it is applied to neutral species through gas puff, if fb_actuator=2 it is applied to the ionized species only, if actuator=3 it is applied only to fb_species. For fb_type=6 and actuator=3 then it is also applied to all ionized species (as old feedback treatment). Overridden by switch 'b2stbc_isfeedback'
Sets the time of reference (in s) for the feedback (IFB).
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 Can be overridden by switches 'b2stbc_nesepm_alpha', 'b2stbc_fnaycore_alpha', 'b2stbc_fheycore_alpha', 'b2stbc_fhiycore_alpha', 'b2stbc_fchycore_alpha', 'b2stbc_volrec_alpha'
Factor by which the rescaling is slowed. See above. Options 4 and 5: Corresponds to parameter V (ffb_rtvn). Cna be overridden by 'b2stbc_volrec_beta'
Options 4 and 5: Corresponds to parameter C. If negative, C is computed as the initial total particle content of the sequence.
Index of boundary condition through which the feedback is being applied.
Eirene strata where the gas puff feedback is applied. Overridden by switch 'eirene_nesepm_istra'
Minimum gas puff being applied. Can be overrdidden by switch 'b2stbc_nesepm_minpuff'
Maximum gas puff being applied. Can be overrdidden by switch 'b2stbc_nesepm_maxpuff'
If the density is larger than target*overshoot, the gas puff is turned off. Can be overridden by switch 'b2stbc_nesepm_overshoot', 'b2stbc_volrec_overshoot'
Pointing for each feedback to the first index in FB_REG, and its number of CVs or faces
Listing for each feedback the cells or faces where controlled quantity is calculated. If different feedbacks use the same domain then this can be defined only ones and simply point to it twice with fb_regP. nMxFbReg is a maximum length set in b2us_feedback
If the two FB_REG and FB_REGP are not defined, then the code will look for this array. For fb_type 6-10-11-12-13, fb_reg_par(ifb,:) indicates the starting (iFb,1) and ending (iFb,2) face region label over which the controlled quantity is calculated. For fb_type 1-2-4-5-8, fb_reg_par(ifb,:) indicates the starting (ifB,1) and ending (iFb,2) volume region label over which the controlled quantity is calculated.
Last value used for the feedback actuator for each feedback.
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. Obsolete for WG.
Starting label index on the B2.5 grid. Similar use as BCSTART from /BOUNDARY/ namelist.
Ending label index on the B2.5 grid. 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. Obsolete for WG.
Contains the Y-coordinate of the cells where recycling boundaries are applied. Similar use as BC_LIST_Y from /BOUNDARY/ namelist. Obsolete for WG.
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. If B2RECYC is not explicitly specified, the values from RECYC will be taken.
Multiplies Eirene recycling fluxes if 'eirene_use_recyceir' is set to 0 (see b2cdci for details). If B2RECYC is not specified, RECYC will also be applied to B2.5 neutral fluid species.
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: 'A' - (automatic) topological mesh boundaries (same use as BCCHAR in /BOUNDARY/ namelist). The boundary mesh faces that have fcLbl.ge.RCSTART and fcLbl.le.RCEND belong to this stratum. '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.
Scaling factor for the fluid neutral B2.5 sources stemming from the stratum. This scaling is required for the spatially hybrid approach, where the hybrid model for hydrogenic strata is typically combined with a fully kinetic model for other species. In that case, STRASCLFL has to made 1.0 for hydrogenic strata and STRASCLFL has to become small for the other strata. STRASCLFL is a stratum-dependent multiplier in addition to b2mndr_rescale_neutrals_sources. b2mndr_rescale_neutrals_sources is used for all strata and has to be made 1.0 for the hybrid model.
Contains the surface material for surface strata, used for AFN neutrals. Examples: 'Be' (beryllium), 'C' (carbon), 'Fe' (iron), 'Mo' (molybdenum) and 'W' (tungsten).
If 1, adding sheath acceleration for the ions hitting the boundary. Used for AFN boundary conditions.
Determines the approximation of the particle-velocity distribution of the incident atoms for the AFN boundary condition. MAXW=0: a diffusion approach is used; MAXW=1: a truncated drifting Maxwellian is used. For some cases, MAXW=0 slightly reduces the discrepancy between the fluid and kinetic solution due to an intrinsic correction for lower collisionality. The diffusion and Maxwellian approach are combined when using MAXW=2, for which the Maxwellian approach is used when the local Knudsen number (Kn) is smaller than b2stbr_kn_b1, the diffusion approach is used for Kn > b2stbr_kn_b2 and a linear interpolation between the two approaches for b2stbr_kn_b1 < Kn < b2stbr_kn_b2. The macroscopic length scale to determine the local Knudsen number is defined by the parameter l_macro_afn.
Option to define a purely molecular stratum. Useful if you want to treat only molecules kinetically, whereas the recycled atoms are treated on the B2.5 fluid side.
When HYB_TYPE='A' (automatic), possibility to have automatic switching between fluid and kinetic recycling for the stratum for the advanced hybrid model. See description of use_auto_spatial_hyb. HYB_TYPE should be 'M' (manual) for all strata for multispecies simulations (ns > 1).
Franck-Condon energy assigned to the thermally released fraction of recycled hydrogen neutrals, used in the AFN boundary condition. It is assumed that thermally released molecules are dissociated immediately at the surface and get an energy E_FC (eV) and are emitted isotropically. 3.0 eV is the value typically used in the EIRENE input file.
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 the B2.5 fluid neutral model. If not specified, B2SPECIES_START = SPECIES_START.
Specifies the index of the last B2.5 species involved in stratum (istra) for the B2.5 fluid neutral model. If not specified, B2SPECIES_END = SPECIES_END.
Specifies the index of the first B2.5 species involved in stratum (istra) for EIRENE. SPECIES_START is also used for the B2.5 fluid neutral model, if B2SPECIES_START is not specified.
Specifies the index of the last B2.5 species involved in stratum (istra) for EIRENE. SPECIES_END is also used for the B2.5 fluid neutral model, if B2SPECIES_END is not specified.
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.
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 fluid neutral heat equation in region (ireg).
List of flux tube indices for which individual timestep multipliers can be applied.
Individual timestep multiplier to all continuity equations applied in the corresponding flux tube from the DTFTS list. Must be positive.
Individual timestep multiplier to all parallel momentum equations applied in the corresponding flux tube from the DTFTS list. Must be positive.
Individual timestep multiplier to the electron heat equation applied in the corresponding flux tube from the DTFTS list. Must be positive.
Individual timestep multiplier to the ion heat equation applied in the corresponding flux tube from the DTFTS list. Must be positive.
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 fluid neutral 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 the 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.
Indicates whether the ExB turbulent kinetic energy equation is to be solved in region (ireg). Subservient to 'no_solve' switch.
Indicates whether the ExB turbulent enstrophy equation is to be solved in region (ireg). Subservient to 'no_solve' switch.
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 cost functions (icf = 1, NCF).
Specifies the type of cost function for index icf. Valid options are: 0 - Sum of weighted least squares cost functions. CFSTART, CFEND indicate over which icf indexes the cost functions are to be summed (valid indexes are 1-3 and 7-9). 1 - Difference, on the desired CVs, of calculated electron density and experimental values read from file, normalized with average experimental value (units must be 10^{19}m^{-3}). 2 - Difference, on the desired CVs, of calculated electron temperature and values read from file, normalized with average experimental value (units must be eV). 3 - Difference, on the desired CVs, of calculated ion temperature and values read from file, normalized with average experimental value (units must be eV). 4 - Maximum electron temperature on the desired CVs (output in eV). 5 - Maximum heat flux on the desired FCs (output in W^{2}/m^{2}). 6 - Bayesian MAP cost function. CFSTART, CFEND indicate over which icf indexes the cost functions are to be summed (valid indexes are 1-3 and 7-9). 7 - Difference, on the desired CVs, of calculated electron density radial gradient and experimental values read from file, normalized with average experimental value (units must be 10^{19}m^{-4}). 8 - Difference, on the desired CVs, of calculated electron temperature radial gradient and values read from file, normalized with average experimental value (units must be eV/m). 9 - Difference, on the desired CVs, of calculated ion temperature radial gradient and values read from file, normalized with average experimental value (units must be eV/m). 10 - Difference, on the desired CVs, of calculated ion saturation current (as ne*cs*eV) and values read from file, normalized with average experimental value (units must be kA/m^{2}). 11 - Heat flux peakedness on the desired FCs (output in W^{2}/m^{2}). 12 - Difference, on the desired FCs, of calculated perpendicular or parallel heat flux and values read from file, normalized with average experimental value (units must be in MW^{2}/m^{2}). Use PARALLEL_HF to specify whehter the parallel or perpendicular heat flux is used.
Indicates for cost function 12 if the heat flux considered is the one parallel to the magnetic field (TRUE) or perpendicular to the surface (FALSE).
Indicates how the domain of cost function ifc is specified. Valid options are: 1 - The cost function is defined on the OMP, using the rzomp array. In this case CFSTART,CFEND are not necessary and ignored. 2 - The cost function is defined using the FC labels specified in CFSTART,CFEND. Cost functions requiring CVs or FCs can both be defined in this way. 3 - The cost function is defined using the CV labels specified in CFSTART,CFEND. Only cost functions requiring CVs can be defined in this way. 4 - The cost function is defined using the CF_REG and CF_REGP arrays. Cost functions requiring CVs or FCs can both be defined in this way.
Specify the starting/ending, face (FC) label fcLbl or control volume (CV) label cvLbl where the cost function icf will be evaluated. For CFTYPE = 0 and CFTYPE = 6 they instead identify which cost functions will be summed together.
Listing for each cost function the corresponding domain CVs or FCs. WARNING! No attempt is made to check if these CVs or FCs are connected, sorted or even available in the current geometry: it is up to the user to make sure the list is correct. mxnCf is a hard coded maximum length defined in b2mod_par_opt.
CF_REGP(icf, 1) points, for cost function icf, to the first index in the CF_REG list. CF_REGP(icf, 2) number of items in CF_REG list for cost function icf.
Multipliers to re-scale each cost function separately.
Indicates if cost function icf is to be re-mapped to the OMP, according to CFTYPE. Warning! The mapping is a rigid translation of the data to the OMP and is safe to work only for non-extended grids and for cost functions defined on the same number of elements as the OMP list.
Number of standard deviation variables used for Bayesian MAP cost function (normally this should be equal to the number of cost functions summed for CFTYPE = 6, i.e. one for each variable on each domain) (isigma = 1, NSIGMA}).
Value of the prediction error standard deviation. For absolute STD the units are the same as the experimental data used for the cost function it refers to.
Indicates if SIGMA(isigma) is intended as relative standard deviation (SCALE_SIGMA(isigma)=.TRUE.) or absolute (SCALE_SIGMA(isigma)=.FALSE.).
(Experimental! Do not use!) Specifies if sigma can be read from the experimental data files cfi.dat.
Number of predictin error mean variables used for Bayesian MAP cost function.
Value of the predictin error mean.
Defines the type of prior distribution for each parameter for the MAP cost function. Possible values are: 0 - Uniform distribution (i.e. $\pi(\theta)=1.0$). 1 - Uninformative Gaussian prior. 2 - Gamma distribution, defined through mean and standard deviation.
Indicates parameters for the prior distributions. For PRIOR_TYPE=0, not used. For PRIOR_TYPE=1 and 2, PRIOR_TYPE(ii,1) indicates the mean and PRIOR_TYPE(ii,2) indicates the standard deviation of the distribution.
Indicates the valid range for the prior distributions, outside of which the program will return an infinite cost function value. PRIOR_RANGE(:,1) sets the lower range and PRIOR_RANGE(:,2) sets the upper range. (p_infty is a parameter currently set to 1e30).
For each cost function where experimental data is read, specifies whether or not to shift the separatrix position in the experimental data spatial coordinate as follows: (r - r_{sep}) - SHIFT_VALUE$. Unique integers larger than zero for each cost function means that each has its specific shift value. Using the same integer for two (or more) cost functions makes to code use the same shift value for those cost functions (the first SHIFT_VALUE among the cost function with the same integer is adopted). Possible usage: Thomson scattering data of density and temperature which should be defined in two separate cost functions but come from the same source and thus have the same shift. NSHIFT is then the number of unique integers and is used to define other parameters below other parameters.
Defines the value of the shift in millimeters for each cost function with SHIFT_CF_DATA .neq. 0. Tis is also used as initial guess in case of optimization.
Same as PRIOR_TYPE but defined for each unique cost function shift.
Same as PRIOR_PAR but defined for each unique cost function shift.
Same as PRIOR_RANGE but defined for each unique cost function shift.
Specifies the type of correlation model for the covariance matrix in the MAP cost function.Possible values are: 0 - No correlation, thus simple diagonal covariance matrix. 1 - Correlation model based on exponential decay with a correlation length specified using CORR_LENGTH.
Defines the value of the correlation length in millimeters for each MAP-type cost function with CORR_MODEL>0. This is also used as initial guess in case of optimization.
For the exponential decay correlation model, off-diagonal elements with value of the exponential below this threshold are set to zero.
Same as PRIOR_TYPE but defined for each MAP cost function where sensitivity for CORR_LENGTH is calculated or this parameter optimized (see CORR_OPT). NCORR_OPT is the number of correlation lengths which are optimized, so for which a prior is required.
Same as PRIOR_PAR but defined for each MAP cost function where sensitivity for CORR_LENGTH is calculated or this parameter optimized.
Same as PRIOR_RANGE but defined for each MAP cost function where sensitivity for CORR_LENGTH is calculated or this parameter optimized.
Number of sensitivity/optimization parameters (radially varying transport coefficients count as 1).
Indicates the type of sensitivity/optimization parameter. Valid options are: -1 - mean, error mean in likelihood function for Bayesian inference (one for each plasma quantity and location used). If present, these must ALWAYS be the last parameters in the vector! 0 - sigma, standard deviation in likelihood function for Bayesian inference (one for each plasma quantity and location used). If present, these must ALWAYS be defined after each physical parameter (PARTYPE>0) and before the mean (PARTYPE=-1). 1 - parm_dna or tdata(:,:,1,:) 2 - parm_dpa or tdata(:,:,2,:) 3 - parm_hci or tdata(:,:,3,:) 4 - parm_hce or tdata(:,:,4,:) 5 - tdata(:,:,5,:) (not yet tested!) 6 - parm_vla or tdata(:,:,6,:) 7 - parm_vsa or tdata(:,:,7,:) 8 - parm_sig or tdata(:,:,8,:) 9 - parm_alf or tdata(:,:,9,:) 10 - enepar 11 - enipar 12 - conpar 13 - mompar 14 - potpar 15 - enkpar 16 - b2recyc 17 - b2tqna_ballooning 18 - b2tqna_ballooning_rescale 19 - keps_cd 20 - keps_heat 21 - keps_heat_i 22 - keps_sig 23 - keps_alf 24 - keps_visc 25 - keps_dkt 26 - keps_dzt 27 - b2sikt_fac_diss 28 - b2sikt_fac_diss_core 29 - b2sikt_fac_sheath 30 - b2sikt_fac_sheath_core 31 - keps_shear 32 - b2tfhi_fsigkt 33 - b2sikt_fac_vis_RS 34 - b2tfhi_fflokt 35 - b2tfhi_fconkt 36 - b2tfhi_fflozt 37 - b2tfhi_fconzt 38 - b2tfhi_fkt_hie 39 - b2tfhe_vis_kt
For each sigma, specifies whether its sensitivity is computed or not, and thus whether it is optimized or fixed.
For each mean, specifies whether its sensitivity is computed or not, and thus whether it is optimized or fixed.
For each SHIFT_CF_VALUE, specifies whether its sensitivity is computed or not, and whether is optimized or fixed.
For each MAP cost function where CORR_MODEL>0, specifies whether CORR_LENGTH sensitivity is computed or not, and whether is optimized or fixed.
For each sensitivity/optimization variables, specifies if radially varying coefficients are optimized. Meaningful only for partype = 1 to 9.
Number of spatial points for sensitivity/optimization parameter ipar (should be the same as ndata in b2.transport.inputfile). The actual number of optimized variables (NPAR_OPT) is then equal to the total number of spatial points + any other additional parameter.
Specifies the species of the sensitivity/optimization parameter. Only meaningful for partype = [1-3,5,7,12,13,16].
Specifies the boundary/strata of the sensitivity/optimization parameter. Only meaningful for partype = [10-16].
Coefficient to rescale optimization parameters (and their gradient). Note dimension is NPAR_OPT, so need to specify it also for each point in case of spatially varying coefficients.
Initial guess for each optimization parameter. Must be defined (p_infty is a parameter currently set to 1e30).
Lower bound for each optimization parameter.
Upper bound for each optimization parameter.
Lower bound for SHIFT_CF_DATA parameters.
Upper bound for SHIFT_CF_DATA parameters.
Lower bound for CORR_LENGTH parameters. Must be specified for each cost function!
Upper bound for CORR_LENGTH parameters. Must be specified for each cost function!
Same as PAR_RESCALE but defined for CORR_LENGTH} parameters. Must be specified for each cost function!
Maximum number of optimization iterations. Can overridden with PETCs/Tao by using command line options for the optimization.
Maximum amount of CPU time in seconds after which the optimization is stopped.
Tolerance below which optimization is stopped. For PETCs/Tao it is used for the absolute gradient norm, the relative gradient norm, and the gradient reduction stopping criteria. Can be overridden with PETCs/Tao by using command line options for the optimization.
Type of Hessian approximation employed. Depends on optimization library. Possible values are currently: 'limited-memory' and 'exact'. When using the 'exact' option the steepest descent algorithm is used.
Specifies which kind of Hessian approximation is employed if HESSIAN_APPROXIMATION is set to 'limited-memory'. Can overridden with PETCs/Tao by using command line options for the optimization.
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 'set_transport_i[xy]ref' 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 'set_transport_i[xy]ref' 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)
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_IP_TIME_MOD), reaches or exceeds TRANSPORT_IP_TIME_SWITCH, reads the new namelist from TRANSPORT_IP_FILENAME. Also switches to the new namelist if the ELM count (here time/TRANSPORT_IP_TIME_MOD) changes. Only active if TRANSPORT_IP_TIME_MOD is greater than 0.
If TRANSPORT_IP_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_IP_FILENAME. Only active if TRANSPORT_IP_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 all flux limit multipliers, setting them to zero turns off the corresponding flux limiter. Otherwise, the multipliers must be positive.
Multiplier to the electron heat flux limit.
Multiplier to the ion heat flux limit.
Multiplier to the viscous heat flux limit.
Multiplier to the thermo-electric coefficient flux limit.
Multiplier to the friction force flux limit.
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.
Coordinates defining a segment. All control volumes intersected by this segment will be part of the outer midplane. rzomp(1,1:2) holds the R-coordinates of the first and second point of the segment. rzomp(2,1:2) holds the Z-coordinates of the first and second point of the segment. Warning! This variable needs always be defined for SOLPS-ITER version 3.2.0 or younger!
Coordinates defining a segment. All control volumes intersected by this segment will be part of the inner midplane. rzimp(1,1:2) holds the R-coordinates of the first and second point of the segment. rzimp(2,1:2) holds the Z-coordinates of the first and second point of the segment.
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.
Switches to turn on/off the tracing files controlled by the ank_tracing switch. The files are defined in b2mod_diag, in order, starting with element 2 of the filedata array: test.trc, residuals.trc, sources.trc, blnn.trc, blne.trc, integral.trc, user.trc, blnm.trc, sepdata.trc. If there is no b2.user.parameters file present, the flags for user.trc are set to .false.. If the tracing files are to be appended but a reading error occurs when opening them, the corresponding filedata element is overwritten to .false..
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.