TALK=T;RUN(1,1)
DISPLAY
Library Case I405: Deposition, downpipe, Re=9.894E3
---------------------------------------------------
The case considered is 2d steady, isothermal, turbulent flow
of air and particles down a pipe with particle deposition on
the wall of the pipe. The Eulerian drift flux model is used to
simulate particle transport with deposition by means of the
3-layer model of Bin & Zhou(2006), which accounts for deposition
via gravity, Brownian and turbulent diffusion, and turbophoresis.
This case has been studied experimentally by Liu & Agrawal(1974),
who for fully-developed flow, measured deposition rates of partic$
les
to the walls of the pipe with nominal sizes ranging from 1.4 to 21
microns in a 0.0127m pipe at air speeds of 11.85m/s and 60 m/s,
corresponding to Reynolds numbers of 9,894 & 50,000, respectively.
The present Q1 file considers the air speed of 11.85 m/s at a
temperature of 22degC so as to replicate the experimental deposit$
ion
of 14 different particles sizes onto the pipe wall in the
fully-developed region of the flow.
ENDDIS
PHOTON USE
p
up -z
view y
con C6 x 1 fi;.1;pa
vec x 1;pa
ENDUSE
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(I405: Deposition, downpipe, Re=9.894E3 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = F
NX = 1
XULAST =0.1
XFRAC(1)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 10
YVLAST =6.35E-03
YFRAC(1)=0.049252 ;YFRAC(2)=0.105892
YFRAC(3)=0.171028 ;YFRAC(4)=0.245934
YFRAC(5)=0.332076 ;YFRAC(6)=0.43114
YFRAC(7)=0.545063 ;YFRAC(8)=0.676074
YFRAC(9)=0.826737 ;YFRAC(10)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 80
ZWLAST =1.02
ZFRAC(1)=0.0125 ;ZFRAC(2)=0.025
ZFRAC(3)=0.0375 ;ZFRAC(4)=0.05
ZFRAC(5)=0.0625 ;ZFRAC(6)=0.075
ZFRAC(7)=0.0875 ;ZFRAC(8)=0.1
ZFRAC(9)=0.1125 ;ZFRAC(10)=0.125
ZFRAC(11)=0.1375 ;ZFRAC(12)=0.15
ZFRAC(13)=0.1625 ;ZFRAC(14)=0.175
ZFRAC(15)=0.1875 ;ZFRAC(16)=0.2
ZFRAC(17)=0.2125 ;ZFRAC(18)=0.225
ZFRAC(19)=0.2375 ;ZFRAC(20)=0.25
ZFRAC(21)=0.2625 ;ZFRAC(22)=0.275
ZFRAC(23)=0.2875 ;ZFRAC(24)=0.3
ZFRAC(25)=0.3125 ;ZFRAC(26)=0.325
ZFRAC(27)=0.3375 ;ZFRAC(28)=0.35
ZFRAC(29)=0.3625 ;ZFRAC(30)=0.375
ZFRAC(31)=0.3875 ;ZFRAC(32)=0.4
ZFRAC(33)=0.4125 ;ZFRAC(34)=0.425
ZFRAC(35)=0.4375 ;ZFRAC(36)=0.45
ZFRAC(37)=0.4625 ;ZFRAC(38)=0.475
ZFRAC(39)=0.4875 ;ZFRAC(40)=0.5
ZFRAC(41)=0.5125 ;ZFRAC(42)=0.525
ZFRAC(43)=0.5375 ;ZFRAC(44)=0.55
ZFRAC(45)=0.5625 ;ZFRAC(46)=0.575
ZFRAC(47)=0.5875 ;ZFRAC(48)=0.6
ZFRAC(49)=0.6125 ;ZFRAC(50)=0.625
ZFRAC(51)=0.6375 ;ZFRAC(52)=0.65
ZFRAC(53)=0.6625 ;ZFRAC(54)=0.675
ZFRAC(55)=0.6875 ;ZFRAC(56)=0.7
ZFRAC(57)=0.7125 ;ZFRAC(58)=0.725
ZFRAC(59)=0.7375 ;ZFRAC(60)=0.75
ZFRAC(61)=0.7625 ;ZFRAC(62)=0.775
ZFRAC(63)=0.7875 ;ZFRAC(64)=0.8
ZFRAC(65)=0.8125 ;ZFRAC(66)=0.825
ZFRAC(67)=0.8375 ;ZFRAC(68)=0.85
ZFRAC(69)=0.8625 ;ZFRAC(70)=0.875
ZFRAC(71)=0.8875 ;ZFRAC(72)=0.9
ZFRAC(73)=0.9125 ;ZFRAC(74)=0.925
ZFRAC(75)=0.9375 ;ZFRAC(76)=0.95
ZFRAC(77)=0.9625 ;ZFRAC(78)=0.975
ZFRAC(79)=0.9875 ;ZFRAC(80)=1.
************************************************************
Group 6. Body-Fitted Coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
ONEPHS = T
NAME(1)=P1 ;NAME(5)=V1
NAME(7)=W1 ;NAME(12)=KE
NAME(13)=EP ;NAME(21)=C6
NAME(22)=C7 ;NAME(23)=C8
NAME(24)=C9 ;NAME(25)=C10
NAME(26)=C11 ;NAME(27)=C12
NAME(28)=C13 ;NAME(29)=C14
NAME(30)=C15 ;NAME(31)=C16
NAME(32)=C17 ;NAME(33)=C18
NAME(34)=C19 ;NAME(35)=C20
NAME(36)=C21 ;NAME(37)=C22
NAME(38)=C23 ;NAME(46)=VP18
NAME(47)=TP18 ;NAME(49)=VP17
NAME(50)=TP17 ;NAME(54)=VP16
NAME(55)=TP16 ;NAME(59)=VP15
NAME(60)=TP15 ;NAME(64)=VP14
NAME(65)=TP14 ;NAME(69)=VP13
NAME(70)=TP13 ;NAME(72)=TR13
NAME(74)=VP12 ;NAME(75)=TP12
NAME(79)=VP11 ;NAME(80)=TP11
NAME(84)=VP10 ;NAME(85)=TP10
NAME(89)=VP9 ;NAME(90)=TP9
NAME(94)=VP8 ;NAME(95)=TP8
NAME(96)=VD8 ;NAME(97)=TR8
NAME(99)=VP7 ;NAME(100)=TP7
NAME(104)=VP6 ;NAME(105)=TP6
NAME(108)=VP5 ;NAME(109)=TP5
NAME(110)=VP4 ;NAME(111)=TP4
NAME(112)=VP3 ;NAME(113)=TP3
NAME(114)=VP2 ;NAME(115)=TP2
NAME(116)=VP1 ;NAME(117)=TP1
NAME(121)=DEP5 ;NAME(125)=DEP4
NAME(129)=DEP3 ;NAME(133)=DEP2
NAME(137)=VSTR ;NAME(140)=DEP1
NAME(144)=ENUL ;NAME(145)=YPLS
NAME(146)=STRS ;NAME(147)=EPKE
NAME(148)=DEN1 ;NAME(149)=EL1
NAME(150)=ENUT
* 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,Y)
SOLUTN(V1,Y,Y,N,N,N,Y)
SOLUTN(W1,Y,Y,N,N,N,Y)
SOLUTN(KE,Y,Y,N,N,N,N)
SOLUTN(EP,Y,Y,N,N,N,N)
SOLUTN(C6,Y,Y,Y,N,N,Y)
SOLUTN(C7,Y,Y,Y,N,N,Y)
SOLUTN(C8,Y,Y,Y,N,N,Y)
SOLUTN(C9,Y,Y,Y,N,N,Y)
SOLUTN(C10,Y,Y,Y,N,N,Y)
SOLUTN(C11,Y,Y,Y,N,N,Y)
SOLUTN(C12,Y,Y,Y,N,N,Y)
SOLUTN(C13,Y,Y,Y,N,N,Y)
SOLUTN(C14,Y,Y,Y,N,N,Y)
SOLUTN(C15,Y,Y,Y,N,N,Y)
SOLUTN(C16,Y,Y,Y,N,N,Y)
SOLUTN(C17,Y,Y,Y,N,N,Y)
SOLUTN(C18,Y,Y,Y,N,N,Y)
SOLUTN(C19,Y,Y,Y,N,N,Y)
SOLUTN(C20,Y,Y,Y,N,N,Y)
SOLUTN(C21,Y,Y,Y,N,N,Y)
SOLUTN(C22,Y,Y,Y,N,N,Y)
SOLUTN(C23,Y,Y,Y,N,N,Y)
SOLUTN(VP18,Y,N,N,N,N,Y)
SOLUTN(TP18,Y,N,N,N,N,Y)
SOLUTN(VP17,Y,N,N,N,N,Y)
SOLUTN(TP17,Y,N,N,N,N,Y)
SOLUTN(VP16,Y,N,N,N,N,Y)
SOLUTN(TP16,Y,N,N,N,N,Y)
SOLUTN(VP15,Y,N,N,N,N,Y)
SOLUTN(TP15,Y,N,N,N,N,Y)
SOLUTN(VP14,Y,N,N,N,N,Y)
SOLUTN(TP14,Y,N,N,N,N,Y)
SOLUTN(VP13,Y,N,N,N,N,Y)
SOLUTN(TP13,Y,N,N,N,N,Y)
SOLUTN(TR13,Y,N,N,N,N,Y)
SOLUTN(VP12,Y,N,N,N,N,Y)
SOLUTN(TP12,Y,N,N,N,N,Y)
SOLUTN(VP11,Y,N,N,N,N,Y)
SOLUTN(TP11,Y,N,N,N,N,Y)
SOLUTN(VP10,Y,N,N,N,N,Y)
SOLUTN(TP10,Y,N,N,N,N,Y)
SOLUTN(VP9,Y,N,N,N,N,Y)
SOLUTN(TP9,Y,N,N,N,N,Y)
SOLUTN(VP8,Y,N,N,N,N,Y)
SOLUTN(TP8,Y,N,N,N,N,Y)
SOLUTN(VD8,Y,N,N,N,N,Y)
SOLUTN(TR8,Y,N,N,N,N,Y)
SOLUTN(VP7,Y,N,N,N,N,Y)
SOLUTN(TP7,Y,N,N,N,N,Y)
SOLUTN(VP6,Y,N,N,N,N,Y)
SOLUTN(TP6,Y,N,N,N,N,Y)
SOLUTN(VP5,Y,N,N,N,N,Y)
SOLUTN(TP5,Y,N,N,N,N,Y)
SOLUTN(VP4,Y,N,N,N,N,Y)
SOLUTN(TP4,Y,N,N,N,N,Y)
SOLUTN(VP3,Y,N,N,N,N,Y)
SOLUTN(TP3,Y,N,N,N,N,Y)
SOLUTN(VP2,Y,N,N,N,N,Y)
SOLUTN(TP2,Y,N,N,N,N,Y)
SOLUTN(VP1,Y,N,N,N,N,Y)
SOLUTN(TP1,Y,N,N,N,N,Y)
SOLUTN(DEP5,Y,N,N,N,N,Y)
SOLUTN(DEP4,Y,N,N,N,N,Y)
SOLUTN(DEP3,Y,N,N,N,N,Y)
SOLUTN(DEP2,Y,N,N,N,N,Y)
SOLUTN(VSTR,Y,N,N,N,N,Y)
SOLUTN(DEP1,Y,N,N,N,N,Y)
SOLUTN(ENUL,Y,N,N,N,N,Y)
SOLUTN(YPLS,Y,N,N,N,N,Y)
SOLUTN(STRS,Y,N,N,N,N,Y)
SOLUTN(EPKE,Y,N,N,N,N,Y)
SOLUTN(DEN1,Y,N,N,N,N,Y)
SOLUTN(EL1,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
DEN1 = 148
VISL = 144
VIST = 150
LEN1 = 149
************************************************************
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(V1,Y,Y,Y,Y,Y,Y)
TERMS(W1,Y,Y,Y,Y,Y,Y)
TERMS(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
TERMS(C6,N,Y,Y,Y,N,Y)
TERMS(C7,N,Y,Y,Y,Y,Y)
TERMS(C8,N,Y,Y,Y,N,Y)
TERMS(C9,N,Y,Y,Y,Y,Y)
TERMS(C10,N,Y,Y,Y,N,Y)
TERMS(C11,N,Y,Y,Y,Y,Y)
TERMS(C12,N,Y,Y,Y,N,Y)
TERMS(C13,N,Y,Y,Y,Y,Y)
TERMS(C14,N,Y,Y,Y,N,Y)
TERMS(C15,N,Y,Y,Y,Y,Y)
TERMS(C16,N,Y,Y,Y,N,Y)
TERMS(C17,N,Y,Y,Y,Y,Y)
TERMS(C18,N,Y,Y,Y,N,Y)
TERMS(C19,N,Y,Y,Y,Y,Y)
TERMS(C20,N,Y,Y,Y,N,Y)
TERMS(C21,N,Y,Y,Y,Y,Y)
TERMS(C22,N,Y,Y,Y,N,Y)
TERMS(C23,N,Y,Y,Y,Y,Y)
DIFCUT =0.5 ;ZDIFAC =1.
GALA = F ;ADDDIF = F
NEWENT = T
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.198 ;TMP1 = GRND1 ;EL1 = GRND4
TSURR =0. ;TEMP0 =273. ;PRESS0 =1.01325E+05
DVO1DT =3.41E-03 ;DRH1DP =0.
TMP1A =22. ;TMP1B =0. ;TMP1C =0.
TMP2A =0. ;TMP2B =0. ;TMP2C =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =1.524E-05 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(V1)=1. ;PRNDTL(W1)=1.
PRNDTL(KE)=1. ;PRNDTL(EP)=1.
PRNDTL(C6)=1. ;PRNDTL(C7)=1.
PRNDTL(C8)=1. ;PRNDTL(C9)=1.
PRNDTL(C10)=1. ;PRNDTL(C11)=1.
PRNDTL(C12)=1. ;PRNDTL(C13)=1.
PRNDTL(C14)=1. ;PRNDTL(C15)=1.
PRNDTL(C16)=1. ;PRNDTL(C17)=1.
PRNDTL(C18)=1. ;PRNDTL(C19)=1.
PRNDTL(C20)=1. ;PRNDTL(C21)=1.
PRNDTL(C22)=1. ;PRNDTL(C23)=1.
PRT(V1)=1. ;PRT(W1)=1.
PRT(KE)=1. ;PRT(EP)=1.314
PRT(C6)=1. ;PRT(C7)=1.
PRT(C8)=1. ;PRT(C9)=1.
PRT(C10)=1. ;PRT(C11)=1.
PRT(C12)=1. ;PRT(C13)=1.
PRT(C14)=1. ;PRT(C15)=1.
PRT(C16)=1. ;PRT(C17)=1.
PRT(C18)=1. ;PRT(C19)=1.
PRT(C20)=1. ;PRT(C21)=1.
PRT(C22)=1. ;PRT(C23)=1.
CP1 =1005. ;CP2 =1.
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initial field variables (PHIs)
FIINIT(P1)=0. ;FIINIT(V1)=1.0E-10
FIINIT(W1)=11.85 ;FIINIT(KE)=0.351056
FIINIT(EP)=53.825245 ;FIINIT(C6)=1.
FIINIT(C7)=1. ;FIINIT(C8)=1.
FIINIT(C9)=1. ;FIINIT(C10)=1.
FIINIT(C11)=1. ;FIINIT(C12)=1.
FIINIT(C13)=1. ;FIINIT(C14)=1.
FIINIT(C15)=1. ;FIINIT(C16)=1.
FIINIT(C17)=1. ;FIINIT(C18)=1.
FIINIT(C19)=1. ;FIINIT(C20)=1.
FIINIT(C21)=1. ;FIINIT(C22)=1.
FIINIT(C23)=1. ;FIINIT(VP18)=1.0E-10
FIINIT(TP18)=1.0E-10 ;FIINIT(VP17)=1.0E-10
FIINIT(TP17)=1.0E-10 ;FIINIT(VP16)=1.0E-10
FIINIT(TP16)=1.0E-10 ;FIINIT(VP15)=1.0E-10
FIINIT(TP15)=1.0E-10 ;FIINIT(VP14)=1.0E-10
FIINIT(TP14)=1.0E-10 ;FIINIT(VP13)=1.0E-10
FIINIT(TP13)=1.0E-10 ;FIINIT(TR13)=1.0E-10
FIINIT(VP12)=1.0E-10 ;FIINIT(TP12)=1.0E-10
FIINIT(VP11)=1.0E-10 ;FIINIT(TP11)=1.0E-10
FIINIT(VP10)=1.0E-10 ;FIINIT(TP10)=1.0E-10
FIINIT(VP9)=1.0E-10 ;FIINIT(TP9)=1.0E-10
FIINIT(VP8)=1.0E-10 ;FIINIT(TP8)=1.0E-10
FIINIT(VD8)=1.0E-10 ;FIINIT(TR8)=1.0E-10
FIINIT(VP7)=1.0E-10 ;FIINIT(TP7)=1.0E-10
FIINIT(VP6)=1.0E-10 ;FIINIT(TP6)=1.0E-10
FIINIT(VP5)=1.0E-10 ;FIINIT(TP5)=1.0E-10
FIINIT(VP4)=1.0E-10 ;FIINIT(TP4)=1.0E-10
FIINIT(VP3)=1.0E-10 ;FIINIT(TP3)=1.0E-10
FIINIT(VP2)=1.0E-10 ;FIINIT(TP2)=1.0E-10
FIINIT(VP1)=1.0E-10 ;FIINIT(TP1)=1.0E-10
FIINIT(DEP5)=1.0E-10 ;FIINIT(DEP4)=1.0E-10
FIINIT(DEP3)=1.0E-10 ;FIINIT(DEP2)=1.0E-10
FIINIT(VSTR)=1.0E-10 ;FIINIT(DEP1)=1.0E-10
FIINIT(ENUL)=1.0E-10 ;FIINIT(YPLS)=1.0E-10
FIINIT(STRS)=1.0E-10 ;FIINIT(EPKE)=1.0E-10
FIINIT(DEN1)=1.0E-10 ;FIINIT(EL1)=1.0E-10
FIINIT(ENUT)=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(DFLUX ,CELL , 1, 1, 1, 10, 1, 80, 1, 1)
COVAL(DFLUX ,C6 , GRND4 , GRND4 )
COVAL(DFLUX ,C7 , GRND4 , GRND4 )
COVAL(DFLUX ,C8 , GRND4 , GRND4 )
COVAL(DFLUX ,C9 , GRND4 , GRND4 )
COVAL(DFLUX ,C10 , GRND4 , GRND4 )
COVAL(DFLUX ,C11 , GRND4 , GRND4 )
COVAL(DFLUX ,C12 , GRND4 , GRND4 )
COVAL(DFLUX ,C13 , GRND4 , GRND4 )
COVAL(DFLUX ,C14 , GRND4 , GRND4 )
COVAL(DFLUX ,C15 , GRND4 , GRND4 )
COVAL(DFLUX ,C16 , GRND4 , GRND4 )
COVAL(DFLUX ,C17 , GRND4 , GRND4 )
COVAL(DFLUX ,C18 , GRND4 , GRND4 )
COVAL(DFLUX ,C19 , GRND4 , GRND4 )
COVAL(DFLUX ,C20 , GRND4 , GRND4 )
COVAL(DFLUX ,C21 , GRND4 , GRND4 )
COVAL(DFLUX ,C22 , GRND4 , GRND4 )
COVAL(DFLUX ,C23 , GRND4 , GRND4 )
PATCH(KESOURCE,PHASEM, 1, 1, 1, 10, 1, 80, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(OB1 ,NWALL , 1, 1, 10, 10, 1, 80, 1, 1)
COVAL(OB1 ,W1 , GRND2 ,0. )
COVAL(OB1 ,KE , GRND2 , GRND2 )
COVAL(OB1 ,EP , GRND2 , GRND2 )
PATCH(PW1 ,NORTH , 1, 1, 10, 10, 1, 80, 1, 1)
COVAL(PW1 ,C6 , GRND3 ,0. )
COVAL(PW1 ,C7 , GRND3 ,0. )
COVAL(PW1 ,C8 , GRND3 ,0. )
COVAL(PW1 ,C9 , GRND3 ,0. )
COVAL(PW1 ,C10 , GRND3 ,0. )
COVAL(PW1 ,C11 , GRND3 ,0. )
COVAL(PW1 ,C12 , GRND3 ,0. )
COVAL(PW1 ,C13 , GRND3 ,0. )
COVAL(PW1 ,C14 , GRND3 ,0. )
COVAL(PW1 ,C15 , GRND3 ,0. )
COVAL(PW1 ,C16 , GRND3 ,0. )
COVAL(PW1 ,C17 , GRND3 ,0. )
COVAL(PW1 ,C18 , GRND3 ,0. )
COVAL(PW1 ,C19 , GRND3 ,0. )
COVAL(PW1 ,C20 , GRND3 ,0. )
COVAL(PW1 ,C21 , GRND3 ,0. )
COVAL(PW1 ,C22 , GRND3 ,0. )
COVAL(PW1 ,C23 , GRND3 ,0. )
PATCH(OB2 ,LOW , 1, 1, 1, 10, 1, 1, 1, 1)
COVAL(OB2 ,P1 , FIXFLU ,14.1963 )
COVAL(OB2 ,V1 ,0. ,0. )
COVAL(OB2 ,W1 ,0. ,11.85 )
COVAL(OB2 ,KE ,0. ,0.351056 )
COVAL(OB2 ,EP ,0. ,53.825245 )
COVAL(OB2 ,C6 ,0. ,1. )
COVAL(OB2 ,C7 ,0. ,1. )
COVAL(OB2 ,C8 ,0. ,1. )
COVAL(OB2 ,C9 ,0. ,1. )
COVAL(OB2 ,C10 ,0. ,1. )
COVAL(OB2 ,C11 ,0. ,1. )
COVAL(OB2 ,C12 ,0. ,1. )
COVAL(OB2 ,C13 ,0. ,1. )
COVAL(OB2 ,C14 ,0. ,1. )
COVAL(OB2 ,C15 ,0. ,1. )
COVAL(OB2 ,C16 ,0. ,1. )
COVAL(OB2 ,C17 ,0. ,1. )
COVAL(OB2 ,C18 ,0. ,1. )
COVAL(OB2 ,C19 ,0. ,1. )
COVAL(OB2 ,C20 ,0. ,1. )
COVAL(OB2 ,C21 ,0. ,1. )
COVAL(OB2 ,C22 ,0. ,1. )
COVAL(OB2 ,C23 ,0. ,1. )
PATCH(OB3 ,HIGH , 1, 1, 1, 10, 80, 80, 1, 1)
COVAL(OB3 ,P1 ,1000. ,0. )
COVAL(OB3 ,V1 ,0. ,0. )
COVAL(OB3 ,W1 ,0. ,0. )
COVAL(OB3 ,KE ,0. , SAME )
COVAL(OB3 ,EP ,0. , SAME )
XCYCLE = F
EGWF = T
WALLCO = GRND2
BUOYA =0. ; BUOYB =0. ; BUOYC =9.81
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 140 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-04
************************************************************
Group 16. Terminate Iterations
LITER(P1)=50 ;LITER(V1)=10
LITER(W1)=10 ;LITER(KE)=20
LITER(EP)=20 ;LITER(C6)=20
LITER(C7)=20 ;LITER(C8)=20
LITER(C9)=20 ;LITER(C10)=20
LITER(C11)=20 ;LITER(C12)=20
LITER(C13)=20 ;LITER(C14)=20
LITER(C15)=20 ;LITER(C16)=20
LITER(C17)=20 ;LITER(C18)=20
LITER(C19)=20 ;LITER(C20)=20
LITER(C21)=20 ;LITER(C22)=20
LITER(C23)=20
ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
ENDIT(EP)=1.0E-03 ;ENDIT(C6)=1.0E-03
ENDIT(C7)=1.0E-03 ;ENDIT(C8)=1.0E-03
ENDIT(C9)=1.0E-03 ;ENDIT(C10)=1.0E-03
ENDIT(C11)=1.0E-03 ;ENDIT(C12)=1.0E-03
ENDIT(C13)=1.0E-03 ;ENDIT(C14)=1.0E-03
ENDIT(C15)=1.0E-03 ;ENDIT(C16)=1.0E-03
ENDIT(C17)=1.0E-03 ;ENDIT(C18)=1.0E-03
ENDIT(C19)=1.0E-03 ;ENDIT(C20)=1.0E-03
ENDIT(C21)=1.0E-03 ;ENDIT(C22)=1.0E-03
ENDIT(C23)=1.0E-03
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(V1,FALSDT,1.)
RELAX(W1,FALSDT,1.)
RELAX(KE,LINRLX,0.5)
RELAX(EP,LINRLX,0.5)
RELAX(C6,FALSDT,10.)
RELAX(C7,FALSDT,1.0E+09)
RELAX(C8,FALSDT,1.0E+09)
RELAX(C9,FALSDT,1.0E+09)
RELAX(C10,FALSDT,1.0E+09)
RELAX(C11,FALSDT,10.)
RELAX(C12,FALSDT,10.)
RELAX(C13,FALSDT,10.)
RELAX(C14,FALSDT,10.)
RELAX(C15,FALSDT,10.)
RELAX(C16,FALSDT,10.)
RELAX(C17,FALSDT,10.)
RELAX(C18,FALSDT,10.)
RELAX(C19,FALSDT,10.)
RELAX(C20,FALSDT,10.)
RELAX(C21,FALSDT,10.)
RELAX(C22,FALSDT,10.)
RELAX(C23,FALSDT,10.)
RELAX(VP18,LINRLX,1.)
RELAX(TP18,LINRLX,1.)
RELAX(VP17,LINRLX,1.)
RELAX(TP17,LINRLX,1.)
RELAX(VP16,LINRLX,1.)
RELAX(TP16,LINRLX,1.)
RELAX(VP15,LINRLX,1.)
RELAX(TP15,LINRLX,1.)
RELAX(VP14,LINRLX,1.)
RELAX(TP14,LINRLX,1.)
RELAX(VP13,LINRLX,1.)
RELAX(TP13,LINRLX,1.)
RELAX(TR13,LINRLX,1.)
RELAX(VP12,LINRLX,1.)
RELAX(TP12,LINRLX,1.)
RELAX(VP11,LINRLX,1.)
RELAX(TP11,LINRLX,1.)
RELAX(VP10,LINRLX,1.)
RELAX(TP10,LINRLX,1.)
RELAX(VP9,LINRLX,1.)
RELAX(TP9,LINRLX,1.)
RELAX(VP8,LINRLX,1.)
RELAX(TP8,LINRLX,1.)
RELAX(VD8,LINRLX,1.)
RELAX(TR8,LINRLX,1.)
RELAX(VP7,LINRLX,1.)
RELAX(TP7,LINRLX,1.)
RELAX(VP6,LINRLX,1.)
RELAX(TP6,LINRLX,1.)
RELAX(VP5,LINRLX,1.)
RELAX(TP5,LINRLX,1.)
RELAX(VP4,LINRLX,1.)
RELAX(TP4,LINRLX,1.)
RELAX(VP3,LINRLX,1.)
RELAX(TP3,LINRLX,1.)
RELAX(VP2,LINRLX,1.)
RELAX(TP2,LINRLX,1.)
RELAX(VP1,LINRLX,1.)
RELAX(TP1,LINRLX,1.)
RELAX(DEP5,LINRLX,1.)
RELAX(DEP4,LINRLX,1.)
RELAX(DEP3,LINRLX,1.)
RELAX(DEP2,LINRLX,1.)
RELAX(VSTR,LINRLX,1.)
RELAX(DEP1,LINRLX,1.)
RELAX(ENUL,LINRLX,1.)
RELAX(YPLS,LINRLX,1.)
RELAX(STRS,LINRLX,1.)
RELAX(EPKE,LINRLX,1.)
RELAX(DEN1,LINRLX,1.)
RELAX(EL1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
KELIN = 3
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
VARMAX(C6)=1. ;VARMIN(C6)=0.
VARMAX(C7)=1.0E+10 ;VARMIN(C7)=-1.0E+10
VARMAX(C8)=1.0E+10 ;VARMIN(C8)=-1.0E+10
VARMAX(C9)=1.0E+10 ;VARMIN(C9)=-1.0E+10
VARMAX(C10)=1.0E+10 ;VARMIN(C10)=-1.0E+10
VARMAX(C11)=1.0E+10 ;VARMIN(C11)=-1.0E+10
VARMAX(C12)=1.0E+10 ;VARMIN(C12)=-1.0E+10
VARMAX(C13)=1.0E+10 ;VARMIN(C13)=-1.0E+10
VARMAX(C14)=1.0E+10 ;VARMIN(C14)=-1.0E+10
VARMAX(C15)=1.0E+10 ;VARMIN(C15)=-1.0E+10
VARMAX(C16)=1.0E+10 ;VARMIN(C16)=-1.0E+10
VARMAX(C17)=1.0E+10 ;VARMIN(C17)=-1.0E+10
VARMAX(C18)=1.0E+10 ;VARMIN(C18)=-1.0E+10
VARMAX(C19)=1.0E+10 ;VARMIN(C19)=-1.0E+10
VARMAX(C20)=1.0E+10 ;VARMIN(C20)=-1.0E+10
VARMAX(C21)=1.0E+10 ;VARMIN(C21)=-1.0E+10
VARMAX(C22)=1.0E+10 ;VARMIN(C22)=-1.0E+10
VARMAX(C23)=1.0E+10 ;VARMIN(C23)=-1.0E+10
VARMAX(VP18)=1.0E+10 ;VARMIN(VP18)=-1.0E+10
VARMAX(TP18)=1.0E+10 ;VARMIN(TP18)=-1.0E+10
VARMAX(VP17)=1.0E+10 ;VARMIN(VP17)=-1.0E+10
VARMAX(TP17)=1.0E+10 ;VARMIN(TP17)=-1.0E+10
VARMAX(VP16)=1.0E+10 ;VARMIN(VP16)=-1.0E+10
VARMAX(TP16)=1.0E+10 ;VARMIN(TP16)=-1.0E+10
VARMAX(VP15)=1.0E+10 ;VARMIN(VP15)=-1.0E+10
VARMAX(TP15)=1.0E+10 ;VARMIN(TP15)=-1.0E+10
VARMAX(VP14)=1.0E+10 ;VARMIN(VP14)=-1.0E+10
VARMAX(TP14)=1.0E+10 ;VARMIN(TP14)=-1.0E+10
VARMAX(VP13)=1.0E+10 ;VARMIN(VP13)=-1.0E+10
VARMAX(TP13)=1.0E+10 ;VARMIN(TP13)=-1.0E+10
VARMAX(TR13)=1.0E+10 ;VARMIN(TR13)=-1.0E+10
VARMAX(VP12)=1.0E+10 ;VARMIN(VP12)=-1.0E+10
VARMAX(TP12)=1.0E+10 ;VARMIN(TP12)=-1.0E+10
VARMAX(VP11)=1.0E+10 ;VARMIN(VP11)=-1.0E+10
VARMAX(TP11)=1.0E+10 ;VARMIN(TP11)=-1.0E+10
VARMAX(VP10)=1.0E+10 ;VARMIN(VP10)=-1.0E+10
VARMAX(TP10)=1.0E+10 ;VARMIN(TP10)=-1.0E+10
VARMAX(VP9)=1.0E+10 ;VARMIN(VP9)=-1.0E+10
VARMAX(TP9)=1.0E+10 ;VARMIN(TP9)=-1.0E+10
VARMAX(VP8)=1.0E+10 ;VARMIN(VP8)=-1.0E+10
VARMAX(TP8)=1.0E+10 ;VARMIN(TP8)=-1.0E+10
VARMAX(VD8)=1.0E+10 ;VARMIN(VD8)=-1.0E+10
VARMAX(TR8)=1.0E+10 ;VARMIN(TR8)=-1.0E+10
VARMAX(VP7)=1.0E+10 ;VARMIN(VP7)=-1.0E+10
VARMAX(TP7)=1.0E+10 ;VARMIN(TP7)=-1.0E+10
VARMAX(VP6)=1.0E+10 ;VARMIN(VP6)=-1.0E+10
VARMAX(TP6)=1.0E+10 ;VARMIN(TP6)=-1.0E+10
VARMAX(VP5)=1.0E+10 ;VARMIN(VP5)=-1.0E+10
VARMAX(TP5)=1.0E+10 ;VARMIN(TP5)=-1.0E+10
VARMAX(VP4)=1.0E+10 ;VARMIN(VP4)=-1.0E+10
VARMAX(TP4)=1.0E+10 ;VARMIN(TP4)=-1.0E+10
VARMAX(VP3)=1.0E+10 ;VARMIN(VP3)=-1.0E+10
VARMAX(TP3)=1.0E+10 ;VARMIN(TP3)=-1.0E+10
VARMAX(VP2)=1.0E+10 ;VARMIN(VP2)=-1.0E+10
VARMAX(TP2)=1.0E+10 ;VARMIN(TP2)=-1.0E+10
VARMAX(VP1)=1.0E+10 ;VARMIN(VP1)=-1.0E+10
VARMAX(TP1)=1.0E+10 ;VARMIN(TP1)=-1.0E+10
VARMAX(DEP5)=1.0E+10 ;VARMIN(DEP5)=-1.0E+10
VARMAX(DEP4)=1.0E+10 ;VARMIN(DEP4)=-1.0E+10
VARMAX(DEP3)=1.0E+10 ;VARMIN(DEP3)=-1.0E+10
VARMAX(DEP2)=1.0E+10 ;VARMIN(DEP2)=-1.0E+10
VARMAX(VSTR)=1.0E+10 ;VARMIN(VSTR)=-1.0E+10
VARMAX(DEP1)=1.0E+10 ;VARMIN(DEP1)=-1.0E+10
VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10
VARMAX(YPLS)=1.0E+10 ;VARMIN(YPLS)=-1.0E+10
VARMAX(STRS)=1.0E+10 ;VARMIN(STRS)=-1.0E+10
VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
VARMAX(EL1)=1.0E+10 ;VARMIN(EL1)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
NAMGRD =FLAR
GENK = T
PARSOL = F
ISG62 = 0
SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,0)
SPEDAT(SET,DFLUX,DIAP6,R,2.5E-06)
SPEDAT(SET,DFLUX,DENP6,R,920.)
SPEDAT(SET,DFLUX,DIAP7,R,2.6E-06)
SPEDAT(SET,DFLUX,DENP7,R,920.)
SPEDAT(SET,DFLUX,DIAP8,R,3.2E-06)
SPEDAT(SET,DFLUX,DENP8,R,920.)
SPEDAT(SET,DFLUX,DIAP9,R,3.6E-06)
SPEDAT(SET,DFLUX,DENP9,R,920.)
SPEDAT(SET,DFLUX,DIAP10,R,3.7E-06)
SPEDAT(SET,DFLUX,DENP10,R,920.)
SPEDAT(SET,DFLUX,DIAP11,R,5.3E-06)
SPEDAT(SET,DFLUX,DENP11,R,920.)
SPEDAT(SET,DFLUX,DIAP12,R,7.2E-06)
SPEDAT(SET,DFLUX,DENP12,R,920.)
SPEDAT(SET,DFLUX,DIAP13,R,8.100001E-06)
SPEDAT(SET,DFLUX,DENP13,R,920.)
SPEDAT(SET,DFLUX,DIAP14,R,1.0E-05)
SPEDAT(SET,DFLUX,DENP14,R,920.)
SPEDAT(SET,DFLUX,DIAP15,R,1.01E-05)
SPEDAT(SET,DFLUX,DENP15,R,920.)
SPEDAT(SET,DFLUX,DIAP16,R,1.4E-05)
SPEDAT(SET,DFLUX,DENP16,R,920.)
SPEDAT(SET,DFLUX,DIAP17,R,1.68E-05)
SPEDAT(SET,DFLUX,DENP17,R,920.)
SPEDAT(SET,DFLUX,DIAP18,R,2.1E-05)
SPEDAT(SET,DFLUX,DENP18,R,920.)
SPEDAT(SET,DFLUX,DFMODL,L,T)
SPEDAT(SET,DFLUX,DEPOMOD,I,3)
SPEDAT(SET,DFLUX,DENP1,R,920.)
SPEDAT(SET,DFLUX,DIAP1,R,5.0E-08)
SPEDAT(SET,DFLUX,DENP2,R,920.)
SPEDAT(SET,DFLUX,DIAP2,R,3.0E-07)
SPEDAT(SET,DFLUX,DENP3,R,920.)
SPEDAT(SET,DFLUX,DIAP3,R,7.5E-07)
SPEDAT(SET,DFLUX,DENP4,R,920.)
SPEDAT(SET,DFLUX,DIAP4,R,1.2E-06)
SPEDAT(SET,DFLUX,DENP5,R,920.)
SPEDAT(SET,DFLUX,DIAP5,R,1.4E-06)
SPEDAT(SET,MAKE,VS_NY,C,=0)
SPEDAT(SET,STORED,VS_NY,C,=VSTR[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,V*,C,=VS_NY)
SPEDAT(SET,MAKE,YP_NY,C,=0)
SPEDAT(SET,STORED,YP_NY,C,=YPLS[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,Y+,C,=YP_NY)
SPEDAT(SET,MAKE,TP_1,C,=0)
SPEDAT(SET,STORED,TP_1,C,=TP1[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_1,C,=TP_1)
SPEDAT(SET,MAKE,VP_1,C,=0)
SPEDAT(SET,STORED,VP_1,C,=VP1[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_1,C,=VP_1)
SPEDAT(SET,MAKE,TP_2,C,=0)
SPEDAT(SET,STORED,TP_2,C,=TP2[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_2,C,=TP_2)
SPEDAT(SET,MAKE,VP_2,C,=0)
SPEDAT(SET,STORED,VP_2,C,=VP2[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_2,C,=VP_2)
SPEDAT(SET,MAKE,TP_3,C,=0)
SPEDAT(SET,STORED,TP_3,C,=TP3[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_3,C,=TP_3)
SPEDAT(SET,MAKE,VP_3,C,=0)
SPEDAT(SET,STORED,VP_3,C,=VP3[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_3,C,=VP_3)
SPEDAT(SET,MAKE,TP_4,C,=0)
SPEDAT(SET,STORED,TP_4,C,=TP4[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_4,C,=TP_4)
SPEDAT(SET,MAKE,VP_4,C,=0)
SPEDAT(SET,STORED,VP_4,C,=VP4[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_4,C,=VP_4)
SPEDAT(SET,MAKE,TP_5,C,=0)
SPEDAT(SET,STORED,TP_5,C,=TP5[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_5,C,=TP_5)
SPEDAT(SET,MAKE,VP_5,C,=0)
SPEDAT(SET,STORED,VP_5,C,=VP5[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_5,C,=VP_5)
SPEDAT(SET,MAKE,TP_6,C,=0)
SPEDAT(SET,STORED,TP_6,C,=TP6[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_6,C,=TP_6)
SPEDAT(SET,MAKE,VP_6,C,=0)
SPEDAT(SET,STORED,VP_6,C,=VP6[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_6,C,=VP_6)
SPEDAT(SET,MAKE,TP_7,C,=0)
SPEDAT(SET,STORED,TP_7,C,=TP7[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_7,C,=TP_7)
SPEDAT(SET,MAKE,VP_7,C,=0)
SPEDAT(SET,STORED,VP_7,C,=VP7[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_7,C,=VP_7)
SPEDAT(SET,MAKE,TP_8,C,=0)
SPEDAT(SET,STORED,TP_8,C,=TP8[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_8,C,=TP_8)
SPEDAT(SET,MAKE,VP_8,C,=0)
SPEDAT(SET,STORED,VP_8,C,=VP8[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_8,C,=VP_8)
SPEDAT(SET,MAKE,TP_9,C,=0)
SPEDAT(SET,STORED,TP_9,C,=TP9[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_9,C,=TP_9)
SPEDAT(SET,MAKE,VP_9,C,=0)
SPEDAT(SET,STORED,VP_9,C,=VP9[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_9,C,=VP_9)
SPEDAT(SET,MAKE,TP_10,C,=0)
SPEDAT(SET,STORED,TP_10,C,=TP10[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_10,C,=TP_10)
SPEDAT(SET,MAKE,VP_10,C,=0)
SPEDAT(SET,STORED,VP_10,C,=VP10[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_10,C,=VP_10)
SPEDAT(SET,MAKE,TP_11,C,=0)
SPEDAT(SET,STORED,TP_11,C,=TP11[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_11,C,=TP_11)
SPEDAT(SET,MAKE,VP_11,C,=0)
SPEDAT(SET,STORED,VP_11,C,=VP11[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_11,C,=VP_11)
SPEDAT(SET,MAKE,TP_12,C,=0)
SPEDAT(SET,STORED,TP_12,C,=TP12[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_12,C,=TP_12)
SPEDAT(SET,MAKE,VP_12,C,=0)
SPEDAT(SET,STORED,VP_12,C,=VP12[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_12,C,=VP_12)
SPEDAT(SET,MAKE,TP_13,C,=0)
SPEDAT(SET,STORED,TP_13,C,=TP13[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_13,C,=TP_13)
SPEDAT(SET,MAKE,VP_13,C,=0)
SPEDAT(SET,STORED,VP_13,C,=VP13[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_13,C,=VP_13)
SPEDAT(SET,MAKE,TP_14,C,=0)
SPEDAT(SET,STORED,TP_14,C,=TP14[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_14,C,=TP_14)
SPEDAT(SET,MAKE,VP_14,C,=0)
SPEDAT(SET,STORED,VP_14,C,=VP14[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_14,C,=VP_14)
SPEDAT(SET,MAKE,TP_15,C,=0)
SPEDAT(SET,STORED,TP_15,C,=TP15[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_15,C,=TP_15)
SPEDAT(SET,MAKE,VP_15,C,=0)
SPEDAT(SET,STORED,VP_15,C,=VP15[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_15,C,=VP_15)
SPEDAT(SET,MAKE,TP_16,C,=0)
SPEDAT(SET,STORED,TP_16,C,=TP16[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_16,C,=TP_16)
SPEDAT(SET,MAKE,VP_16,C,=0)
SPEDAT(SET,STORED,VP_16,C,=VP16[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_16,C,=VP_16)
SPEDAT(SET,MAKE,TP_17,C,=0)
SPEDAT(SET,STORED,TP_17,C,=TP17[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_17,C,=TP_17)
SPEDAT(SET,MAKE,VP_17,C,=0)
SPEDAT(SET,STORED,VP_17,C,=VP17[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_17,C,=VP_17)
SPEDAT(SET,MAKE,TP_18,C,=0)
SPEDAT(SET,STORED,TP_18,C,=TP18[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,TR+_18,C,=TP_18)
SPEDAT(SET,MAKE,VP_18,C,=0)
SPEDAT(SET,STORED,VP_18,C,=VP18[1&NY&75]!ZSLFIN)
SPEDAT(SET,PRINT,VD+_18,C,=VP_18)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,OBJNAM,!OB1,C,WALLN)
SPEDAT(SET,OBJTYP,!OB1,C,PLATE)
SPEDAT(SET,OBJNAM,!PW1,C,WALLN)
SPEDAT(SET,OBJTYP,!PW1,C,USER_DEFINED)
SPEDAT(SET,OBJNAM,!OB2,C,INL)
SPEDAT(SET,OBJTYP,!OB2,C,INLET)
SPEDAT(SET,OBJNAM,!OB3,C,OUT)
SPEDAT(SET,OBJTYP,!OB3,C,OPENING)
SPEDAT(SET,ARATIO,!OB3,R,1.)
SPEDAT(SET,FACETDAT,NUMOBJ,I,3)
************************************************************
Group 20. Preliminary Printout
DISTIL = T ;NULLPR = F
NDST = 0
DSTTOL =1.0E-02
EX(P1)=102.099998 ;EX(V1)=2.702E-03
EX(W1)=13.17 ;EX(KE)=0.8443
EX(EP)=339.399994 ;EX(C6)=0.9999
EX(C7)=0.9998 ;EX(C8)=0.9998
EX(C9)=0.9996 ;EX(C10)=0.9993
EX(C11)=0.9992 ;EX(C12)=0.9916
EX(C13)=0.9996 ;EX(C14)=0.9993
EX(C15)=0.9992 ;EX(C16)=0.9916
EX(C17)=0.8103 ;EX(C18)=0.5956
EX(C19)=0.4797 ;EX(C20)=0.4788
EX(C21)=0.4527 ;EX(C22)=0.4378
EX(C23)=0.4238 ;EX(VP18)=0.03118
EX(TP18)=4.429 ;EX(VP17)=0.02722
EX(TP17)=2.84 ;EX(VP16)=0.02393
EX(TP16)=1.976 ;EX(VP15)=0.01967
EX(TP15)=1.033 ;EX(VP14)=0.01953
EX(TP14)=1.013 ;EX(VP13)=9.044E-03
EX(TP13)=0.6672 ;EX(TR13)=1.872E-04
EX(VP12)=2.487E-03 ;EX(TP12)=0.5285
EX(VP11)=8.919E-05 ;EX(TP11)=0.2887
EX(VP10)=8.692E-06 ;EX(TP10)=0.1425
EX(VP9)=7.68E-06 ;EX(TP9)=0.1351
EX(VP8)=4.832E-06 ;EX(TP8)=0.1073
EX(VD8)=3.572E-06 ;EX(TR8)=3.011E-05
EX(VP7)=2.688E-06 ;EX(TP7)=0.07166
EX(VP6)=2.474E-06 ;EX(TP6)=0.06641
EX(VP5)=1.398E-06 ;EX(TP5)=0.02185
EX(VP4)=1.368E-06 ;EX(TP4)=0.01633
EX(VP3)=1.496E-06 ;EX(TP3)=6.846E-03
EX(VP2)=2.379E-06 ;EX(TP2)=1.406E-03
EX(VP1)=7.507E-06 ;EX(TP1)=1.249E-04
EX(DEP5)=1.235E-06 ;EX(DEP4)=1.208E-06
EX(DEP3)=1.32E-06 ;EX(DEP2)=2.099E-06
EX(VSTR)=0.07367 ;EX(DEP1)=6.619E-06
EX(ENUL)=1.524E-05 ;EX(YPLS)=2.659
EX(STRS)=0.05432 ;EX(EPKE)=277.100006
EX(DEN1)=1.198 ;EX(EL1)=7.267E-04
EX(ENUT)=3.3E-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(P1,Y,N,Y,Y,Y,Y)
OUTPUT(V1,Y,N,Y,Y,Y,Y)
OUTPUT(W1,Y,N,Y,Y,Y,Y)
OUTPUT(KE,Y,N,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(C6,Y,N,Y,Y,Y,Y)
OUTPUT(C7,Y,N,Y,Y,Y,Y)
OUTPUT(C8,Y,N,Y,Y,Y,Y)
OUTPUT(C9,Y,N,Y,Y,Y,Y)
OUTPUT(C10,Y,N,Y,Y,Y,Y)
OUTPUT(C11,Y,N,Y,Y,Y,Y)
OUTPUT(C12,Y,N,Y,Y,Y,Y)
OUTPUT(C13,Y,N,Y,Y,Y,Y)
OUTPUT(C14,Y,N,Y,Y,Y,Y)
OUTPUT(C15,Y,N,Y,Y,Y,Y)
OUTPUT(C16,Y,N,Y,Y,Y,Y)
OUTPUT(C17,Y,N,Y,Y,Y,Y)
OUTPUT(C18,Y,N,Y,Y,Y,Y)
OUTPUT(C19,Y,N,Y,Y,Y,Y)
OUTPUT(C20,Y,N,Y,Y,Y,Y)
OUTPUT(C21,Y,N,Y,Y,Y,Y)
OUTPUT(C22,Y,N,Y,Y,Y,Y)
OUTPUT(C23,Y,N,Y,Y,Y,Y)
OUTPUT(VP18,Y,N,Y,N,N,N)
OUTPUT(TP18,Y,N,Y,N,N,N)
OUTPUT(VP17,Y,N,Y,N,N,N)
OUTPUT(TP17,Y,N,Y,N,N,N)
OUTPUT(VP16,Y,N,Y,N,N,N)
OUTPUT(TP16,Y,N,Y,N,N,N)
OUTPUT(VP15,Y,N,Y,N,N,N)
OUTPUT(TP15,Y,N,Y,N,N,N)
OUTPUT(VP14,Y,N,Y,N,N,N)
OUTPUT(TP14,Y,N,Y,N,N,N)
OUTPUT(VP13,Y,N,Y,N,N,N)
OUTPUT(TP13,Y,N,Y,N,N,N)
OUTPUT(TR13,Y,N,Y,N,N,N)
OUTPUT(VP12,Y,N,Y,N,N,N)
OUTPUT(TP12,Y,N,Y,N,N,N)
OUTPUT(VP11,Y,N,Y,N,N,N)
OUTPUT(TP11,Y,N,Y,N,N,N)
OUTPUT(VP10,Y,N,Y,N,N,N)
OUTPUT(TP10,Y,N,Y,N,N,N)
OUTPUT(VP9,Y,N,Y,N,N,N)
OUTPUT(TP9,Y,N,Y,N,N,N)
OUTPUT(VP8,Y,N,Y,N,N,N)
OUTPUT(TP8,Y,N,Y,N,N,N)
OUTPUT(VD8,Y,N,Y,N,N,N)
OUTPUT(TR8,Y,N,Y,N,N,N)
OUTPUT(VP7,Y,N,Y,N,N,N)
OUTPUT(TP7,Y,N,Y,N,N,N)
OUTPUT(VP6,Y,N,Y,N,N,N)
OUTPUT(TP6,Y,N,Y,N,N,N)
OUTPUT(VP5,Y,N,Y,N,N,N)
OUTPUT(TP5,Y,N,Y,N,N,N)
OUTPUT(VP4,Y,N,Y,N,N,N)
OUTPUT(TP4,Y,N,Y,N,N,N)
OUTPUT(VP3,Y,N,Y,N,N,N)
OUTPUT(TP3,Y,N,Y,N,N,N)
OUTPUT(VP2,Y,N,Y,N,N,N)
OUTPUT(TP2,Y,N,Y,N,N,N)
OUTPUT(VP1,Y,N,Y,N,N,N)
OUTPUT(TP1,Y,N,Y,N,N,N)
OUTPUT(DEP5,Y,N,Y,N,N,N)
OUTPUT(DEP4,Y,N,Y,N,N,N)
OUTPUT(DEP3,Y,N,Y,N,N,N)
OUTPUT(DEP2,Y,N,Y,N,N,N)
OUTPUT(VSTR,Y,N,Y,N,N,N)
OUTPUT(DEP1,Y,N,Y,N,N,N)
OUTPUT(ENUL,Y,N,Y,N,N,N)
OUTPUT(YPLS,Y,N,Y,N,N,N)
OUTPUT(STRS,Y,N,Y,N,N,N)
OUTPUT(EPKE,Y,N,Y,N,N,N)
OUTPUT(DEN1,Y,N,Y,N,N,N)
OUTPUT(EL1,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 10 ;IZMON = 75
NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
UWATCH = T ;USTEER = T
HIGHLO = F
************************************************************
Group 23.Field Print-Out & Plot Control
NPRINT = 100000 ;NUMCLS = 5
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