TALK=T;RUN(1,1)
ENDMAIN
save26end
save26end
save26end
save26end
save26end
save26end
save26end
EOF:
VRV USE
* Start of frame
* Setting display switches
TEXT CLEAR
AXIS ON
CELPOS OFF
CONTOUR SCALE ON
POSITION CONTOURKEY 2.252252E-02 8.389262E-02
TEXT ON
POSITION TITLE 2.590090E-01 8.993289E-01
POSITION PROBE 7.972973E-01 8.389262E-02
GRID OFF
* Domain scaling factors
SCALE 1.000000E+00 1.000000E+00 1.000000E+00
* Settings for current slice
PROBE 5.468750E+00 1.450000E+01 5.000000E-02; PROBE ON
SLICE Y
SLICE OUTLINE ON
* View and up directions
VIEW -6.841546E-02 -9.966097E-01 -4.569499E-02
UP -4.894052E-02 -4.239444E-02 9.979016E-01
* View centre
VIEW CENTRE 2.900758E+01 1.717579E+01 3.993689E+00
* View size
VIEW SIZE 1.704600E+01
* View perspective
VIEW DEPTH 3.000000E+00;VIEW TILT 0.8
DOMAIN ON
* Setting object visibility and painting status
OBJECT SHOW TYPE BLOCKAGE
OBJECT PAINT TYPE BLOCKAGE OFF
OBJECT WIREFRAME TYPE BLOCKAGE OFF
OBJECT SHOW TYPE INLET
OBJECT PAINT TYPE INLET OFF
OBJECT WIREFRAME TYPE INLET OFF
OBJECT SHOW TYPE OUTLET
OBJECT PAINT TYPE OUTLET OFF
OBJECT WIREFRAME TYPE OUTLET OFF
OBJECT SHOW TYPE USER_DEFINED
OBJECT PAINT TYPE USER_DEFINED OFF
OBJECT WIREFRAME TYPE USER_DEFINED OFF
VARIABLE u1
VECTOR ON
VECTOR COLOUR MULTI
VECTOR TYPE TOTAL
CONTOUR ON
CONTOUR BLANK ON
CONTOUR AVERAGE ON
SURFACE OFF
MINMAX OFF
CONTOUR OPAQUENESS 100
PAUSE
VRVEND
************************************************************
Group 1. Run Title and Number
************************************************************
************************************************************
TEXT(gas-release with wind angle = 0.785398 )
************************************************************
************************************************************
IRUNN = 1 ;LIBREF = 162
************************************************************
Group 2. Time dependence
STEADY = T
************************************************************
Group 3. X-Direction Grid Spacing
CARTES = T
NX = 50
XULAST =60.
XFRAC(1)=0.025083 ;XFRAC(2)=0.046894
XFRAC(3)=0.065861 ;XFRAC(4)=0.082353
XFRAC(5)=0.096694 ;XFRAC(6)=0.100661
XFRAC(7)=0.104628 ;XFRAC(8)=0.108595
XFRAC(9)=0.112562 ;XFRAC(10)=0.116529
XFRAC(11)=0.120496 ;XFRAC(12)=0.124463
XFRAC(13)=0.12843 ;XFRAC(14)=0.132397
XFRAC(15)=0.136364 ;XFRAC(16)=0.156282
XFRAC(17)=0.179188 ;XFRAC(18)=0.20553
XFRAC(19)=0.235824 ;XFRAC(20)=0.270661
XFRAC(21)=0.305499 ;XFRAC(22)=0.335792
XFRAC(23)=0.362134 ;XFRAC(24)=0.38504
XFRAC(25)=0.404959 ;XFRAC(26)=0.414876
XFRAC(27)=0.424793 ;XFRAC(28)=0.434711
XFRAC(29)=0.444628 ;XFRAC(30)=0.454545
XFRAC(31)=0.471041 ;XFRAC(32)=0.488362
XFRAC(33)=0.506549 ;XFRAC(34)=0.525645
XFRAC(35)=0.545696 ;XFRAC(36)=0.56675
XFRAC(37)=0.588856 ;XFRAC(38)=0.612067
XFRAC(39)=0.636439 ;XFRAC(40)=0.66203
XFRAC(41)=0.6889 ;XFRAC(42)=0.717114
XFRAC(43)=0.746738 ;XFRAC(44)=0.777844
XFRAC(45)=0.810504 ;XFRAC(46)=0.844798
XFRAC(47)=0.880807 ;XFRAC(48)=0.918616
XFRAC(49)=0.958315 ;XFRAC(50)=1.
************************************************************
Group 4. Y-Direction Grid Spacing
NY = 30
YVLAST =30.
YFRAC(1)=0.062964 ;YFRAC(2)=0.119631
YFRAC(3)=0.170631 ;YFRAC(4)=0.216532
YFRAC(5)=0.257842 ;YFRAC(6)=0.295022
YFRAC(7)=0.328483 ;YFRAC(8)=0.358599
YFRAC(9)=0.385702 ;YFRAC(10)=0.410096
YFRAC(11)=0.43205 ;YFRAC(12)=0.451809
YFRAC(13)=0.469591 ;YFRAC(14)=0.485596
YFRAC(15)=0.5 ;YFRAC(16)=0.514404
YFRAC(17)=0.530409 ;YFRAC(18)=0.548191
YFRAC(19)=0.56795 ;YFRAC(20)=0.589904
YFRAC(21)=0.614298 ;YFRAC(22)=0.641401
YFRAC(23)=0.671517 ;YFRAC(24)=0.704978
YFRAC(25)=0.742157 ;YFRAC(26)=0.783468
YFRAC(27)=0.829368 ;YFRAC(28)=0.880369
YFRAC(29)=0.937036 ;YFRAC(30)=1.
************************************************************
Group 5. Z-Direction Grid Spacing
PARAB = F
NZ = 35
ZWLAST =8.75
ZFRAC(1)=0.02 ;ZFRAC(2)=0.04
ZFRAC(3)=0.06 ;ZFRAC(4)=0.08
ZFRAC(5)=0.1 ;ZFRAC(6)=0.12
ZFRAC(7)=0.14 ;ZFRAC(8)=0.16
ZFRAC(9)=0.18 ;ZFRAC(10)=0.2
ZFRAC(11)=0.22 ;ZFRAC(12)=0.24
ZFRAC(13)=0.26 ;ZFRAC(14)=0.28
ZFRAC(15)=0.3 ;ZFRAC(16)=0.32
ZFRAC(17)=0.34 ;ZFRAC(18)=0.36
ZFRAC(19)=0.38 ;ZFRAC(20)=0.4
ZFRAC(21)=0.42 ;ZFRAC(22)=0.44
ZFRAC(23)=0.46 ;ZFRAC(24)=0.48
ZFRAC(25)=0.5 ;ZFRAC(26)=0.531373
ZFRAC(27)=0.565883 ;ZFRAC(28)=0.603844
ZFRAC(29)=0.645601 ;ZFRAC(30)=0.691533
ZFRAC(31)=0.74206 ;ZFRAC(32)=0.797638
ZFRAC(33)=0.858775 ;ZFRAC(34)=0.926025
ZFRAC(35)=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(12)=KE ;NAME(13)=EP
NAME(143)=DEN1 ;NAME(144)=ENUT
NAME(145)=EPOR ;NAME(146)=MARK
NAME(147)=G.25 ;NAME(148)=GAS
NAME(149)=VABS ;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(P1,Y,Y,Y,N,N,Y)
SOLUTN(U1,Y,Y,Y,N,N,Y)
SOLUTN(V1,Y,Y,Y,N,N,Y)
SOLUTN(W1,Y,Y,Y,N,N,Y)
SOLUTN(KE,Y,Y,Y,N,N,Y)
SOLUTN(EP,Y,Y,Y,N,N,Y)
SOLUTN(DEN1,Y,N,N,N,N,Y)
SOLUTN(ENUT,Y,N,N,N,N,Y)
SOLUTN(EPOR,Y,N,N,N,N,Y)
SOLUTN(MARK,Y,N,N,N,N,Y)
SOLUTN(G.25,Y,N,N,N,N,Y)
SOLUTN(GAS,Y,Y,Y,N,N,Y)
SOLUTN(VABS,Y,N,N,N,N,Y)
SOLUTN(PRPS,Y,N,N,N,N,Y)
DEN1 = 143
VIST = 144
EPOR = 145 ;HPOR = 0 ;NPOR = 0 ;VPOR = 0
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
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(KE,N,Y,Y,Y,Y,N)
TERMS(EP,N,Y,Y,Y,Y,N)
TERMS(GAS,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. ;TMP1 =0. ;EL1 = GRND4
TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
DVO1DT =0. ;DRH1DP =0.
EMISS =0. ;SCATT =0.
RADIA =0. ;RADIB =0.
EL1A =0. ;EL1B =0. ;EL1C =0.
ENUL =1.0E-05 ;ENUT = GRND3
ENUTA =0. ;ENUTB =0. ;ENUTC =0.
IENUTA = 0
PRNDTL(U1)=1. ;PRNDTL(V1)=1.
PRNDTL(W1)=1. ;PRNDTL(KE)=1.
PRNDTL(EP)=1. ;PRNDTL(GAS)=1.
PRT(U1)=1. ;PRT(V1)=1.
PRT(W1)=1. ;PRT(KE)=1.
PRT(EP)=1.314 ;PRT(GAS)=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(KE)=1.0E-10 ;FIINIT(EP)=1.0E-10
FIINIT(DEN1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
FIINIT(EPOR)=1. ;FIINIT(MARK)=1.0E-10
FIINIT(G.25)=1.0E-10 ;FIINIT(GAS)=1.0E-10
FIINIT(VABS)=1.0E-10 ;FIINIT(PRPS)=-1.
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(KESOURCE,PHASEM, 1, 50, 1, 30, 1, 35, 1, 1)
COVAL(KESOURCE,KE , GRND4 , GRND4 )
COVAL(KESOURCE,EP , GRND4 , GRND4 )
PATCH(WHOLE ,CELL , 1, 50, 1, 30, 1, 35, 1, 1)
PATCH(GRAVITY ,PHASEM, 1, 50, 1, 30, 1, 35, 1, 1)
COVAL(GRAVITY ,W1 ,In-Form:source - see Grp 19)
PATCH(CEILING ,HIGH , 1, 50, 1, 30, 35, 35, 1, 1)
COVAL(CEILING ,P1 ,1. ,0. )
COVAL(CEILING ,U1 ,In-Form:source - see Grp 19)
COVAL(CEILING ,V1 ,In-Form:source - see Grp 19)
COVAL(CEILING ,GAS ,In-Form:source - see Grp 19)
PATCH(CEILING2,HIGH , 1, 50, 1, 30, 34, 34, 1, 1)
COVAL(CEILING2,W1 ,In-Form:source - see Grp 19)
COVAL(CEILING2,GAS ,In-Form:source - see Grp 19)
PATCH(FLOOR ,LOW , 1, 50, 1, 30, 1, 1, 1, 1)
COVAL(FLOOR ,U1 ,In-Form:source - see Grp 19)
COVAL(FLOOR ,V1 ,In-Form:source - see Grp 19)
PATCH(INFOB1 ,WEST , 1, 1, 1, 30, 1, 35, 1, 1)
COVAL(INFOB1 ,P1 ,In-Form:source - see Grp 19)
COVAL(INFOB1 ,U1 ,In-Form:source - see Grp 19)
COVAL(INFOB1 ,V1 ,In-Form:source - see Grp 19)
COVAL(INFOB1 ,W1 ,In-Form:source - see Grp 19)
COVAL(INFOB1 ,GAS ,In-Form:source - see Grp 19)
PATCH(INFOB2 ,EAST , 48, 50, 1, 30, 1, 35, 1, 1)
COVAL(INFOB2 ,P1 ,In-Form:source - see Grp 19)
COVAL(INFOB2 ,U1 ,In-Form:source - see Grp 19)
COVAL(INFOB2 ,V1 ,In-Form:source - see Grp 19)
COVAL(INFOB2 ,W1 ,In-Form:source - see Grp 19)
COVAL(INFOB2 ,GAS ,In-Form:source - see Grp 19)
PATCH(INFOB3 ,SOUTH , 1, 50, 1, 1, 1, 35, 1, 1)
COVAL(INFOB3 ,P1 ,In-Form:source - see Grp 19)
COVAL(INFOB3 ,U1 ,In-Form:source - see Grp 19)
COVAL(INFOB3 ,V1 ,In-Form:source - see Grp 19)
COVAL(INFOB3 ,W1 ,In-Form:source - see Grp 19)
COVAL(INFOB3 ,GAS ,In-Form:source - see Grp 19)
PATCH(INFOB4 ,NORTH , 1, 50, 29, 30, 1, 35, 1, 1)
COVAL(INFOB4 ,P1 ,In-Form:source - see Grp 19)
COVAL(INFOB4 ,U1 ,In-Form:source - see Grp 19)
COVAL(INFOB4 ,V1 ,In-Form:source - see Grp 19)
COVAL(INFOB4 ,W1 ,In-Form:source - see Grp 19)
COVAL(INFOB4 ,GAS ,In-Form:source - see Grp 19)
PATCH(INFOB5 ,CELL , 24, 31, 1, 30, 1, 17, 1, 1)
COVAL(INFOB5 ,U1 ,In-Form:source - see Grp 19)
COVAL(INFOB5 ,V1 ,In-Form:source - see Grp 19)
COVAL(INFOB5 ,W1 ,In-Form:source - see Grp 19)
PATCH(INFOB6 ,CELL , 1, 25, 7, 23, 1, 17, 1, 1)
COVAL(INFOB6 ,P1 ,In-Form:source - see Grp 19)
COVAL(INFOB6 ,U1 ,In-Form:source - see Grp 19)
COVAL(INFOB6 ,W1 ,In-Form:source - see Grp 19)
COVAL(INFOB6 ,GAS ,In-Form:source - see Grp 19)
PATCH(INFOB7 ,PHASEM, 1, 25, 1, 13, 1, 17, 1, 1)
COVAL(INFOB7 ,U1 ,In-Form:source - see Grp 19)
COVAL(INFOB7 ,V1 ,In-Form:source - see Grp 19)
COVAL(INFOB7 ,W1 ,In-Form:source - see Grp 19)
PATCH(INFOB8 ,PHASEM, 1, 25, 16, 30, 1, 17, 1, 1)
COVAL(INFOB8 ,U1 ,In-Form:source - see Grp 19)
COVAL(INFOB8 ,V1 ,In-Form:source - see Grp 19)
COVAL(INFOB8 ,W1 ,In-Form:source - see Grp 19)
XCYCLE = F
EGWF = T
WALLCO = GRND2
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 50 ;ISWC1 = 1
LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
SELREF = T
RESFAC =1.0E-03
************************************************************
Group 16. Terminate Iterations
LITER(P1)=200 ;LITER(U1)=10
LITER(V1)=10 ;LITER(W1)=10
LITER(KE)=20 ;LITER(EP)=20
LITER(GAS)=200
ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03
ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
ENDIT(GAS)=0.
************************************************************
Group 17. Relaxation
RELAX(P1,LINRLX,1.)
RELAX(U1,LINRLX,1.)
RELAX(V1,LINRLX,1.)
RELAX(W1,LINRLX,1.)
RELAX(KE,LINRLX,0.5)
RELAX(EP,LINRLX,0.5)
RELAX(DEN1,LINRLX,1.)
RELAX(ENUT,LINRLX,1.)
RELAX(MARK,LINRLX,1.)
RELAX(G.25,LINRLX,1.)
RELAX(GAS,LINRLX,1.)
RELAX(VABS,LINRLX,1.)
RELAX(PRPS,LINRLX,1.)
KELIN = 0
OVRRLX =0.
EXPERT = F ;NNORSL = F
************************************************************
Group 18. Limits
VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
VARMAX(U1)=1. ;VARMIN(U1)=-1.0E+11
VARMAX(V1)=1. ;VARMIN(V1)=-1.0E+11
VARMAX(W1)=1. ;VARMIN(W1)=-1.0E+11
VARMAX(KE)=1.0E+10 ;VARMIN(KE)=0.02
VARMAX(EP)=1.0E+10 ;VARMIN(EP)=0.028284
VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
VARMAX(EPOR)=1.0E+10 ;VARMIN(EPOR)=-1.0E+10
VARMAX(MARK)=1.0E+10 ;VARMIN(MARK)=-1.0E+10
VARMAX(G.25)=1.0E+10 ;VARMIN(G.25)=-1.0E+10
VARMAX(GAS)=1.0E+10 ;VARMIN(GAS)=0.
VARMAX(VABS)=1.0E+10 ;VARMIN(VABS)=-1.0E+10
VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
************************************************************
Group 19. Data transmitted to GROUND
GENK = T
PARSOL = F
CONWIZ = T
CALFOR = T
ISG50 = 1
ISG52 = 1
ISG62 = 1
SPEDAT(SET,STORED,G.25,C,=GAS^0.25!ZSLFIN)
SPEDAT(SET,PROPERTY,RHO1,C,=1.189*(1.+GAS*(34.080002/28.9-1.)))
SPEDAT(SET,PROPERTY,ENUL,C,=1.544E-5)
SPEDAT(SET,INITIAL,U1!WHOLE,C,=1.414214*(ZG/8.75)^(1./7.))
SPEDAT(SET,INITIAL,V1!WHOLE,C,=-1.414213*(ZG/8.75)^(1./7.))
SPEDAT(SET,INITIAL,GAS!WHOLE,C,=0.0)
SPEDAT(SET,SOURCE,W1!GRAVITY,C,=COVAL(2.0E-10&-9.81*(DEN1-1.189)))
SPEDAT(SET,SOURCE,U1!CEILING,C,=COVAL(2.0E+10&1.414214*(ZG/8.75)^$)
SPEDAT(SET,SOURCE,U1!CEILING,C,(1./7.)))
SPEDAT(SET,SOURCE,V1!CEILING,C,=COVAL(2.0E+10&-1.414213*(ZG/8.75)$)
SPEDAT(SET,SOURCE,V1!CEILING,C,^(1./7.)))
SPEDAT(SET,SOURCE,GAS!CEILING,C,=COVAL(2.0E+10&0.0))
SPEDAT(SET,SOURCE,W1!CEILING2,C,=COVAL(1.E3&0.0))
SPEDAT(SET,SOURCE,GAS!CEILING2,C,=COVAL(2.0E+10&0.0))
SPEDAT(SET,SOURCE,U1!FLOOR,C,=COVAL(DEN1*VABS*5.0E-03&0.0))
SPEDAT(SET,SOURCE,V1!FLOOR,C,=COVAL(DEN1*VABS*5.0E-03&0.0))
SPEDAT(SET,INFOB,*!INFOB1,C,=BOX(0.&0.&0.&1.504977&30.&8.75&0.&0.$)
SPEDAT(SET,INFOB,*!INFOB1,C,&0.)!INFOB_1)
SPEDAT(SET,SOURCE,R1!INFOB1,C,=-1.E10*P1!INFOB_1!LINE)
SPEDAT(SET,SOURCE,U1!INFOB1,C,=1.414214*(ZG/8.75)^(1./7.)!INFOB_1$)
SPEDAT(SET,SOURCE,U1!INFOB1,C,!FIXV)
SPEDAT(SET,SOURCE,V1!INFOB1,C,=-1.414213*(ZG/8.75)^(1./7.)!INFOB_$)
SPEDAT(SET,SOURCE,V1!INFOB1,C,1!FIXV)
SPEDAT(SET,SOURCE,W1!INFOB1,C,=0.!INFOB_1!FIXV)
SPEDAT(SET,SOURCE,GAS!INFOB1,C,=0.0!INFOB_1!ONLY)
SPEDAT(SET,INITIAL,MARK!INFOB1,C,=1+1!INFOB_1)
SPEDAT(SET,INFOB,*!INFOB2,C,=BOX(55.116959&0.&0.&4.883041&30.&8.7$)
SPEDAT(SET,INFOB,*!INFOB2,C,5&0.&0.&0.)!INFOB_2)
SPEDAT(SET,SOURCE,R1!INFOB2,C,=-1.E10*P1!INFOB_2!LINE)
SPEDAT(SET,SOURCE,U1!INFOB2,C,=1.414214*(ZG/8.75)^(1./7.)!INFOB_2$)
SPEDAT(SET,SOURCE,U1!INFOB2,C,!FIXV)
SPEDAT(SET,SOURCE,V1!INFOB2,C,=-1.414213*(ZG/8.75)^(1./7.)!INFOB_$)
SPEDAT(SET,SOURCE,V1!INFOB2,C,2!FIXV)
SPEDAT(SET,SOURCE,W1!INFOB2,C,=0.!INFOB_2!FIXV)
SPEDAT(SET,SOURCE,GAS!INFOB2,C,=0.0!INFOB_2!ONLY)
SPEDAT(SET,INITIAL,MARK!INFOB2,C,=2+1!INFOB_2)
SPEDAT(SET,INFOB,*!INFOB3,C,=BOX(0.&0.&0.&60.&1.88891&8.75&0.&0.&$)
SPEDAT(SET,INFOB,*!INFOB3,C,0.)-BOX(24.&-30.&0.&2.5&90.&2.5&0.&0.$)
SPEDAT(SET,INFOB,*!INFOB3,C,&0.)!INFOB_3)
SPEDAT(SET,SOURCE,R1!INFOB3,C,=-1.E10*P1!INFOB_3!LINE)
SPEDAT(SET,SOURCE,U1!INFOB3,C,=1.414214*(ZG/8.75)^(1./7.)!INFOB_3$)
SPEDAT(SET,SOURCE,U1!INFOB3,C,!FIXV)
SPEDAT(SET,SOURCE,V1!INFOB3,C,=-1.414213*(ZG/8.75)^(1./7.)!INFOB_$)
SPEDAT(SET,SOURCE,V1!INFOB3,C,3!FIXV)
SPEDAT(SET,SOURCE,W1!INFOB3,C,=0.!INFOB_3!FIXV)
SPEDAT(SET,SOURCE,GAS!INFOB3,C,=0.0!INFOB_3!ONLY)
SPEDAT(SET,INITIAL,MARK!INFOB3,C,=3+1!INFOB_3)
SPEDAT(SET,INFOB,*!INFOB4,C,=BOX(0.&26.41107&0.&60.&3.58893&8.75&$)
SPEDAT(SET,INFOB,*!INFOB4,C,0.&0.&0.)-BOX(24.&-30.&0.&2.5&90.&2.5$)
SPEDAT(SET,INFOB,*!INFOB4,C,&0.&0.&0.)!INFOB_4)
SPEDAT(SET,SOURCE,R1!INFOB4,C,=-1.E10*P1!INFOB_4!LINE)
SPEDAT(SET,SOURCE,U1!INFOB4,C,=1.414214*(ZG/8.75)^(1./7.)!INFOB_4$)
SPEDAT(SET,SOURCE,U1!INFOB4,C,!FIXV)
SPEDAT(SET,SOURCE,V1!INFOB4,C,=-1.414213*(ZG/8.75)^(1./7.)!INFOB_$)
SPEDAT(SET,SOURCE,V1!INFOB4,C,4!FIXV)
SPEDAT(SET,SOURCE,W1!INFOB4,C,=0.!INFOB_4!FIXV)
SPEDAT(SET,SOURCE,GAS!INFOB4,C,=0.0!INFOB_4!ONLY)
SPEDAT(SET,INITIAL,MARK!INFOB4,C,=4+1!INFOB_4)
SPEDAT(SET,INFOB,*!INFOB5,C,=BOX(25.&-30.&0.&0.5&90.&2.5&0.&0.&0.$)
SPEDAT(SET,INFOB,*!INFOB5,C,)!INFOB_5)
SPEDAT(SET,SOURCE,U1!INFOB5,C,=1.E9*(0.0-U1)!INFOB_5!LINE)
SPEDAT(SET,SOURCE,V1!INFOB5,C,=1.E9*(0.0-V1)!INFOB_5!LINE)
SPEDAT(SET,SOURCE,W1!INFOB5,C,=1.E9*(0.0-W1)!INFOB_5!LINE)
SPEDAT(SET,INITIAL,MARK!INFOB5,C,=1.0!INFOB_5)
SPEDAT(SET,INFOB,*!INFOB6,C,=POINT(6.25&15.&1.2&0.01)!INFOB_6)
SPEDAT(SET,SOURCE,R1!INFOB6,C,=(2*1.0E+05/1.189*34.080002/28.9)^.$)
SPEDAT(SET,SOURCE,R1!INFOB6,C,5*1.189*34.080002/28.9*0.0508^2*3.1$)
SPEDAT(SET,SOURCE,R1!INFOB6,C,4159/4!INFOB_6)
SPEDAT(SET,SOURCE,U1!INFOB6,C,=(2*1.0E+05/1.189*34.080002/28.9)^.$)
SPEDAT(SET,SOURCE,U1!INFOB6,C,5*1.267591E-06!INFOB_6!ONLY)
SPEDAT(SET,SOURCE,W1!INFOB6,C,=(2*1.0E+05/1.189*34.080002/28.9)^.$)
SPEDAT(SET,SOURCE,W1!INFOB6,C,5*1!INFOB_6!ONLY)
SPEDAT(SET,SOURCE,GAS!INFOB6,C,=1.0!INFOB_6!ONLY)
SPEDAT(SET,INFOB,*!INFOB7,C,=ELLPSD(6.25&15.&1.+0.2&0.2&100.0&0.2$)
SPEDAT(SET,INFOB,*!INFOB7,C,&0.&0.&0.)!INFOB_7)
SPEDAT(SET,SOURCE,U1!INFOB7,C,=1.E9*(0.0-U1)!INFOB_7!LINE)
SPEDAT(SET,SOURCE,V1!INFOB7,C,=1.E9*(0.0-V1)!INFOB_7!LINE)
SPEDAT(SET,SOURCE,W1!INFOB7,C,=1.E9*(0.0-W1)!INFOB_7!LINE)
SPEDAT(SET,INITIAL,MARK!INFOB7,C,=1.0!INFOB_7)
SPEDAT(SET,INFOB,*!INFOB8,C,=ELLPSD(6.25&15.&1.+0.2&0.2&100.0&0.2$)
SPEDAT(SET,INFOB,*!INFOB8,C,&0.&0.&0.)!INFOB_8)
SPEDAT(SET,SOURCE,U1!INFOB8,C,=1.E9*(0.0-U1)!INFOB_8!LINE)
SPEDAT(SET,SOURCE,V1!INFOB8,C,=1.E9*(0.0-V1)!INFOB_8!LINE)
SPEDAT(SET,SOURCE,W1!INFOB8,C,=1.E9*(0.0-W1)!INFOB_8!LINE)
SPEDAT(SET,INITIAL,MARK!INFOB8,C,=1.0!INFOB_8)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,DRAG,MOMCENX,R,0.)
SPEDAT(SET,DRAG,MOMCENY,R,0.)
SPEDAT(SET,DRAG,MOMCENZ,R,0.)
SPEDAT(SET,FACETDAT,NUMOBJ,I,0)
************************************************************
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(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(KE,Y,N,Y,Y,Y,Y)
OUTPUT(EP,Y,N,Y,Y,Y,Y)
OUTPUT(DEN1,Y,N,Y,N,N,N)
OUTPUT(ENUT,Y,N,Y,N,N,N)
OUTPUT(EPOR,Y,N,Y,N,N,N)
OUTPUT(MARK,Y,N,Y,N,N,N)
OUTPUT(G.25,Y,N,Y,Y,Y,Y)
OUTPUT(GAS,Y,N,Y,Y,Y,Y)
OUTPUT(VABS,Y,N,Y,N,N,N)
OUTPUT(PRPS,Y,N,Y,N,N,N)
************************************************************
Group 22. Monitor Print-Out
IXMON = 5 ;IYMON = 16 ;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 = -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