TALK=T;RUN( 1, 1)

 ************************************************************
   Q1 created by VDI menu, Version 3.4, Date 04/09/01
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  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

  
   RADIATING PARTICLES IN A FLAME SYSTEM

  Wall-fired furnace: 24 burners

  Dry coal, no water in solid phase .

  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.

  ENDDIS
 ************************************************************
 IRUNN   =       1 ;LIBREF =       0
 ************************************************************
  Group 1. Run Title
 TEXT(Wall-fired model furnace: P-1 radiation model)
 ************************************************************
  Group 2. Transience
 STEADY  =    T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,5,5,20)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS  =    F
    * Non-default variable names
 NAME(105) =CRAD ; NAME(106) =MCCO
 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
    * Solved variables list
 SOLVE(P1  ,U1  ,U2  ,V1  ,V2  ,W1  ,W2  ,R1  )
 SOLVE(R2  ,RS  ,H1  ,H2  ,CRAD,YH2O,YCO2,YCO )
 SOLVE(YO2 ,YCH4,COL2,CHA2,YHCN,YNOX)
    * Stored variables list
 STORE(EPKE,ENUT,EL1 ,CMDO,VRAT,CFIP,TMP1,TMP2)
 STORE(NUSS,VREL,CD  ,REYN,APRJ,DEN1,SIZE,KD  )
 STORE(KC  ,HCOF,VCOK,ASH2,YN2 ,VCOM,CCOM,WMIX)
 STORE(BURN,VCOD,CCOD,CCOK,FK1 ,MOX ,RNOX,BRN1)
 STORE(BRN2,XXX ,MCCO)
    * Additional solver options
 SOLUTN(P1  ,Y,Y,Y,N,N,Y)
 SOLUTN(R1  ,Y,Y,N,Y,N,Y)
 SOLUTN(R2  ,Y,Y,N,Y,N,Y)
 TURMOD(KEMODL)

 ************************************************************
  Group 8. Terms & Devices
 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)
 NEWRH1  =    T
 ************************************************************
  Group 9. Properties
 RHO1    = GRND5
 RHO1A   = 0.000000E+00 ;RHO1B  = 3.496503E-03
 RHO1C   = 7.142857E-01
 PRESS0  = 1.000000E+05
 TMP1    = GRND
 TMP1A   = 5.000000E-01 ;TMP1B  = 1.300000E-01
 TMP1C   = 0.000000E+00
 CP1     = 1.050000E+03
 RHO2    = 1.000000E+03
 RHO2A   = 3.000000E-01 ;RHO2B  = 0.000000E+00
 RHO2C   = 0.000000E+00
 TMP2    = GRND
 TMP2A   = 5.000000E-01 ;TMP2B  = 5.000000E-01
 TMP2C   = 0.000000E+00
 CP2     = 1.800000E+03
 ENUL    = 1.000000E-05
 DRH1DP  = GRND5
 PRNDTL(R1  ) =  1.000000E+09 ;PRNDTL(R2  ) =  1.000000E+09
 PRNDTL(CRAD) =  GRND
 PRT   (R1  ) =  1.000000E+09 ;PRT   (R2  ) =  1.000000E+09
 PRT   (EP  ) =  1.314000E+00 ;PRT   (CRAD) =  1.000000E+09
 EMISS   = 0.000000E+00 ;SCATT  = 0.000000E+00
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 CFIPS   = GRND7
 RLOLIM  = 1.000000E-09 ;CMDOT  = GRND
 CFIPA   = 1.000000E-09 ;CFIPB  = 2.500000E-05
 CFIPC   = 1.000000E+02 ;CFIPD  = 0.000000E+00
 CINT  (H1  ) =  1.000000E-20 ;CINT  (H2  ) =  1.000000E-20
 ************************************************************
  Echo PLANT settings
  PLANTBEGIN

    TMP1A is gas absorption coefficient, 1/m
    TMP1B is gas scattering coefficient, 1/m
    TMP2A is particle emmisivity, dimensionless
    TMP2B is particle reflectivity, 1-TMP2A

NAMSAT=MOSG
PRNDTL(CRAD)=GRND
  LAMPR(CRAD)=1e-5*DEN1*3*(TMP1A+TMP1B+APRJ*TMP2A*(2-TMP2A))
    * 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=TMP1A*(CRAD-4.*5.68e-8*TMP1**4)+$
          APRJ*TMP2A*(CRAD-4.*5.68e-8*TMP2**4)
COVAL(RADSOR,H1,FIXFLU,GRND)
   VAL=TMP1A*4*APRJ*(CRAD/4.-5.68e-8*TMP2**4)
COVAL(RADSOR,H2,FIXFLU,GRND)

    * Radiosity sources
      * from gas
PATCH(RDSTGAS,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
   VAL=4.*5.68e-8*TMP1**4
COVAL(RDSTGAS,CRAD,TMP1A,GRND)
      * from fines
PATCH(RDSTFINE,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)
   CO =APRJ*TMP2A
   VAL=4.*5.68e-8*TMP2**4
COVAL(RDSTFINE,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/               $
              EXP(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 11.Initialise Var/Porosity Fields
 FIINIT(R1  ) =  5.000000E-01 ;FIINIT(R2  ) =  5.000000E-01
 FIINIT(RS  ) =  5.000000E-01
   No PATCHes used for this Group


 INIADD  =    F
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =    1250
 RESFAC  = 1.000000E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER (R1  ) =    1 ;LITER (R2  ) =    1
 LITER (RS  ) =    1
 ************************************************************
  Group 17. Relaxation
 RELAX(P1  ,LINRLX, 5.000000E-01)
 RELAX(U1  ,FALSDT, 1.000000E-02)
 RELAX(U2  ,FALSDT, 1.000000E-02)
 RELAX(V1  ,FALSDT, 1.000000E-02)
 RELAX(V2  ,FALSDT, 1.000000E-02)
 RELAX(W1  ,FALSDT, 1.000000E-02)
 RELAX(W2  ,FALSDT, 1.000000E-02)
 RELAX(R1  ,LINRLX, 5.000000E-01)
 RELAX(R2  ,LINRLX, 5.000000E-01)
 RELAX(KE  ,LINRLX, 2.000000E-01)
 RELAX(EP  ,LINRLX, 2.000000E-01)
 RELAX(H1  ,FALSDT, 1.000000E-02)
 RELAX(H2  ,FALSDT, 1.000000E-02)
 RELAX(CRAD,FALSDT, 1.000000E-02)
 RELAX(YH2O,FALSDT, 1.000000E-02)
 RELAX(YCO2,FALSDT, 1.000000E-02)
 RELAX(YCO ,FALSDT, 1.000000E-02)
 RELAX(YO2 ,FALSDT, 1.000000E-02)
 RELAX(YCH4,FALSDT, 1.000000E-02)
 RELAX(DEN1,LINRLX, 1.000000E-01)
 RELAX(COL2,FALSDT, 1.000000E-02)
 RELAX(CHA2,FALSDT, 1.000000E-02)
 RELAX(CMDO,LINRLX, 3.000000E-01)
 RELAX(YHCN,FALSDT, 1.000000E-02)
 RELAX(YNOX,FALSDT, 1.000000E-02)
 KELIN   =       3
 ************************************************************
  Group 18. Limits
 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06
 VARMAX(U2  ) = 1.000000E+06 ;VARMIN(U2  ) =-1.000000E+06
 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06
 VARMAX(V2  ) = 1.000000E+06 ;VARMIN(V2  ) =-1.000000E+06
 VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06
 VARMAX(W2  ) = 1.000000E+06 ;VARMIN(W2  ) =-1.000000E+06
 VARMAX(R1  ) = 1.000000E+00 ;VARMIN(R1  ) = 1.000000E-06
 VARMAX(R2  ) = 1.000000E+00 ;VARMIN(R2  ) = 1.000000E-09
 VARMAX(H1  ) = 1.000000E+10 ;VARMIN(H1  ) = 0.000000E+00
 VARMAX(H2  ) = 1.000000E+10 ;VARMIN(H2  ) = 0.000000E+00
 VARMAX(CCOM) = 1.000000E+10 ;VARMIN(CCOM) = 0.000000E+00
 VARMAX(VCOM) = 1.000000E+10 ;VARMIN(VCOM) = 0.000000E+00
 VARMAX(YN2 ) = 1.000000E+00 ;VARMIN(YN2 ) = 0.000000E+00
 VARMAX(ASH2) = 1.000000E+00 ;VARMIN(ASH2) = 0.000000E+00
 VARMAX(YH2O) = 1.000000E+00 ;VARMIN(YH2O) = 0.000000E+00
 VARMAX(YCO2) = 1.000000E+00 ;VARMIN(YCO2) = 0.000000E+00
 VARMAX(YCO ) = 1.000000E+00 ;VARMIN(YCO ) = 0.000000E+00
 VARMAX(YO2 ) = 1.000000E+00 ;VARMIN(YO2 ) = 0.000000E+00
 VARMAX(YCH4) = 1.000000E+00 ;VARMIN(YCH4) = 0.000000E+00
 VARMAX(COL2) = 1.000000E+00 ;VARMIN(COL2) = 0.000000E+00
 VARMAX(CHA2) = 1.000000E+00 ;VARMIN(CHA2) = 0.000000E+00
 VARMAX(TMP2) = 1.000000E+10 ;VARMIN(TMP2) = 3.000000E+02
 VARMAX(TMP1) = 1.000000E+10 ;VARMIN(TMP1) = 3.000000E+02
 VARMAX(CMDO) = 1.000000E+10 ;VARMIN(CMDO) = 0.000000E+00
 VARMAX(YHCN) = 1.000000E+00 ;VARMIN(YHCN) = 0.000000E+00
 VARMAX(YNOX) = 1.000000E+00 ;VARMIN(YNOX) = 0.000000E+00
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 NAMSAT  =MOSG
 GENK    =    T
 ASAP    =    T
 ************************************************************
  Group 20. Preliminary Printout
 ECHO    =    T
 ************************************************************
  Group 21. Print-out of Variables
 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)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON   =       3 ;IYMON  =       4 ;IZMON  =       9
 NPRMON  =  100000
 NPRMNT  =       1
 TSTSWP  =      -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT  =  100000
 NZPRIN  =       1
 ISWPRF  =       1 ;ISWPRL =  100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 NOWIPE  =    T

 GVIEW(P,-8.292519E-01,3.366233E-01,4.461233E-01)
 GVIEW(UP,4.125013E-01,-1.699202E-01,8.949692E-01)

> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+01
> DOM,    MONIT,       5.000000E-01, 7.000000E-01, 4.250000E+00
> DOM,    SCALE,       3.000000E+00, 3.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> DOM,    RELAX,       3.000000E-01

> OBJ,    NAME,        OUTLET
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 1.000000E+01
> OBJ,    SIZE,        1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ,    CLIPART,     cube12t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        OUTLET
> OBJ,    PRESSURE,      0.000000E+00
> OBJ,    TEMPERATURE,  -1.026000E+04
> OBJ,    COEFFICIENT,   1.000000E+00
> OBJ,    TEMPERATURE-2,-1.026000E+04
> OBJ,    COEFFICIENT-2, 1.000000E+04
> OBJ,    TURBULENCE,   -1.026000E+04,-1.026000E+04

> OBJ,    NAME,        B4
> OBJ,    POSITION,    6.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        B6
> OBJ,    POSITION,    2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        B7
> OBJ,    POSITION,    2.000000E-01, 0.000000E+00, 1.500000E+00
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        B9
> OBJ,    POSITION,    6.000000E-01, 0.000000E+00, 1.500000E+00
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        B8
> OBJ,    POSITION,    6.000000E-01, 0.000000E+00, 2.500000E+00
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        B10
> OBJ,    POSITION,    2.000000E-01, 0.000000E+00, 2.500000E+00
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        B11
> OBJ,    POSITION,    6.000000E-01, 0.000000E+00, 3.500000E+00
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        B12
> OBJ,    POSITION,    2.000000E-01, 0.000000E+00, 3.500000E+00
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        B13
> OBJ,    POSITION,    6.000000E-01, 0.000000E+00, 4.500000E+00
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        B14
> OBJ,    POSITION,    2.000000E-01, 0.000000E+00, 4.500000E+00
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        B15
> OBJ,    POSITION,    6.000000E-01, 0.000000E+00, 5.500000E+00
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        B16
> OBJ,    POSITION,    2.000000E-01, 0.000000E+00, 5.500000E+00
> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    GRID,              2
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       6.930000E-01
> OBJ,    DENSITY-2,     7.000000E-02
> OBJ,    VELOCITY,      0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    VELOCITY-2,    0.000000E+00, 1.000000E+01, 0.000000E+00
> OBJ,    TEMPERATURE,   5.250000E+05
> OBJ,    TEMPERATURE-2, 9.000000E+05
> OBJ,    INLET_YO2,     3.500000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00
STOP