TALK=T;RUN(1,1)
DISPLAY
PLANT is used for specifying two time-step sizes.
The first is set equal to 5 sec.
The second time step is computed by dividing the smallest cell
size by the largest velocity.
The solved-for variable is H1; but, since the actual solution is
of not of interest, no initial or boundary conditions are
introduced.
The non-uniform grid and velocity field are initialised
for expected size of the second time step to be 0.1 sec.
PLANT information :
* Data input groups used: 2, 19
* Ground groups planted : 1, 2, 19-3
* Headings used : SCTS??
* Functions used : SUM
* Commands used : IF, REGION, TEXT.
ENDDIS
PLANTBEGIN
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** Ask PLANT to introduce first time step, 5 seconds.
DT=5.
REGION(1,1,1,1,1,1,1,1)
REGION command is used to control this simple setting.
Its first six arguments limit DO loop to perform single
action over first time step defined by the last argument
pair.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** Choose the smallest cell size at the start of each iz-slab
for the last sweep of the first time step.
RG(1)=GREAT
DXG2D - cell size in X-direction,
DYG2D - cell size in Y-direction,
DZWNZ - cell size in Z-direction,
RG(1)=AMIN1(RG(1),DXG2D,DYG2D,DZWNZ)
IF(ISTEP.EQ.1.AND.ISWEEP.EQ.LSWEEP)
Command IF restricts the whole-domain-default extents of
PLANT statement. The arguments of IF command instruct
PLANT to introduce the logical conditions of the
bracketed expression. The execution of the statement
returns the auxilliary variable, RG(1), with the
smallest cell size over the whole domain.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** Output of smallest cell size using summation in one
cell at 1st time step, last sweep for checking.
SIZMIN=SUM(RG(1))
TEXT(Smallest cell size)
REGION(1,1,1,1,1,1,1,1) /ISWEEP.EQ.LSWEEP
The function SUM is used in above three lines just to
dump into globcalc file the value of RG(1) headed by the
character arguments of TEXT command.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
** Choose the largest velocity value just before the
second time step calculation
RG(2)=AMAX1(RG(2),U1,V1,W1)
IF(ISTEP.EQ.2)
** Output the largest velocity using summation in one
cell at 2nd time step for checking.
VELMAX=SUM(RG(2))
TEXT(Largest velocity)
REGION(1,1,1,1,1,1,2,2)
Auxiliary variable, RG(2), will be the largest velocity
value as result of the above statements at the second
time moment. It will be dumped in globcalc file and
appropriately headed.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
DT=SIZMIN/VELMAX
REGION(1,1,1,1,1,1,2,2)
The final stage of time step settings: second time step
size, DT, is set as division of smallest size by largest
velocity.
<<<<<<<<<<<<<<<<<<<<<<< Comment ends <<<<<<<<<<<<<<<<<<<<<
PLANTEND
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(TIME STEP CALCULATIONS )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 619
************************************************************
Group 2. Time dependence
STEADY = F
* Set overall time and no. of steps
TFIRST =0. ;TLAST = GRND
FSTEP = 1 ;LSTEP = 2
TFRAC(1)=1.
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 2
XULAST =2.
XFRAC(1)=0.5 ;XFRAC(2)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 2
YVLAST =4.
YFRAC(1)=0.5 ;YFRAC(2)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 2
ZWLAST =6.
ZFRAC(1)=0.5 ;ZFRAC(2)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(3)=U1 ;NAME(5)=V1
NAME(7)=W1 ;NAME(14)=H1
* 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(W1,Y,N,N,N,N,Y)
SOLUTN(H1,Y,Y,N,N,N,Y)
************************************************************
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
TERMS(H1,Y,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
HUNIT =1.
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.
PRNDTL(H1)=1.
PRT(H1)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(U1)=1. ;FIINIT(V1)=0.5
FIINIT(W1)=0.1 ;FIINIT(H1)=1.0E-10
PATCH(MAXVEL ,INIVAL, 1, 1, 2, 2, 2, 2, 1, 1)
INIT(MAXVEL ,U1 ,0. ,10. )
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 = 2 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(H1)=20
ENDIT(H1)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(U1,LINRLX,1.)
RELAX(V1,LINRLX,1.)
RELAX(W1,LINRLX,1.)
RELAX(H1,FALSDT,1.0E+09)
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(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
NAMSAT =MOSG
PARSOL = F
ISG62 = 1
SPEDAT(SET,GXMONI,TRANSIENT,L,F)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
RG( 1) =1.0E+20
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(U1)=3.25 ;EX(V1)=0.5
EX(W1)=0.1 ;EX(H1)=5.445E-04
************************************************************
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(W1,Y,N,Y,N,N,N)
OUTPUT(H1,Y,N,Y,Y,Y,Y)
************************************************************
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
NTPRIN = 100000 ;ISTPRF = 1 ;ISTPRL = 100000
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