Encyclopaedia Index
Reserved Names of Patches or Variables
Contents
- Introduction
- Patch names
- Variable names
- Fortran variables used for properties
1. INTRODUCTION
The use of special names for patches or variables is a powerful and
flexible technique which users can employ for defining and controlling
the flow simulations performed by PHOENICS.
Many such names are employed in the versions of PHOENICS which are delivered
to users. They are here called 'reserved names'; and users are warned that,
when introducing new variables for their own use, they should always
avoid using those which are reserved. In particular, although they may
rename any of the first 17 reserved variables, they should not seek to give them
any non-standard physical significance. For example, to re-name H1
(first-phase sepcific enthalpy) as EPOT (electric potential) could have
undesired consequences.
This Encyclopaedia article provides lists, and brief explanations of the
names which have been reserved so far.
Numerous new reserved names of variables have been introduced into PHOENICS
in recent years, especially in order to facilitate:
The recently-introduced features are listed in the 'included' COMMON file
lbnamer where explanatory comments
are to be found.
Click here for reserved
variable names
or
here for reserved patch names
2. PATCH NAMES
Contents
First character
. as the first character. see PHENC entry: DOT-PATCH
& as the first character. see PHENC entry: Group 12
* as the first character. see PHENC entry: STAR-NAME
$ as the first character. see PHENC entry: DOLLAR-NAME
> as the first character. see PHENC entry: GT-NAME
+ as the first character. see PHENC entry: LINK
- as the first character. see PHENC entry: none
+ as the second character. see PHENC entry: LINK
- as the second character. see PHENC entry: LINEAR LINKS
Second character
% as the second character. followed by:-
SE as the 2nd and 3rd characters. See PHENC entry: MASK
SW as the 2nd and 3rd characters. See PHENC entry: MASK
SH as the 2nd and 3rd characters. See PHENC entry: MASK
SL as the 2nd and 3rd characters. See PHENC entry: MASK
WN as the 2nd and 3rd characters. See PHENC entry: MASK
WH as the 2nd and 3rd characters. See PHENC entry: MASK
WL as the 2nd and 3rd characters. See PHENC entry: MASK
EX as the 2nd and 3rd characters. See PHENC entry: MASK
EY as the 2nd and 3rd characters. See PHENC entry: MASK
EZ as the 2nd and 3rd characters. See PHENC entry: MASK
DX as the 2nd and 3rd characters. See PHENC entry: MASK
PX as the 2nd and 3rd characters. See PHENC entry: MASK
DY as the 2nd and 3rd characters. See PHENC entry: MASK
OY as the 2nd and 3rd characters. See PHENC entry: MASK
DZ as the 2nd and 3rd characters. See PHENC entry: MASK
OZ as the 2nd and 3rd characters. See PHENC entry: MASK
NE as the 2nd and 3rd characters
RAD as the 2nd, 3rd and 4th characters,
see the PHENC entry: STAR-NAME
POL as the 2nd, 3rd and 4th characters,
TIM as the 2nd, 3rd and 4th characters, see
First four characters
BUOY as the 1st, 2nd, 3rd and 4th characters,
BLIN as the 1st, 2nd, 3rd and 4th characters,
CHSO as the 1st, 2nd, 3rd and 4th characters,
CLDA as the 1st, 2nd, 3rd and 4th characters,
GP12 as the 1st, 2nd, 3rd and 4th characters, see GROUP 12 features
IBFC as the 1st, 2nd, 3rd and 4th characters,
IMBL as the 1st, 2nd, 3rd and 4th characters, see IMBAL
IPST as the 1st, 2nd, 3rd and 4th characters,
KEDI as the 1st, 2nd, 3rd and 4th characters,
KESO as the 1st, 2nd, 3rd and 4th characters,
LATG as the 1st, 2nd, 3rd and 4th characters,
LESO as the 1st, 2nd, 3rd and 4th characters,
OUTL as the 1st, 2nd, 3rd and 4th characters,
PROF as the 1st, 2nd, 3rd and 4th characters,
RADI as the 1st, 2nd, 3rd and 4th characters,
ROTA as the 1st, 2nd, 3rd and 4th characters,
SHSO as the 1st, 2nd, 3rd and 4th characters,
NOCP as the 1st, 2nd, 3rd and 4th characters,
First five characters
INIPOL as the 1st 5 characters, see
INIPOL
3. RESERVED NAMES OF SOLVED & STORED VARIABLES
3.1 The first 17 variables: P1......C2
These variables are:
number | name | meaning | typical units |
1 | P1 | first-phase pressure | N/m2 |
2 | P2 | second-phase pressure | N/m2 |
3 | U1 | first-phase x-direction velocity | m/s |
4 | U2 | second-phase x-direction velocity | m/s |
5 | V1 | first-phase y-direction velocity | m/s |
6 | V2 | second-phase y-direction velocity | m/s |
7 | W1 | first-phase z-direction velocity | m/s |
8 | W2 | second-phase z-direction velocity | m/s |
9 | R1 | first-phase volume fraction | none |
10 | R2 | second-phase volume fraction | none |
11 | RS | second-phase shadow volume fraction | none |
12 | KE | kinetic energy of turbulencs | J/kg |
13 | EP | dissipation rate ofturbulence energy | J/(kg s) |
14 | H1 | first-phase epecific enthalpy | J/kg |
15 | H2 | second-phase epecific enthalpy | J/kg |
16 | C1 | first-phase component concentration | none |
17 | C2 | first-phase component concentration | none |
3.2 Variables which may be placed in any order
When the command: STORE(reserved name)
appears in a Q1 file, the PHOENICS solver module, EARTH, computes, and makes
available for printing or graphical display, the corresponding quantity,
as explained below.
It should however be understood that the STORE statement may not be
sufficient. Thus:
- STORE(AREE) will always produce the east areas, because
they always are computed, whereas
- STORE(CFIP) will produce non-default values only if other
commands (viz ONEPHS = F, etc) have been supplied to ensure that
the inter-phase friction coefficient is actually being computed.
The discussion deals with complete names first and partial names second.
Complete names
the name : the corresponding quantity
-
- APRJ : projected area of particles per unit
volume in an inter-phase-friction model
- AREE : east cell-face area
- AREN : north cell-face area
- AREH : high cell-face area
- ARRT : flame-arrival time
- BLOK
-
BTAU : The magnitude of the shear-stress tensor for the Non-Newtonian fluid.
-
CD : non-dimensional drag coefficient in an inter-phase-friction model
- CFIP : inter-phase friction coefficient
- CHAR : char in wood combustion model
- CMDO : mass-transfer rate from phase 2 to phase 2 for the cell
- CMI1 : mass of phase 1 entering the cell from outside the domain MIN1, divided by
the pressure difference which drives it
- CMI2 : mass of phase 2 entering the cell from outside the domain MIN2, divided by
the pressure difference which drives it
- CNE1 : outward convection flux (kg/s) of phase 1 through east cell face
- CNE2 : outward convection flux (kg/s) of phase 2 through east cell face
- CNN1 : outward convection flux (kg/s) of phase 1 through north cell face
- CNN2 : outward convection flux (kg/s) of phase 2 through north cell face
- CNH1 : outward convection flux (kg/s) of phase 1 through high cell face
- CNH2 : outward convection flux (kg/s) of phase 2 through high cell face
- CONI : as for CNE1 for (possibly moving) body-fitted coordinates
- CONJ : as for CNN1 for (possibly moving) body-fitted coordinates
- CONK : as for CNH1 for (possibly moving) body-fitted coordinates
- CP1 : phase-1 specific-heat at constant pressure
- CP2 : phase-2 specific-heat at constant pressure
- CREY : local cell Reynolds number = VABS*VOL1./3/(ENUT+ENUL)
- CTI1 : inward convection flux (kg/s) of phase 1 through 'time face'
- CTI2 : inward convection flux (kg/s) of phase 2 through 'time face'
- CTO1 : outward convection flux (kg/s) of phase 1 through 'time face'
- CTO2 : outward convection flux (kg/s) of phase 2 through 'time face'
- CUTC : indicator of whether a cell has 'cut-cell' status
- DEN1 : phase-1 density
- DEN2 : : phase-2 density
- DIAM : diameter of particle in various two-phase models
- DP1X : X-direction pressure correction difference
- DP1Y : Y-direction pressure correction difference
- DP1Z : Z-direction pressure correction difference
- DRH1 : rate of change of log(phase-1 density) with pressure
- DRH2 : rate of change of log(phase-1 density) with pressure
- DTDX : X-direction temperature gradient (TEM1 if TEM1 is SOLVEd, TEMP1 if TEMP1 is STOREd)
- DTDY : Y-direction temperature gradient (TEM1 if TEM1 is SOLVEd, TEMP1 if TEMP1 is STOREd)
- DTDZ : Z-direction temperature gradient (TEM1 if TEM1 is SOLVEd, TEMP1 if TEMP1 is STOREd)
- DTX2 : X-direction temperature gradient (TEM2 if TEM2 is SOLVEd, TEMP2 if TEMP2 is STOREd)
- DTY2 : Y-direction temperature gradient (TEM2 if TEM2 is SOLVEd, TEMP2 if TEMP2 is STOREd)
- DTZ2 : Z-direction temperature gradient (TEM2 if TEM2 is SOLVEd, TEMP2 if TEMP2 is STOREd)
- DU1P : rate of change of U1 with pressure difference
- DV1P : rate of change of V1 with pressure difference
- DVO1 : rate of change of log(phase-1 specific volume) with temperature
- DVO2 : rate of change of log(phase-2 specific volume) with temperature
- DW1P : rate of change of W1 with pressure difference
- DU2P: rate of change of U2 with pressure difference
- DV2P: rate of change of V2 with pressure difference
- DW2P: rate of change of W2 with pressure difference
- DUDX: 1st Phase velocity gradient = ∂U1/∂x
- DUDY: 1st Phase velocity gradient = ∂U1/∂y
- DUDZ: 1st Phase velocity gradient = ∂U1/∂z
- DVDX: 1st Phase velocity gradient = ∂V1/∂x
- DVDY: 1st Phase velocity gradient = ∂V1/∂y
- DVDZ: 1st Phase velocity gradient = ∂V1/∂z
- DWDX: 1st Phase velocity gradient = ∂W1/∂x
- DWDY: 1st Phase velocity gradient = ∂W1/∂y
- DWDZ: 1st Phase velocity gradient = ∂W1/∂z
- DU2X: 2nd Phase velocity gradient = ∂U2/∂x
- DU2Y: 2nd Phase velocity gradient = ∂U2/∂y
- DU2Z: 2nd Phase velocity gradient = ∂U2/∂z
- DV2X: 2nd Phase velocity gradient = ∂V2/∂x
- DV2Y: 2nd Phase velocity gradient = ∂V2/∂y
- DV2Z: 2nd Phase velocity gradient = ∂V2/∂z
- DW2X: 2nd Phase velocity gradient = ∂W2/∂x
- DW2Y: 2nd Phase velocity gradient = ∂W2/∂y
- DW2Z: 2nd Phase velocity gradient = ∂W2/∂z
- EFEH: Parsol edge fraction
- EFNE: Parsol edge fraction
- EFNH: Parsol edge fraction
- EL1 : mixing-length of turbulence, for phase 1
- EL2 : mixing-length of turbulence, for phase 2
- EMIS: emissivity (also absorptivity) of radiant energy per unit volume
- ENUL: laminar kinematic viscosity
- ENUT : turbulent contribution to the
effective kinematic viscosity
- EOTV: Eotvos number used in Ellipsoidal-Bubble "Clean-water" Drag correlation
- EPKE: dissipation rate of turbulence energy per unit volume, epsilon, divided by turbulence energy, k
- EPOT: potential - no convection or transient terms, no turbulent diffusion
- FONE: Damping function in Lam-Bremhorst
Low-Reynolds Number turbulence model
- FSQ : root-mean square of composition variable F in
Presumed-PDF combustion model.
- FTWO: as FONE
- FMU : as FONE
- FOMG: multiplier FΩ in
MMK turbulence model
- FUEL : mass fraction of unburned fuel in a
simple chemically-reacting system (SCRS)
- FWD : Mass fraction Wood derivative
- FC : Carbon content for wood combustion model
- FH : Hydrogen content for wood combustion model
- FO : Oxygen content for wood combustion model
- FN : Nitrogen content for wood combustion model
- GEN1: sum of squares of phase-1 velocity gradients where (GEN1)1/2=√(2.*Sij*Sij)
where Sij= 0.5*(∂Ui/∂xj + ∂Uj/∂xi)
- GEN2: sum of squares of phase-2 velocity gradients
- GENK: effective viscosity times GEN1
- GNK2: effective viscosity times GEN2
- H0_1: enthalpy of phase 1 at zero temperature
- H0_2: enthalpy of phase 2 at zero temperature
- HTCO: heat-transfer coefficient derived from wall function
- HFLX: heat flux
- IMB1: continuity imbalances of phase 1
- IMB2: continuity imbalances of phase 2
- INTF: as for CFIP
- INTM: as for CMDO
- ISVR : for debugging cut cell masking
- No variables starting 'J' at present.
-
KOND: the thermal conductivity (for phase 1)
- KND2: the thermal conductivity for phase 2
- LEN1: as for EL1
- LEN2: as for EL2
- LIMB: similar to object-identifier, OBID; used for MOFOR
- LTLS: variable solved in order to deduce
values of WDIS and WGAP
- MACH: MACH Number of phase 1 fluid
- MACZ: MACH Number based on Z velocity used for parabolic option
IPARAB=5.
- MAC2: Mach Number of phase 2 fluid
- MARK: variable used in InForm to 'mark' cells
- MAS1: mass of phase 1 in the cell
- MAS2: mass of phase 2 in the cell
- MDOT: as for CMDO
- MIXF: 'mixture fraction', i.e. mass fraction
of (burned or unburned) material emanating from the fuel-bearing stream in a
Simple Chemically-Reacting System.
- MIN1: mass of phase 1 entering the cell from outside the domain
- MIN2: mass of phase 2 entering the cell from outside the domain
- MUL: dynamic laminar viscosity (in Pa.s), i.e. μl=ρ.νl
- MXF1 - MXF7: mass fraction of material from inlet 1 - 7 in wood combustion model (7 gasses)
- OBID: object identifier, used for MOFOR
- OXID: oxidant mass fraction when a Simple Chemically Reacting System is being modelled
- PDCX: cell-to-cell pressure decrement in x-direction
- PDCY: cell-to-cell pressure decrement in y-direction
- PDCZ: cell-to-cell pressure decrement in z-direction
- PHDE: phase 2 diffusion flux in X direction
- PHDN: phase 2 diffusion flux in Y direction
- PHDH: phase 2 diffusion flux in Z direction
- POT : velocity potential
- PROD: product mass fraction when a Simple Chemically Reacting System is being modelled
- PRPS: material marker
- PRL : Prandtl Number for variables with
PRNDTL(phi)=GRND1
- PSOX: x-direction momentum source resulting from the pressure gradient
- PSOY: y-direction momentum source resulting from the pressure gradient
- PSOZ: z-direction momentum source resulting from the pressure gradient
- PTOT: total pressure
(P1+0.5*DEN1*VABS^2) based on phase 1 velocity. Requires MACH to be STOREd
- PTO2 : total pressure based on phase 2 velocity. Requires MAC2 to be STOREd
- QCR1: Q-criterion where Q=0.5*(Ω2-S2) where Ω is the
vorticity magnitude and S the strain rate magnitude
- QCR2: Q-criterion for phase 2.
- QDX : x-direction diffusive/conductive heat flux derived from gradients of TEM1 or H1
- QDY : y-direction diffusive/conductive heat flux derived from gradients of TEM1 or H1
- QDZ : z-direction diffusive/conductive heat flux derived from gradients of TEM1 or H1
- QRX : x-direction radiation flux computed by
immersol
- QRY : y-direction radiation flux computed by
immersol
- QRZ : z-direction radiation flux computed by
immersol
- RADX: net X-direction radiative heat flux for
6-Flux model
- RADY: net Y-direction radiative heat flux for 6-Flux model
- RADZ: net Z-direction radiative heat flux for 6-Flux model
- REYD: Reynolds Number based on interphase slip velocity used in interphase drag laws
- REYN: Reynolds Number in Lam-Bremhorst Low-
Reynolds Number turbulence model
- REYT: Reynolds Number in Lam-Bremhorst Low-
Reynolds Number turbulence model
- RHO1: as for DEN1
- RHO2: as for DEN2
- SCAT: scattering coefficient of radiant energy per unit volume
- SHRX: X-direction friction force from wall functions
- SHRY: Y-direction friction force from wall functions
- SHRZ: Z-direction friction force from wall functions
- SIZE: particle/droplet diameter in twophase flow
- SKIN: skin friction coefficient
- SLPU: interphase slip velocity in X direction
(U1-U2)
- SLPV: interphase slip velocity in Y direction
(V1-V2)
- SLPW: interphase slip velocity in Z direction
(W1-W2)
- SPH1: as for CP1
- SPH2: as for CP2
- STAN: Stanton Number, i.e. dimensionless heat-transfer coefficient
- STRS: Wall shear stress divided by the local density
- SRM1: Strain rate magnitude S=(GEN1)1/2=√(2.*Sij*Sij)where
Sij= 0.5*(∂Ui/∂xj + ∂Uj/∂xi)
- SRM2: Strain rate magnitude of 2nd phase S=(GEN2)1/2.
- SURN: current-time volume fraction of liquid in
free-surface problems, solved with VOF,
SEM or HOL
- SRN2: current-time volume fraction of second liquid in
3-phase free-surface problems, solved with VOF
- TFNN: temperature function fT used in
Non-Newtonian viscosity function.
- TREY: local turbulent Reynolds number, = ENUT/ENUL
- TWAL: wall surface temperature associated with the THINPLT object
- T3 : radiation temperature, used in IMMERSOL radiation model
- TEMP: phase-1 temperature, when enthalpy is computed
directly
- TEM1: phase-1 temperature, when computed directly
- TEM2: phase-2 temperature, when computed directly
- TMP1: phase-1 temperature, when enthalpy is
computed directly
- U1AD: the addition to U1 resulting from pressure
correction
- U1SL: U1 after solution but before pressure correction
- U2AD: the addition to U2 resulting from pressure correction
- U2SL: U2 after solution but before pressure correction
- U2CM: x-direction phase 2 velocity component for BFC
- U2CR: x-direction cartesian phase 2 velocity component for BFC
- UCMP: x-direction velocity component for BFC
- UCRT: x-direction cartesian velocity component for BFC
- V1AD: the addition to V1 resulting from
pressure correction
- V1SL: V1 after solution but before pressure correction
- V2AD: the addition to V2 resulting from pressure correction
- V2CM: y-direction phase 2 velocity component for BFC
- V2CR: y-direction cartesian phase 2velocity component for BFC
- V2SL: V2 after solution but before pressure correction
- VABS: Absolute velocity of phase 1
- VAB2: Absolute velocity of phase 2
- VCMP: y-direction velocity component for BFC
- VCRT: y-direction cartesian
velocity component for BFC
- VFOL: old-time volume fraction of liquid in
free-surface problems, solved with VOF,
SEM or HOL
- VFL2: old-time volume fraction of second liquid in
3-phase free-surface problems, solved with VOF
- VISL: as for ENUL
- VIST: as for ENUT
- VLSQ: sum of squares of phase 1 velocities, i.e. twice the kinetic energy
- VOLU: cell volume
- VOR1: Magnitude of the vorticity vector Ω=√(2.*ωijωij)1/2 where
ωij= 0.5*(∂U1,i/∂xj - ∂U1,j/∂xi)
- VOR2: Magnitude of the 2nd-phase vorticity vector.
- VOX1: vorticity component in YZ plane = ∂W1/∂y - ∂V1/∂z
- VOY1: vorticity component in XZ plane = ∂U1/∂z - ∂W1/∂x
- VOZ1: vorticity component in XY plane = ∂V1/∂x - ∂U1/∂y
- VOX2: phase-2 vorticity component in YZ plane = ∂W2/∂y - ∂V2/∂z
- VOY2: phase-2 vorticity component in XZ plane = ∂U2/∂z - ∂W2/∂x
- VOZ2: phase-2 vorticity component in XY plane = ∂V2/∂x - ∂U2/∂y
- VOXY: vorticity in XY plane based on circulation theory
- VOYZ: vorticity in YZ plane based on circulation theory
- VOZX: vorticity in ZX plane based on circulation theory
- VREL: relative absolute velocity between the
two phases when ONEPHS=F
- W1AD: the addition to W1 resulting from
pressure correction
- W1SL: W1 after solution but before pressure correction
- W2AD: the addition to absolute velocity resulting from pressure correction
- W2CM: z-direction phase 2 velocity component for BFC
- W2CR: z-direction cartesian phase 2velocity component for BFC
- W2SL: W2 after solution but before pressure correction
- WAVE: phase-average z-direction velocity
- WCMP: z-direction velocity component for BFC
- WCRT: z-direction cartesian velocity component for BFC:
- WDIS: distance from the wall, calculated by the
LTLS method
- WEB : Weber number in interphase drag law for CFIPD=4 or 5
- WGAP: distance between walls, calculated by the
LTLS method
- XCEN: x-direction cartesian coordinate of
cell centre in a BFC grid
- YCEN: y-direction cartesian coordinate of cell
centre in a BFC grid
- YPLS: 'y plus', the non-dimensional wall distance
- ZCEN: z-direction cartesian coordinate of cell
centre in a BFC grid
hash-name variables
For a reason explained under
'partial names' below, many reserved names have # as the first character.
The following allow the
'status' of the cells
in the domain to be recorded, for
printing or graphical display, in 3D-variable manner.
They correspond to the names of logical variables which users who engage in
the addition of
Fortran programming to
PHOENICS can exploit.
- #EBO, indicating that the East boundary is
either a domain boundary or adjoins a VAC cell
or a zero face porosity
- #EF , indicating that the East face is a phase boundary
- #EF0 which East face is blocked by zero face porosity;
- #EWP, indicating that the East wall type patch is present
- #HBO, indicating that the High boundary is
either a domain boundary or adjoins a VAC cell
or a zero face porosity
- #HF , indicating that the High face is a phase boundary
- #HF0, indicating that the High face is blocked by zero face porosity;
- #HWP, indicating that the High wall type patch is present
- #LBO, indicating that the Low boundary is
either a domain boundary or adjoins a VAC cell
or a zero face porosity
- #LF , indicating that the Low face is a phase boundary
- #LF0, indicating that the Low face is blocked by zero face porosity;
- #LWP, indicating that the Low wall type patch is present
- #MSFP, indicating that the that there is a mass-flow patch present in a cell;
- #NBO, indicating that the Low boundary is
either a domain boundary or adjoins a VAC cell
or a zero face porosity
- #NF, indicating that the North face is a phase boundary
- #NF0, indicating that the North face is blocked by zero face porosity;
- #NWP, indicating that the North wall type patch is present
- #POR, indicating that the there is a non-participating solid with friction;
- #SBO, indicating that the North boundary is
either a domain boundary or adjoins a VAC cell
or a zero face porosity
- #SF, indicating that the South face is a phase boundary
- #SF0, indicating that the South face is blocked by zero face porosity;
- #SLD which SLD means the cell is occupied by a solid material;
- #SWP, indicating that the South wall type patch is present
- #WBO, indicating that the West boundary is
either a domain boundary or adjoins a VAC cell
or a zero face porosity
- #WF0, indicating that the West boundary is
either a domain boundary or adjoins a VAC cell
or a zero face porosity
- #WF, indicating that the West face is blocked by zero face porosity;
- #WWP, indicating that the West wall type patch is present
- #VAC, indicating that the cell is empty, or 'vacant', and therefore
plays no part in the flow simulation.
Partial names.
Additionally, significance is attached to certain characters when they
appear in particular positions in the names of variables, as follows:
- # as the first character of the name will prevent
EARTH from automatically setting its value to
zero within cells occupied by solids, as is done
for all variables except:
R1, R2, RS, H1, H2, TEM1, TEM2, T3,EPOR, NPOR, HPOR, BLOK, and PRPS.
- % as the fourth character allows the
RESIDUALS of a
solved-for variable to be stored and printed. However the
In-Form 'RESI'
feature does the same job more neatly.
- C as the fourth character allows the values of the 'corrections'
of a solved-for variable to be stored and printed. However the
In-Form 'CORR'
feature does the same job more neatly.
4. FORTRAN VARIABLES USED FOR PROPERTIES
PHOENICS users who make use of the Ground-coding facility should also be
aware of the integer variables used for accessing material properties, namely those contained in COMMON/IPROPS/ to be found in the file prpcmn in D_INCLUD.
They are:
COMMON/IPROPS/DENST1,CMPRS1,DENST2,CMPRS2,VISTRB,VISCLM,PRNLAM,
1 INTVL1,INTVL2,TEMPR1,TEMPR2,MIXLN1,MIXLN2,SPEHT1,
1 SPEHT2,THRME1,THRME2,EMISSV,SCATTR,NAMSP3,INTFCO,
1 INTMTR,INTRC1,INTRC2,IVRMCO,NAMSP4(5)
The meanings of all of these quantities are best understood by reading comments at the top of the Fortran file gxprutil.for ,
as follows
It should be noted that each has a prescribed value, which should on
no account be changed.