TALK=T;RUN(1,1)
PHOTON USE
p;;;;
vi y
msg(vectors and contours at symmetry plane
vec y 1 sh
pause
msg(shell-side velocity vectors
vec off
msg(shell-side pressure contours
con ps y 1 fi;.001
pause
msg(shell-side temperature contours
con tems y 1 fi;.001
pause
msg(tube-side temperature contours
con temt y 1 fi;.001
pause
ENDUSE
DISPLAY
A 3D shell-and-tube heat-exchanger, with one tube-side pass,
and two shell-side baffles, in steady flow.
The shell-side fluid is saturated water; the tube-side fluid
is SAE_5W-30_engine_oil.
In-Form is used for many purposes.
Note by dbs on 24.08.08.
This file, dated 20.05.05, was probably that which was later
used for the 2005 ASME paper. The latter however was developed
further, not least by the use of MUSES to enable the stresses in
the tubes to be calculated.
The diagram has today been changed because it still referred to
a five-pass exchanger!
Also the nomenclature has been streamlined so thats rather than
1 refers to the shell side and t to the tube side.
ENDDIS
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(3D Shell and (1-pass) Tube Heat Exchange)
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 30
XULAST =1.
XFRAC(1)=0.033333 ;XFRAC(2)=0.066667
XFRAC(3)=0.1 ;XFRAC(4)=0.133333
XFRAC(5)=0.166667 ;XFRAC(6)=0.2
XFRAC(7)=0.233333 ;XFRAC(8)=0.266667
XFRAC(9)=0.3 ;XFRAC(10)=0.333333
XFRAC(11)=0.366667 ;XFRAC(12)=0.4
XFRAC(13)=0.433333 ;XFRAC(14)=0.466667
XFRAC(15)=0.5 ;XFRAC(16)=0.533333
XFRAC(17)=0.566667 ;XFRAC(18)=0.6
XFRAC(19)=0.633333 ;XFRAC(20)=0.666667
XFRAC(21)=0.7 ;XFRAC(22)=0.733333
XFRAC(23)=0.766667 ;XFRAC(24)=0.8
XFRAC(25)=0.833333 ;XFRAC(26)=0.866667
XFRAC(27)=0.9 ;XFRAC(28)=0.933333
XFRAC(29)=0.966667 ;XFRAC(30)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 5
YVLAST =0.5
YFRAC(1)=0.2 ;YFRAC(2)=0.4
YFRAC(3)=0.6 ;YFRAC(4)=0.8
YFRAC(5)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 30
ZWLAST =4.
ZFRAC(1)=0.033333 ;ZFRAC(2)=0.066667
ZFRAC(3)=0.1 ;ZFRAC(4)=0.133333
ZFRAC(5)=0.166667 ;ZFRAC(6)=0.2
ZFRAC(7)=0.233333 ;ZFRAC(8)=0.266667
ZFRAC(9)=0.3 ;ZFRAC(10)=0.333333
ZFRAC(11)=0.366667 ;ZFRAC(12)=0.4
ZFRAC(13)=0.433333 ;ZFRAC(14)=0.466667
ZFRAC(15)=0.5 ;ZFRAC(16)=0.533333
ZFRAC(17)=0.566667 ;ZFRAC(18)=0.6
ZFRAC(19)=0.633333 ;ZFRAC(20)=0.666667
ZFRAC(21)=0.7 ;ZFRAC(22)=0.733333
ZFRAC(23)=0.766667 ;ZFRAC(24)=0.8
ZFRAC(25)=0.833333 ;ZFRAC(26)=0.866667
ZFRAC(27)=0.9 ;ZFRAC(28)=0.933333
ZFRAC(29)=0.966667 ;ZFRAC(30)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=PS ;NAME(3)=U1
NAME(5)=V1 ;NAME(7)=W1
NAME(14)=HS ;NAME(124)=ENUM
NAME(125)=FLUX ;NAME(126)=CONT
NAME(127)=CPT ;NAME(128)=RHOT
NAME(129)=ENUP ;NAME(130)=COND
NAME(131)=CPS ;NAME(132)=RHO1
NAME(133)=ENUL ;NAME(134)=HPOR
NAME(135)=NPOR ;NAME(136)=EPOR
NAME(137)=TEMM ;NAME(138)=COEU
NAME(139)=COET ;NAME(140)=PRNT
NAME(141)=NUST ;NAME(142)=REYT
NAME(143)=COES ;NAME(144)=PRNS
NAME(145)=NUSS ;NAME(146)=REYS
NAME(147)=VABS ;NAME(148)=TEMT
NAME(149)=TEMS ;NAME(150)=HT
* Y in SOLUTN argument list denotes:
* 1-stored 2-solved 3-whole-field
* 4-point-by-point 5-explicit 6-harmonic averaging
SOLUTN(PS,Y,Y,Y,N,N,Y)
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(HS,Y,Y,N,N,N,Y)
SOLUTN(ENUM,Y,N,N,N,N,N)
SOLUTN(FLUX,Y,N,N,N,N,N)
SOLUTN(CONT,Y,N,N,N,N,Y)
SOLUTN(CPT,Y,N,N,N,N,Y)
SOLUTN(RHOT,Y,N,N,N,N,Y)
SOLUTN(ENUP,Y,N,N,N,N,Y)
SOLUTN(COND,Y,N,N,N,N,Y)
SOLUTN(CPS,Y,N,N,N,N,Y)
SOLUTN(RHO1,Y,N,N,N,N,Y)
SOLUTN(ENUL,Y,N,N,N,N,Y)
SOLUTN(HPOR,Y,N,N,N,N,Y)
SOLUTN(NPOR,Y,N,N,N,N,Y)
SOLUTN(EPOR,Y,N,N,N,N,Y)
SOLUTN(TEMM,Y,N,N,N,N,Y)
SOLUTN(COEU,Y,N,N,N,N,Y)
SOLUTN(COET,Y,N,N,N,N,Y)
SOLUTN(PRNT,Y,N,N,N,N,Y)
SOLUTN(NUST,Y,N,N,N,N,Y)
SOLUTN(REYT,Y,N,N,N,N,Y)
SOLUTN(COES,Y,N,N,N,N,Y)
SOLUTN(PRNS,Y,N,N,N,N,Y)
SOLUTN(NUSS,Y,N,N,N,N,Y)
SOLUTN(REYS,Y,N,N,N,N,Y)
SOLUTN(VABS,Y,N,N,N,N,Y)
SOLUTN(TEMT,Y,N,N,N,N,Y)
SOLUTN(TEMS,Y,N,N,N,N,Y)
SOLUTN(HT,Y,Y,N,N,N,Y)
DEN1 = 132
VISL = 133
EPOR = 136 ;HPOR = 134 ;NPOR = 135 ;VPOR = 0
************************************************************
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(PS,Y,Y,Y,N,Y,Y)
TERMS(U1,Y,Y,N,Y,Y,Y)
TERMS(V1,Y,Y,N,Y,Y,Y)
TERMS(W1,Y,Y,N,Y,Y,Y)
TERMS(HS,N,Y,N,Y,Y,Y)
TERMS(HT,N,N,N,Y,N,N)
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(PS)=1.0E-10 ;FIINIT(U1)=1.0E-10
FIINIT(V1)=1.0E-10 ;FIINIT(W1)=1.0E-10
FIINIT(HS)=1.171054E+06 ;FIINIT(ENUM)=1.0E-10
FIINIT(FLUX)=1.0E-10 ;FIINIT(CONT)=1.0E-10
FIINIT(CPT)=1.0E-10 ;FIINIT(RHOT)=1.0E-10
FIINIT(ENUP)=1.0E-10 ;FIINIT(COND)=1.0E-10
FIINIT(CPS)=1.0E-10 ;FIINIT(RHO1)=1.0E-10
FIINIT(ENUL)=1.0E-10 ;FIINIT(HPOR)=0.5
FIINIT(NPOR)=0.5 ;FIINIT(EPOR)=0.5
FIINIT(TEMM)=1.0E-10 ;FIINIT(COEU)=1.0E-10
FIINIT(COET)=1.0E-10 ;FIINIT(PRNT)=1.0E-10
FIINIT(NUST)=1.0E-10 ;FIINIT(REYT)=1.0E-10
FIINIT(COES)=1.0E-10 ;FIINIT(PRNS)=1.0E-10
FIINIT(NUSS)=1.0E-10 ;FIINIT(REYS)=1.0E-10
FIINIT(VABS)=1.0E-10 ;FIINIT(TEMT)=353.
FIINIT(TEMS)=283. ;FIINIT(HT)=6.7776E+05
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(INLET1 ,CELL , 1, 1, 3, 5, 1, 3, 1, 1)
COVAL(INLET1 ,PS , FIXFLU ,5.555555 )
COVAL(INLET1 ,HS ,0. ,1.171054E+06 )
PATCH(OUTLET1 ,EAST , 30, 30, 3, 5, 28, 30, 1, 1)
COVAL(OUTLET1 ,PS ,1. ,0. )
PATCH(INLET2 ,CELL , 1, 30, 1, 5, 30, 30, 1, 1000)
COVAL(INLET2 ,HT ,0.166667 ,6.7776E+05 )
PATCH(NEHS ,CELL , 1, 30, 1, 5, 1, 29, 1, 1)
COVAL(NEHS ,HT ,0.166667 , GRND8 )
PATCH(HEATEXCH,VOLUME, 1, 30, 1, 5, 1, 30, 1, 1000)
COVAL(HEATEXCH,HS ,In-Form:source - see Grp 19)
COVAL(HEATEXCH,HT ,In-Form:source - see Grp 19)
PATCH(BAFFLE1 ,HIGH , 1, 18, 1, 5, 10, 10, 1, 1000)
COVAL(BAFFLE1 ,W1 , FIXVAL ,0. )
PATCH(BAFFLE2 ,HIGH , 13, 30, 1, 5, 20, 20, 1, 1000)
COVAL(BAFFLE2 ,W1 , FIXVAL ,0. )
PATCH(RESIST ,PHASEM, 1, 30, 1, 5, 1, 30, 1, 1000)
COVAL(RESIST ,U1 ,100. ,0. )
COVAL(RESIST ,V1 ,100. ,0. )
COVAL(RESIST ,W1 ,50. ,0. )
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 200 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(PS)=100 ;LITER(U1)=10
LITER(V1)=10 ;LITER(W1)=10
LITER(HS)=20 ;LITER(HT)=20
ENDIT(PS)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
ENDIT(HS)=1.0E-03 ;ENDIT(HT)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(PS,LINRLX,0.5)
RELAX(U1,FALSDT,0.1)
RELAX(V1,FALSDT,0.1)
RELAX(W1,FALSDT,0.1)
RELAX(HS,FALSDT,1.0E+09)
RELAX(ENUM,LINRLX,1.)
RELAX(FLUX,LINRLX,1.)
RELAX(CONT,LINRLX,1.)
RELAX(CPT,LINRLX,1.)
RELAX(RHOT,LINRLX,1.)
RELAX(ENUP,LINRLX,1.)
RELAX(COND,LINRLX,1.)
RELAX(CPS,LINRLX,1.)
RELAX(RHO1,LINRLX,1.)
RELAX(ENUL,LINRLX,1.)
RELAX(TEMM,LINRLX,1.)
RELAX(COEU,LINRLX,1.)
RELAX(COET,LINRLX,1.)
RELAX(PRNT,LINRLX,1.)
RELAX(NUST,LINRLX,1.)
RELAX(REYT,LINRLX,1.)
RELAX(COES,LINRLX,1.)
RELAX(PRNS,LINRLX,1.)
RELAX(NUSS,LINRLX,1.)
RELAX(REYS,LINRLX,1.)
RELAX(VABS,LINRLX,1.)
RELAX(TEMT,LINRLX,1.)
RELAX(TEMS,LINRLX,1.)
RELAX(HT,FALSDT,1.0E+09)
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(PS)=1.0E+10 ;VARMIN(PS)=-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(HS)=1.0E+10 ;VARMIN(HS)=-1.0E+10
VARMAX(ENUM)=1.0E+10 ;VARMIN(ENUM)=-1.0E+10
VARMAX(FLUX)=1.0E+10 ;VARMIN(FLUX)=-1.0E+10
VARMAX(CONT)=1.0E+10 ;VARMIN(CONT)=-1.0E+10
VARMAX(CPT)=1.0E+10 ;VARMIN(CPT)=-1.0E+10
VARMAX(RHOT)=5000. ;VARMIN(RHOT)=500.
VARMAX(ENUP)=1.0E+10 ;VARMIN(ENUP)=-1.0E+10
VARMAX(COND)=1.0E+10 ;VARMIN(COND)=-1.0E+10
VARMAX(CPS)=1.0E+10 ;VARMIN(CPS)=-1.0E+10
VARMAX(RHO1)=5000. ;VARMIN(RHO1)=500.
VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10
VARMAX(HPOR)=1.0E+10 ;VARMIN(HPOR)=-1.0E+10
VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10
VARMAX(EPOR)=1.0E+10 ;VARMIN(EPOR)=-1.0E+10
VARMAX(TEMM)=1.0E+10 ;VARMIN(TEMM)=-1.0E+10
VARMAX(COEU)=1.0E+10 ;VARMIN(COEU)=-1.0E+10
VARMAX(COET)=1.0E+10 ;VARMIN(COET)=-1.0E+10
VARMAX(PRNT)=1.0E+10 ;VARMIN(PRNT)=-1.0E+10
VARMAX(NUST)=1.0E+10 ;VARMIN(NUST)=-1.0E+10
VARMAX(REYT)=1.0E+10 ;VARMIN(REYT)=-1.0E+10
VARMAX(COES)=1.0E+10 ;VARMIN(COES)=-1.0E+10
VARMAX(PRNS)=1.0E+10 ;VARMIN(PRNS)=-1.0E+10
VARMAX(NUSS)=1.0E+10 ;VARMIN(NUSS)=-1.0E+10
VARMAX(REYS)=1.0E+10 ;VARMIN(REYS)=-1.0E+10
VARMAX(VABS)=1.0E+10 ;VARMIN(VABS)=-1.0E+10
VARMAX(TEMT)=1.0E+10 ;VARMIN(TEMT)=-1.0E+10
VARMAX(TEMS)=1.0E+10 ;VARMIN(TEMS)=-1.0E+10
VARMAX(HT)=1.0E+10 ;VARMIN(HT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
PARSOL = F
ISG50 = 1
ISG52 = 2
ISG62 = 1
SPEDAT(SET,STORED,TEMS,C,=HS/4138.!ZSLFIN)
SPEDAT(SET,STORED,TEMT,C,=HT/1920.!ZSLFIN)
SPEDAT(SET,STORED,TEMM,C,=(COES*TEMT+COET*TEMS)/(COES+COET)!ZSLFIN)
SPEDAT(SET,STORED,REYS,C,=1.0E-02*VABS/ENUL!ZSLFIN)
SPEDAT(SET,STORED,PRNS,C,=CPS*RHO1*ENUL/COND!ZSLFIN)
SPEDAT(SET,STORED,NUSS,C,=0.2*REYS^0.6*PRNS^0.33!ZSLFIN)
SPEDAT(SET,STORED,COES,C,=201.062378*NUSS*COND/1.0E-02!ZSLFIN)
SPEDAT(SET,STORED,REYT,C,=1.0E-02*0.046393/ENUP!ZSLFIN)
SPEDAT(SET,STORED,PRNT,C,=CPT*RHOT*ENUP/CONT!ZSLFIN)
SPEDAT(SET,STORED,NUST,C,=MAX(2.0&0.328*(REYT*PRNT)^0.33)!ZSLFIN)
SPEDAT(SET,STORED,COET,C,=(ENUM/ENUP)^0.14*201.062378*NUST*CONT/1$)
SPEDAT(SET,STORED,COET,C,.0E-02!ZSLFIN)
SPEDAT(SET,STORED,COEU,C,=COES*COET/(COES+COET)!ZSLFIN)
SPEDAT(SET,STORED,FLUX,C,=COEU*(TEMT-TEMS)!ZSLFIN)
SPEDAT(SET,STORED,COND,C,=.001*POL5(TEMS&62.282&-1.768417&.03499&$)
SPEDAT(SET,STORED,COND,C,-1.15706E-4&1.53599E-7&-7.7477E-11)!ZSLF$)
SPEDAT(SET,STORED,COND,C,IN)
SPEDAT(SET,PROPERTY,RHO1,C,=POL5(TEMS&2446.&-20.6741&.11576&-3.12$)
SPEDAT(SET,PROPERTY,RHO1,C,895E-4&4.0505E-7&-2.054E-10))
SPEDAT(SET,PROPERTY,ENUL,C,=1.E-7*EXP((1.12646-.039638*TEMS)/(1.-$)
SPEDAT(SET,PROPERTY,ENUL,C,7.29769E-3*TEMS))/(RHO1))
SPEDAT(SET,STORED,CPS,C,=EXP((8.29041-.012557*TEMS)/(1.-1.52373E-$)
SPEDAT(SET,STORED,CPS,C,3*TEMS))!ZSLFIN)
SPEDAT(SET,STORED,RHOT,C,=1052.3-0.6420*TEMT!ZSLFIN)
SPEDAT(SET,STORED,ENUP,C,=10.^(POL4(TEMT&58.2987&-.53817&1.92827E$)
SPEDAT(SET,STORED,ENUP,C,-3&-3.16448E-6&1.97922E-9)-2)/(RHOT)!ZSL$)
SPEDAT(SET,STORED,ENUP,C,FIN)
SPEDAT(SET,STORED,CPT,C,=753.7+3.65*TEMT!ZSLFIN)
SPEDAT(SET,STORED,CONT,C,=0.1447-2.3073E-5*TEMT!ZSLFIN)
SPEDAT(SET,STORED,ENUM,C,=10.^(POL4(TEMM&58.2987&-.53817&1.92827E$)
SPEDAT(SET,STORED,ENUM,C,-3&-3.16448E-6&1.97922E-9)-2)/(RHOT)!ZSL$)
SPEDAT(SET,STORED,ENUM,C,FIN)
SPEDAT(SET,SOURCE,HS!HEATEXCH,C,=COEU*(TEMT-TEMS))
SPEDAT(SET,SOURCE,HT!HEATEXCH,C,=COEU*(TEMS-TEMT))
SPEDAT(SET,LONGNAME,HS,C,shell-side_fluid_enthalpy)
SPEDAT(SET,LONGNAME,TEMS,C,shell-side_fluid_temperature)
SPEDAT(SET,LONGNAME,RHO1,C,shell-side_fluid_density)
SPEDAT(SET,LONGNAME,CPS,C,shell-side_fluid_specific_heat)
SPEDAT(SET,LONGNAME,COND,C,shell-side_fluid_thermal_conductivity)
SPEDAT(SET,LONGNAME,ENUL,C,shell-side_fluid_kinematic_viscosity)
SPEDAT(SET,LONGNAME,HT,C,tube-side_fluid_enthalpy)
SPEDAT(SET,LONGNAME,TEMT,C,tube-side_fluid_temperature)
SPEDAT(SET,LONGNAME,RHOT,C,tube-side_fluid_density)
SPEDAT(SET,LONGNAME,CPT,C,tube-side_fluid_specific_heat)
SPEDAT(SET,LONGNAME,CONT,C,tube-side_fluid_thermal_conductivity)
SPEDAT(SET,LONGNAME,ENUP,C,tube-side_fluid_kinematic_viscosity)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
************************************************************
Group 20. Preliminary Printout
************************************************************
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(PS,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(HS,Y,N,Y,Y,Y,Y)
OUTPUT(ENUM,Y,N,Y,N,N,N)
OUTPUT(FLUX,Y,N,Y,N,N,N)
OUTPUT(CONT,Y,N,Y,N,N,N)
OUTPUT(CPT,Y,N,Y,N,N,N)
OUTPUT(RHOT,Y,N,Y,N,N,N)
OUTPUT(ENUP,Y,N,Y,N,N,N)
OUTPUT(COND,Y,N,Y,N,N,N)
OUTPUT(CPS,Y,N,Y,N,N,N)
OUTPUT(RHO1,Y,N,Y,N,N,N)
OUTPUT(ENUL,Y,N,Y,N,N,N)
OUTPUT(HPOR,N,N,N,N,N,N)
OUTPUT(NPOR,N,N,N,N,N,N)
OUTPUT(EPOR,N,N,N,N,N,N)
OUTPUT(TEMM,Y,N,Y,N,N,N)
OUTPUT(COEU,Y,N,Y,N,N,N)
OUTPUT(COET,Y,N,Y,N,N,N)
OUTPUT(PRNT,Y,N,Y,N,N,N)
OUTPUT(NUST,Y,N,Y,N,N,N)
OUTPUT(REYT,Y,N,Y,N,N,N)
OUTPUT(COES,Y,N,Y,N,N,N)
OUTPUT(PRNS,Y,N,Y,N,N,N)
OUTPUT(NUSS,Y,N,Y,N,N,N)
OUTPUT(REYS,Y,N,Y,N,N,N)
OUTPUT(VABS,Y,N,Y,N,N,N)
OUTPUT(TEMT,Y,N,Y,N,N,N)
OUTPUT(TEMS,Y,N,Y,N,N,N)
OUTPUT(HT,Y,N,Y,Y,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 5 ;IYMON = 5 ;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 = T ;YZPR = F
IPLTF = 1 ;IPLTL = 200 ;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(TABMAP ,CONTUR, 1, 30, 1, 5, 15, 15, 1, 1000)
PLOT(TABMAP ,HS ,0. ,10. )
PLOT(TABMAP ,HT ,0. ,10. )
PATCH(TABYEQ3 ,CONTUR, 1, 30, 3, 3, 1, 30, 1, 1000)
PLOT(TABYEQ3 ,PS ,0. ,10. )
PLOT(TABYEQ3 ,HS ,0. ,10. )
PLOT(TABYEQ3 ,HT ,0. ,10. )
************************************************************
Group 24. Dumps For Restarts
SAVE = T ;NOWIPE = F
NSAVE =CHAM
STOP