TALK=T;RUN(1,1)
  PHOTON USE
  p




  con vpor x 1 fil;.001
  set vec comp - v1 gvel
  vec x 1
  con vpor x 1;0.5 0.501;int 1
  msg   Porosity distribution, gas velocity vectors and
  msg   cavity contour
  pause
  vec cl;red
  con vpor x 1;0.5 0.501;int 1
  set vec comp - v2 cvel
  vec x 1
  msg   Fines velocity vectors
  pause
  cl;
  con rcar x 1 fil;.001
  con vpor x 1;0.5 0.501;int 1
  msg Burning rate of coal fines
  pause
  cl;
  con rbed x 1 fil;.001
  con vpor x 1;0.5 0.501;int 1
  msg Burning rate of coke lumps
  pause
  cl;
  con yo2 x 1 fil;.001
  con vpor x 1;0.5 0.501;int 1
  msg Oxygen mass fraction
  pause
  cl;
  con yn2 x 1 fil;.001
  con vpor x 1;0.5 0.501;int 1
  msg Nitrogen mass fraction
  pause
  cl;
  con yh2 x 1 fil;.001
  con vpor x 1;0.5 0.501;int 1
  msg Hydrogen mass fraction
  pause
  cl;
  con yh2o x 1 fil;.001
  con vpor x 1;0.5 0.501;int 1
  msg Water vapour mass fraction
  pause
  cl;
  con yco2 x 1 fil;.001
  con vpor x 1;0.5 0.501;int 1
  msg Carbon dioxide mass fraction
  pause
  cl;
  con yco x 1 fil;.001
  con vpor x 1;0.5 0.501;int 1
  msg Carbon monoxide mass fraction
  pause
  cl;
  con tmp1 x 1 fil;.001
  con vpor x 1;0.5 0.501;int 1
  msg  Gas temperature distribution
  pause
  cl;
  con tmp2 x 1 fil;.001
  con vpor x 1;0.5 0.501;int 1
  msg  Fines temperature distribution
  pause
  cl;
  con fuel x 1 fil;.001
  con vpor x 1;0.5 0.501;int 1
  msg Coal fines volume fraction
  pause
  ENDUSE
  DISPLAY

  Combustion-driven raceway: a pulverised air/coal flame
  in a coke bed.

  Two-phase, 2D, one space, shadow method for volume
  porosity of the bed.

  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Coal fines flame in a coke bed: raceway )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 1
 XULAST =1.
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 20
 YVLAST =15.
 YFRAC(1)=0.05 ;YFRAC(2)=0.1
 YFRAC(3)=0.15 ;YFRAC(4)=0.2
 YFRAC(5)=0.25 ;YFRAC(6)=0.3
 YFRAC(7)=0.35 ;YFRAC(8)=0.4
 YFRAC(9)=0.45 ;YFRAC(10)=0.5
 YFRAC(11)=0.55 ;YFRAC(12)=0.6
 YFRAC(13)=0.65 ;YFRAC(14)=0.7
 YFRAC(15)=0.75 ;YFRAC(16)=0.8
 YFRAC(17)=0.85 ;YFRAC(18)=0.9
 YFRAC(19)=0.95 ;YFRAC(20)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 20
 ZWLAST =10.
 ZFRAC(1)=0.05 ;ZFRAC(2)=0.1
 ZFRAC(3)=0.15 ;ZFRAC(4)=0.2
 ZFRAC(5)=0.25 ;ZFRAC(6)=0.3
 ZFRAC(7)=0.35 ;ZFRAC(8)=0.4
 ZFRAC(9)=0.45 ;ZFRAC(10)=0.5
 ZFRAC(11)=0.55 ;ZFRAC(12)=0.6
 ZFRAC(13)=0.65 ;ZFRAC(14)=0.7
 ZFRAC(15)=0.75 ;ZFRAC(16)=0.8
 ZFRAC(17)=0.85 ;ZFRAC(18)=0.9
 ZFRAC(19)=0.95 ;ZFRAC(20)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = F
 EQUVEL = F
 NAME(1)=P1 ;NAME(5)=V1
 NAME(6)=V2 ;NAME(7)=GVEL
 NAME(8)=CVEL ;NAME(9)=GAS
 NAME(10)=FUEL ;NAME(11)=SHAD
 NAME(14)=H1 ;NAME(15)=H2
 NAME(16)=MIXF ;NAME(18)=POR
 NAME(121)=TOTB ;NAME(122)=RCAR
 NAME(123)=GOPA ;NAME(124)=GOFU
 NAME(125)=GH ;NAME(126)=GC
 NAME(127)=GO ;NAME(128)=FRAC
 NAME(129)=FLIM ;NAME(130)=YSUM
 NAME(131)=HSUB ;NAME(132)=RMIX
 NAME(133)=VVPO ;NAME(134)=SPOR
 NAME(135)=RBED ;NAME(136)=TMP2
 NAME(137)=TMP1 ;NAME(138)=NPOR
 NAME(139)=HPOR ;NAME(140)=VPOR
 NAME(141)=RHO2 ;NAME(142)=RHO1
 NAME(143)=YH2O ;NAME(144)=YH2
 NAME(145)=YN2 ;NAME(146)=YCO2
 NAME(147)=YO2 ;NAME(148)=YCO
 NAME(149)=CFIP ;NAME(150)=MDOT
    * 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(V2,Y,Y,N,N,N,Y)
 SOLUTN(GVEL,Y,Y,N,N,N,Y)
 SOLUTN(CVEL,Y,Y,N,N,N,Y)
 SOLUTN(GAS,Y,Y,N,Y,N,Y)
 SOLUTN(FUEL,Y,Y,N,Y,N,Y)
 SOLUTN(SHAD,Y,Y,N,Y,N,Y)
 SOLUTN(H1,Y,Y,N,N,N,Y)
 SOLUTN(H2,Y,Y,N,N,N,Y)
 SOLUTN(MIXF,Y,Y,Y,N,N,Y)
 SOLUTN(POR,Y,Y,N,N,N,Y)
 SOLUTN(TOTB,Y,N,N,N,N,Y)
 SOLUTN(RCAR,Y,N,N,N,N,Y)
 SOLUTN(GOPA,Y,N,N,N,N,Y)
 SOLUTN(GOFU,Y,N,N,N,N,Y)
 SOLUTN(GH,Y,N,N,N,N,Y)
 SOLUTN(GC,Y,N,N,N,N,Y)
 SOLUTN(GO,Y,N,N,N,N,Y)
 SOLUTN(FRAC,Y,N,N,N,N,Y)
 SOLUTN(FLIM,Y,N,N,N,N,Y)
 SOLUTN(YSUM,Y,N,N,N,N,Y)
 SOLUTN(HSUB,Y,N,N,N,N,Y)
 SOLUTN(RMIX,Y,N,N,N,N,Y)
 SOLUTN(VVPO,Y,N,N,N,N,Y)
 SOLUTN(SPOR,Y,N,N,N,N,Y)
 SOLUTN(RBED,Y,N,N,N,N,Y)
 SOLUTN(TMP2,Y,N,N,N,N,Y)
 SOLUTN(TMP1,Y,N,N,N,N,Y)
 SOLUTN(NPOR,Y,N,N,N,N,Y)
 SOLUTN(HPOR,Y,N,N,N,N,Y)
 SOLUTN(VPOR,Y,N,N,N,N,Y)
 SOLUTN(RHO2,Y,N,N,N,N,Y)
 SOLUTN(RHO1,Y,N,N,N,N,Y)
 SOLUTN(YH2O,Y,N,N,N,N,Y)
 SOLUTN(YH2,Y,N,N,N,N,Y)
 SOLUTN(YN2,Y,N,N,N,N,Y)
 SOLUTN(YCO2,Y,N,N,N,N,Y)
 SOLUTN(YO2,Y,N,N,N,N,Y)
 SOLUTN(YCO,Y,N,N,N,N,Y)
 SOLUTN(CFIP,Y,N,N,N,N,Y)
 SOLUTN(MDOT,Y,N,N,N,N,Y)
 DEN1 = 142
 DEN2 = 141
 INTMDT = 150
 EPOR = 0 ;HPOR = 139 ;NPOR = 138 ;VPOR = 140
 TEMP1 = 137
 TEMP2 = 136
 ************************************************************
  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(V2,Y,Y,Y,Y,N,Y)
 TERMS(GVEL,Y,Y,Y,Y,Y,Y)
 TERMS(CVEL,Y,Y,Y,Y,N,Y)
 TERMS(GAS,Y,Y,Y,Y,Y,Y)
 TERMS(FUEL,Y,Y,Y,Y,N,Y)
 TERMS(SHAD,Y,Y,Y,Y,N,Y)
 TERMS(H1,N,Y,Y,Y,Y,Y)
 TERMS(H2,N,Y,Y,Y,N,Y)
 TERMS(MIXF,N,Y,Y,Y,Y,Y)
 TERMS(POR,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 EQDVDP = 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 = GRND ;TMP1 =0. ;EL1 =0.
 TSURR =0. ;TEMP0 =0. ;PRESS0 =1.0E+05
 DVO1DT =0. ;DRH1DP =0.
 RHO1A =0. ;RHO1B =0. ;RHO1C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 RHO2 =1000. ;TMP2 =0. ;EL2 =0.
 DVO2DT =0. ;DRH2DP =0.
 ENUL =1.0E-05 ;ENUT =0.
 PHINT(V1)=-2.022E+04 ;PHINT(V2)=-2.022E+04
 PHINT(GVEL)=-2.022E+04 ;PHINT(CVEL)=-2.022E+04
 PHINT(GAS)=-2.022E+04 ;PHINT(FUEL)=-2.022E+04
 PHINT(SHAD)=-2.022E+04 ;PHINT(H1)=GRND
 PHINT(H2)=GRND ;PHINT(MIXF)=-2.022E+04
 PHINT(POR)=-2.022E+04
 PHNH1A =0. ;PHNH1B =0. ;PHNH1C =0.
 PHNH2A =0. ;PHNH2B =0. ;PHNH2C =0.
 PRNDTL(V1)=1. ;PRNDTL(V2)=1.
 PRNDTL(GVEL)=1. ;PRNDTL(CVEL)=1.
 PRNDTL(GAS)=1. ;PRNDTL(FUEL)=1.
 PRNDTL(SHAD)=1. ;PRNDTL(H1)=1.
 PRNDTL(H2)=1. ;PRNDTL(MIXF)=1.
 PRNDTL(POR)=1.
 PRT(V1)=1. ;PRT(V2)=1.
 PRT(GVEL)=1. ;PRT(CVEL)=1.
 PRT(GAS)=1. ;PRT(FUEL)=1.
 PRT(SHAD)=1. ;PRT(H1)=1.
 PRT(H2)=1. ;PRT(MIXF)=1.
 PRT(POR)=1.
 CP1 =1100. ;CP2 =1100.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 CFIPS = GRND1
 RLOLIM =0. ;CMDOT = GRND
 CFIPA =0. ;CFIPB =0.
 CFIPC =1.0E+09 ;CFIPD =0.
 CMDTA =0. ;CMDTB =0.
 CMDTC =0. ;CMDTD =0.
 CINT(P1)=1. ;CINT(V1)=1.
 CINT(V2)=1. ;CINT(GVEL)=1.
 CINT(CVEL)=1. ;CINT(GAS)=1.
 CINT(FUEL)=1. ;CINT(SHAD)=1.
 CINT(H1)=GRND ;CINT(H2)=GRND
 CINT(MIXF)=1. ;CINT(POR)=1.
 PHINT(P1)=-2.022E+04 ;PHINT(V1)=-2.022E+04
 PHINT(V2)=-2.022E+04 ;PHINT(GVEL)=-2.022E+04
 PHINT(CVEL)=-2.022E+04 ;PHINT(GAS)=-2.022E+04
 PHINT(FUEL)=-2.022E+04 ;PHINT(SHAD)=-2.022E+04
 PHINT(H1)=GRND ;PHINT(H2)=GRND
 PHINT(MIXF)=-2.022E+04 ;PHINT(POR)=-2.022E+04
 CINH1A =0. ;CINH1B =0. ;CINH1C =0.
 CINH2A =0. ;CINH2B =0. ;CINH2C =0.
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(V1)=1.0E-10
 FIINIT(V2)=1.0E-10 ;FIINIT(GVEL)=85.590065
 FIINIT(CVEL)=10. ;FIINIT(GAS)=0.999
 FIINIT(FUEL)=1.0E-03 ;FIINIT(SHAD)=1.0E-03
 FIINIT(H1)=7.7E+05 ;FIINIT(H2)=3.79674E+07
 FIINIT(MIXF)=0. ;FIINIT(POR)=1.0E-10
 FIINIT(TOTB)=1.0E-10 ;FIINIT(RCAR)=1.0E-10
 FIINIT(GOPA)=1.0E-10 ;FIINIT(GOFU)=1.0E-10
 FIINIT(GH)=1.0E-10 ;FIINIT(GC)=1.0E-10
 FIINIT(GO)=1.0E-10 ;FIINIT(FRAC)=1.0E-10
 FIINIT(FLIM)=1.0E-10 ;FIINIT(YSUM)=1.0E-10
 FIINIT(HSUB)=1.0E-10 ;FIINIT(RMIX)=1.0E-10
 FIINIT(VVPO)=1.0E-10 ;FIINIT(SPOR)=0.5
 FIINIT(RBED)=1.0E-10 ;FIINIT(TMP2)=700.
 FIINIT(TMP1)=700. ;FIINIT(NPOR)=0.5
 FIINIT(HPOR)=0.5 ;FIINIT(VPOR)=0.5
 FIINIT(RHO2)=1.0E-10 ;FIINIT(RHO1)=0.49705
 FIINIT(YH2O)=1.0E-10 ;FIINIT(YH2)=1.0E-10
 FIINIT(YN2)=1.0E-10 ;FIINIT(YCO2)=1.0E-10
 FIINIT(YO2)=1.0E-10 ;FIINIT(YCO)=1.0E-10
 FIINIT(CFIP)=1.0E-10 ;FIINIT(MDOT)=0.1
 
 PATCH(WALS ,INIVAL, 1, 1, 4, 4, 1, 4, 1, 1)
 INIT(WALS ,NPOR,0. ,0. )
 
 PATCH(WALN ,INIVAL, 1, 1, 5, 5, 1, 4, 1, 1)
 INIT(WALN ,NPOR,0. ,0. )
 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(INLETS ,LOW , 1, 1, 5, 5, 1, 1, 1, 1)
 COVAL(INLETS ,CVEL,0. ,10. )
 COVAL(INLETS ,FUEL, FIXFLU ,10. )
 COVAL(INLETS ,H2 ,0. ,3.79674E+07 )
 
 PATCH(INLETG ,LOW , 1, 1, 5, 5, 1, 1, 1, 1)
 COVAL(INLETG ,GVEL,0. ,85.590065 )
 COVAL(INLETG ,GAS , FIXFLU ,42.5 )
 COVAL(INLETG ,H1 ,0. ,7.7E+05 )
 COVAL(INLETG ,MIXF,0. ,0. )
 COVAL(INLETG ,POR ,0. ,0. )
 
 PATCH(OUTLET ,CELL , 1, 1, 20, 20, 1, 20, 1, 1)
 COVAL(OUTLET ,GAS ,100. ,0. )
 COVAL(OUTLET ,FUEL,1.0E+06 ,0. )
 
 PATCH(BED2GAS ,VOLUME, 1, 1, 1, 20, 1, 20, 1, 1)
 COVAL(BED2GAS ,GAS , FIXFLU , GRND )
 COVAL(BED2GAS ,H1 ,0. , GRND )
 COVAL(BED2GAS ,MIXF, GRND ,1. )
 COVAL(BED2GAS ,POR , GRND , GRND )
 
 PATCH(FRIC ,VOLUME, 1, 1, 1, 20, 1, 20, 1, 1)
 COVAL(FRIC ,V1 , GRND ,0. )
 COVAL(FRIC ,V2 , GRND ,0. )
 COVAL(FRIC ,GVEL, GRND ,0. )
 COVAL(FRIC ,CVEL, GRND ,0. )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 500 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 ISWR1 = 1 ;ISWR2 = 10000
 SELREF = T
 RESFAC =1.0E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(V1)=10
 LITER(V2)=10 ;LITER(GVEL)=10
 LITER(CVEL)=10 ;LITER(GAS)=1
 LITER(FUEL)=1 ;LITER(SHAD)=1
 LITER(H1)=20 ;LITER(H2)=20
 LITER(MIXF)=20 ;LITER(POR)=20
 ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
 ENDIT(V2)=1.0E-03 ;ENDIT(GVEL)=1.0E-03
 ENDIT(CVEL)=1.0E-03 ;ENDIT(GAS)=1.0E-03
 ENDIT(FUEL)=1.0E-03 ;ENDIT(SHAD)=1.0E-03
 ENDIT(H1)=1.0E-03 ;ENDIT(H2)=1.0E-03
 ENDIT(MIXF)=1.0E-03 ;ENDIT(POR)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.3)
 RELAX(V1,FALSDT,1.0E-04)
 RELAX(V2,FALSDT,1.0E-04)
 RELAX(GVEL,FALSDT,1.0E-04)
 RELAX(CVEL,FALSDT,1.0E-04)
 RELAX(GAS,LINRLX,0.1)
 RELAX(FUEL,LINRLX,0.1)
 RELAX(SHAD,LINRLX,0.3)
 RELAX(H1,FALSDT,1.0E-04)
 RELAX(H2,FALSDT,1.0E-04)
 RELAX(MIXF,FALSDT,1.0E-04)
 RELAX(POR,FALSDT,0.1)
 RELAX(TOTB,LINRLX,1.)
 RELAX(RCAR,LINRLX,1.)
 RELAX(GOPA,LINRLX,1.)
 RELAX(GOFU,LINRLX,1.)
 RELAX(GH,LINRLX,1.)
 RELAX(GC,LINRLX,1.)
 RELAX(GO,LINRLX,1.)
 RELAX(FRAC,LINRLX,1.)
 RELAX(FLIM,LINRLX,1.)
 RELAX(YSUM,LINRLX,1.)
 RELAX(HSUB,LINRLX,1.)
 RELAX(RMIX,LINRLX,1.)
 RELAX(VVPO,LINRLX,1.)
 RELAX(SPOR,LINRLX,1.)
 RELAX(RBED,LINRLX,1.)
 RELAX(TMP2,LINRLX,1.)
 RELAX(TMP1,LINRLX,1.)
 RELAX(RHO2,LINRLX,1.)
 RELAX(RHO1,LINRLX,1.)
 RELAX(YH2O,LINRLX,1.)
 RELAX(YH2,LINRLX,1.)
 RELAX(YN2,LINRLX,1.)
 RELAX(YCO2,LINRLX,1.)
 RELAX(YO2,LINRLX,1.)
 RELAX(YCO,LINRLX,1.)
 RELAX(CFIP,LINRLX,1.)
 RELAX(MDOT,LINRLX,0.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(V2)=1.0E+06 ;VARMIN(V2)=-1.0E+06
 VARMAX(GVEL)=1.0E+06 ;VARMIN(GVEL)=-1.0E+06
 VARMAX(CVEL)=1.0E+06 ;VARMIN(CVEL)=-1.0E+06
 VARMAX(GAS)=1. ;VARMIN(GAS)=1.0E-06
 VARMAX(FUEL)=1. ;VARMIN(FUEL)=1.0E-09
 VARMAX(SHAD)=1. ;VARMIN(SHAD)=1.0E-06
 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10
 VARMAX(H2)=1.0E+10 ;VARMIN(H2)=-1.0E+10
 VARMAX(MIXF)=0.154804 ;VARMIN(MIXF)=0.
 VARMAX(POR)=1. ;VARMIN(POR)=0.
 VARMAX(TOTB)=1.0E+10 ;VARMIN(TOTB)=-1.0E+10
 VARMAX(RCAR)=1.0E+10 ;VARMIN(RCAR)=-1.0E+10
 VARMAX(GOPA)=1.0E+10 ;VARMIN(GOPA)=-1.0E+10
 VARMAX(GOFU)=1.0E+10 ;VARMIN(GOFU)=-1.0E+10
 VARMAX(GH)=1.0E+10 ;VARMIN(GH)=-1.0E+10
 VARMAX(GC)=1.0E+10 ;VARMIN(GC)=-1.0E+10
 VARMAX(GO)=1.0E+10 ;VARMIN(GO)=-1.0E+10
 VARMAX(FRAC)=1.0E+10 ;VARMIN(FRAC)=-1.0E+10
 VARMAX(FLIM)=1.0E+10 ;VARMIN(FLIM)=-1.0E+10
 VARMAX(YSUM)=1.0E+10 ;VARMIN(YSUM)=-1.0E+10
 VARMAX(HSUB)=1.0E+10 ;VARMIN(HSUB)=-1.0E+10
 VARMAX(RMIX)=1.0E+10 ;VARMIN(RMIX)=-1.0E+10
 VARMAX(VVPO)=1.0E+10 ;VARMIN(VVPO)=-1.0E+10
 VARMAX(SPOR)=1.0E+10 ;VARMIN(SPOR)=-1.0E+10
 VARMAX(RBED)=1.0E+10 ;VARMIN(RBED)=-1.0E+10
 VARMAX(TMP2)=1.0E+10 ;VARMIN(TMP2)=-1.0E+10
 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10
 VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10
 VARMAX(HPOR)=1.0E+10 ;VARMIN(HPOR)=-1.0E+10
 VARMAX(VPOR)=1.0E+10 ;VARMIN(VPOR)=-1.0E+10
 VARMAX(RHO2)=1.0E+10 ;VARMIN(RHO2)=-1.0E+10
 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10
 VARMAX(YH2O)=1. ;VARMIN(YH2O)=0.
 VARMAX(YH2)=1. ;VARMIN(YH2)=0.
 VARMAX(YN2)=1. ;VARMIN(YN2)=0.
 VARMAX(YCO2)=1. ;VARMIN(YCO2)=0.
 VARMAX(YO2)=1. ;VARMIN(YO2)=0.
 VARMAX(YCO)=1. ;VARMIN(YCO)=0.
 VARMAX(CFIP)=1.0E+10 ;VARMIN(CFIP)=-1.0E+10
 VARMAX(MDOT)=1.0E+10 ;VARMIN(MDOT)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMSAT =MOSG
 PARSOL = F
 ISG62 = 1
 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(P1,Y,Y,Y,Y,Y,Y)
 OUTPUT(V1,Y,Y,Y,Y,Y,Y)
 OUTPUT(V2,Y,Y,Y,Y,Y,Y)
 OUTPUT(GVEL,Y,Y,Y,Y,Y,Y)
 OUTPUT(CVEL,Y,Y,Y,Y,Y,Y)
 OUTPUT(GAS,Y,Y,Y,Y,Y,Y)
 OUTPUT(FUEL,Y,Y,Y,Y,Y,Y)
 OUTPUT(SHAD,Y,Y,Y,Y,Y,Y)
 OUTPUT(H1,Y,Y,Y,Y,Y,Y)
 OUTPUT(H2,Y,Y,Y,Y,Y,Y)
 OUTPUT(MIXF,Y,Y,Y,Y,Y,Y)
 OUTPUT(POR,Y,N,Y,Y,Y,Y)
 OUTPUT(TOTB,Y,N,Y,N,N,N)
 OUTPUT(RCAR,Y,N,Y,N,N,N)
 OUTPUT(GOPA,Y,N,Y,N,N,N)
 OUTPUT(GOFU,Y,N,Y,N,N,N)
 OUTPUT(GH,Y,N,Y,N,N,N)
 OUTPUT(GC,Y,N,Y,N,N,N)
 OUTPUT(GO,Y,N,Y,N,N,N)
 OUTPUT(FRAC,Y,N,Y,N,N,N)
 OUTPUT(FLIM,Y,N,Y,N,N,N)
 OUTPUT(YSUM,Y,N,Y,N,N,N)
 OUTPUT(HSUB,Y,N,Y,N,N,N)
 OUTPUT(RMIX,Y,N,Y,N,N,N)
 OUTPUT(VVPO,Y,N,Y,N,N,N)
 OUTPUT(SPOR,Y,N,Y,N,N,N)
 OUTPUT(RBED,Y,N,Y,N,N,N)
 OUTPUT(TMP2,Y,N,Y,N,N,N)
 OUTPUT(TMP1,Y,N,Y,N,N,N)
 OUTPUT(NPOR,Y,N,Y,N,N,N)
 OUTPUT(HPOR,Y,N,Y,N,N,N)
 OUTPUT(VPOR,Y,N,Y,N,N,N)
 OUTPUT(RHO2,Y,N,Y,N,N,N)
 OUTPUT(RHO1,Y,N,Y,N,N,N)
 OUTPUT(YH2O,Y,N,Y,N,N,N)
 OUTPUT(YH2,Y,N,Y,N,N,N)
 OUTPUT(YN2,Y,N,Y,N,N,N)
 OUTPUT(YCO2,Y,N,Y,N,N,N)
 OUTPUT(YO2,Y,N,Y,N,N,N)
 OUTPUT(YCO,Y,N,Y,N,N,N)
 OUTPUT(CFIP,Y,N,Y,N,N,N)
 OUTPUT(MDOT,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 18 ;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
 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 = 1 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.2
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
 
 PATCH(PROFIL1 ,PROFIL, 1, 1, 1, 20, 1, 20, 1, 1)
 PLOT(PROFIL1 ,FUEL,0. ,0. )
 PLOT(PROFIL1 ,MIXF,0. ,0. )
 PLOT(PROFIL1 ,YCO2,0. ,0. )
 PLOT(PROFIL1 ,YCO ,0. ,0. )
 
 PATCH(PROFIL2 ,PROFIL, 1, 1, 1, 1, 20, 20, 1, 1)
 PLOT(PROFIL2 ,P1 ,0. ,0. )
 PLOT(PROFIL2 ,V1 ,0. ,0. )
 PLOT(PROFIL2 ,GVEL,0. ,0. )
 PLOT(PROFIL2 ,RHO1,0. ,0. )
 
 PATCH(PROFIL3 ,PROFIL, 1, 1, 1, 20, 1, 20, 1, 1)
 PLOT(PROFIL3 ,TMP2,0. ,0. )
 PLOT(PROFIL3 ,TMP1,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP