TALK=T;RUN(1,1)
PHOTON USE
p;;;;;
vi z
con temp z 5 fi;0.0001
pause; cl; con cp1 z 5 fi;0.0001
msg specific heat formulae
pause; cl; con s5cp z 5 fi;0.0001
msg specific heat formulae
pause; cl; con l3cp z 5 fi;0.0001
msg specific heat formulae
pause; cl; con lfcp z 5 fi;0.0001
msg specific heat formulae
pause; cl; con enul z 5 fi;0.0001
msg kinematic viscosity formulae
pause; cl; con s5en z 5 fi;0.0001
msg kinematic viscosity formulae
pause; cl; con l3en z 5 fi;0.0001
msg kinematic viscosity formulae
pause; cl; con lfen z 5 fi;0.0001
msg kinematic viscosity formulae
pause; cl; cl; con rho1 z 5 fi;0.0001
msg density formulae
pause; cl; con s5rh z 5 fi;0.0001
msg density formulae
pause; cl; con l3rh z 5 fi;0.0001
msg density formulae
pause; cl; con lfrh z 5 fi;0.0001
msg density formulae
pause; cl; con cond z 5 fi;0.0001
msg thermal conductity formulae
pause; cl; con s5cn z 5 fi;0.0001
msg thermal conductity formulae
pause; cl; con l3cn z 5 fi;0.0001
msg thermal conductity formulae
pause; cl; con lfcn z 5 fi;0.0001
msg thermal conductity formulae
pause; cl; con ent1 z 5 fi;0.0001
msg enthalpy
ENDUSE
PHOTON USE
p
red
msg u1 contours. Press return for w1
con u1 z 1 fi;0.001;con u1 z 5 fi;0.001;con u1 z 10 fi;0.001
gr ou z m;gr ou y m;gr ou x m
pause
con off;red
msg w1 contours, and velocity vectors on the top.
msg Press return for h1
con w1 y 1 fi;0.001;con w1 y 5 fi;0.001;con w1 y m fi;0.001
vec y m
pause
vec off;con off;red
msg temperature contours. Press e to end
con h1 z 1 fi;0.001;con h1 z 5 fi;0.001;con h1 z 10 fi;0.001
enduse
DISPLAY
This simulation is similar to that of case 249; but the
cavity is limited in the third direction by adiabatic
end walls, at rest. Economy could be affected by noting
that symmetry must prevail about the NZ/2 plane; but
this is not done here. The symmetry of the solution
can therefore be checked..
_ _ _ _ _ _ _
/ /|
/ north wall / |
/ at H1=1. --/--|--->moving
/_ _ _ _ _ _ _/ | y
| | | ^ ^
| |east wall | /z
| low wall | /H1=0. | /
| adiabatic | / |/----> x
|_ _ _ _ _ _ _|/
EXPERT is in use
The Q1 contains PHOTON USE commands
ENDDIS
DISPLAY
Here the configuration of case 277 is used, for the purpose of
showing the various ways in which In-Form allows fluid properties
to be set.
The fluid is ethylene glycol.
Specifically illustrated are:
* properties from the polynomials in case 089
(these are the properties used in the solution for velocity
and temperature)
* a 3-part piece-wise linear function
* a spline function
* a multi-part piece-wise linear function from files, here called:
DENPRP, ENUPRP, CPPRP and CNDPRP,
These files are provided in /phoenics/d_earth/d_core/inplib
To make comparison easy, the properties computed in these four
different ways have been divided by the values actually used.
Their closeness to unity therefore measures their agreement.
The Q1 contains PHOTON USE commands
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(3D SQUARE CAVITY; various fluids: 763 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 763
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 10
XULAST =0.1
XFRAC(1)=0.02 ;XFRAC(2)=0.08
XFRAC(3)=0.18 ;XFRAC(4)=0.32
XFRAC(5)=0.5 ;XFRAC(6)=0.68
XFRAC(7)=0.82 ;XFRAC(8)=0.92
XFRAC(9)=0.98 ;XFRAC(10)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 10
YVLAST =0.1
YFRAC(1)=0.02 ;YFRAC(2)=0.08
YFRAC(3)=0.18 ;YFRAC(4)=0.32
YFRAC(5)=0.5 ;YFRAC(6)=0.68
YFRAC(7)=0.82 ;YFRAC(8)=0.92
YFRAC(9)=0.98 ;YFRAC(10)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 10
ZWLAST =0.1
ZFRAC(1)=0.02 ;ZFRAC(2)=0.08
ZFRAC(3)=0.18 ;ZFRAC(4)=0.32
ZFRAC(5)=0.5 ;ZFRAC(6)=0.68
ZFRAC(7)=0.82 ;ZFRAC(8)=0.92
ZFRAC(9)=0.98 ;ZFRAC(10)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(3)=U1
NAME(5)=V1 ;NAME(7)=W1
NAME(14)=TEMP ;NAME(134)=ENT1
NAME(135)=LFCN ;NAME(136)=S5CN
NAME(137)=L3CN ;NAME(138)=COND
NAME(139)=LFCP ;NAME(140)=S5CP
NAME(141)=L3CP ;NAME(142)=CP1
NAME(143)=LFEN ;NAME(144)=S5EN
NAME(145)=L3EN ;NAME(146)=ENUL
NAME(147)=LFRH ;NAME(148)=S5RH
NAME(149)=L3RH ;NAME(150)=RHO1
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(P1,Y,Y,Y,N,N,N)
SOLUTN(U1,Y,Y,N,N,N,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(TEMP,Y,Y,N,N,N,Y)
SOLUTN(ENT1,Y,N,N,N,N,Y)
SOLUTN(LFCN,Y,N,N,N,N,Y)
SOLUTN(S5CN,Y,N,N,N,N,Y)
SOLUTN(L3CN,Y,N,N,N,N,Y)
SOLUTN(COND,Y,N,N,N,N,Y)
SOLUTN(LFCP,Y,N,N,N,N,Y)
SOLUTN(S5CP,Y,N,N,N,N,Y)
SOLUTN(L3CP,Y,N,N,N,N,Y)
SOLUTN(CP1,Y,N,N,N,N,Y)
SOLUTN(LFEN,Y,N,N,N,N,Y)
SOLUTN(S5EN,Y,N,N,N,N,Y)
SOLUTN(L3EN,Y,N,N,N,N,Y)
SOLUTN(ENUL,Y,N,N,N,N,Y)
SOLUTN(LFRH,Y,N,N,N,N,Y)
SOLUTN(S5RH,Y,N,N,N,N,Y)
SOLUTN(L3RH,Y,N,N,N,N,Y)
SOLUTN(RHO1,Y,N,N,N,N,Y)
DEN1 = 150
VISL = 146
************************************************************
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(P1,Y,Y,Y,N,Y,Y)
TERMS(U1,Y,Y,Y,Y,Y,Y)
TERMS(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(TEMP,N,Y,Y,Y,Y,Y)
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.
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(W1)=1. ;PRNDTL(TEMP)=0.7
PRT(U1)=1. ;PRT(V1)=1.
PRT(W1)=1. ;PRT(TEMP)=1.
CP1 =1. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10
FIINIT(TEMP)=350. ;FIINIT(ENT1)=1.0E-10
FIINIT(LFCN)=1.0E-10 ;FIINIT(S5CN)=1.0E-10
FIINIT(L3CN)=1.0E-10 ;FIINIT(COND)=1.0E-10
FIINIT(LFCP)=1.0E-10 ;FIINIT(S5CP)=1.0E-10
FIINIT(L3CP)=1.0E-10 ;FIINIT(CP1)=1.0E-10
FIINIT(LFEN)=1.0E-10 ;FIINIT(S5EN)=1.0E-10
FIINIT(L3EN)=1.0E-10 ;FIINIT(ENUL)=1.0E-10
FIINIT(LFRH)=1.0E-10 ;FIINIT(S5RH)=1.0E-10
FIINIT(L3RH)=1.0E-10 ;FIINIT(RHO1)=1.0E-10
No PATCHes yet used for this Group
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
PATCH(SOUTH ,SWALL , 1, 10, 1, 1, 1, 10, 1, 1)
COVAL(SOUTH ,U1 ,1. ,0. )
COVAL(SOUTH ,W1 ,1. ,0. )
COVAL(SOUTH ,TEMP,1. ,280. )
PATCH(MOVING ,NWALL , 1, 10, 10, 10, 1, 10, 1, 1)
COVAL(MOVING ,U1 ,1. ,1.0E-02 )
COVAL(MOVING ,W1 ,1. ,0. )
COVAL(MOVING ,TEMP,1. ,420. )
PATCH(WEST ,WWALL , 1, 1, 1, 10, 1, 10, 1, 1)
COVAL(WEST ,V1 ,1. ,0. )
COVAL(WEST ,W1 ,1. ,0. )
COVAL(WEST ,TEMP,1. ,280. )
PATCH(EAST ,EWALL , 10, 10, 1, 10, 1, 10, 1, 1)
COVAL(EAST ,V1 ,1. ,0. )
COVAL(EAST ,W1 ,1. ,0. )
COVAL(EAST ,TEMP,1. ,280. )
PATCH(HIGH ,HWALL , 1, 10, 1, 10, 10, 10, 1, 1)
COVAL(HIGH ,U1 ,1. ,0. )
COVAL(HIGH ,V1 ,1. ,0. )
PATCH(LOW ,LWALL , 1, 10, 1, 10, 1, 1, 1, 1)
COVAL(LOW ,U1 ,1. ,0. )
COVAL(LOW ,V1 ,1. ,0. )
PATCH(RELIEF ,CELL , 5, 5, 5, 5, 1, 1, 1, 1)
COVAL(RELIEF ,P1 ,1. ,0. )
COVAL(RELIEF ,TEMP,0. , SAME )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 250 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-06
************************************************************
Group 16. Terminate Iterations
LITER(P1)=100 ;LITER(U1)=20
LITER(V1)=20 ;LITER(W1)=20
LITER(TEMP)=20
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
ENDIT(TEMP)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,FALSDT,0.1)
RELAX(V1,FALSDT,0.1)
RELAX(W1,LINRLX,0.1)
RELAX(TEMP,FALSDT,1.0E+09)
RELAX(ENT1,LINRLX,1.)
RELAX(LFCN,LINRLX,1.)
RELAX(S5CN,LINRLX,1.)
RELAX(L3CN,LINRLX,1.)
RELAX(COND,LINRLX,1.)
RELAX(LFCP,LINRLX,1.)
RELAX(S5CP,LINRLX,1.)
RELAX(L3CP,LINRLX,1.)
RELAX(CP1,LINRLX,1.)
RELAX(LFEN,LINRLX,1.)
RELAX(S5EN,LINRLX,1.)
RELAX(L3EN,LINRLX,1.)
RELAX(ENUL,LINRLX,1.)
RELAX(LFRH,LINRLX,1.)
RELAX(S5RH,LINRLX,1.)
RELAX(L3RH,LINRLX,1.)
RELAX(RHO1,LINRLX,1.)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
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(TEMP)=1.0E+10 ;VARMIN(TEMP)=-1.0E+10
VARMAX(ENT1)=1.0E+10 ;VARMIN(ENT1)=-1.0E+10
VARMAX(LFCN)=1.0E+10 ;VARMIN(LFCN)=-1.0E+10
VARMAX(S5CN)=1.0E+10 ;VARMIN(S5CN)=-1.0E+10
VARMAX(L3CN)=1.0E+10 ;VARMIN(L3CN)=-1.0E+10
VARMAX(COND)=1.0E+10 ;VARMIN(COND)=-1.0E+10
VARMAX(LFCP)=1.0E+10 ;VARMIN(LFCP)=-1.0E+10
VARMAX(S5CP)=1.0E+10 ;VARMIN(S5CP)=-1.0E+10
VARMAX(L3CP)=1.0E+10 ;VARMIN(L3CP)=-1.0E+10
VARMAX(CP1)=1.0E+10 ;VARMIN(CP1)=-1.0E+10
VARMAX(LFEN)=1.0E+10 ;VARMIN(LFEN)=-1.0E+10
VARMAX(S5EN)=1.0E+10 ;VARMIN(S5EN)=-1.0E+10
VARMAX(L3EN)=1.0E+10 ;VARMIN(L3EN)=-1.0E+10
VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10
VARMAX(LFRH)=1.0E+10 ;VARMIN(LFRH)=-1.0E+10
VARMAX(S5RH)=1.0E+10 ;VARMIN(S5RH)=-1.0E+10
VARMAX(L3RH)=1.0E+10 ;VARMIN(L3RH)=-1.0E+10
VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
CONWIZ = T
ISG50 = 1
ISG52 = 2
ISG62 = 1
CSG3 ='CNGR'
SPEDAT(SET,LONGNAME,L3RH,C,3-Piece-Wise_Linear_Density/rho1)
SPEDAT(SET,LONGNAME,L3EN,C,3-Piece-Wise_Linear_Viscosity/enul)
SPEDAT(SET,LONGNAME,L3CP,C,3-Piece-Wise_Linear_Specific_Heat/cp1)
SPEDAT(SET,LONGNAME,L3CN,C,3-Piece-Wise_Linear_Conductivity/cond)
SPEDAT(SET,LONGNAME,S5RH,C,Spline_Curve_Density/rho1)
SPEDAT(SET,LONGNAME,S5EN,C,Spline_Curve_Viscosity/enul)
SPEDAT(SET,LONGNAME,S5CP,C,Spline_Curve_Specific_Heat/cp1)
SPEDAT(SET,LONGNAME,S5CN,C,Spline_Curve_Conductivity/cond)
SPEDAT(SET,LONGNAME,LFRH,C,Piece-Wise_Linear_Density_From_File/rho1)
SPEDAT(SET,LONGNAME,LFEN,C,Piece-Wise_Linear_Visc_From_File/enul)
SPEDAT(SET,LONGNAME,LFCP,C,Piece-Wise_Linear_Spe_Heat_From_File/cp1)
SPEDAT(SET,LONGNAME,LFCN,C,Piece-Wise_Linear_Cond_From_File/cond)
SPEDAT(SET,LONGNAME,TEMP,C,The_H1_Store_Is_Used_For_This)
SPEDAT(SET,LONGNAME,ENT1,C,The_Enthalpy_I.E._CP*temp)
SPEDAT(SET,LONGNAME,ENUL,C,The_Viscosity_Actually_Used)
SPEDAT(SET,LONGNAME,RHO1,C,The_Density_Actually_Used)
SPEDAT(SET,LONGNAME,CP1,C,The_Specific_Heat_Actually_Used)
SPEDAT(SET,LONGNAME,COND,C,The_Conductivity_Actually_Used)
SPEDAT(SET,STORED,ENT1,C,=TEMP*CP1!ZSLFIN)
SPEDAT(SET,PROPERTY,RHO1,C,=POL3(TEMP&1082.22&1.34031&-5.63885E-3$)
SPEDAT(SET,PROPERTY,RHO1,C,&4.97135E-6))
SPEDAT(SET,PROPERTY,ENUL,C,=10.^(POL3(TEMP&21.5824&-.149&3.40615E$)
SPEDAT(SET,PROPERTY,ENUL,C,-4&-2.71336E-7)-2)/RHO1)
SPEDAT(SET,PROPERTY,CP1,C,=POL3(TEMP&889.59&5.27817&-9.16533E-4&8$)
SPEDAT(SET,PROPERTY,CP1,C,.80859E-7))
SPEDAT(SET,STORED,COND,C,=.001*POL3(TEMP&526.92&-.88347&3.1519E-4$)
SPEDAT(SET,STORED,COND,C,&-2.42873E-7)!ZSLFIN)
SPEDAT(SET,PROPERTY,CN014,C,=COND/CP1)
SPEDAT(SET,STORED,L3RH,C,=PWL3(TEMP&280.&1125.&320.&1096.&370.&10$)
SPEDAT(SET,STORED,L3RH,C,60.&420.&1019.)/RHO1!ZSLFIN)
SPEDAT(SET,STORED,L3EN,C,=PWL3(TEMP&280.&.0409&320.&.00781&370.&.$)
SPEDAT(SET,STORED,L3EN,C,00268&420.&.00097)/(RHO1*ENUL)!ZSLFIN)
SPEDAT(SET,STORED,L3CP,C,=PWL3(TEMP&280.&2315.&320.&2514.&370.&27$)
SPEDAT(SET,STORED,L3CP,C,62.&420.&3010.)/CP1!ZSLFIN)
SPEDAT(SET,STORED,L3CN,C,=PWL3(TEMP&280.&.2989&320.&.2685&370.&.2$)
SPEDAT(SET,STORED,L3CN,C,309&420.&.1935)/COND!ZSLFIN)
SPEDAT(SET,STORED,S5RH,C,=SPL5(TEMP&280.&1125.&310.&1104.&340.&10$)
SPEDAT(SET,STORED,S5RH,C,82.&380.&1051.&420.&1019.)/RHO1!ZSLFIN)
SPEDAT(SET,STORED,S5EN,C,=SPL5(TEMP&280.&.0409&310.&.0107&340.&.0$)
SPEDAT(SET,STORED,S5EN,C,0428&380.&.00182&420.&.00097)/(RHO1*ENUL$)
SPEDAT(SET,STORED,S5EN,C,)!ZSLFIN)
SPEDAT(SET,STORED,S5CP,C,=SPL5(TEMP&280.&2315.&310.&2464.&340.&26$)
SPEDAT(SET,STORED,S5CP,C,13.&380.&2811.&420.&3010.)/CP1!ZSLFIN)
SPEDAT(SET,STORED,S5CN,C,=SPL5(TEMP&280.&.2989&310.&.2761&340.&.2$)
SPEDAT(SET,STORED,S5CN,C,534&380.&.2234&420.&.1935)/COND!ZSLFIN)
SPEDAT(SET,STORED,LFRH,C,=PWLF(/PHOENICS/D_EARTH/D_CORE/INPLIB/DE$)
SPEDAT(SET,STORED,LFRH,C,NPRP&TEMP)/RHO1!ZSLFIN)
SPEDAT(SET,STORED,LFEN,C,=PWLF(/PHOENICS/D_EARTH/D_CORE/INPLIB/EN$)
SPEDAT(SET,STORED,LFEN,C,UPRP&TEMP)/(RHO1*ENUL)!ZSLFIN)
SPEDAT(SET,STORED,LFCP,C,=PWLF(/PHOENICS/D_EARTH/D_CORE/INPLIB/CP$)
SPEDAT(SET,STORED,LFCP,C,PRP&TEMP)/CP1!ZSLFIN)
SPEDAT(SET,STORED,LFCN,C,=PWLF(/PHOENICS/D_EARTH/D_CORE/INPLIB/CN$)
SPEDAT(SET,STORED,LFCN,C,DPRP&TEMP)/COND!ZSLFIN)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=3.798E-03 ;EX(U1)=6.416E-04
EX(V1)=1.697E-04 ;EX(W1)=1.482E-05
EX(TEMP)=301.299988 ;EX(ENT1)=7.36E+05
EX(LFCN)=0.9999 ;EX(S5CN)=1.
EX(L3CN)=1. ;EX(COND)=0.2828
EX(LFCP)=1. ;EX(S5CP)=1.
EX(L3CP)=1. ;EX(CP1)=2421.
EX(LFEN)=1.013 ;EX(S5EN)=1.02
EX(L3EN)=1.149 ;EX(ENUL)=2.479E-05
EX(LFRH)=1. ;EX(S5RH)=1.
EX(L3RH)=1. ;EX(RHO1)=1109.
************************************************************
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(P1,Y,N,Y,Y,Y,Y)
OUTPUT(U1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(TEMP,Y,N,Y,Y,Y,Y)
OUTPUT(ENT1,Y,N,Y,N,N,N)
OUTPUT(LFCN,Y,N,Y,N,N,N)
OUTPUT(S5CN,Y,N,Y,N,N,N)
OUTPUT(L3CN,Y,N,Y,N,N,N)
OUTPUT(COND,Y,N,Y,N,N,N)
OUTPUT(LFCP,Y,N,Y,N,N,N)
OUTPUT(S5CP,Y,N,Y,N,N,N)
OUTPUT(L3CP,Y,N,Y,N,N,N)
OUTPUT(CP1,Y,N,Y,N,N,N)
OUTPUT(LFEN,Y,N,Y,N,N,N)
OUTPUT(S5EN,Y,N,Y,N,N,N)
OUTPUT(L3EN,Y,N,Y,N,N,N)
OUTPUT(ENUL,Y,N,Y,N,N,N)
OUTPUT(LFRH,Y,N,Y,N,N,N)
OUTPUT(S5RH,Y,N,Y,N,N,N)
OUTPUT(L3RH,Y,N,Y,N,N,N)
OUTPUT(RHO1,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 5 ;IYMON = 9 ;IZMON = 5
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 = 5 ;IZPRL = 5
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
PATCH(MIDX ,PROFIL, 5, 5, 1, 10, 5, 5, 1, 1)
PLOT(MIDX ,U1 ,0. ,0. )
PLOT(MIDX ,TEMP,0. ,0. )
PATCH(MAP ,CONTUR, 1, 10, 1, 10, 5, 5, 1, 1)
PLOT(MAP ,U1 ,0. ,10. )
PLOT(MAP ,V1 ,0. ,10. )
PLOT(MAP ,W1 ,0. ,10. )
PLOT(MAP ,TEMP,0. ,10. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP