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.170632 ;YFRAC(4)=0.216532
 YFRAC(5)=0.257843 ;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.742158 ;YFRAC(26)=0.783468
 YFRAC(27)=0.829369 ;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.411074&0.&60.&3.588926&8.7$)
 SPEDAT(SET,INFOB,*!INFOB4,C,5&0.&0.&0.)-BOX(24.&-30.&0.&2.5&90.&2$)
 SPEDAT(SET,INFOB,*!INFOB4,C,.5&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,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 = 15 ;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