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

  set con fi dep 3
  con yo2 x 1 fi;0.05
  gr ou x 1
  gr ou x 1 y 25 32 z 1 5
  te
  1

  oxygen mass-fraction contours
   0.13265E+04 0.19215E+04 CR
  msg Press  to continue
  pause
  cl
  con fue x 1 fi;0.05
  gr ou x 1
  gr ou x 1 y 25 32 z 1 5
  te
  1

  contours of volume fraction of solid
   0.13265E+04 0.21215E+04 CR
  msg Press  to continue
  pause
  cl
  con yco2 x 1 fi;0.05
  gr ou x 1
  gr ou x 1 y 25 32 z 1 5
  te
  1

  carbon dioxide contours
   0.13265E+04 0.21215E+04 CR
  msg Press  to continue
  pause
  cl
  con tmp1 x 1 fi;0.05
  gr ou x 1
  gr ou x 1 y 25 32 z 1 5
  te
  1

  gas temperature contours
   0.13265E+04 0.23215E+04 CR
  msg Press  to continue
  pause
  cl
  con tmp2 x 1 fi;0.05
  gr ou x 1
  gr ou x 1 y 25 32 z 1 5
  te
  1

  coal temperature contours
   0.13265E+04 0.23215E+04 CR
  msg Press  to continue
  pause
  cl
  con mixf x 1 fi;0.05
  gr ou x 1
  gr ou x 1 y 25 32 z 1 5
  te
  1

  mixture fraction contours
   0.13265E+04 0.23215E+04 CR
  msg Press  to continue
  pause
  enduse
  DISPLAY
   The case considered is the 2d turbulent flow of combusting
   pulverised coal particles in the near-field region of a
   power-station burner. The burner comprises 3 air streams,
   tertiary air enters through a central air stream, swirling
   primary air carrying coal particles enters through an inner
   annulus located above the tertiary-air stream, and swirling
   secondary air enters through an outer annulus located above the
   primary air stream. The main combustion process is simulated
   using an IPSA-based 'equilibrium' two-phase combustion model. A
   2nd post-processing run can be performed which computes thermal
   NOx based on CHO equilibrium as determined by the CREK module.
   ===================================================
   Coal-combustion model; solid phase - carbon only
   reactions: C (s) + 0.5 O2  > CO       (exothermic )
            CO    + 0.5 O2  > CO2      (exothermic )
            C(s)  + CO2     > 2CO      (endothermic)
            C(s)  + H2O     > CO  + H2 (endothermic)
            H2    + 0.5 O2  > H2O      (exothermic )
   ===================================================
  ENDDIS
  noxdbegin
    ns   7
    nr   3
  noxdend
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(2-D Coal Combustion With NOX            )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 1
 XULAST =1.0E-02
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 32
 YVLAST =1.
 YFRAC(1)=0.03993 ;YFRAC(2)=0.078285
 YFRAC(3)=0.114752 ;YFRAC(4)=0.148828
 YFRAC(5)=0.179533 ;YFRAC(6)=0.2032
 YFRAC(7)=0.219075 ;YFRAC(8)=0.23495
 YFRAC(9)=0.250825 ;YFRAC(10)=0.2667
 YFRAC(11)=0.282575 ;YFRAC(12)=0.29845
 YFRAC(13)=0.314325 ;YFRAC(14)=0.3302
 YFRAC(15)=0.351836 ;YFRAC(16)=0.379905
 YFRAC(17)=0.411056 ;YFRAC(18)=0.444393
 YFRAC(19)=0.479456 ;YFRAC(20)=0.515959
 YFRAC(21)=0.553704 ;YFRAC(22)=0.592547
 YFRAC(23)=0.632375 ;YFRAC(24)=0.6731
 YFRAC(25)=0.700059 ;YFRAC(26)=0.735036
 YFRAC(27)=0.773852 ;YFRAC(28)=0.815391
 YFRAC(29)=0.859082 ;YFRAC(30)=0.904567
 YFRAC(31)=0.9516 ;YFRAC(32)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 25
 ZWLAST =5.28735
 ZFRAC(1)=0.010869 ;ZFRAC(2)=0.021739
 ZFRAC(3)=0.032608 ;ZFRAC(4)=0.043477
 ZFRAC(5)=0.054347 ;ZFRAC(6)=0.060154
 ZFRAC(7)=0.073215 ;ZFRAC(8)=0.091938
 ZFRAC(9)=0.11565 ;ZFRAC(10)=0.143931
 ZFRAC(11)=0.176481 ;ZFRAC(12)=0.213073
 ZFRAC(13)=0.253521 ;ZFRAC(14)=0.297675
 ZFRAC(15)=0.345406 ;ZFRAC(16)=0.396601
 ZFRAC(17)=0.451163 ;ZFRAC(18)=0.509004
 ZFRAC(19)=0.570048 ;ZFRAC(20)=0.634224
 ZFRAC(21)=0.701467 ;ZFRAC(22)=0.771718
 ZFRAC(23)=0.844924 ;ZFRAC(24)=0.921033
 ZFRAC(25)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = F
 EQUVEL = F
 NAME(1)=P1 ;NAME(3)=U1
 NAME(4)=U2 ;NAME(5)=V1
 NAME(6)=V2 ;NAME(7)=W1
 NAME(8)=W2 ;NAME(9)=GAS
 NAME(10)=FUE ;NAME(11)=SHAD
 NAME(12)=KE ;NAME(13)=EP
 NAME(14)=H1 ;NAME(15)=H2
 NAME(16)=MIXF ;NAME(17)=C2
 NAME(133)=PRPS ;NAME(134)=NPOR
 NAME(136)=EPKE ;NAME(137)=RADY
 NAME(138)=RADZ ;NAME(139)=ENUT
 NAME(140)=RHO1 ;NAME(141)=TMP2
 NAME(142)=TMP1 ;NAME(143)=YH2O
 NAME(144)=YH2 ;NAME(145)=YN2
 NAME(146)=YCO2 ;NAME(147)=YCO
 NAME(148)=YO2 ;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(U1,Y,Y,N,N,N,Y)
 SOLUTN(U2,Y,Y,N,N,N,Y)
 SOLUTN(V1,Y,Y,N,N,N,Y)
 SOLUTN(V2,Y,Y,N,N,N,Y)
 SOLUTN(W1,Y,Y,N,N,N,Y)
 SOLUTN(W2,Y,Y,N,N,N,Y)
 SOLUTN(GAS,Y,Y,N,N,N,Y)
 SOLUTN(FUE,Y,Y,N,N,N,Y)
 SOLUTN(SHAD,Y,Y,N,N,N,Y)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,Y,N,N,N,N)
 SOLUTN(H1,Y,Y,Y,N,N,Y)
 SOLUTN(H2,Y,Y,Y,N,N,Y)
 SOLUTN(MIXF,Y,Y,N,N,N,Y)
 SOLUTN(C2,Y,Y,N,N,N,Y)
 SOLUTN(PRPS,Y,N,N,N,N,N)
 SOLUTN(NPOR,Y,N,N,N,N,N)
 SOLUTN(EPKE,Y,N,N,N,N,N)
 SOLUTN(RADY,Y,Y,N,N,N,N)
 SOLUTN(RADZ,Y,Y,Y,N,N,N)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 SOLUTN(RHO1,Y,N,N,N,N,Y)
 SOLUTN(TMP2,Y,N,N,N,N,Y)
 SOLUTN(TMP1,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(YCO,Y,N,N,N,N,Y)
 SOLUTN(YO2,Y,N,N,N,N,Y)
 SOLUTN(CFIP,Y,N,N,N,N,Y)
 SOLUTN(MDOT,Y,N,N,N,N,Y)
 DEN1 = 140
 VIST = 139
 INTMDT = 150
 EPOR = 0 ;HPOR = 0 ;NPOR = 134 ;VPOR = 0
 TEMP1 = 142
 TEMP2 = 141
 PRPS = 133
 ************************************************************
  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(U2,Y,Y,Y,Y,N,Y)
 TERMS(V1,Y,Y,Y,Y,Y,Y)
 TERMS(V2,Y,Y,Y,Y,N,Y)
 TERMS(W1,Y,Y,Y,Y,Y,Y)
 TERMS(W2,Y,Y,Y,Y,N,Y)
 TERMS(GAS,Y,Y,Y,Y,Y,Y)
 TERMS(FUE,Y,Y,Y,Y,N,Y)
 TERMS(SHAD,Y,Y,Y,Y,N,Y)
 TERMS(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 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(C2,N,Y,Y,Y,N,Y)
 TERMS(RADY,N,N,Y,N,Y,N)
 TERMS(RADZ,N,N,Y,N,Y,N)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 NEWENT = 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 = GRND7 ;TMP1 =0.
 EL1 = GRND4
 TSURR =0. ;TEMP0 =0.
 PRESS0 =1.0E+05
 DVO1DT =0. ;DRH1DP =0.
 RHO1A =0.95 ;RHO1B =0.05
 RHO1C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0.5 ;RADIB =0.02
 EL1A =0. ;EL1B =0.
 EL1C =0.
 RHO2 =1350. ;TMP2 =0.
 EL2 =0.
 DVO2DT =0. ;DRH2DP =0.
 ENUL =1.0E-05 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0.
 ENUTC =0.
 IENUTA = 0
 PHINT(U1)=-2.022E+04 ;PHINT(U2)=-2.022E+04
 PHINT(V1)=-2.022E+04 ;PHINT(V2)=-2.022E+04
 PHINT(W1)=-2.022E+04 ;PHINT(W2)=-2.022E+04
 PHINT(GAS)=-2.022E+04 ;PHINT(FUE)=-2.022E+04
 PHINT(SHAD)=-2.022E+04 ;PHINT(H1)=GRND7
 PHINT(H2)=GRND7 ;PHINT(MIXF)=1.
 PHINT(C2)=-2.022E+04
 PHNH1A =0. ;PHNH1B =0.
 PHNH1C =0.
 PHNH2A =0. ;PHNH2B =0.
 PHNH2C =0.
 PRNDTL(U1)=1. ;PRNDTL(U2)=1.
 PRNDTL(V1)=1. ;PRNDTL(V2)=1.
 PRNDTL(W1)=1. ;PRNDTL(W2)=1.
 PRNDTL(GAS)=1. ;PRNDTL(FUE)=1.
 PRNDTL(SHAD)=1. ;PRNDTL(KE)=1.
 PRNDTL(EP)=1. ;PRNDTL(H1)=1.
 PRNDTL(H2)=1. ;PRNDTL(MIXF)=1.
 PRNDTL(C2)=1. ;PRNDTL(RADY)=-GRND1
 PRNDTL(RADZ)=-GRND1
 PRT(U1)=1. ;PRT(U2)=1.
 PRT(V1)=1. ;PRT(V2)=1.
 PRT(W1)=1. ;PRT(W2)=1.
 PRT(GAS)=1. ;PRT(FUE)=1.
 PRT(SHAD)=1. ;PRT(KE)=1.
 PRT(EP)=1.314 ;PRT(H1)=1.
 PRT(H2)=1. ;PRT(MIXF)=1.
 PRT(C2)=1. ;PRT(RADY)=1.0E+10
 PRT(RADZ)=1.0E+10
 CP1 =1100. ;CP2 =1100.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 CFIPS = GRND1
 RLOLIM =1.0E-06 ;CMDOT = GRND3
 CFIPA =1.0E-03 ;CFIPB =0.
 CFIPC =1.0E+05 ;CFIPD =0.
 CMDTA =1.0E+04 ;CMDTB =0.
 CMDTC =0.154804 ;CMDTD =0.
 CINT(P1)=1. ;CINT(U1)=1.
 CINT(U2)=1. ;CINT(V1)=1.
 CINT(V2)=1. ;CINT(W1)=1.
 CINT(W2)=1. ;CINT(GAS)=1.
 CINT(FUE)=1. ;CINT(SHAD)=1.
 CINT(H1)=1. ;CINT(H2)=1.
 CINT(MIXF)=0. ;CINT(C2)=0.
 PHINT(P1)=-2.022E+04 ;PHINT(U1)=-2.022E+04
 PHINT(U2)=-2.022E+04 ;PHINT(V1)=-2.022E+04
 PHINT(V2)=-2.022E+04 ;PHINT(W1)=-2.022E+04
 PHINT(W2)=-2.022E+04 ;PHINT(GAS)=-2.022E+04
 PHINT(FUE)=-2.022E+04 ;PHINT(SHAD)=-2.022E+04
 PHINT(H1)=GRND7 ;PHINT(H2)=GRND7
 PHINT(MIXF)=1. ;PHINT(C2)=-2.022E+04
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=10.
 FIINIT(U2)=10. ;FIINIT(V1)=0.
 FIINIT(V2)=0. ;FIINIT(W1)=28.07303
 FIINIT(W2)=21.354288 ;FIINIT(GAS)=0.99999
 FIINIT(FUE)=1.0E-05 ;FIINIT(SHAD)=1.0E-05
 FIINIT(KE)=1.970238 ;FIINIT(EP)=13.250978
 FIINIT(H1)=5.863E+05 ;FIINIT(H2)=3.756981E+07
 FIINIT(MIXF)=0. ;FIINIT(C2)=1.
 FIINIT(PRPS)=-1. ;FIINIT(NPOR)=1.
 FIINIT(EPKE)=1.0E-10 ;FIINIT(RADY)=1.0E-10
 FIINIT(RADZ)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
 FIINIT(RHO1)=0.652786 ;FIINIT(TMP2)=338.555542
 FIINIT(TMP1)=1.0E-10 ;FIINIT(YH2O)=1.0E-10
 FIINIT(YH2)=1.0E-10 ;FIINIT(YN2)=1.0E-10
 FIINIT(YCO2)=1.0E-10 ;FIINIT(YCO)=1.0E-10
 FIINIT(YO2)=1.0E-10 ;FIINIT(CFIP)=1.0E-10
 FIINIT(MDOT)=4.309906E-05
 
 PATCH(INIT1 ,INIVAL, 1, 1, 1, 6, 1, 25, 1, 1)
 INIT(INIT1 ,U1 ,0. ,0. )
 INIT(INIT1 ,U2 ,0. ,0. )
 INIT(INIT1 ,W1 ,0. ,41.934875 )
 INIT(INIT1 ,KE ,0. ,4.396334 )
 INIT(INIT1 ,EP ,0. ,74.533218 )
 INIT(INIT1 ,H1 ,0. ,5.863E+05 )
 INIT(INIT1 ,RHO1,0. ,0.652786 )
 
 PATCH(INIT2 ,INIVAL, 1, 1, 7, 14, 1, 25, 1, 1)
 INIT(INIT2 ,U1 ,0. ,10. )
 INIT(INIT2 ,U2 ,0. ,10. )
 INIT(INIT2 ,W1 ,0. ,21.355698 )
 INIT(INIT2 ,KE ,0. ,1.140164 )
 INIT(INIT2 ,EP ,0. ,15.750166 )
 INIT(INIT2 ,H1 ,0. ,3.724111E+05 )
 INIT(INIT2 ,RHO1,0. ,1.027704 )
 
 PATCH(BLOCK1 ,INIVAL, 1, 1, 25, 32, 1, 5, 1, 1)
 INIT(BLOCK1 ,PRPS,0. ,199. )
 
 PATCH(WALL1 ,INIVAL, 1, 1, 24, 24, 1, 5, 1, 1)
 INIT(WALL1 ,V1 ,0. ,0. )
 INIT(WALL1 ,V2 ,0. ,0. )
 INIT(WALL1 ,NPOR,0. ,0. )
 INIADD = F
 FSWEEP = 1
 IURINI = -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(RADISO ,FREEVL, 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(RADISO ,H1 , GRND1 , GRND1 )
 COVAL(RADISO ,RADY, GRND1 , GRND1 )
 COVAL(RADISO ,RADZ, GRND1 , GRND1 )
 
 PATCH(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(WALL1-NW,NWALL , 1, 1, 24, 24, 1, 5, 1, 1)
 COVAL(WALL1-NW,U1 , GRND2 ,0. )
 COVAL(WALL1-NW,W1 , GRND2 ,0. )
 COVAL(WALL1-NW,KE , GRND2 , GRND2 )
 COVAL(WALL1-NW,EP , GRND2 , GRND2 )
 
 PATCH(WALL1-SW,SWALL , 1, 1, 25, 25, 1, 5, 1, 1)
 COVAL(WALL1-SW,U1 , GRND2 ,0. )
 COVAL(WALL1-SW,W1 , GRND2 ,0. )
 COVAL(WALL1-SW,KE , GRND2 , GRND2 )
 COVAL(WALL1-SW,EP , GRND2 , GRND2 )
 
 PATCH(INPRMA ,LOW , 1, 1, 7, 14, 1, 1, 1, 1)
 COVAL(INPRMA ,U1 ,0. ,10. )
 COVAL(INPRMA ,V1 ,0. ,0. )
 COVAL(INPRMA ,W1 ,0. ,21.355698 )
 COVAL(INPRMA ,GAS , FIXFLU ,22.33717 )
 COVAL(INPRMA ,KE ,0. ,1.140164 )
 COVAL(INPRMA ,EP ,0. ,15.750166 )
 COVAL(INPRMA ,H1 ,0. ,3.724111E+05 )
 COVAL(INPRMA ,MIXF,0. ,0. )
 
 PATCH(INPRMC ,LOW , 1, 1, 7, 14, 1, 1, 1, 1)
 COVAL(INPRMC ,U2 ,0. ,10. )
 COVAL(INPRMC ,V2 ,0. ,0. )
 COVAL(INPRMC ,W2 ,0. ,21.355698 )
 COVAL(INPRMC ,FUE , FIXFLU ,12.724513 )
 COVAL(INPRMC ,H2 ,0. ,3.756981E+07 )
 COVAL(INPRMC ,C2 ,0. ,1. )
 
 PATCH(INSEC ,LOW , 1, 1, 15, 24, 1, 1, 1, 1)
 COVAL(INSEC ,U1 ,0. ,10. )
 COVAL(INSEC ,V1 ,0. ,0. )
 COVAL(INSEC ,W1 ,0. ,28.07303 )
 COVAL(INSEC ,GAS , FIXFLU ,18.663469 )
 COVAL(INSEC ,KE ,0. ,1.970238 )
 COVAL(INSEC ,EP ,0. ,13.250978 )
 COVAL(INSEC ,H1 ,0. ,5.863E+05 )
 COVAL(INSEC ,MIXF,0. ,0. )
 
 PATCH(INTER ,LOW , 1, 1, 1, 6, 1, 1, 1, 1)
 COVAL(INTER ,U1 ,0. ,0. )
 COVAL(INTER ,V1 ,0. ,0. )
 COVAL(INTER ,W1 ,0. ,41.934875 )
 COVAL(INTER ,GAS , FIXFLU ,27.869837 )
 COVAL(INTER ,KE ,0. ,4.396334 )
 COVAL(INTER ,EP ,0. ,74.533218 )
 COVAL(INTER ,H1 ,0. ,5.863E+05 )
 COVAL(INTER ,MIXF,0. ,0. )
 
 PATCH(OUTLET ,HIGH , 1, 1, 1, 32, 25, 25, 1, 1)
 COVAL(OUTLET ,U1 ,0. ,0. )
 COVAL(OUTLET ,U2 ,0. ,0. )
 COVAL(OUTLET ,V1 ,0. ,0. )
 COVAL(OUTLET ,V2 ,0. ,0. )
 COVAL(OUTLET ,W1 ,0. ,0. )
 COVAL(OUTLET ,W2 ,0. ,0. )
 COVAL(OUTLET ,GAS ,100. ,0. )
 COVAL(OUTLET ,FUE ,1.35E+06 ,0. )
 COVAL(OUTLET ,KE ,0. , SAME )
 COVAL(OUTLET ,EP ,0. , SAME )
 COVAL(OUTLET ,H1 ,0. , SAME )
 COVAL(OUTLET ,H2 ,0. , SAME )
 COVAL(OUTLET ,MIXF,0. , SAME )
 COVAL(OUTLET ,C2 ,0. ,0. )
 COVAL(OUTLET ,RADY,0. , SAME )
 COVAL(OUTLET ,RADZ,0. , SAME )
 
 PATCH(WALNR ,NORTH , 1, 1, 24, 24, 1, 5, 1, 1)
 COVAL(WALNR ,RADY,1. ,1.21226E+04 )
 
 PATCH(WALLR ,LOW , 1, 1, 25, 32, 6, 6, 1, 1)
 COVAL(WALLR ,RADZ,1. ,1.21226E+04 )
 IURVAL = -1
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 350 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 ISWR1 = 1 ;ISWR2 = 10000
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(U1)=10
 LITER(U2)=10 ;LITER(V1)=10
 LITER(V2)=10 ;LITER(W1)=10
 LITER(W2)=10 ;LITER(GAS)=1
 LITER(FUE)=1 ;LITER(SHAD)=1
 LITER(KE)=20 ;LITER(EP)=20
 LITER(H1)=20 ;LITER(H2)=20
 LITER(MIXF)=20 ;LITER(C2)=20
 LITER(RADY)=2 ;LITER(RADZ)=2
 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03
 ENDIT(U2)=1.0E-03 ;ENDIT(V1)=1.0E-03
 ENDIT(V2)=1.0E-03 ;ENDIT(W1)=1.0E-03
 ENDIT(W2)=1.0E-03 ;ENDIT(GAS)=1.0E-03
 ENDIT(FUE)=1.0E-03 ;ENDIT(SHAD)=1.0E-03
 ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
 ENDIT(H1)=1.0E-03 ;ENDIT(H2)=1.0E-03
 ENDIT(MIXF)=1.0E-03 ;ENDIT(C2)=1.0E-03
 ENDIT(RADY)=1.0E-03 ;ENDIT(RADZ)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,1.)
 RELAX(U2,FALSDT,1.)
 RELAX(V1,FALSDT,1.)
 RELAX(V2,FALSDT,1.)
 RELAX(W1,FALSDT,1.)
 RELAX(W2,FALSDT,1.)
 RELAX(GAS,LINRLX,0.4)
 RELAX(FUE,LINRLX,0.4)
 RELAX(SHAD,LINRLX,0.4)
 RELAX(KE,LINRLX,0.5)
 RELAX(EP,LINRLX,0.5)
 RELAX(H1,FALSDT,1.0E+09)
 RELAX(H2,FALSDT,1.0E+09)
 RELAX(MIXF,FALSDT,1.0E+09)
 RELAX(C2,FALSDT,1.0E+09)
 RELAX(PRPS,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(RADY,FALSDT,1.0E+09)
 RELAX(RADZ,FALSDT,1.0E+09)
 RELAX(ENUT,LINRLX,1.)
 RELAX(RHO1,LINRLX,1.)
 RELAX(TMP2,LINRLX,1.)
 RELAX(TMP1,LINRLX,1.)
 RELAX(YH2O,LINRLX,1.)
 RELAX(YH2,LINRLX,1.)
 RELAX(YN2,LINRLX,1.)
 RELAX(YCO2,LINRLX,1.)
 RELAX(YCO,LINRLX,1.)
 RELAX(YO2,LINRLX,1.)
 RELAX(CFIP,LINRLX,1.)
 RELAX(MDOT,LINRLX,1.)
 KELIN = 3
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(U1)=1.0E+06 ;VARMIN(U1)=-1.0E+06
 VARMAX(U2)=1.0E+06 ;VARMIN(U2)=-1.0E+06
 VARMAX(V1)=1.0E+06 ;VARMIN(V1)=-1.0E+06
 VARMAX(V2)=1.0E+06 ;VARMIN(V2)=-1.0E+06
 VARMAX(W1)=1.0E+06 ;VARMIN(W1)=-1.0E+06
 VARMAX(W2)=1.0E+06 ;VARMIN(W2)=-1.0E+06
 VARMAX(GAS)=1. ;VARMIN(GAS)=1.0E-06
 VARMAX(FUE)=1. ;VARMIN(FUE)=1.0E-12
 VARMAX(SHAD)=1. ;VARMIN(SHAD)=1.0E-12
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 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(C2)=1. ;VARMIN(C2)=0.
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(RADY)=1.0E+10 ;VARMIN(RADY)=-1.0E+10
 VARMAX(RADZ)=1.0E+10 ;VARMIN(RADZ)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=0.069686
 VARMAX(TMP2)=1.0E+10 ;VARMIN(TMP2)=-1.0E+10
 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10
 VARMAX(YH2O)=1.0E+10 ;VARMIN(YH2O)=-1.0E+10
 VARMAX(YH2)=1.0E+10 ;VARMIN(YH2)=-1.0E+10
 VARMAX(YN2)=1.0E+10 ;VARMIN(YN2)=-1.0E+10
 VARMAX(YCO2)=1.0E+10 ;VARMIN(YCO2)=-1.0E+10
 VARMAX(YCO)=1.0E+10 ;VARMIN(YCO)=-1.0E+10
 VARMAX(YO2)=1.0E+10 ;VARMIN(YO2)=-1.0E+10
 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
 GENK = T
 PARSOL = F
 CONWIZ = T
 ISG62 = 1
 SPEDAT(SET,LONGNAME,epke,C,NGNAME,epke,c,micro-mixing_rate)
 SPEDAT(SET,LONGNAME,rady,C,NGNAME,rady,c,y-direction_flux_sum)
 SPEDAT(SET,LONGNAME,radz,C,NGNAME,radz,c,z-direction_flux_sum)
 SPEDAT(SET,MATERIAL,199,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=57.52 ;EX(U1)=2.407
 EX(U2)=2.496 ;EX(V1)=1.866
 EX(V2)=1.714 ;EX(W1)=35.029999
 EX(W2)=29.549999 ;EX(GAS)=0.95
 EX(FUE)=0. ;EX(SHAD)=0.
 EX(KE)=10.85 ;EX(EP)=231.100006
 EX(H1)=1.147E+06 ;EX(H2)=3.604E+07
 EX(MIXF)=0.02047 ;EX(C2)=0.
 EX(PRPS)=0. ;EX(NPOR)=0.9913
 EX(EPKE)=0. ;EX(RADY)=1.106E+05
 EX(RADZ)=2.517E+05 ;EX(ENUT)=0.05748
 EX(RHO1)=0.3856 ;EX(TMP2)=635.400024
 EX(TMP1)=1042. ;EX(YH2O)=0.
 EX(YH2)=0. ;EX(YN2)=0.7139
 EX(YCO2)=0.0713 ;EX(YCO)=0.
 EX(YO2)=0.1556 ;EX(CFIP)=0.
 EX(MDOT)=0.
 ************************************************************
  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(U2,Y,N,Y,Y,Y,Y)
 OUTPUT(V1,Y,N,Y,Y,Y,Y)
 OUTPUT(V2,Y,N,Y,Y,Y,Y)
 OUTPUT(W1,Y,N,Y,Y,Y,Y)
 OUTPUT(W2,Y,N,Y,Y,Y,Y)
 OUTPUT(GAS,Y,N,Y,Y,Y,Y)
 OUTPUT(FUE,Y,N,Y,Y,Y,Y)
 OUTPUT(SHAD,Y,N,Y,Y,Y,Y)
 OUTPUT(KE,Y,N,Y,Y,Y,Y)
 OUTPUT(EP,Y,N,Y,Y,Y,Y)
 OUTPUT(H1,Y,N,Y,Y,Y,Y)
 OUTPUT(H2,Y,N,Y,Y,Y,Y)
 OUTPUT(MIXF,Y,N,Y,Y,Y,Y)
 OUTPUT(C2,N,N,N,N,N,N)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 OUTPUT(NPOR,Y,N,Y,N,N,N)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(RADY,Y,N,Y,Y,Y,Y)
 OUTPUT(RADZ,Y,N,Y,Y,Y,Y)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 OUTPUT(RHO1,Y,Y,Y,Y,Y,Y)
 OUTPUT(TMP2,Y,Y,Y,Y,Y,Y)
 OUTPUT(TMP1,Y,Y,Y,Y,Y,Y)
 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(YCO,Y,N,Y,N,N,N)
 OUTPUT(YO2,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 = 16 ;IZMON = 18
 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 = 10
 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