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

  2 2 0.9

  up z
  use patgeo
  vec x 3 sh
  pause;cl
  use patgeo
  con tmp1 z 1 fil;.001
  con tmp1 x 2 fil;.001
  con tmp1 x 4 fil;.001
  con tmp1 z m fil;.001
  con tmp1 z 10 fil;.001
  vec x 3
  msg(Gas temperature contours
  pause;cl

  use patgeo
  con yhcn z 1 fil;.001
  con yhcn x 2 fil;.001
  con yhcn x 4 fil;.001
  con yhcn z m fil;.001
  con yhcn z 10 fil;.001
  vec x 3
  msg(HCN mass fractions
  pause;cl

  use patgeo
  con ynox z 1 fil;.001
  con ynox x 2 fil;.001
  con ynox x 4 fil;.001
  con ynox z m fil;.001
  con ynox z 10 fil;.001
  vec x 3
  msg(NOX mass fractions
  pause;cl

  use patgeo
  con tmp2 z 1 fil;.001
  con tmp2 x 2 fil;.001
  con tmp2 x 4 fil;.001
  con tmp2 z m fil;.001
  con tmp2 z 10 fil;.001
  vec x 3
  msg(Solid temperature contours
  pause;cl

  use patgeo
  con VRAT z 1 fil;.001
  con VRAT x 2 fil;.001
  con VRAT x 4 fil;.001
  con VRAT z m fil;.001
  con VRAT z 10 fil;.001
  vec x 3
  msg(RATe of Volatalization, kg/s
  pause;cl

  use patgeo
  con VCOM z 1 fil;.001
  con VCOM x 2 fil;.001
  con VCOM x 4 fil;.001
  con VCOM z m fil;.001
  con VCOM z 10 fil;.001
  vec x 3
  msg(Rate of Volatiles, CH4, COMbustion, kg/s
  pause;cl

  use patgeo
  con CCOM z 1 fil;.001
  con CCOM x 2 fil;.001
  con CCOM x 4 fil;.001
  con CCOM z m fil;.001
  con CCOM z 10 fil;.001
  vec x 3
  msg(Rate of CO COMbustion, kg/s
  pause;cl

  use patgeo
  con BURN z 1 fil;.001
  con BURN x 2 fil;.001
  con BURN x 4 fil;.001
  con BURN z m fil;.001
  con BURN z 10 fil;.001
  vec x 3
  msg(Rate of char BURNing, kg/s
  pause;cl

  use patgeo
  con YCH4 z 1 fil;.001
  con YCH4 x 2 fil;.001
  con YCH4 x 4 fil;.001
  con YCH4 z m fil;.001
  con YCH4 z 10 fil;.001
  vec x 3
  msg(Volatiles, CH4, mass fractions
  pause;cl
  use patgeo
  con YCO2 z 1 fil;.001
  con YCO2 x 2 fil;.001
  con YCO2 x 4 fil;.001
  con YCO2 z m fil;.001
  con YCO2 z 10 fil;.001
  vec x 3
  msg(CO2  mass fractions
  pause;cl
  use patgeo
  con YCO z 1 fil;.001
  con YCO x 2 fil;.001
  con YCO x 4 fil;.001
  con YCO z m fil;.001
  con YCO z 10 fil;.001
  vec x 3
  msg(CO  mass fractions
  pause;cl
  use patgeo
  con YH2O z 1 fil;.001
  con YH2O x 2 fil;.001
  con YH2O x 4 fil;.001
  con YH2O z m fil;.001
  con YH2O z 10 fil;.001
  vec x 3
  msg(H2O  mass fractions
  pause;cl
  use patgeo
  con YO2 z 1 fil;.001
  con YO2 x 2 fil;.001
  con YO2 x 4 fil;.001
  con YO2 z m fil;.001
  con YO2 z 10 fil;.001
  vec x 3
  msg(O2  mass fractions
  pause;cl
  use patgeo
  con YN2 z 1 fil;.001
  con YN2 x 2 fil;.001
  con YN2 x 4 fil;.001
  con YN2 z m fil;.001
  con YN2 z 10 fil;.001
  vec x 3
  msg(N2  mass fractions

  pause;cl
  use patgeo
  con COL2 z 1 fil;.001
  con COL2 x 2 fil;.001
  con COL2 x 4 fil;.001
  con COL2 z m fil;.001
  con COL2 z 10 fil;.001
  vec x 3
  msg(COAL particle content
  pause;cl
  use patgeo
  con CHA2 z 1 fil;.001
  con CHA2 x 2 fil;.001
  con CHA2 x 4 fil;.001
  con CHA2 z m fil;.001
  con CHA2 z 10 fil;.001
  vec x 3
  msg(CHAR particle content
  pause;cl
  use patgeo
  con ASH2 z 1 fil;.001
  con ASH2 x 2 fil;.001
  con ASH2 x 4 fil;.001
  con ASH2 z m fil;.001
  con ASH2 z 10 fil;.001
  vec x 3
  msg(ASH particle content
  pause;cl
  use patgeo
  con R2 z 1 fil;.001
  con R2 x 2 fil;.001
  con R2 x 4 fil;.001
  con R2 z m fil;.001
  con R2 z 10 fil;.001
  vec x 3
  msg(Particle volume fractions
  pause;cl
  use patgeo
  con VREL z 1 fil;.001
  con VREL x 2 fil;.001
  con VREL x 4 fil;.001
  con VREL z m fil;.001
  con VREL z 10 fil;.001
  vec x 3
  msg(Particle relative velocity
  pause;cl
  use patgeo
  con SIZE z 1 fil;.001
  con SIZE x 2 fil;.001
  con SIZE x 4 fil;.001
  con SIZE z m fil;.001
  con SIZE z 10 fil;.001
  vec x 3
  msg(Particle diameter
  pause;cl
  use patgeo
  con REYN z 1 fil;.001
  con REYN x 2 fil;.001
  con REYN x 4 fil;.001
  con REYN z m fil;.001
  con REYN z 10 fil;.001
  vec x 3
  msg(Particle Reynolds number
  pause;cl
  use patgeo
  con CD z 1 fil;.001
  con CD x 2 fil;.001
  con CD x 4 fil;.001
  con CD z m fil;.001
  con CD z 10 fil;.001
  vec x 3
  msg(Particle drag coefficient
  pause;cl
  use patgeo
  con NUSS z 1 fil;.001
  con NUSS x 2 fil;.001
  con NUSS x 4 fil;.001
  con NUSS z m fil;.001
  con NUSS z 10 fil;.001
  vec x 3
  msg(Particle Nusselt number
  pause;cl
  use patgeo
  con HCOF z 1 fil;.001
  con HCOF x 2 fil;.001
  con HCOF x 4 fil;.001
  con HCOF z m fil;.001
  con HCOF z 10 fil;.001
  vec x 3
  msg(Interphase heat transfer coefficient
  pause;cl
  use patgeo
  con DEN1 z 1 fil;.001
  con DEN1 x 2 fil;.001
  con DEN1 x 4 fil;.001
  con DEN1 z m fil;.001
  con DEN1 z 10 fil;.001
  vec x 3
  msg(Gas density
  pause;cl
  use patgeo
  con EPKE z 1 fil;.001
  con EPKE x 2 fil;.001
  con EPKE x 4 fil;.001
  con EPKE z m fil;.001
  con EPKE z 10 fil;.001
  vec x 3
  msg(Reciprocal time scale. 1/s
  pause;cl
  ENDUSE
  DISPLAY

  
  FUEL-DUST FLAMES in a FURNACE


  Wall-fired furnace: 24 burners

  Dry coal.

  1 Kinetically released volatiles;
  2 Combustion of:

  * Volatiles, CH4,  by EBU with a provison for kinetic ;
  * Char, C(s), by harmonic blend of kinetic and diffusion and
  * Carbon monoxide, CO, by harmonic blend of kinetic and EBU.

  3 Thermal NOX and
  4 Fuel NOX formations.
  5 Radiosity model for radiating particles

  ENDDIS
  PLANTBEGIN
NAMSAT=MOSG
PRNDTL(CRAD)=GRND
   LAMPR(CRAD)=:ENUL:*DEN1*3./4.*(0.5+5.e-1*6.*R2/SIZE)
    * Interphase mass transfer
       VRAT-volatile release; BURN-char burning
   INTMDT=VRAT+BURN
    * Phase temperatures
   TEMP1=H1/(1059.+0.25*(TMP1-300.))
   TEMP2=H2/1800.

   =========================================
    * Devolatilization related sources/sinks

      # Reaction: 1kg Coal -> Ykg Vol + (1-Y)kg Char
PATCH(PLSR1,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
      ** Coal
   VAL=-VRAT/RHO2A
COVAL(PLSR1,COL2,FIXFLU,GRND)
      ** Char
   VAL=VRAT*(1.-RHO2A)/RHO2A
COVAL(PLSR1,CHA2,FIXFLU,GRND)

    ======================================
           CH4-volatile oxidation
    * YCH4 supply via volatile release
PATCH(PLSR5,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
   VAL=VRAT*1.0
COVAL(PLSR5,YCH4,FIXFLU,GRND)

   Reaction:   2CH4  + 3O2 ->  2CO +4H2O
PATCH(PLSR6,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
  CO=96./32*VCOM/(YO2+tiny)
COVAL(PLSR6,YO2,GRND,0.0)
   CO=VCOM/(YCH4+tiny)
COVAL(PLSR6,YCH4,GRND,0.0)
  VAL=72./32.*VCOM
COVAL(PLSR6,YH2O,FIXFLU,GRND)
  VAL=56./32.*VCOM
COVAL(PLSR6,YCO,FIXFLU,GRND)
     ==========================================
         Carbon-Char oxidation
   Reaction:   C(s)  + (x+1)/2*O2 ->  xCO2 + (1-x)CO
    * Sink of char,C(s), in a 2nd phase
PATCH(PLSR8,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
   VAL=-BURN
COVAL(PLSR8,CHA2,FIXFLU,GRND)
    * Sink/sources of 1st phase components
PATCH(PLSR9,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
   VAL=44./12.*XXX*BURN
COVAL(PLSR9,YCO2,FIXFLU,GRND)
   VAL=28./12.*(1.-XXX)*BURN
COVAL(PLSR9,YCO,FIXFLU,GRND)
   VAL=-32./12.*(1.+XXX)/2.*BURN
COVAL(PLSR9,YO2,FIXFLU,GRND)
     ===========================================
         Carbon monoxide oxidation
   Reaction:   2CO +  O2 -> 2CO2
PATCH(PLSR7,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
   CO=CCOM/(YCO+tiny)
COVAL(PLSR7,YCO,GRND,0.0)
   CO=CCOM*32./56./(YO2+tiny)
COVAL(PLSR7,YO2,GRND,0.0)
  VAL=CCOM*88./56.
COVAL(PLSR7,YCO2,FIXFLU,GRND)
     ==========================================
    * Enthalpy sources due to heats of combustion
       # Volatiles
PATCH(PLSR3,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
   VAL=VCOM*5.5e7+CCOM*7.45e6
COVAL(PLSR3,H1,FIXFLU,GRND)
       # Char
PATCH(PLSR10,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
   VAL=BURN*(XXX*3.28e7+(1.-XXX)*9.21e6)
COVAL(PLSR10,H1,FIXFLU,GRND)

    * Enthalpy sources due to interphase heat transfer
PATCH(PLSR4,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
   CO=HCOF/(1059.+0.25*(TMP1-300.))
   VAL=H2*(1059.+0.25*(TMP1-300.)) /1800.
COVAL(PLSR4,H1,GRND,GRND)
   CO=HCOF/1800.
   VAL=1800./(1059.+0.25*(TMP1-300.)) *H1
COVAL(PLSR4,H2,GRND,GRND)

    * Enthalpy sources due to radiation
PATCH(RADSOR,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
   VAL=4.*0.5*(CRAD-5.68e-8*TMP1**4)
COVAL(RADSOR,H1,FIXFLU,GRND)
   VAL=4.*5.e-1*6.*R2/SIZE*(CRAD-5.68e-8*TMP2**4)
COVAL(RADSOR,H2,FIXFLU,GRND)
    * Radiosity source
PATCH(RDIOSITY,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
   CO=4.*(5.e-1*6.*R2/SIZE+0.5)
   VAL=5.68e-8*(5.e-1*6.*R2/SIZE*TMP2**4+0.5*TMP1**4 )/$
       (5.e-1*6.*R2/SIZE+0.5)
COVAL(RDIOSITY,CRAD,GRND,GRND)
     ============================================
   * Sources to counter the transfer to phase 1
PATCH(PLSR90,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
   VAL=CMDO*COL2
COVAL(PLSR90,COL2,FIXFLU,GRND)
   VAL=CMDO*CHA2
COVAL(PLSR90,CHA2,FIXFLU,GRND)
     ===========================================
PATCH(THERMNOX,VOLUME,1,NX,1,NY,1,NZ,1,1)
     VAL=RNOX
COVAL(THERMNOX,YNOX,FIXFLU,GRND)

PATCH(FUELHCN,CELL,1,NX,1,NY,1,NZ,1,LSTEP)
   VAL=27./18.*(1.e-2*VRAT+1.e-2*BURN)
COVAL(FUELHCN,YHCN,FIXFLU,GRND)

PATCH(FUELNOX,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
   CO=3.e12*YHCN/EXP(2.5e8/(8314.*AMAX1(300.,TMP1)))
   VAL=0.0333*YO2
COVAL(FUELNOX,YNOX,GRND,GRND)

PATCH(NOXRED,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
   CO=6.*R2/SIZE*40./EXP(1.5e8/(8314.*AMAX1(300.,TMP2)))
COVAL(NOXRED,YNOX,GRND,0.0)
     ===========================================
store(HCOF,VRAT,VCOM,VCOD,VCOK,CCOM,CCOD,CCOK,WMIX,KD,KC,KT,BURN)
store(BRN1,BRN2)
    # Char burning rate constants:
           ** Mixture molecular mass
   WMIX=28.8
           ** Diffusion
   KD=1.55*NUSS*1.78e-5*12./8314./AMAX1(300.,TMP2)/SIZE
           ** Kinetic
   KC=0.1309/EXP(26850./4./AMAX1(300.,TMP2))
           ** Harmonic blending
   KT=1./(1./KD+1./KC)

    # Rates:

        # Volatile release, kg/s
   VRAT=2.e3*RHO2A*COL2/EXP(2829./AMAX1(300.,TMP2))*$
  RHO2*R2*VOL*RHO2A

        # Volatile combustion, kg/s
           ** EBU
   VCOD=4.*EPKE*AMIN1(YCH4,32.*YO2/96.)
           ** Kinetic
   VCOK=1.15e9*DEN1*YO2**1.3/EXP(24444./AMAX1(300.,TMP1))/$
               (YCH4**0.3+tiny)
           ** Harmonic blending
   VCOM=1./(1./VCOD+1./VCOK)*VOL*R1*DEN1

        # Carbon monoxide,CO, combustion, kg/s
           ** EBU
   CCOD=4.*EPKE*AMIN1(YCO,32.*YO2/56.)
           ** Kinetic
   CCOK=5.42e9*DEN1*YCO*YO2**0.25*YH2O**0.5/              E$
XP(15152./AMAX1(300.,TMP1))
           ** Harmonic blending
   CCOM=1./(1./CCOD+1./CCOK)*VOL*R1*DEN1

        # Char,C, burning out, kg/s
   BRN1=6./SIZE*KT*(P1+PRESS0)*YO2*WMIX/32.*R2*VOL
   BRN2=4.*DEN1*R1*VOL*EPKE*AMIN1(CHA2*MASS2/MASS1,$
12./(32.*(1.+XXX)/2.)*YO2)
   BURN=1./(1./BRN1+1./BRN2)
        # CO/CO2 of C-combustion
store(MCCO)
   MCCO=2500.*EXP(-6249./AMAX1(300.,TMP2))
   XXX=1./(1.+44./28.*MCCO)

    # Interphase heat transfer coefficient
   NUSS=2.+0.65*REYN**0.5
   HCOF=6.*R2*0.0458*NUSS/SIZE**2.

    # Ash contents of 2nd phase
   ASH2=1.-CHA2-COL2

    # Nitrogen contents of 1st phase
   YN2=1.-YCH4-YO2-YCO-YCO2-YH2O

store(FK1,MOX,RNOX)
               *** Reaction constant
   FK1=1.8e8*EXP(-38370./AMAX1(TMP1,300.))
               *** Mole O atom concentration
   MOX=3.97e5*(YO2*DEN1/32./AMAX1(TMP1,300.))**.5     *EXP($
-31090./AMAX1(TMP1,300.))
               *** NOX reaction rate
   RNOX=30.*2.*FK1*MOX*YN2*DEN1/28.

  PLANTEND
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Wall-fired model furnace                )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 5
 XULAST =1.
 XFRAC(1)=0.2 ;XFRAC(2)=0.4
 XFRAC(3)=0.6 ;XFRAC(4)=0.8
 XFRAC(5)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 5
 YVLAST =1.
 YFRAC(1)=0.2 ;YFRAC(2)=0.4
 YFRAC(3)=0.6 ;YFRAC(4)=0.8
 YFRAC(5)=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(3)=U1
 NAME(4)=U2 ;NAME(5)=V1
 NAME(6)=V2 ;NAME(7)=W1
 NAME(8)=W2 ;NAME(9)=R1
 NAME(10)=R2 ;NAME(11)=RS
 NAME(12)=KE ;NAME(13)=EP
 NAME(14)=H1 ;NAME(15)=H2
 NAME(105)=CRAD ;NAME(106)=MCCO
 NAME(107)=KT ;NAME(108)=XXX
 NAME(109)=BRN2 ;NAME(110)=BRN1
 NAME(111)=RNOX ;NAME(112)=MOX
 NAME(113)=FK1 ;NAME(114)=CCOK
 NAME(115)=CCOD ;NAME(116)=VCOD
 NAME(117)=BURN ;NAME(118)=WMIX
 NAME(119)=CCOM ;NAME(120)=VCOM
 NAME(121)=YN2 ;NAME(122)=ASH2
 NAME(123)=VCOK ;NAME(124)=YH2O
 NAME(125)=YCO2 ;NAME(126)=YCO
 NAME(127)=YO2 ;NAME(128)=YCH4
 NAME(129)=HCOF ;NAME(130)=KC
 NAME(131)=KD ;NAME(132)=SIZE
 NAME(133)=DEN1 ;NAME(134)=APRJ
 NAME(135)=COL2 ;NAME(136)=CHA2
 NAME(137)=REYN ;NAME(138)=CD
 NAME(139)=VREL ;NAME(140)=NUSS
 NAME(141)=TMP2 ;NAME(142)=TMP1
 NAME(143)=CFIP ;NAME(144)=VRAT
 NAME(145)=CMDO ;NAME(146)=EL1
 NAME(147)=ENUT ;NAME(148)=EPKE
 NAME(149)=YHCN ;NAME(150)=YNOX
    * 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(R1,Y,Y,N,Y,N,Y)
 SOLUTN(R2,Y,Y,N,Y,N,Y)
 SOLUTN(RS,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,N,N,N,Y)
 SOLUTN(H2,Y,Y,N,N,N,Y)
 SOLUTN(CRAD,Y,Y,N,N,N,Y)
 SOLUTN(MCCO,Y,N,N,N,N,Y)
 SOLUTN(KT,Y,N,N,N,N,Y)
 SOLUTN(XXX,Y,N,N,N,N,Y)
 SOLUTN(BRN2,Y,N,N,N,N,Y)
 SOLUTN(BRN1,Y,N,N,N,N,Y)
 SOLUTN(RNOX,Y,N,N,N,N,Y)
 SOLUTN(MOX,Y,N,N,N,N,Y)
 SOLUTN(FK1,Y,N,N,N,N,Y)
 SOLUTN(CCOK,Y,N,N,N,N,Y)
 SOLUTN(CCOD,Y,N,N,N,N,Y)
 SOLUTN(VCOD,Y,N,N,N,N,Y)
 SOLUTN(BURN,Y,N,N,N,N,Y)
 SOLUTN(WMIX,Y,N,N,N,N,Y)
 SOLUTN(CCOM,Y,N,N,N,N,Y)
 SOLUTN(VCOM,Y,N,N,N,N,Y)
 SOLUTN(YN2,Y,N,N,N,N,Y)
 SOLUTN(ASH2,Y,N,N,N,N,Y)
 SOLUTN(VCOK,Y,N,N,N,N,Y)
 SOLUTN(YH2O,Y,Y,N,N,N,Y)
 SOLUTN(YCO2,Y,Y,N,N,N,Y)
 SOLUTN(YCO,Y,Y,N,N,N,Y)
 SOLUTN(YO2,Y,Y,N,N,N,Y)
 SOLUTN(YCH4,Y,Y,N,N,N,Y)
 SOLUTN(HCOF,Y,N,N,N,N,Y)
 SOLUTN(KC,Y,N,N,N,N,Y)
 SOLUTN(KD,Y,N,N,N,N,Y)
 SOLUTN(SIZE,Y,N,N,N,N,Y)
 SOLUTN(DEN1,Y,N,N,N,N,Y)
 SOLUTN(APRJ,Y,N,N,N,N,Y)
 SOLUTN(COL2,Y,Y,N,N,N,Y)
 SOLUTN(CHA2,Y,Y,N,N,N,Y)
 SOLUTN(REYN,Y,N,N,N,N,Y)
 SOLUTN(CD,Y,N,N,N,N,Y)
 SOLUTN(VREL,Y,N,N,N,N,Y)
 SOLUTN(NUSS,Y,N,N,N,N,Y)
 SOLUTN(TMP2,Y,N,N,N,N,Y)
 SOLUTN(TMP1,Y,N,N,N,N,Y)
 SOLUTN(CFIP,Y,N,N,N,N,Y)
 SOLUTN(VRAT,Y,N,N,N,N,Y)
 SOLUTN(CMDO,Y,N,N,N,N,Y)
 SOLUTN(EL1,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 SOLUTN(EPKE,Y,N,N,N,N,Y)
 SOLUTN(YHCN,Y,Y,N,N,N,Y)
 SOLUTN(YNOX,Y,Y,N,N,N,Y)
 DEN1 = 133
 VIST = 147
 TEMP1 = 142
 TEMP2 = 141
 LEN1 = 146
 ************************************************************
  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(R1,Y,Y,Y,Y,Y,Y)
 TERMS(R2,Y,Y,Y,Y,N,Y)
 TERMS(RS,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(CRAD,N,N,Y,Y,Y,N)
 TERMS(YH2O,N,Y,Y,Y,Y,N)
 TERMS(YCO2,N,Y,Y,Y,Y,N)
 TERMS(YCO,N,Y,Y,Y,Y,N)
 TERMS(YO2,N,Y,Y,Y,Y,N)
 TERMS(YCH4,N,Y,Y,Y,Y,N)
 TERMS(COL2,N,Y,Y,Y,N,N)
 TERMS(CHA2,N,Y,Y,Y,N,N)
 TERMS(YHCN,N,Y,Y,Y,Y,N)
 TERMS(YNOX,N,Y,Y,Y,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 = GRND5 ;TMP1 = GRND
 EL1 = GRND4
 TSURR =0. ;TEMP0 =0.
 PRESS0 =1.0E+05
 DVO1DT =0. ;DRH1DP = GRND5
 RHO1A =0. ;RHO1B =3.496503E-03
 RHO1C =0.714286
 TMP1A =0. ;TMP1B =0.
 TMP1C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0.
 EL1C =0.
 RHO2 =1000. ;TMP2 = GRND
 EL2 =0.
 DVO2DT =0. ;DRH2DP =0.
 TMP2A =0. ;TMP2B =0.
 TMP2C =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(R1)=-2.022E+04 ;PHINT(R2)=-2.022E+04
 PHINT(RS)=-2.022E+04 ;PHINT(H1)=-2.022E+04
 PHINT(H2)=-2.022E+04
 PRNDTL(U1)=1. ;PRNDTL(U2)=1.
 PRNDTL(V1)=1. ;PRNDTL(V2)=1.
 PRNDTL(W1)=1. ;PRNDTL(W2)=1.
 PRNDTL(R1)=1.0E+09 ;PRNDTL(R2)=1.0E+09
 PRNDTL(RS)=1. ;PRNDTL(KE)=1.
 PRNDTL(EP)=1. ;PRNDTL(H1)=1.
 PRNDTL(H2)=1. ;PRNDTL(CRAD)=GRND
 PRNDTL(YH2O)=1. ;PRNDTL(YCO2)=1.
 PRNDTL(YCO)=1. ;PRNDTL(YO2)=1.
 PRNDTL(YCH4)=1. ;PRNDTL(COL2)=1.
 PRNDTL(CHA2)=1. ;PRNDTL(YHCN)=1.
 PRNDTL(YNOX)=1.
 PRT(U1)=1. ;PRT(U2)=1.
 PRT(V1)=1. ;PRT(V2)=1.
 PRT(W1)=1. ;PRT(W2)=1.
 PRT(R1)=1.0E+09 ;PRT(R2)=1.0E+09
 PRT(RS)=1. ;PRT(KE)=1.
 PRT(EP)=1.314 ;PRT(H1)=1.
 PRT(H2)=1. ;PRT(CRAD)=1.0E+09
 PRT(YH2O)=1. ;PRT(YCO2)=1.
 PRT(YCO)=1. ;PRT(YO2)=1.
 PRT(YCH4)=1. ;PRT(COL2)=1.
 PRT(CHA2)=1. ;PRT(YHCN)=1.
 PRT(YNOX)=1.
 CP1 =1050. ;CP2 =1800.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 CFIPS = GRND7
 RLOLIM =1.0E-09 ;CMDOT = GRND
 CFIPA =1.0E-09 ;CFIPB =2.5E-05
 CFIPC =100. ;CFIPD =0.
 CMDTA =0. ;CMDTB =0.
 CMDTC =0. ;CMDTD =0.
 CINT(P1)=1. ;CINT(U1)=1.
 CINT(U2)=1. ;CINT(V1)=1.
 CINT(V2)=1. ;CINT(W1)=1.
 CINT(W2)=1. ;CINT(R1)=1.
 CINT(R2)=1. ;CINT(RS)=1.
 CINT(H1)=1.0E-20 ;CINT(H2)=1.0E-20
 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(R1)=-2.022E+04
 PHINT(R2)=-2.022E+04 ;PHINT(RS)=-2.022E+04
 PHINT(H1)=-2.022E+04 ;PHINT(H2)=-2.022E+04
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10
 FIINIT(U2)=1.0E-10 ;FIINIT(V1)=1.0E-10
 FIINIT(V2)=1.0E-10 ;FIINIT(W1)=1.0E-10
 FIINIT(W2)=1.0E-10 ;FIINIT(R1)=0.5
 FIINIT(R2)=0.5 ;FIINIT(RS)=0.5
 FIINIT(KE)=0.25 ;FIINIT(EP)=1.437814
 FIINIT(H1)=5.25E+05 ;FIINIT(H2)=9.0E+05
 FIINIT(CRAD)=0. ;FIINIT(MCCO)=1.0E-10
 FIINIT(KT)=1.0E-10 ;FIINIT(XXX)=1.0E-10
 FIINIT(BRN2)=1.0E-10 ;FIINIT(BRN1)=1.0E-10
 FIINIT(RNOX)=1.0E-10 ;FIINIT(MOX)=1.0E-10
 FIINIT(FK1)=1.0E-10 ;FIINIT(CCOK)=1.0E-10
 FIINIT(CCOD)=1.0E-10 ;FIINIT(VCOD)=1.0E-10
 FIINIT(BURN)=1.0E-10 ;FIINIT(WMIX)=1.0E-10
 FIINIT(CCOM)=1.0E-10 ;FIINIT(VCOM)=1.0E-10
 FIINIT(YN2)=1.0E-10 ;FIINIT(ASH2)=1.0E-10
 FIINIT(VCOK)=1.0E-10 ;FIINIT(YH2O)=0.
 FIINIT(YCO2)=0. ;FIINIT(YCO)=0.
 FIINIT(YO2)=0.35 ;FIINIT(YCH4)=0.
 FIINIT(HCOF)=1.0E-10 ;FIINIT(KC)=1.0E-10
 FIINIT(KD)=1.0E-10 ;FIINIT(SIZE)=1.0E-10
 FIINIT(DEN1)=1.0E-10 ;FIINIT(APRJ)=1.0E-10
 FIINIT(COL2)=0.8 ;FIINIT(CHA2)=0.
 FIINIT(REYN)=1.0E-10 ;FIINIT(CD)=1.0E-10
 FIINIT(VREL)=1.0E-10 ;FIINIT(NUSS)=1.0E-10
 FIINIT(TMP2)=1.0E-10 ;FIINIT(TMP1)=1.0E-10
 FIINIT(CFIP)=1.0E-10 ;FIINIT(VRAT)=1.0E-10
 FIINIT(CMDO)=1.0E-10 ;FIINIT(EL1)=1.0E-10
 FIINIT(ENUT)=1.0E-10 ;FIINIT(EPKE)=1.0E-10
 FIINIT(YHCN)=0. ;FIINIT(YNOX)=0.
   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(IN1 ,SOUTH , 4, 4, 1, 1, 2, 2, 1, 1)
 COVAL(IN1 ,CRAD,1. ,142. )
 
 PATCH(IN2 ,SOUTH , 2, 2, 1, 1, 2, 2, 1, 1)
 COVAL(IN2 ,CRAD,1. ,142. )
 
 PATCH(IN3 ,SOUTH , 2, 2, 1, 1, 4, 4, 1, 1)
 COVAL(IN3 ,CRAD,1. ,142. )
 
 PATCH(IN4 ,SOUTH , 4, 4, 1, 1, 4, 4, 1, 1)
 COVAL(IN4 ,CRAD,1. ,142. )
 
 PATCH(IN5 ,SOUTH , 4, 4, 1, 1, 6, 6, 1, 1)
 COVAL(IN5 ,CRAD,1. ,142. )
 
 PATCH(IN6 ,SOUTH , 2, 2, 1, 1, 6, 6, 1, 1)
 COVAL(IN6 ,CRAD,1. ,142. )
 
 PATCH(IN7 ,SOUTH , 4, 4, 1, 1, 8, 8, 1, 1)
 COVAL(IN7 ,CRAD,1. ,142. )
 
 PATCH(IN8 ,SOUTH , 2, 2, 1, 1, 8, 8, 1, 1)
 COVAL(IN8 ,CRAD,1. ,142. )
 
 PATCH(IN9 ,SOUTH , 4, 4, 1, 1, 10, 10, 1, 1)
 COVAL(IN9 ,CRAD,1. ,142. )
 
 PATCH(IN10 ,SOUTH , 2, 2, 1, 1, 10, 10, 1, 1)
 COVAL(IN10 ,CRAD,1. ,142. )
 
 PATCH(IN11 ,SOUTH , 4, 4, 1, 1, 12, 12, 1, 1)
 COVAL(IN11 ,CRAD,1. ,142. )
 
 PATCH(IN12 ,SOUTH , 2, 2, 1, 1, 12, 12, 1, 1)
 COVAL(IN12 ,CRAD,1. ,142. )
 
 PATCH(KESOURCE,PHASEM, 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 PATCH(PLSR1 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(PLSR1 ,COL2, FIXFLU , GRND )
 COVAL(PLSR1 ,CHA2, FIXFLU , GRND )
 
 PATCH(PLSR5 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(PLSR5 ,YCH4, FIXFLU , GRND )
 
 PATCH(PLSR6 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(PLSR6 ,YH2O, FIXFLU , GRND )
 COVAL(PLSR6 ,YCO , FIXFLU , GRND )
 COVAL(PLSR6 ,YO2 , GRND ,0. )
 COVAL(PLSR6 ,YCH4, GRND ,0. )
 
 PATCH(PLSR8 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(PLSR8 ,CHA2, FIXFLU , GRND )
 
 PATCH(PLSR9 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(PLSR9 ,YCO2, FIXFLU , GRND )
 COVAL(PLSR9 ,YCO , FIXFLU , GRND )
 COVAL(PLSR9 ,YO2 , FIXFLU , GRND )
 
 PATCH(PLSR7 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(PLSR7 ,YCO2, FIXFLU , GRND )
 COVAL(PLSR7 ,YCO , GRND ,0. )
 COVAL(PLSR7 ,YO2 , GRND ,0. )
 
 PATCH(PLSR3 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(PLSR3 ,H1 , FIXFLU , GRND )
 
 PATCH(PLSR10 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(PLSR10 ,H1 , FIXFLU , GRND )
 
 PATCH(PLSR4 ,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(PLSR4 ,H1 , GRND , GRND )
 COVAL(PLSR4 ,H2 , GRND , GRND )
 
 PATCH(RADSOR ,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(RADSOR ,H1 , FIXFLU , GRND )
 COVAL(RADSOR ,H2 , FIXFLU , GRND )
 
 PATCH(RDIOSITY,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(RDIOSITY,CRAD, GRND , GRND )
 
 PATCH(PLSR90 ,CELL , 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(PLSR90 ,COL2, FIXFLU , GRND )
 COVAL(PLSR90 ,CHA2, FIXFLU , GRND )
 
 PATCH(THERMNOX,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(THERMNOX,YNOX, FIXFLU , GRND )
 
 PATCH(FUELHCN ,CELL , 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(FUELHCN ,YHCN, FIXFLU , GRND )
 
 PATCH(FUELNOX ,PHASEM, 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(FUELNOX ,YNOX, GRND , GRND )
 
 PATCH(NOXRED ,VOLUME, 1, 5, 1, 5, 1, 20, 1, 1)
 COVAL(NOXRED ,YNOX, GRND ,0. )
 
 PATCH(OB1 ,HIGH , 1, 5, 1, 5, 20, 20, 1, 1)
 COVAL(OB1 ,U1 ,0. ,0. )
 COVAL(OB1 ,U2 ,0. ,0. )
 COVAL(OB1 ,V1 ,0. ,0. )
 COVAL(OB1 ,V2 ,0. ,0. )
 COVAL(OB1 ,W1 ,0. ,0. )
 COVAL(OB1 ,W2 ,0. ,0. )
 COVAL(OB1 ,R1 ,1. ,0. )
 COVAL(OB1 ,R2 ,1.0E+04 ,0. )
 COVAL(OB1 ,KE ,0. , SAME )
 COVAL(OB1 ,EP ,0. , SAME )
 COVAL(OB1 ,H1 ,0. , SAME )
 COVAL(OB1 ,H2 ,0. , SAME )
 
 PATCH(^OB2 ,SOUTH , 2, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OB2 ,U1 ,0. ,0. )
 COVAL(^OB2 ,U2 ,0. ,0. )
 COVAL(^OB2 ,V1 ,0. ,10. )
 COVAL(^OB2 ,V2 ,0. ,10. )
 COVAL(^OB2 ,W1 ,0. ,0. )
 COVAL(^OB2 ,W2 ,0. ,0. )
 COVAL(^OB2 ,R1 , FIXFLU ,6.93 )
 COVAL(^OB2 ,R2 , FIXFLU ,0.7 )
 COVAL(^OB2 ,KE ,0. ,0.25 )
 COVAL(^OB2 ,EP ,0. ,1.437814 )
 COVAL(^OB2 ,H1 ,0. ,5.25E+05 )
 COVAL(^OB2 ,H2 ,0. ,9.0E+05 )
 COVAL(^OB2 ,YO2 ,0. ,0.35 )
 COVAL(^OB2 ,COL2,0. ,0.8 )
 
 PATCH(^OB3 ,SOUTH , 3, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OB3 ,U1 ,0. ,0. )
 COVAL(^OB3 ,U2 ,0. ,0. )
 COVAL(^OB3 ,V1 ,0. ,10. )
 COVAL(^OB3 ,V2 ,0. ,10. )
 COVAL(^OB3 ,W1 ,0. ,0. )
 COVAL(^OB3 ,W2 ,0. ,0. )
 COVAL(^OB3 ,R1 , FIXFLU ,6.93 )
 COVAL(^OB3 ,R2 , FIXFLU ,0.7 )
 COVAL(^OB3 ,KE ,0. ,0.25 )
 COVAL(^OB3 ,EP ,0. ,1.437814 )
 COVAL(^OB3 ,H1 ,0. ,5.25E+05 )
 COVAL(^OB3 ,H2 ,0. ,9.0E+05 )
 COVAL(^OB3 ,YO2 ,0. ,0.35 )
 COVAL(^OB3 ,COL2,0. ,0.8 )
 
 PATCH(^OB4 ,SOUTH , 4, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OB4 ,U1 ,0. ,0. )
 COVAL(^OB4 ,U2 ,0. ,0. )
 COVAL(^OB4 ,V1 ,0. ,10. )
 COVAL(^OB4 ,V2 ,0. ,10. )
 COVAL(^OB4 ,W1 ,0. ,0. )
 COVAL(^OB4 ,W2 ,0. ,0. )
 COVAL(^OB4 ,R1 , FIXFLU ,6.93 )
 COVAL(^OB4 ,R2 , FIXFLU ,0.7 )
 COVAL(^OB4 ,KE ,0. ,0.25 )
 COVAL(^OB4 ,EP ,0. ,1.437814 )
 COVAL(^OB4 ,H1 ,0. ,5.25E+05 )
 COVAL(^OB4 ,H2 ,0. ,9.0E+05 )
 COVAL(^OB4 ,YO2 ,0. ,0.35 )
 COVAL(^OB4 ,COL2,0. ,0.8 )
 
 PATCH(^OB5 ,SOUTH , 5, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OB5 ,U1 ,0. ,0. )
 COVAL(^OB5 ,U2 ,0. ,0. )
 COVAL(^OB5 ,V1 ,0. ,10. )
 COVAL(^OB5 ,V2 ,0. ,10. )
 COVAL(^OB5 ,W1 ,0. ,0. )
 COVAL(^OB5 ,W2 ,0. ,0. )
 COVAL(^OB5 ,R1 , FIXFLU ,6.93 )
 COVAL(^OB5 ,R2 , FIXFLU ,0.7 )
 COVAL(^OB5 ,KE ,0. ,0.25 )
 COVAL(^OB5 ,EP ,0. ,1.437814 )
 COVAL(^OB5 ,H1 ,0. ,5.25E+05 )
 COVAL(^OB5 ,H2 ,0. ,9.0E+05 )
 COVAL(^OB5 ,YO2 ,0. ,0.35 )
 COVAL(^OB5 ,COL2,0. ,0.8 )
 
 PATCH(^OB6 ,SOUTH , 6, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OB6 ,U1 ,0. ,0. )
 COVAL(^OB6 ,U2 ,0. ,0. )
 COVAL(^OB6 ,V1 ,0. ,10. )
 COVAL(^OB6 ,V2 ,0. ,10. )
 COVAL(^OB6 ,W1 ,0. ,0. )
 COVAL(^OB6 ,W2 ,0. ,0. )
 COVAL(^OB6 ,R1 , FIXFLU ,6.93 )
 COVAL(^OB6 ,R2 , FIXFLU ,0.7 )
 COVAL(^OB6 ,KE ,0. ,0.25 )
 COVAL(^OB6 ,EP ,0. ,1.437814 )
 COVAL(^OB6 ,H1 ,0. ,5.25E+05 )
 COVAL(^OB6 ,H2 ,0. ,9.0E+05 )
 COVAL(^OB6 ,YO2 ,0. ,0.35 )
 COVAL(^OB6 ,COL2,0. ,0.8 )
 
 PATCH(^OB7 ,SOUTH , 7, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OB7 ,U1 ,0. ,0. )
 COVAL(^OB7 ,U2 ,0. ,0. )
 COVAL(^OB7 ,V1 ,0. ,10. )
 COVAL(^OB7 ,V2 ,0. ,10. )
 COVAL(^OB7 ,W1 ,0. ,0. )
 COVAL(^OB7 ,W2 ,0. ,0. )
 COVAL(^OB7 ,R1 , FIXFLU ,6.93 )
 COVAL(^OB7 ,R2 , FIXFLU ,0.7 )
 COVAL(^OB7 ,KE ,0. ,0.25 )
 COVAL(^OB7 ,EP ,0. ,1.437814 )
 COVAL(^OB7 ,H1 ,0. ,5.25E+05 )
 COVAL(^OB7 ,H2 ,0. ,9.0E+05 )
 COVAL(^OB7 ,YO2 ,0. ,0.35 )
 COVAL(^OB7 ,COL2,0. ,0.8 )
 
 PATCH(^OB8 ,SOUTH , 8, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OB8 ,U1 ,0. ,0. )
 COVAL(^OB8 ,U2 ,0. ,0. )
 COVAL(^OB8 ,V1 ,0. ,10. )
 COVAL(^OB8 ,V2 ,0. ,10. )
 COVAL(^OB8 ,W1 ,0. ,0. )
 COVAL(^OB8 ,W2 ,0. ,0. )
 COVAL(^OB8 ,R1 , FIXFLU ,6.93 )
 COVAL(^OB8 ,R2 , FIXFLU ,0.7 )
 COVAL(^OB8 ,KE ,0. ,0.25 )
 COVAL(^OB8 ,EP ,0. ,1.437814 )
 COVAL(^OB8 ,H1 ,0. ,5.25E+05 )
 COVAL(^OB8 ,H2 ,0. ,9.0E+05 )
 COVAL(^OB8 ,YO2 ,0. ,0.35 )
 COVAL(^OB8 ,COL2,0. ,0.8 )
 
 PATCH(^OB9 ,SOUTH , 9, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OB9 ,U1 ,0. ,0. )
 COVAL(^OB9 ,U2 ,0. ,0. )
 COVAL(^OB9 ,V1 ,0. ,10. )
 COVAL(^OB9 ,V2 ,0. ,10. )
 COVAL(^OB9 ,W1 ,0. ,0. )
 COVAL(^OB9 ,W2 ,0. ,0. )
 COVAL(^OB9 ,R1 , FIXFLU ,6.93 )
 COVAL(^OB9 ,R2 , FIXFLU ,0.7 )
 COVAL(^OB9 ,KE ,0. ,0.25 )
 COVAL(^OB9 ,EP ,0. ,1.437814 )
 COVAL(^OB9 ,H1 ,0. ,5.25E+05 )
 COVAL(^OB9 ,H2 ,0. ,9.0E+05 )
 COVAL(^OB9 ,YO2 ,0. ,0.35 )
 COVAL(^OB9 ,COL2,0. ,0.8 )
 
 PATCH(^OBA ,SOUTH , 10, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OBA ,U1 ,0. ,0. )
 COVAL(^OBA ,U2 ,0. ,0. )
 COVAL(^OBA ,V1 ,0. ,10. )
 COVAL(^OBA ,V2 ,0. ,10. )
 COVAL(^OBA ,W1 ,0. ,0. )
 COVAL(^OBA ,W2 ,0. ,0. )
 COVAL(^OBA ,R1 , FIXFLU ,6.93 )
 COVAL(^OBA ,R2 , FIXFLU ,0.7 )
 COVAL(^OBA ,KE ,0. ,0.25 )
 COVAL(^OBA ,EP ,0. ,1.437814 )
 COVAL(^OBA ,H1 ,0. ,5.25E+05 )
 COVAL(^OBA ,H2 ,0. ,9.0E+05 )
 COVAL(^OBA ,YO2 ,0. ,0.35 )
 COVAL(^OBA ,COL2,0. ,0.8 )
 
 PATCH(^OBB ,SOUTH , 11, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OBB ,U1 ,0. ,0. )
 COVAL(^OBB ,U2 ,0. ,0. )
 COVAL(^OBB ,V1 ,0. ,10. )
 COVAL(^OBB ,V2 ,0. ,10. )
 COVAL(^OBB ,W1 ,0. ,0. )
 COVAL(^OBB ,W2 ,0. ,0. )
 COVAL(^OBB ,R1 , FIXFLU ,6.93 )
 COVAL(^OBB ,R2 , FIXFLU ,0.7 )
 COVAL(^OBB ,KE ,0. ,0.25 )
 COVAL(^OBB ,EP ,0. ,1.437814 )
 COVAL(^OBB ,H1 ,0. ,5.25E+05 )
 COVAL(^OBB ,H2 ,0. ,9.0E+05 )
 COVAL(^OBB ,YO2 ,0. ,0.35 )
 COVAL(^OBB ,COL2,0. ,0.8 )
 
 PATCH(^OBC ,SOUTH , 12, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OBC ,U1 ,0. ,0. )
 COVAL(^OBC ,U2 ,0. ,0. )
 COVAL(^OBC ,V1 ,0. ,10. )
 COVAL(^OBC ,V2 ,0. ,10. )
 COVAL(^OBC ,W1 ,0. ,0. )
 COVAL(^OBC ,W2 ,0. ,0. )
 COVAL(^OBC ,R1 , FIXFLU ,6.93 )
 COVAL(^OBC ,R2 , FIXFLU ,0.7 )
 COVAL(^OBC ,KE ,0. ,0.25 )
 COVAL(^OBC ,EP ,0. ,1.437814 )
 COVAL(^OBC ,H1 ,0. ,5.25E+05 )
 COVAL(^OBC ,H2 ,0. ,9.0E+05 )
 COVAL(^OBC ,YO2 ,0. ,0.35 )
 COVAL(^OBC ,COL2,0. ,0.8 )
 
 PATCH(^OBD ,SOUTH , 13, 0, 0, 0, 0, 0, 1, 1)
 COVAL(^OBD ,U1 ,0. ,0. )
 COVAL(^OBD ,U2 ,0. ,0. )
 COVAL(^OBD ,V1 ,0. ,10. )
 COVAL(^OBD ,V2 ,0. ,10. )
 COVAL(^OBD ,W1 ,0. ,0. )
 COVAL(^OBD ,W2 ,0. ,0. )
 COVAL(^OBD ,R1 , FIXFLU ,6.93 )
 COVAL(^OBD ,R2 , FIXFLU ,0.7 )
 COVAL(^OBD ,KE ,0. ,0.25 )
 COVAL(^OBD ,EP ,0. ,1.437814 )
 COVAL(^OBD ,H1 ,0. ,5.25E+05 )
 COVAL(^OBD ,H2 ,0. ,9.0E+05 )
 COVAL(^OBD ,YO2 ,0. ,0.35 )
 COVAL(^OBD ,COL2,0. ,0.8 )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 1250 ;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(U1)=10
 LITER(U2)=10 ;LITER(V1)=10
 LITER(V2)=10 ;LITER(W1)=10
 LITER(W2)=10 ;LITER(R1)=1
 LITER(R2)=1 ;LITER(RS)=1
 LITER(KE)=20 ;LITER(EP)=20
 LITER(H1)=20 ;LITER(H2)=20
 LITER(CRAD)=20 ;LITER(YH2O)=20
 LITER(YCO2)=20 ;LITER(YCO)=20
 LITER(YO2)=20 ;LITER(YCH4)=20
 LITER(COL2)=20 ;LITER(CHA2)=20
 LITER(YHCN)=20 ;LITER(YNOX)=20
 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(R1)=1.0E-03
 ENDIT(R2)=1.0E-03 ;ENDIT(RS)=1.0E-03
 ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03
 ENDIT(H1)=1.0E-03 ;ENDIT(H2)=1.0E-03
 ENDIT(CRAD)=1.0E-03 ;ENDIT(YH2O)=1.0E-03
 ENDIT(YCO2)=1.0E-03 ;ENDIT(YCO)=1.0E-03
 ENDIT(YO2)=1.0E-03 ;ENDIT(YCH4)=1.0E-03
 ENDIT(COL2)=1.0E-03 ;ENDIT(CHA2)=1.0E-03
 ENDIT(YHCN)=1.0E-03 ;ENDIT(YNOX)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.5)
 RELAX(U1,FALSDT,1.0E-02)
 RELAX(U2,FALSDT,1.0E-02)
 RELAX(V1,FALSDT,1.0E-02)
 RELAX(V2,FALSDT,1.0E-02)
 RELAX(W1,FALSDT,1.0E-02)
 RELAX(W2,FALSDT,1.0E-02)
 RELAX(R1,LINRLX,0.5)
 RELAX(R2,LINRLX,0.5)
 RELAX(RS,LINRLX,0.4)
 RELAX(KE,LINRLX,0.2)
 RELAX(EP,LINRLX,0.2)
 RELAX(H1,FALSDT,1.0E-02)
 RELAX(H2,FALSDT,1.0E-02)
 RELAX(CRAD,FALSDT,1.0E-02)
 RELAX(MCCO,LINRLX,1.)
 RELAX(KT,LINRLX,1.)
 RELAX(XXX,LINRLX,1.)
 RELAX(BRN2,LINRLX,1.)
 RELAX(BRN1,LINRLX,1.)
 RELAX(RNOX,LINRLX,1.)
 RELAX(MOX,LINRLX,1.)
 RELAX(FK1,LINRLX,1.)
 RELAX(CCOK,LINRLX,1.)
 RELAX(CCOD,LINRLX,1.)
 RELAX(VCOD,LINRLX,1.)
 RELAX(BURN,LINRLX,1.)
 RELAX(WMIX,LINRLX,1.)
 RELAX(CCOM,LINRLX,1.)
 RELAX(VCOM,LINRLX,1.)
 RELAX(YN2,LINRLX,1.)
 RELAX(ASH2,LINRLX,1.)
 RELAX(VCOK,LINRLX,1.)
 RELAX(YH2O,FALSDT,1.0E-02)
 RELAX(YCO2,FALSDT,1.0E-02)
 RELAX(YCO,FALSDT,1.0E-02)
 RELAX(YO2,FALSDT,1.0E-02)
 RELAX(YCH4,FALSDT,1.0E-02)
 RELAX(HCOF,LINRLX,1.)
 RELAX(KC,LINRLX,1.)
 RELAX(KD,LINRLX,1.)
 RELAX(SIZE,LINRLX,1.)
 RELAX(DEN1,LINRLX,0.1)
 RELAX(APRJ,LINRLX,1.)
 RELAX(COL2,FALSDT,1.0E-02)
 RELAX(CHA2,FALSDT,1.0E-02)
 RELAX(REYN,LINRLX,1.)
 RELAX(CD,LINRLX,1.)
 RELAX(VREL,LINRLX,1.)
 RELAX(NUSS,LINRLX,1.)
 RELAX(TMP2,LINRLX,1.)
 RELAX(TMP1,LINRLX,1.)
 RELAX(CFIP,LINRLX,1.)
 RELAX(VRAT,LINRLX,1.)
 RELAX(CMDO,LINRLX,0.3)
 RELAX(EL1,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(YHCN,FALSDT,1.0E-02)
 RELAX(YNOX,FALSDT,1.0E-02)
 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(R1)=1. ;VARMIN(R1)=1.0E-06
 VARMAX(R2)=1. ;VARMIN(R2)=1.0E-09
 VARMAX(RS)=1. ;VARMIN(RS)=1.0E-06
 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)=0.
 VARMAX(H2)=1.0E+10 ;VARMIN(H2)=0.
 VARMAX(CRAD)=1.0E+10 ;VARMIN(CRAD)=-1.0E+10
 VARMAX(MCCO)=1.0E+10 ;VARMIN(MCCO)=-1.0E+10
 VARMAX(KT)=1.0E+10 ;VARMIN(KT)=-1.0E+10
 VARMAX(XXX)=1.0E+10 ;VARMIN(XXX)=-1.0E+10
 VARMAX(BRN2)=1.0E+10 ;VARMIN(BRN2)=-1.0E+10
 VARMAX(BRN1)=1.0E+10 ;VARMIN(BRN1)=-1.0E+10
 VARMAX(RNOX)=1.0E+10 ;VARMIN(RNOX)=-1.0E+10
 VARMAX(MOX)=1.0E+10 ;VARMIN(MOX)=-1.0E+10
 VARMAX(FK1)=1.0E+10 ;VARMIN(FK1)=-1.0E+10
 VARMAX(CCOK)=1.0E+10 ;VARMIN(CCOK)=-1.0E+10
 VARMAX(CCOD)=1.0E+10 ;VARMIN(CCOD)=-1.0E+10
 VARMAX(VCOD)=1.0E+10 ;VARMIN(VCOD)=-1.0E+10
 VARMAX(BURN)=1.0E+10 ;VARMIN(BURN)=-1.0E+10
 VARMAX(WMIX)=1.0E+10 ;VARMIN(WMIX)=-1.0E+10
 VARMAX(CCOM)=1.0E+10 ;VARMIN(CCOM)=0.
 VARMAX(VCOM)=1.0E+10 ;VARMIN(VCOM)=0.
 VARMAX(YN2)=1. ;VARMIN(YN2)=0.
 VARMAX(ASH2)=1. ;VARMIN(ASH2)=0.
 VARMAX(VCOK)=1.0E+10 ;VARMIN(VCOK)=-1.0E+10
 VARMAX(YH2O)=1. ;VARMIN(YH2O)=0.
 VARMAX(YCO2)=1. ;VARMIN(YCO2)=0.
 VARMAX(YCO)=1. ;VARMIN(YCO)=0.
 VARMAX(YO2)=1. ;VARMIN(YO2)=0.
 VARMAX(YCH4)=1. ;VARMIN(YCH4)=0.
 VARMAX(HCOF)=1.0E+10 ;VARMIN(HCOF)=-1.0E+10
 VARMAX(KC)=1.0E+10 ;VARMIN(KC)=-1.0E+10
 VARMAX(KD)=1.0E+10 ;VARMIN(KD)=-1.0E+10
 VARMAX(SIZE)=1.0E+10 ;VARMIN(SIZE)=-1.0E+10
 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
 VARMAX(APRJ)=1.0E+10 ;VARMIN(APRJ)=-1.0E+10
 VARMAX(COL2)=1. ;VARMIN(COL2)=0.
 VARMAX(CHA2)=1. ;VARMIN(CHA2)=0.
 VARMAX(REYN)=1.0E+10 ;VARMIN(REYN)=-1.0E+10
 VARMAX(CD)=1.0E+10 ;VARMIN(CD)=-1.0E+10
 VARMAX(VREL)=1.0E+10 ;VARMIN(VREL)=-1.0E+10
 VARMAX(NUSS)=1.0E+10 ;VARMIN(NUSS)=-1.0E+10
 VARMAX(TMP2)=1.0E+10 ;VARMIN(TMP2)=300.
 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=300.
 VARMAX(CFIP)=1.0E+10 ;VARMIN(CFIP)=-1.0E+10
 VARMAX(VRAT)=1.0E+10 ;VARMIN(VRAT)=-1.0E+10
 VARMAX(CMDO)=1.0E+10 ;VARMIN(CMDO)=0.
 VARMAX(EL1)=1.0E+10 ;VARMIN(EL1)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(YHCN)=1. ;VARMIN(YHCN)=0.
 VARMAX(YNOX)=1. ;VARMIN(YNOX)=0.
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMSAT =MOSG
 GENK = T
 ISG62 = 1
 SPEDAT(SET,OBJNAM,!OB1,C,OUTLET)
 SPEDAT(SET,OBJTYP,!OB1,C,OUTLET)
 SPEDAT(SET,ARATIO,!OB1,R,1.)
 SPEDAT(SET,OBJNAM,^OB2,C,B4)
 SPEDAT(SET,OBJTYP,^OB2,C,INLET)
 SPEDAT(SET,OBJNAM,^OB3,C,B6)
 SPEDAT(SET,OBJTYP,^OB3,C,INLET)
 SPEDAT(SET,OBJNAM,^OB4,C,B7)
 SPEDAT(SET,OBJTYP,^OB4,C,INLET)
 SPEDAT(SET,OBJNAM,^OB5,C,B9)
 SPEDAT(SET,OBJTYP,^OB5,C,INLET)
 SPEDAT(SET,OBJNAM,^OB6,C,B8)
 SPEDAT(SET,OBJTYP,^OB6,C,INLET)
 SPEDAT(SET,OBJNAM,^OB7,C,B10)
 SPEDAT(SET,OBJTYP,^OB7,C,INLET)
 SPEDAT(SET,OBJNAM,^OB8,C,B11)
 SPEDAT(SET,OBJTYP,^OB8,C,INLET)
 SPEDAT(SET,OBJNAM,^OB9,C,B12)
 SPEDAT(SET,OBJTYP,^OB9,C,INLET)
 SPEDAT(SET,OBJNAM,^OBA,C,B13)
 SPEDAT(SET,OBJTYP,^OBA,C,INLET)
 SPEDAT(SET,OBJNAM,^OBB,C,B14)
 SPEDAT(SET,OBJTYP,^OBB,C,INLET)
 SPEDAT(SET,OBJNAM,^OBC,C,B15)
 SPEDAT(SET,OBJTYP,^OBC,C,INLET)
 SPEDAT(SET,OBJNAM,^OBD,C,B16)
 SPEDAT(SET,OBJTYP,^OBD,C,INLET)
 SPEDAT(SET,FACETDAT,NUMOBJ,I,13)
 ************************************************************
  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(U1,Y,Y,Y,N,Y,Y)
 OUTPUT(U2,Y,Y,Y,N,Y,Y)
 OUTPUT(V1,Y,Y,Y,N,Y,Y)
 OUTPUT(V2,Y,Y,Y,N,Y,Y)
 OUTPUT(W1,Y,Y,Y,Y,Y,Y)
 OUTPUT(W2,Y,Y,Y,Y,Y,Y)
 OUTPUT(R1,Y,Y,Y,Y,Y,Y)
 OUTPUT(R2,Y,Y,Y,Y,Y,Y)
 OUTPUT(RS,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(CRAD,Y,N,Y,Y,Y,Y)
 OUTPUT(MCCO,Y,N,Y,N,N,N)
 OUTPUT(KT,Y,N,Y,N,N,N)
 OUTPUT(XXX,Y,N,Y,N,N,N)
 OUTPUT(BRN2,Y,N,Y,N,N,N)
 OUTPUT(BRN1,Y,N,Y,N,N,N)
 OUTPUT(RNOX,Y,N,Y,N,N,N)
 OUTPUT(MOX,Y,N,Y,N,N,N)
 OUTPUT(FK1,Y,N,Y,N,N,N)
 OUTPUT(CCOK,Y,N,Y,N,N,N)
 OUTPUT(CCOD,Y,N,Y,N,N,N)
 OUTPUT(VCOD,Y,N,Y,N,N,N)
 OUTPUT(BURN,Y,N,Y,N,N,N)
 OUTPUT(WMIX,Y,N,Y,N,N,N)
 OUTPUT(CCOM,Y,N,Y,N,N,N)
 OUTPUT(VCOM,Y,N,Y,N,N,N)
 OUTPUT(YN2,Y,N,Y,N,N,N)
 OUTPUT(ASH2,Y,N,Y,N,N,N)
 OUTPUT(VCOK,Y,N,Y,N,N,N)
 OUTPUT(YH2O,Y,N,Y,Y,Y,Y)
 OUTPUT(YCO2,Y,N,Y,Y,Y,Y)
 OUTPUT(YCO,Y,N,Y,Y,Y,Y)
 OUTPUT(YO2,Y,N,Y,Y,Y,Y)
 OUTPUT(YCH4,Y,N,Y,Y,Y,Y)
 OUTPUT(HCOF,Y,N,Y,N,N,N)
 OUTPUT(KC,Y,N,Y,N,N,N)
 OUTPUT(KD,Y,N,Y,N,N,N)
 OUTPUT(SIZE,Y,N,Y,N,N,N)
 OUTPUT(DEN1,Y,N,Y,N,N,N)
 OUTPUT(APRJ,Y,N,Y,N,N,N)
 OUTPUT(COL2,Y,N,Y,Y,Y,Y)
 OUTPUT(CHA2,Y,N,Y,Y,Y,Y)
 OUTPUT(REYN,Y,N,Y,N,N,N)
 OUTPUT(CD,Y,N,Y,N,N,N)
 OUTPUT(VREL,Y,N,Y,N,N,N)
 OUTPUT(NUSS,Y,N,Y,N,N,N)
 OUTPUT(TMP2,Y,N,Y,N,N,N)
 OUTPUT(TMP1,Y,N,Y,N,N,N)
 OUTPUT(CFIP,Y,N,Y,N,N,N)
 OUTPUT(VRAT,Y,N,Y,N,N,N)
 OUTPUT(CMDO,Y,N,Y,N,N,N)
 OUTPUT(EL1,Y,N,Y,N,N,N)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(YHCN,Y,N,Y,Y,Y,Y)
 OUTPUT(YNOX,Y,N,Y,Y,Y,Y)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 3 ;IYMON = 4 ;IZMON = 9
 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 = 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