TALK=T;RUN(1,1)
PHOTON USE
p
msg
msg
msg All examples which follow are assembled as different slabs
msg and number of slabs of one 3D domain called
msg
msg INITBOX.
msg
msg The series of PHOTON plots will follow.
msg
msg They illustrate some of the existing initialisation
msg techniques available in PLANT.
msg
msg The presentation covers:
msg
msg * Flow field initialistions by parametric analytics
msg * Manipulating with initial fields
msg * MARKing sub-domains to create arbitrary initial fields
msg * Geometry initialisations
msg
msg Hit enter to proceed to 1st slab.
msg
pause
clear
vi z
set vec ref 25
vec z 1 sh
msg
msg
msg The 1st slab of the box is initialised by stagnation
msg point flow
msg Hit enter to proceed to 2nd slab.
pause
vec off;red
vec z 2 sh
msg
msg
msg The 2nd slab is initialised by solid body rotation flow
msg
msg Hit enter to proceed to 3rd slab.
pause
vec off;red
vec z 3 sh
msg
msg The flow in the 3rd slab is constructed by adding the
msg velocity distributions of 1st and 2nd slabs.
msg
msg Hit enter to proceed to 4th slab.
pause
vec off;red
con mark z 4 fil;.001
set vec aver off
vec z 4
msg
msg
msg
msg
msg
msg
msg
msg
msg
msg
msg
msg
msg
msg
msg
msg
msg
msg
msg The flow in the 4th slab is constructed by making use
msg the different distributions across specifically MARKed
msg regions so as to have better accord with imposed geometry.
msg
msg The latter is generated by PLANTed geometry functions.
msg
msg Hit enter to proceed to 5th section of the box.
pause
* gr z 4
gr off;vec off;con off;red
vi 1 1 1
con prps z 24 x 21 40 y 1 m fil;.01
con prps x 20 y 1 m z 25 44 fil;.01
msg
msg The 5th section of the box occupies 40 slabs.
msg
msg It illustrates the construction of 3d geometry visualised
msg by PRPS surface.
msg
msg The PRPS values are initialised to follow the distribution
msg of MARK values.
msg
msg PRPS surface is being generated.
msg Please, wait ....
surf mark x 5.99
msg
msg Hit E to finish.
ENDUSE
DISPLAY
This file exemplifies some of the basic techniques
available in PLANT for field initialisations.
Five examples are assembled in a single 3D box. Each
example occupies one or more slabs and is wholly
focussed on initialization procedures. No field
distribution are calculated.
The examples include:
* Flow field initializations by parametric analytics;
* Manipulating with initial fields;
* MARKing sub-domains to create arbitrary initial fields
and
* Make complex 2D/3D geometry initialisations.
ENDDIS
PLANTBEGIN
Example 1: Initialization of stagnation point flow
=======================================
PATCH(INI1,INIVAL,1,NX,1,NY,1,1,1,1)
VAL=XU2D-10.
COVAL(INI1,U1,zero,GRND)
VAL=-(YV2D-10.)
COVAL(INI1,V1,zero,GRND)
The 1st slab of the box is initialized by stagnation
point flow with the cartesian components as follows:
U1 = X - 10 and
V1 = 10- Y.
In above settings:
XU2D stands for the distances of the east faces of
continuity cells from the x=0.0 plane and
YV2D means the distances of the north faces of
continuity cells from the y=0.0 plane.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
Example 2: Solid-body rotation flow
========================
PATCH(INI2,INIVAL,1,NX,1,NY,2,2,1,1)
VAL=YG2D-10.
COVAL(INI2,U1,zero,GRND)
VAL=-(XG2D-10.)
COVAL(INI2,V1,zero,GRND)
The 2nd slab is initialized by solid body rotation flow
about domain centre with cartesian components as follows:
U1 = Y - 10 and
V1 = 10- X.
In above settings:
XG2D stands for the distances of the centres of
continuity cells from the x=0.0 plane and
YG2D means the distances of centres of continuity
cells from the y=0.0 plane.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
Example 3: Flow superposition
==================
PATCH(INI3,INIVAL,1,NX,1,NY,3,3,1,1)
VAL=-U1[,,-1]+U1[,,1]
COVAL(INI3,U1,zero,GRND)
VAL=-V1[,,2]+V1[,,-2]
COVAL(INI3,V1,zero,GRND)
The flow in the 3rd slab is initialized by adding the
velocity distributions of 1st and 2nd slabs.
In above settings:
a) Indeces [,,-1] and [,,2] refer to the cell velocities
at 2nd slab and
b) [,,1] and [,,-2] refer to the 1st slab.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
Example 4: : 2D geometry and velocities initializations
==========================================
The flow in the 4th slab is initialized by making use
the different distributions over specifically MARKed
regions so as to have better accord with imposed
geometry. The latter is the circle chamber with two
tangentiall passages.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** 2D geometry, X-Y plane
* Half-circle marked by 2.
FIINIT(MARK)=1.0
PATCH(INIT70,INIVAL,1,NX/2,1,NY,4,4,1,1)
VAL=XYELLP(2.,10.,10.,8.,8.,0.,0.)
INIT (INIT70,MARK,0.,GRND)
In above statement, XYELLP function is used to make the
half-circle of 16 m diameter as follows:
a) In the west half of 4th slab,
b) place the ellipse marked by 2 (1st argument), with the
centre at XC=10 m (2nd argument) and YC=10 m (3rd
argument) and both half-axes equal to 8 m ( 4th and
5th arguments).
c) The 6th and 7th function arguments are insignificant
for the circle shape and CARTES=T.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
* Half-circle marked by 3.
PATCH(INIT71,INIVAL,NX/2+1,NX,1,NY,4,4,1,1)
VAL=XYELLP(3.,10.,10.,8.,8.,0.,0.)
INIT (INIT71,MARK,0.,GRND)
In above statement, XYELLP function is used to make the
half-circle of 16 m diameter in the east half of the
domain. The cells inside ellipse are marked by 3.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
* Rectangular marked by 4.
PATCH(INIT08,INIVAL,1,NX/2,1,NY,4,4,1,1)
VAL=XYELLP(4.,10.,17.,100.,1.,0.,0.)
INIT (INIT08,MARK,0.,GRND)
In above statement, XYELLP function is used to make the
rectangular of 2 m width and 10 m length, tangential to
the circle, as follows:
a) In the west half of 4th slab,
b) place the ellipse marked by 4 (1st argument), with the
centre at XC=10 m (2nd argument) and YC=17 m (3rd
argument), large (100.m) X-direction half-axis (4th
argument) and Y-direction half axis equals to 1.m (5th
argument).
c) The 6th and 7th function arguments are insignificant
for the circle shape and CARTES=T.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
* Rectangular marked by 5.
PATCH(INIT09,INIVAL,NX/2+1,NX,1,NY,4,4,1,1)
VAL=XYELLP(5.,10.,3.,100.,1.,0.,0.)
INIT (INIT09,MARK,0.,GRND)
In above statement, XYELLP function is used to make the
rectangular of 2 m width and 10 m length, tangential to
the circle, in the east half of the domain. The cells
inside ellipse are marked by 5.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** Velocities
The MARKed region can not be used be used to initialize
the velocities in Group 11. Two altenatives are usefull:
1) Make the initializations at the start of IZ slab for
the first sweep as shown below:
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
* Uniform negative velocity of 10 m/s
in sub-domain MARked by 4
U1=-10.
REGION() 4 /ISWEEP.EQ.1
V1=0.
REGION() 4 /ISWEEP.EQ.1
* Uniform negative velocity of 10 m/s
in sub-domain MARked by 5
U1=-10.
REGION() 5 /ISWEEP.EQ.1
V1=0.
REGION() 5 /ISWEEP.EQ.1
* Upward velocity field of second slab
in sub-domain MARked by 2
U1=U1[,,2]
REGION() 2 /ISWEEP.EQ.1
V1=V1[,,2]
REGION() 2 /ISWEEP.EQ.1
* Downward velocity field of second slab
in sub-domain MARked by 3
U1=-U1[,,2]
REGION() 3 /ISWEEP.EQ.1
V1=-V1[,,2]
REGION() 3 /ISWEEP.EQ.1
The crucial features of the above settings are:
1 Using dummy REGION command with MARK as parameter
which ensures that the velocities are initialized in
accord with the statement relationship for all cells
marked by parameter value;
2 The setting of the switch to first sweep and
3 Using indicial expressions to refer to the values
required.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
U1AD=GRND
VELAD=-10.
REGION() 4 /ISWEEP.EQ.1
2) Add initial velocities to the default ( zeros) values
by means of VELAD option in Group 8 for the first
sweep as exemplified for region marked 4 above.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
Example 5: 3D geometry and PRPS initializations
====================================
Unlike for previous examples, the section of 40 IZ-slabs
is used here to represent the generation of 3D geometry.
The geometry in question is a spherical chamber of 16 m
diameter with two cylindrical passages tangentially
joined the sphere.
The generation procedure is to set the MARK values equal
to 6.0 inside the surface of two-cylinders-sphere
assembly.
Then the PRPS are initialized to follow the distribution
of MARK values.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** 3D geometry, X-Y-Z space
* Provide the 90 degree angle as a parameter
REAL(ANG)
ANG=3.14159/2. ; RG(2)=ANG
* Sphere marked by 6.
PATCH(INIT10,INIVAL,1,NX,1,NY,5,NZ,1,1)
VAL=SPHERE(6.,10.,10.,12.,8.)
INIT (INIT10,MARK,0.,GRND)
In above statement, SPHERE function is used to make the
sphere of 16 m diameter as follows:
a) In the section of 40 IZ-slabs,
b) place the sphere marked by 6 (1st argument), with the
centre at XC=10 m (2nd argument), YC=10 m (3rd
argument) and ZC=12 m (4th argument). The 5th argument
is the sphere radius, 8 m.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
* Cylinder marked by 6.
PATCH(INIT11,INIVAL,1,NX/2,1,NY,5,NZ,1,1)
VAL=ELLPSD(6.,10.,16.,12.,100.,2.,2., 0.,RG(2$
),RG(2),RG(2),0.,RG(2),RG(2),RG(2),0.,0.)
INIT (INIT11,MARK,0.,GRND)
In above statement, ELLPSD function is used to make the
cylinder of 2 m radius, tangentially joining the sphere,
as follows:
a) In the east half of the section,
b) place the ellipsoid marked by 6 (1st argument), with
the centre at XC=10 m (2nd argument), YC=16 m (3rd
argument) and ZC=12 m (4th argument).
c) X-direction ellipsoid half-axis (5th argument) is set
to large (100.m) value to ensure that the shape of
ellipsoid is virtually cylindrical.
d) The Y- (6th argument) and Z- (7th argument) directions
half-axis are set to be cylinder radia, 2 m.
e) The arguments 8.9 and 10 set zero angle between
ellipsoid X-axis and corresponding cartesian frame
direction, while the angles between the former and
other frame directions are kept 90 degrees.
f) The arguments 11.12 and 13 set 90 degree angle between
ellipsoid Y-axis and X-cartesian frame direction, the
angle between the former and the Y-cartesian frame
direction is nulified and angle between ellipsoid
Y-axis and Z-frame direction is kept 90 degree.
g) The arguments 14.15 and 16 set 90 degree angles
between ellipsoid Z-axis and X-, Y-cartesian frame
directions, while the angle between the former and the
Z-cartesian frame direction is nulified.
h) The 17th argument, set to zero, is insignificant for
CARTES=T.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
* Cylinder marked by 6.
PATCH(INIT12,INIVAL,NX/2+1,NX,1,NY,5,NZ,1,1)
VAL=ELLPSD(6.,10.,4.,12.,100.,2.,2., 0.,RG(2)$
,RG(2),RG(2),0.,RG(2),RG(2),RG(2),0.,0.)
INIT (INIT12,MARK,0.,GRND)
In above statement, ELLPSD function is used to make the
cylinder of 2 m radius, tangentially joining the sphere
for the west half of the section as explained earlier.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** Properties
PATCH(SS001EXT,INIVAL,1,NX,1,NY,5,NZ,1,1)
VAL=40.
INIT(SS001EXT,PRPS,0.,GRND)
Initialization of PRPS field by 40 over all cells having
MARK value equal unity. The crucial feature of above
settings is that special PATCH name, SS001???, is to be
used, in which first two characters specify the
particular action, Set Space source, and number 001 is
the MARK value required.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PATCH(SS006INT,INIVAL,1,NX,1,NY,5,NZ,1,1)
INIT16 VAL=1.
INIT(SS006INT,PRPS,0.,GRND)
Initialization of PRPS field by unity over all cells
having MARK value equal six.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PLANTEND
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT( Five examples initialization box )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 605
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 40
XULAST =20.
XFRAC(1)=0.025 ;XFRAC(2)=0.05
XFRAC(3)=0.075 ;XFRAC(4)=0.1
XFRAC(5)=0.125 ;XFRAC(6)=0.15
XFRAC(7)=0.175 ;XFRAC(8)=0.2
XFRAC(9)=0.225 ;XFRAC(10)=0.25
XFRAC(11)=0.275 ;XFRAC(12)=0.3
XFRAC(13)=0.325 ;XFRAC(14)=0.35
XFRAC(15)=0.375 ;XFRAC(16)=0.4
XFRAC(17)=0.425 ;XFRAC(18)=0.45
XFRAC(19)=0.475 ;XFRAC(20)=0.5
XFRAC(21)=0.525 ;XFRAC(22)=0.55
XFRAC(23)=0.575 ;XFRAC(24)=0.6
XFRAC(25)=0.625 ;XFRAC(26)=0.65
XFRAC(27)=0.675 ;XFRAC(28)=0.7
XFRAC(29)=0.725 ;XFRAC(30)=0.75
XFRAC(31)=0.775 ;XFRAC(32)=0.8
XFRAC(33)=0.825 ;XFRAC(34)=0.85
XFRAC(35)=0.875 ;XFRAC(36)=0.9
XFRAC(37)=0.925 ;XFRAC(38)=0.95
XFRAC(39)=0.975 ;XFRAC(40)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 40
YVLAST =20.
YFRAC(1)=0.025 ;YFRAC(2)=0.05
YFRAC(3)=0.075 ;YFRAC(4)=0.1
YFRAC(5)=0.125 ;YFRAC(6)=0.15
YFRAC(7)=0.175 ;YFRAC(8)=0.2
YFRAC(9)=0.225 ;YFRAC(10)=0.25
YFRAC(11)=0.275 ;YFRAC(12)=0.3
YFRAC(13)=0.325 ;YFRAC(14)=0.35
YFRAC(15)=0.375 ;YFRAC(16)=0.4
YFRAC(17)=0.425 ;YFRAC(18)=0.45
YFRAC(19)=0.475 ;YFRAC(20)=0.5
YFRAC(21)=0.525 ;YFRAC(22)=0.55
YFRAC(23)=0.575 ;YFRAC(24)=0.6
YFRAC(25)=0.625 ;YFRAC(26)=0.65
YFRAC(27)=0.675 ;YFRAC(28)=0.7
YFRAC(29)=0.725 ;YFRAC(30)=0.75
YFRAC(31)=0.775 ;YFRAC(32)=0.8
YFRAC(33)=0.825 ;YFRAC(34)=0.85
YFRAC(35)=0.875 ;YFRAC(36)=0.9
YFRAC(37)=0.925 ;YFRAC(38)=0.95
YFRAC(39)=0.975 ;YFRAC(40)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 44
ZWLAST =22.
ZFRAC(1)=0.022727 ;ZFRAC(2)=0.045455
ZFRAC(3)=0.068182 ;ZFRAC(4)=0.090909
ZFRAC(5)=0.113636 ;ZFRAC(6)=0.136364
ZFRAC(7)=0.159091 ;ZFRAC(8)=0.181818
ZFRAC(9)=0.204545 ;ZFRAC(10)=0.227273
ZFRAC(11)=0.25 ;ZFRAC(12)=0.272727
ZFRAC(13)=0.295455 ;ZFRAC(14)=0.318182
ZFRAC(15)=0.340909 ;ZFRAC(16)=0.363636
ZFRAC(17)=0.386364 ;ZFRAC(18)=0.409091
ZFRAC(19)=0.431818 ;ZFRAC(20)=0.454545
ZFRAC(21)=0.477273 ;ZFRAC(22)=0.5
ZFRAC(23)=0.522727 ;ZFRAC(24)=0.545455
ZFRAC(25)=0.568182 ;ZFRAC(26)=0.590909
ZFRAC(27)=0.613636 ;ZFRAC(28)=0.636364
ZFRAC(29)=0.659091 ;ZFRAC(30)=0.681818
ZFRAC(31)=0.704545 ;ZFRAC(32)=0.727273
ZFRAC(33)=0.75 ;ZFRAC(34)=0.772727
ZFRAC(35)=0.795455 ;ZFRAC(36)=0.818182
ZFRAC(37)=0.840909 ;ZFRAC(38)=0.863636
ZFRAC(39)=0.886364 ;ZFRAC(40)=0.909091
ZFRAC(41)=0.931818 ;ZFRAC(42)=0.954545
ZFRAC(43)=0.977273 ;ZFRAC(44)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(3)=U1 ;NAME(5)=V1
NAME(149)=MARK ;NAME(150)=PRPS
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(U1,Y,N,N,N,N,Y)
SOLUTN(V1,Y,N,N,N,N,Y)
SOLUTN(MARK,Y,N,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,Y)
PRPS = 150
************************************************************
Group 8. Terms & Devices
* Y in TERMS argument list denotes:
* 1-built-in source 2-convection 3-diffusion 4-transient
* 5-first phase variable 6-interphase transport
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
************************************************************
Group 9. Properties used if PRPS is not
stored, and where PRPS = -1.0 if it is!
RHO1 =1. ;TMP1 =0. ;EL1 =0.
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
ENUL =1.0E-05 ;ENUT =0.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(U1)=1.0E-10 ;FIINIT(V1)=1.0E-10
FIINIT(MARK)=1. ;FIINIT(PRPS)=-1.
PATCH(INI1 ,INIVAL, 1, 40, 1, 40, 1, 1, 1, 1)
INIT(INI1 ,U1 ,0. , GRND )
INIT(INI1 ,V1 ,0. , GRND )
PATCH(INI2 ,INIVAL, 1, 40, 1, 40, 2, 2, 1, 1)
INIT(INI2 ,U1 ,0. , GRND )
INIT(INI2 ,V1 ,0. , GRND )
PATCH(INI3 ,INIVAL, 1, 40, 1, 40, 3, 3, 1, 1)
INIT(INI3 ,U1 ,0. , GRND )
INIT(INI3 ,V1 ,0. , GRND )
PATCH(INIT70 ,INIVAL, 1, 20, 1, 40, 4, 4, 1, 1)
INIT(INIT70 ,MARK,0. , GRND )
PATCH(INIT71 ,INIVAL, 21, 40, 1, 40, 4, 4, 1, 1)
INIT(INIT71 ,MARK,0. , GRND )
PATCH(INIT08 ,INIVAL, 1, 20, 1, 40, 4, 4, 1, 1)
INIT(INIT08 ,MARK,0. , GRND )
PATCH(INIT09 ,INIVAL, 21, 40, 1, 40, 4, 4, 1, 1)
INIT(INIT09 ,MARK,0. , GRND )
PATCH(INIT10 ,INIVAL, 1, 40, 1, 40, 5, 44, 1, 1)
INIT(INIT10 ,MARK,0. , GRND )
PATCH(INIT11 ,INIVAL, 1, 20, 1, 40, 5, 44, 1, 1)
INIT(INIT11 ,MARK,0. , GRND )
PATCH(INIT12 ,INIVAL, 21, 40, 1, 40, 5, 44, 1, 1)
INIT(INIT12 ,MARK,0. , GRND )
PATCH(SS001EXT,INIVAL, 1, 40, 1, 40, 5, 44, 1, 1)
INIT(SS001EXT,PRPS,0. , GRND )
INIADD = F
FSWEEP = 1
NAMFI =CHAM
************************************************************
Group 12. Patchwise adjustment of terms
Patches for this group are printed with those
for Group 13.
Their names begin either with GP12 or &
************************************************************
Group 13. Boundary & Special Sources
No PATCHes yet used for this Group
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
************************************************************
Group 17. Relaxation
RELAX(U1,LINRLX,1.)
RELAX(V1,LINRLX,1.)
RELAX(MARK,LINRLX,1.)
RELAX(PRPS,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(MARK)=1.0E+10 ;VARMIN(MARK)=-1.0E+10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
NAMSAT =MOSG
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,PLOTALL,L,T)
RG( 2) =1.570795
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(U1)=0.4311 ;EX(V1)=0.4126
EX(MARK)=2.324 ;EX(PRPS)=27.610001
************************************************************
Group 21. Print-out of Variables
INIFLD = F ;SUBWGR = F
* Y in OUTPUT argument list denotes:
* 1-field 2-correction-eq. monitor 3-selective dumping
* 4-whole-field residual 5-spot-value table 6-residual table
OUTPUT(U1,Y,N,Y,N,N,N)
OUTPUT(V1,Y,N,Y,N,N,N)
OUTPUT(MARK,Y,N,Y,N,N,N)
OUTPUT(PRPS,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 1 ;IZMON = 1
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
XZPR = F ;YZPR = F
IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
ISWPRF = 1 ;ISWPRL = 100000
ITABL = 3 ;IPROF = 1
ABSIZ =0.5 ;ORSIZ =0.4
NTZPRF = 1 ;NCOLPF = 50
ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
No PATCHes yet used for this Group
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP