___________________________________________________________________ AEAST.....name used to denote an array containing the 'free' east-face area after blockages are deducted. ___________________________________________________________________

AHIGH.....name used to denote an array containing the 'free' high-face area after blockages are deducted.

___________________________________________________________________

ANORTH....name used to denote an array containing the 'free' north-face area after blockages are deducted.

___________________________________________________________________

AZDZ....determines the variation of DZ with Z distance. It must be set to GRND as a part of block for a setting the formulae for direct setting DZ.

____________________________________________________________________

AZXU....determines the variation of XULAST with Z distance. It must be set to GRND as a part of block for a setting the formulae to calculate the XRAT as the factor which will multiply XULAST.

____________________________________________________________________

AZYV....determines the variation of YVLAST with Z distance. It must be set to GRND as a part of block for a setting the formulae to calculate the YRAT as the factor which will multiply YVLAST.

____________________________________________________________________

BOX....function used as part of PLANT statement.

The syntax is :

MARK=BOX(MRK,XPP,YPP,ZPP,XSIZE,YSIZE,ZSIZE, ANGXX,ANGXY,ANGXZ,ANGYX,ANGYY,ANGYZ,ANGZX,ANGZY,ANGZZ, ROTANG)

It sets the MARK values equal to MRK inside the surface of three dimensional box placed on the grid.

The function has 17 arguments, namely (in order):

MRK - MARK value which is to be associated with the object;

XPP - X - coordinate of west south low vertice of the box;

YPP - Y - coordinate of west south low vertice of the box;

ZPP - Z - coordinate of west south low vertice of the box;

XSIZE - box size in X-direction;

YSIZE - box size in Y-direction;

ZSIZE - box size in Z-direction;

ANGXX - angle between X-axis of object coordinate system and

X-axis of reference cartesian frame, rad;

ANGXY - angle between X-axis of object coordinate system and

Y-axis of reference cartesian frame, rad;

ANGXZ - angle between X-axis of object coordinate system and

Z-axis of reference cartesian frame, rad;

ANGYX - angle between Y-axis of object coordinate system and

X-axis of reference cartesian frame, rad;

ANGYY - angle between Y-axis of object coordinate system and

Y-axis of reference cartesian frame, rad;

ANGYZ - angle between Y-axis of object coordinate system and

Z-axis of reference cartesian frame, rad;

ANGZX - angle between Z-axis of object coordinate system and

X-axis of reference cartesian frame, rad;

ANGZY - angle between Z-axis of object coordinate system and

Y-axis of reference cartesian frame, rad;

ANGZZ - angle between Z-axis of object coordinate system and

Y-axis of reference cartesian frame, rad;

ROTANG- angle of shifting reference cartesian frame in the

X-direction ( significant for CARTES=F), rad.

___________________________________________________________________

CFIPS .... is the reference interphase-transport coefficient.

It must be set to GRND as a part of block for a setting the formulae to calculate the filed values at each Z-slab.

____________________________________________________________________

CINT(PHI) .... is the ratio of the bulk-to-interface transfer coefficient, for the variable indicated.

It must be set to GRND as a part of block for a setting the formulae to calculate the filed values at each Z-slab.

____________________________________________________________________

CMDOT .... is the interphase mass transfer parameter.

It must be set to GRND as a part of block for a setting the formulae to calculate the filed values at each Z-slab.

____________________________________________________________________

CO....name used to denote the 'coefficient' array, and which is used to set the coefficients of linearised-source expressions.

___________________________________________________________________

COI1(PHI)....name used to denote the inter-phase transfer coefficient from phase 1 to the interface in response to the command, CINT(PHI)=GRND .

___________________________________________________________________

COI2(PHI)....name used to denote the inter-phase transfer coefficient from phase 2 to the interface in response to the command, CINT(PHI)=GRND .

___________________________________________________________________

CP1 .... is the heat capacity at the constant pressure of the phase-1.

____________________________________________________________________

CP2 .... is the heat capacity at the constant pressure of the phase-2.

____________________________________________________________________

C1, C3, C5, C7, C9, C11, C13, C15, C17, C19, C21, C23, C25, C27, C29, C31, C33, C35....standard names used to denote first-phase concentration full-field variables.

___________________________________________________________________

C2, C4, C6, C8, C10, C12, C14, C16, C18, C20, C22, C24, C26, C28, C30, C32, C34....standard names used to denote second-phase concentration full-field variables.

___________________________________________________________________

D1DP....name used to denote the values for the derivative of the natural log of first-phase density with respect to pressure in response to the command, DRH1DP=GRND .

___________________________________________________________________

D2DP....name used to denote the values for the derivative of the natural log of second-phase density with respect to pressure in response to the command, DRH2DP=GRND .

___________________________________________________________________

DEN1....name used to denote the values for the densities of the first phase in response of command, RHO1=GRND

___________________________________________________________________

DEN2....name used to denote the values for the densities of the second phase in response of command, RHO2=GRND

___________________________________________________________________

DISPLC .... function used to add the displacements of solids to cell corner coordinates ( BFC=T) so that the distortion of the solid block is computed. Then it can be visualised by PHOTON to get the image of the possible solid shapes.

The syntax is

XC=DISPLC(scale1) YC=DISPLC(scale2) ZC=DISPLC(scale3)

It has one argument - the scale factor to magnify the contribution of displacements to corresponding cell corner coordinate. By default, all scales are set to unity.

=> Example:

{MXYZ30> XC=DISPLC(0.25) {MXYZ30> YC=DISPLC(RG(1)) {MXYZ30> ZC=DISPLC(0.25*RG(2)) IF(ISWEEP.EQ.LSWEEP.AND.STORE(W1)) REGION(2,5,4,8,1,NZ)

___________________________________________________________________

DISPLS....function used to calculate the location and shape of solid object affected by displacements on XY plane of cartesian (CARTES=T) grids.

It does so by re-calculation of resulting locations for each solid cell marked by certain PRPS value. Then PHOTON can be used to display the new solid shape.

The syntax is

PRPS=DISPLS(PRPM, XSCALE, YSCALE )

It has three arguments, namely (in order):

PRPM - PRPS value which is to be associated with solid cells;

XSCALE - the scale factor to magnify the contribution of

displacements to the shifting solid cells in

X-direction.

YSCALE - the scale factor to magnify the contribution of

displacements to the shifting solid cells in

Y-direction.

By default, all scales are set to unity.

For example, the statement

{SC0601> PRPS=DISPLS(111., SCAL[1,1,1],RG(1)*SCAL[1,3.1])

provides the calculation of displacement affected steel (PRPS=111.) cell movements at the end of each z-slab with scale factors taken from field distribution of SCAL variable.

___________________________________________________________________

DRH1DP .... is the compressibility of the phase-1.

___________________________________________________________________

DRH2DP .... is the compressibility of the phase-2.

____________________________________________________________________

DT....is the name used to denote the values for the current time-step sizein response of command, TLAST=GRND

___________________________________________________________________

DTFALS(PHI)....name used to denote value of false-time-step of under-relaxation for the variable indicated, and which is used to set the specific expressions.

___________________________________________________________________

DXG2D....name used to denote the values of the distances in the IX-direction between centres of continuity cells.

___________________________________________________________________

DXU2D....name used to denote the values of the distances in the IX-direction between walls of continuity cells, where are stored the u velocity resolutes.

___________________________________________________________________

DYG2D....name used to denote values of the distances in the IY-direction between centres of continuity cells.

___________________________________________________________________

DYV2D....name used to denote the values of the distances in the IY-direction between walls of continuity cells, where are stored the v velocity resolutes.

___________________________________________________________________

DZ....name used to denote the the current width of the slab, and the forward-step size in parabolic calculations in response of the command, AZDZ=GRND.

___________________________________________________________________

DZG....is the distance between the cell centres of the current-slab cell and the high slab-cell.

___________________________________________________________________

DZGL....is the DZG of the low-slab.

___________________________________________________________________

DZGNZ....name used to denote values of the distances in the IZ-direction between centres of continuity cells.

___________________________________________________________________

DZL....is the z-thickness of the previous slab.

___________________________________________________________________

DZRAT....is set in EARTH to DZ / DZL.

___________________________________________________________________

DZWNZ....name used to denote values of the distances in the IZ-direction between cell faces.

___________________________________________________________________

EAST(PHI)....name used to denote the value of PHI on the east side of the cell in question.

___________________________________________________________________

EL1 .... is the mixing length scale of the first-phase fluid.

___________________________________________________________________

EL2 .... is the mixing length scale of the second-phase fluid.

___________________________________________________________________

ELLPSD...function used as part of PLANT statement.

The syntax is:

MARK=ELLPSD(MRK,XPP,YPP,ZPP,XAX,YAX,ZAX, ANGXX,ANGXY,ANGXZ,ANGYX,ANGYY,ANGYZ,ANGZX,ANGZY,ANGZZ, ROTANG)

It sets the MARK values equal to MRK inside the surface of three dimensional ellipsoid placed on grid.

The function has 16 arguments, namely (in order):

MRK - MARK value which is to be associated with the object;

XPP - X - coordinate of ellipsoid centre;

YPP - Y - coordinate of ellipsoid centre;

ZPP - Z - coordinate of ellipsoid centre;

XAX - half axis size in X-direction;

YAX - half axis size in Y-direction;

ZAX - half axis size in Z-direction;

ANGXX - angle between X-axis of object coordinate system and

X-axis of reference cartesian frame, rad;

ANGXY - angle between X-axis of object coordinate system and

Y-axis of reference cartesian frame, rad;

ANGXZ - angle between X-axis of object coordinate system and

Z-axis of reference cartesian frame, rad;

ANGYX - angle between Y-axis of object coordinate system and

X-axis of reference cartesian frame, rad;

ANGYY - angle between Y-axis of object coordinate system and

Y-axis of reference cartesian frame, rad;

ANGYZ - angle between Y-axis of object coordinate system and

Z-axis of reference cartesian frame, rad;

ANGZX - angle between Z-axis of object coordinate system and

X-axis of reference cartesian frame, rad;

ANGZY - angle between Z-axis of object coordinate system and

Y-axis of reference cartesian frame, rad;

ANGZZ - angle between Z-axis of object coordinate system and

Y-axis of reference cartesian frame, rad;

ROTANG- angle of shifting reference cartesian frame in the
X-direction ( significant for CARTES=F), rad.

___________________________________________________________________

ENUL .... is the laminar kinematic viscosity.

___________________________________________________________________

ENUT .... is the turbulent kinematic viscosity.

___________________________________________________________________

FII1(PHI)....name used to denote the values of the current phase-1 variable at the phase interface in response to the command, PHINT(PHI)=GRND.

___________________________________________________________________

FII2(PHI)....name used to denote the values of the current phase-2 variable at the phase interface in response to the command, PHINT(PHI)=GRND.

___________________________________________________________________

GREAT....A large number, used in EARTH. Do not re-set. The default value may be re-set during installation to suit the precision of the computer.

___________________________________________________________________

H1....standard name used to denote first-phase enthalpy. The H1 term may be conveniently used to denote first-phase temperature, in which case the term may be reNAMEd (eg to 'TEMP').

___________________________________________________________________

H2....standard name used to denote second-phase enthalpy.

___________________________________________________________________

HIGH(PHI)....name used to denote the value of PHI on the high side of the cell in question.

___________________________________________________________________

INDVAR....is the current dependent- or auxiliary- variable number.

___________________________________________________________________

INTFRC....name used to denote interphase friction coefficient in response of the command, CFIPS=GRND.

___________________________________________________________________

INTMDT....name used to denote interphase mass-transfer rate in response of the command, CMDOT=GRND.

___________________________________________________________________

ISPH1....name used to denote specific heat capacity at constant pressure of the first phase in response of the command, CP1=GRND.

___________________________________________________________________

ISPH2....name used to denote specific heat capacity at constant pressure of the second phase in response of the command, CP2=GRND.

___________________________________________________________________

ISTEP....is the time-step number.

___________________________________________________________________

ISWEEP....is the current sweep number.

___________________________________________________________________

IX....ranges from 1 to NX, denoting cell location in X-direction of PHOENICS.

___________________________________________________________________

IXF....is the west extent of the current patch within the slab.

___________________________________________________________________

IXL....is the east extent of the current patch within the slab.

___________________________________________________________________

IY....ranges from 1 to NY, denoting grid-point location in y-direction of PHOENICS, but not necessarily in cartesian-y direction for BFC=T. In the latter IY=J-1 (see BODY-F for indicial nomenclature in use in body-fitted coordinates).

___________________________________________________________________

IYF....is the south extent of the current patch within the slab.

___________________________________________________________________

IYL....is the north extent of the current patch within the slab.

___________________________________________________________________

IZ....ranges from 1 to NZ, denoting grid-point location in z-direction of PHOENICS, but not necessarily in cartesian-z direction for BFC=T. In the later IZ=K-1 (see BODY-F for indicial nomenclature in use in body-fitted coordinates).

___________________________________________________________________

IZSTEP....is the current z-slab.

___________________________________________________________________

LAMPR(PHI)....name used to denote the of values of the PRANDTL number in response of command, PRNDTL(PHI)=GRND. If PRNDTL(phi) = -GRND this name is used instead to set the values for diffusivities.

___________________________________________________________________

LD12....name used to denote the of values of the virtual mass in response of command, CVM=GRND.

___________________________________________________________________

LEN1....name used to denote the values of the length scale of phase 1 in response of command, EL1=GRND.

___________________________________________________________________

LEN2....name used to denote the values of the length scale of phase 2 in response of command, EL2=GRND.

___________________________________________________________________

LOW(PHI)....name used to denote the value of PHI on the low side of the cell in question.

___________________________________________________________________

MASS1....name used to denote the values, pertaining to the mass of phase-1 in each cell (i.e the product of phase-1 volume fraction, phase-1 density, and the volume of the cell). If the flow is single phase (ONEPHS=T) then MASS1 refers to the mass of fluid in the cells.

___________________________________________________________________

MASS2....name used to denote the values, pertaining to the mass of phase-2 in each cell (i.e the product of phase-2 volume fraction, phase-2 density, and the volume of the cell). If the flow is single phase (ONEPHS=T) then MASS2 has no significance.

___________________________________________________________________

MXYZ?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 19, Section 2 of GROUND to set the cartesian coordinates of cell corners for BFC grids.

___________________________________________________________________

NORTH(PHI)....name used to denote the value of PHI on the north side of the cell in question.

___________________________________________________________________

NZSTEP....is the number of z-slabs into which the domain is sub-divided.

___________________________________________________________________

OLD(PHI)....name used to denote the value of PHI appropriate to the previous instant of time.

___________________________________________________________________

P1....standard name used to denote first-phase pressure. P1 is also used to specify the first-phase mass-flow.

___________________________________________________________________

P2....standard name used to denote second-phase mass-flow.

___________________________________________________________________

PHINT(PHI) .... is the interface value of the variable indicated.

___________________________________________________________________

PRNDTL(PHI) .... is the laminar Prandtl number of the variable indicated.

It must be set to GRND as a part of block for a setting the formulae to calculate the filed values at each Z-slab. <> ___________________________________________________________________

PRPT?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 9 or 10 of GROUND to set properties of the medium ( media) including inter-phase-transfer processes.

____________________________________________________________________

R1....standard name used to denote the first-phase volume fraction.

___________________________________________________________________

R2....standard name used to denote the second-phase volume fraction.

___________________________________________________________________

RESREF(PHI)....name used to denote reference value of residual for the variable indicated, and which is used to set the specific expressions.

___________________________________________________________________

RHO1 .... is the density of the first phase.

___________________________________________________________________

RHO2 .... is the density of the second phase.

____________________________________________________________________

RS....standard name used to denote the second-phase shadow volume fraction, i.e. that which would be present in the absence of inter-phase-mass-transfer.

___________________________________________________________________

RG2D....name used to denote the array of values of distances from the axis of symmetry to the cell centres. RG2D is only of significance in cylindrical-polar coordinates, ie. when CARTES=F and BFC=F.

___________________________________________________________________

RV2D....name used to denote the array of values of distances from the axis of symmetry to the north wall of the continuity cells, at which the v velocities are stored. RV2D is only of significance in cylindrical-polar coordinates, ie. when CARTES=F and BFC=F.

___________________________________________________________________

SC01??....statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Section 1 of GROUND Group 19. This section is called at the start of each time step.

____________________________________________________________________

SC02??....statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Section 2 of GROUND Group 19. This section is called at the start of each sweep. ____________________________________________________________________

SC03??....statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Section 3 of GROUND Group 19. This section is called at the start of each iz slab.

____________________________________________________________________

SC04??....statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Section 4 of GROUND Group 19. This section is called at the start of each iteration.

____________________________________________________________________

SC05??....statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Section 5 of GROUND Group 19. This section is called at the finish of iteration.

____________________________________________________________________

SC06??....statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Section 6 of GROUND Group 19. This section is called at the finish of each iz slab.

____________________________________________________________________

SC07??....statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Section 7 of GROUND Group 19. This section is called at the finish of each sweep.

____________________________________________________________________

SC08??....statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Section 8 of GROUND Group 19. This section is called at the finish of each time step.

____________________________________________________________________

SC09??....statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Section 9 of GROUND Group 19. This section is called at the start of solution sequence for a variable.

____________________________________________________________________

SC10??....statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Section 10 of GROUND Group 19. This section is called at the finish of solution sequence for a variable.

____________________________________________________________________

SCTS?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 2 of GROUND to set the time step, DT.

____________________________________________________________________

SCUF?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 8, section 1, of GROUND to set the formulae for the values to add to the field of U1.

____________________________________________________________________ SCVF?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 8, section 3, of GROUND to set the formulae for the values to add to the field of V1.

____________________________________________________________________

SCWF?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 8, section 5, of GROUND to set the formulae for the values to add to the field of W1.

____________________________________________________________________

SCUS?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 8, section 2, of GROUND to set the formulae for the values to add to the field of U2.

____________________________________________________________________

SCVS?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 8, section 4, of GROUND to set the formulae for the values to add to the field of V2.

____________________________________________________________________

SCWS?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 8, section 6, of GROUND to set the formulae for the values to add to the field of W2.

____________________________________________________________________

SCXS?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 3 of GROUND to set the variable XRAT.

____________________________________________________________________

SCYS?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 4 of GROUND to set the variable YRAT.

____________________________________________________________________

SCZS?? .... statement heading used to instruct PLANT that it must write its FORTRAN block of number ?? in Group 5 of GROUND to set the variable DZ.

____________________________________________________________________

SISBC...function used to introduce boundary conditions for the displacements at the solid fluid interface on non-BFC grids.

The syntax is

{SORC01} CO =SISBC(expression) and/or {SORC01} VAL=SISBC(expression)

Notes:

- The function SISBC is especially usefull when irregular solid shapes are approximated by series of "blocked-off" regions.
- By SISBC, the pressure load boundary conditions can be easily introduced, if appropriate.
- When used, SISBC must be applied both to set CO and VAL;
- It is used along with grid-free planting of COefficient and VAlue of linearized sources as the following settings exemplified.

**=> Example 1.**

The next three lines

PATCH(SS111NOR,NORTH,1,NX,1,NY,1,NZ,1,LSTEP) {SORC01} CO =SISBC(FIXFLU) {SORC01} VAL=SISBC(0.0/FIXFLU) COVAL(SS111NOR,V1,GRND,GRND)

set zero-normal stress boundary conditions at the north faces of all cells adjacent to interface between fluid and solid with PRPS=111. (steel).

**=> Example 2.**

The constrains are introduced in similar manner, so that

PATCH(SS111NOR,NORTH,1,NX,1,NY,1,NZ,1,LSTEP) {SORC01} CO =SISBC(FIXVAL) {SORC01} VAL=SISBC(0.0) COVAL(SS111NOR,V1,GRND,GRND)

would introduce the constrains at all north faces of the cells adjacent to interface between fluid and solid with PRPS=111.

The similar blocks are used for WEST-EAST and LOW-HIGH directions.

____________________________________________________________________

SORC?? .... statement heading used to instruct PLANT that it must
write its FORTRAN block of number ?? in Group 13
- Section 1 when VAR is CO and/or

- Section 12 when VAR is VAL

of GROUND to set the formulae for the linearized sources.

____________________________________________________________________

SOUTH(PHI)....name used to denote the value of PHI on the south side of the cell in question.

___________________________________________________________________

SPHERE...function used as part of PLANT statement.

The syntax is:

MARK=SPHERE(MRK,XPP,YPP,ZPP,RADIUS)

It sets the MARK values equal to MRK inside the surface of spherical object placed on grid.

The function has 5 arguments, namely (in order):

MRK - MARK value which is to be associated with the object;

XPP - X - coordinate of the centre;

YPP - Y - coordinate of the centre;

ZPP - Z - coordinate of the centre;

RADIUS- sphere radius;

___________________________________________________________________

SUM...function which plants the GROUND coding in above Groups of GROUND to calculate the single value or the number of single values as a result of execution of the summation over the whole domain ( or over the specified part of flow domain ) for the expression specified as an argument.

The syntax is:

VAR = SUM (expression),

where VAR stands for automatically declared real variable name.

The function is followed by the mean to dump the resulting value(s) in special file called GLOBCALC with the facility for labelling them by command TEXT.

=> Example:

{SC0601> R=sum(SQRT(U1**2+V1**2)) REGION(1,NX-1,1,NY-1,1,1)

Here the line marked as {SC0601> contains the value to be computed, R, and the mathematical expression to be summed over the extent of the flow domain specified by REGION command. For the particular case shown, R will be the sum of velocity magnitudes at XY plane.

____________________________________________________________________

TEMP1....name used to denote first phase temperature in response of the command, TMP1=GRND.

___________________________________________________________________

TEMP2....name used to denote second phase temperature in response of the command, TMP2=GRND.

___________________________________________________________________

TEXT .... command used to transmit up to 30 characters of text which will be printed in GLOBCALC file. It is used inconjunction with SUM function only.

The syntax is:

TEXT ( Character string )

Note: it can be used in conjunction with SUM function only.

=> Example:

{SC0601> R=sum(SQRT(U1**2+V1**2)) TEXT (Variable R) REGION(1,NX-1,1,NY-1,1,1)

Here the line marked as {SC0601> contains the value to be computed, R, and the mathematical expression to be summed over the extent of the flow domain specified by REGION command. For the particular case shown, R will be the sum of velocity magnitudes at XY plane.

The globcalc file of working directory will contain the following lines:

Global calculations: Variable R = 4.70132

____________________________________________________________________

TIM....is the current time.

___________________________________________________________________

TINY....a small real number, defaulted to 1.E-20 which is used in EARTH to guard against division by zero. This value may be too small for machines with short word length, and would normally be set to an appropriate value at the time of installation.

___________________________________________________________________

TLAST....multiplier of TFRAC.

It must be set to GRND as a part of block for a setting the formulae to calculate time-step size DT.

____________________________________________________________________

TMP1 .... is the temperature of the first phase.

____________________________________________________________________

TMP2 .... is the temperature of the second phase.

___________________________________________________________________

U1....standard name used to denote the first-phase velocity component in the x-direction.

___________________________________________________________________

U2....standard name used to denote the second-phase velocity component in the x-direction.

___________________________________________________________________

U1AD .... extra velocity to be added to U1. It must be set to GRND as a part of block for a setting the formulae to calculate the values to add to the field of U1 at each Z-slab.

___________________________________________________________________

U2AD .... extra velocity to be added to U2. It must be set to GRND as a part of block for a setting the formulae to calculate the values to add to the field of U2 at each Z-slab.

___________________________________________________________________

V1....standard name used to denote the first-phase velocity component in the y-direction.

___________________________________________________________________

V2....standard name used to denote the second-phase velocity component in the y-direction.

___________________________________________________________________

V1AD .... extra velocity to be added to V1. It must be set to GRND as a part of block for a setting the formulae to calculate the values to add to the field of V1 at each Z-slab.

___________________________________________________________________

V2AD .... extra velocity to be added to V2. It must be set to GRND as a part of block for a setting the formulae to calculate the values to add to the field of V2 at each Z-slab.

___________________________________________________________________

VAL....name which refers to the 'value' array, and which is used to set the values of linearised source expressions.

___________________________________________________________________

VARMAX(PHI) .... name used to denote the maximum value allowed to variable indicated, and which is used to set the specific relationship. If VARMAX(PHI) is less than -1.e-10, it will be treated as the maximum absolute value of the increment of the value at any updating stage.

___________________________________________________________________

VARMIN(PHI) .... name used to denote the minimum value allowed to variable indicated, and which is used to set the specific relationship. If VARMIN(PHI) is less than -1.e-10, it will be assumed that the absolute change made to the variable at all updating stages is limited to VARMAX.

___________________________________________________________________

VELAD....name used to denote extra addition to the current slab store of velocities, at the point in the computation at which the convection fluxes are assembled.

___________________________________________________________________

VISL....name used to denote the values of laminar kinematic viscosity of phase 1 in response of command, ENUL=GRND.

___________________________________________________________________

VIST....name used to denote the values of the turbulent kinematic viscosity of phase 1, in response of command, ENUT=GRND.

___________________________________________________________________

VOL....name used to denote an array containing the 'free' cell volumes available for fluid after blockages are deducted.

___________________________________________________________________

W1....standard name used to denote the first-phase velocity component in the z-direction.

___________________________________________________________________

W2....standard name used to denote the second-phase velocity component in the z-direction.

___________________________________________________________________

W1AD .... extra velocity to be added to W1. It must be set to GRND as a part of block for a setting the formulae to calculate the values to add to the field of W1 at each Z-slab.

___________________________________________________________________

W2AD .... extra velocity to be added to W2. It must be set to GRND as a part of block for a setting the formulae to calculate the values to add to the field of W2 at each Z-slab.

____________________________________________________________________

WALDIS....name used to denote the distance from the wall to the first internal node. It is used in connection with the wall functions for turbulence models.

___________________________________________________________________

WEST(PHI)....name used to denote the value of PHI on the west side of the cell in question.

___________________________________________________________________

WGRDNZ....name used to denote the values of the w-velocity of the grid at the high face of each z- slab when the moving grid option is active.

___________________________________________________________________

XC....is the symbol used for body-fitted coordinate grids to represent the cartesian coordinate x of the low-south-west corner of the continuity cell I, J, K.

___________________________________________________________________

XCO....is the symbol used for body-fitted coordinate grids to represent the "old" cartesian coordinate x of the low-south-west corner of the continuity cell I, J, K.

___________________________________________________________________

XG2D....name used to denote the array of values of: either (for CARTES=T or BFC=T) distances of the centres of the continuity cells from the x=0.0 plane in a Cartesian coordinate system, or (for CARTES=F and BFC=F) angular coordinates of the centres of the continuity cells in a cylindrical-polar system.

___________________________________________________________________

XRAT....is the name used to denote the multiplier of the previous XULAST which gives the current XULAST in response of the command, AZXU = GRND.

___________________________________________________________________

XU2D....name used to denote the values, as for XG2D except that the east faces of continuity cells are involved, not cell centres. Unlike XG2D, XU2D has no significance when BFC=T.

___________________________________________________________________

XYBOX...function used as part of PLANT statement.

The syntax is:

MARK=XYBOX(MRK,XPP,YPP,XSIZE,YSIZE,ANG,ROTANG)

It sets the MARK values equal to MRK inside the surface of rectangular object placed on XY plane.

The function has 6 arguments, namely (in order):

MRK - MARK value which is to be associated with the object;

XPP - X - coordinate of south low corner of rectangle;

YPP - Y - coordinate of south low corner of rectangle;

XSIZE - object size in X-direction;

YSIZE - object size in Y-direction;

ANG - angle of counter clockwise rotation of object local

coordinate system about Z-axis , rad.

ROTANG- angle of shifting reference cartesian frame in the

X-direction ( significant for CARTES=F), rad.

___________________________________________________________________

XYCIRC...function used as part of PLANT statement.

The syntax is:

MARK=XYCIRC(MRK,XPP,YPP,RADIUS)

It sets the MARK values equal to MRK inside the surface of circle object placed on XY plane.

The function has 4 arguments, namely (in order):

MRK - MARK value which is to be associated with the object;

XPP - X - coordinate of the centre;

YPP - Y - coordinate of the centre;

RADIUS- circle radius;

___________________________________________________________________

XYELLP...function used as part of PLANT statement.

The syntax is:

MARK=XYELLP(MRK,XPP,YPP,AX1,AX2,ANG,ROTANG)

It sets the MARK values equal to MRK inside the surface of ellipse object placed on XY plane.

The function has 7 arguments, namely (in order):

MRK - MARK value which is to be associated with the object;

XPP - X - coordinate of the centre;

YPP - Y - coordinate of the centre;

AX1 - half axis size in X-direction;

AX2 - half axis size in Y-direction;

ANG - angle of counter clockwise rotation of object local

coordinate system about Z-axis , rad.

ROTANG- angle of shifting reference cartesian frame in the

X-direction ( significant for CARTES=F), rad.

___________________________________________________________________

XYWEDG...function used as part of PLANT statement.

The syntax is:

MARK=XYWEDG(MRK,XPP,YPP,XSIZE,YSIZE,ANG,ROTANG)

It sets the MARK values equal to MRK inside the surface of right angle wedge object placed on XY plane.

The function has 7 arguments, namely (in order):

MRK - MARK value which is to be associated with the object;

XPP - X - coordinate of right angle vertice;

YPP - Y - coordinate of right angle vertice;

XSIZE - wedge size in X-direction;

YSIZE - wedge size in Y-direction;

ANG - angle of counter clockwise rotation of object local

coordinate system about Z-axis , rad.

ROTANG- angle of shifting reference cartesian frame in the

X-direction ( significant for CARTES=F), rad.

___________________________________________________________________

YC....is the symbol used for body-fitted coordinate grids to represent the cartesian y-coordinate of the low-south-west corner for the continuity cells I, J, K.

___________________________________________________________________

YCO....is the symbol used for body-fitted coordinate grids to represent the "old" cartesian y-coordinate of the low-south-west corner of the continuity cell I, J, K.

___________________________________________________________________

YG2D....name used to denote the values of the distances of the centres of the continuity cells from the y=0.0 plane (for CARTES=T or BFC=T) or cylinder (for CARTES=F and BFC=F).

___________________________________________________________________

YRAT....is the name used to denote the multiplier of the previous YVLAST which gives the current YVLAST in response of the command, AZYV = GRND.

___________________________________________________________________

YV2D....name used to denote the values, as for YG2D, except that the north faces of cells are involved, not cell centres. Unlike YG2D, YV2D has no significance when BFC=T.

___________________________________________________________________

YZBOX...function used as part of PLANT statement.

The syntax is:

MARK=YZBOX(MRK,YPP,ZPP,YSIZE,ZSIZE,ANG,ROTANG)

It sets the MARK values equal to MRK inside the surface of rectangular object placed on YZ plane.

The function has 6 arguments, namely (in order):

MRK - MARK value which is to be associated with the object;

YPP - Y - coordinate of south low corner of rectangle;

ZPP - Z - coordinate of south low corner of rectangle;

YSIZE - object size in Y-direction;

ZSIZE - object size in Z-direction;

ANG - angle of counter clockwise rotation of object local

coordinate system about X-axis , rad.

ROTANG- angle of shifting reference cartesian frame in the

X-direction ( significant for CARTES=F), rad.

___________________________________________________________________

ZC....is the symbol used for body-fitted coordinate grids to represent the cartesian z-coordinate of the low-south-west corner for the continuity cell I, J, K.

___________________________________________________________________

ZCO....is the symbol used for body-fitted coordinate grids to represent the "old" cartesian z-coordinate of the low-south-west corner of the continuity cell I, J, K.

___________________________________________________________________

ZGNZ....name used to denote the values for the distances of the centres of the continuity cells from the z=0.0 plane.

___________________________________________________________________

ZW....is the z coordinate of the high face of the current slab.

___________________________________________________________________

ZWNZ...name used to denote the values for distances of the HIGH faces of the continuity cells from the z=0.0 plane.

___________________________________________________________________

**
**