************************************************************
   Q1 created by VDI menu, Version 3.4, Date 04/09/01
 CPVNAM=VDI;SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  PHOTON USE
  p



  use patgeo
  con h1 z 1 fil;.001
  con h1 x 3 fil;.001
  con h1 z m fil;.001
  con h1 z 10 fil;.001
  vec x 3
  msg(Single-fluid gas enthalpy contours
  pause;cl
  use patgeo
  con hfav z 1 fil;.001
  con hfav x 3 fil;.001
  con hfav z m fil;.001
  con hfav z 10 fil;.001
  vec x 3
  msg(Fluid-averaged gas enthalpy contours
  pause;cl
  use patgeo
  con tmp1 z 1 fil;.001
  con tmp1 x 3 fil;.001
  con tmp1 z m fil;.001
  con tmp1 z 10 fil;.001
  vec x 3
  msg(Single-fliuid gas temperature contours
  pause;cl
  use patgeo
  con tfav z 1 fil;.001
  con tfav x 3 fil;.001
  con tfav z m fil;.001
  con tfav z 10 fil;.001
  vec x 3
  msg(Fluid-averaged gas temperature contours
  pause;cl
  use patgeo
  con tf1 z 1 fil;.001
  con tf1 x 3 fil;.001
  con tf1 z m fil;.001
  con tf1 z 10 fil;.001
  vec x 3
  msg(1st fluid temperature contours
  pause;cl
  use patgeo
  con tf2 z 1 fil;.001
  con tf2 x 3 fil;.001
  con tf2 z m fil;.001
  con tf2 z 10 fil;.001
  vec x 3
  msg(2nd fluid temperature contours
  pause;cl
  use patgeo
  con tf3 z 1 fil;.001
  con tf3 x 3 fil;.001
  con tf3 z m fil;.001
  con tf3 z 10 fil;.001
  vec x 3
  msg(3rd fluid temperature contours
  pause;cl
  use patgeo
  con tf4 z 1 fil;.001
  con tf4 x 3 fil;.001
  con tf4 z m fil;.001
  con tf4 z 10 fil;.001
  vec x 3
  msg(4th fluid temperature contours
  pause;cl
  use patgeo
  con tf5 z 1 fil;.001
  con tf5 x 3 fil;.001
  con tf5 z m fil;.001
  con tf5 z 10 fil;.001
  vec x 3
  msg(5th fluid temperature contours
  pause;cl

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

  use patgeo
  con YNOX z 1 fil;.001
  con ynox x 3 fil;.001
  con ynox z m fil;.001
  con ynox z 10 fil;.001
  vec x 3
  msg(Single-fluid thermal NOX mass fraction, kg/kg
  pause;cl
  use patgeo
  con noav z 1 fil;.001
  con noav x 3 fil;.001
  con noav z m fil;.001
  con noav z 10 fil;.001
  vec x 3
  msg(Fluid-averaged thermal NOX mass fraction, kg/kg
  pause;cl
  use patgeo
  con no1 z 1 fil;.001
  con no1 x 3 fil;.001
  con no1 z m fil;.001
  con no1 z 10 fil;.001
  vec x 3
  msg(1st fluid thermal NOX mass fraction, kg/kg
  pause;cl
  use patgeo
  con no2 z 1 fil;.001
  con no2 x 3 fil;.001
  con no2 z m fil;.001
  con no2 z 10 fil;.001
  vec x 3
  msg(2nd fluid thermal NOX mass fraction, kg/kg
  pause;cl
  use patgeo
  con no3 z 1 fil;.001
  con no3 x 3 fil;.001
  con no3 z m fil;.001
  con no3 z 10 fil;.001
  vec x 3
  msg(3rd fluid thermal NOX mass fraction, kg/kg
  pause;cl
  use patgeo
  con no4 z 1 fil;.001
  con no4 x 3 fil;.001
  con no4 z m fil;.001
  con no4 z 10 fil;.001
  vec x 3
  msg(4th fluid thermal NOX mass fraction, kg/kg
  pause;cl
  use patgeo
  con no5 z 1 fil;.001
  con no5 x 3 fil;.001
  con no5 z m fil;.001
  con no5 z 10 fil;.001
  vec x 3
  msg(5th fluid thermal NOX mass fraction, kg/kg
  pause;cl

  use patgeo
  con vrat z 1 fil;.001
  con vrat x 3 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 3 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 3 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 3 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 3 fil;.001
  con YCH4 z m fil;.001
  con YCH4 z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(Single-fluid volatiles, CH4, mass fractions
  pause;cl
  use patgeo
  con YCO2 z 1 fil;.001
  con YCO2 x 3 fil;.001
  con YCO2 z m fil;.001
  con YCO2 z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(Single-fluid CO2  mass fractions
  pause;cl
  use patgeo
  con YCO z 1 fil;.001
  con YCO x 3 fil;.001
  con YCO z m fil;.001
  con YCO z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(Single-fluid CO  mass fractions
  pause;cl
  use patgeo
  con YH2O z 1 fil;.001
  con YH2O x 3 fil;.001
  con YH2O z m fil;.001
  con YH2O z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(Single-fluid H2O  mass fractions
  pause;cl
  use patgeo
  con YO2 z 1 fil;.001
  con YO2 x 3 fil;.001
  con YO2 z m fil;.001
  con YO2 z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(Single-fluid O2  mass fractions
  pause;cl
  use patgeo
  con oxav z 1 fil;.001
  con oxav x 3 fil;.001
  con oxav z m fil;.001
  con oxav z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(Fluid-averaged O2  mass fractions
  pause;cl
  use patgeo
  con ox1 z 1 fil;.001
  con ox1 x 3 fil;.001
  con ox1 z m fil;.001
  con ox1 z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(1st fluid O2  mass fractions
  pause;cl
  use patgeo
  con ox2 z 1 fil;.001
  con ox2 x 3 fil;.001
  con ox2 z m fil;.001
  con ox2 z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(2nd fluid O2  mass fractions
  pause;cl
  use patgeo
  con ox3 z 1 fil;.001
  con ox3 x 3 fil;.001
  con ox3 z m fil;.001
  con ox3 z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(3rd fluid O2  mass fractions
  pause;cl
  use patgeo
  con ox4 z 1 fil;.001
  con ox4 x 3 fil;.001
  con ox4 z m fil;.001
  con ox4 z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(4th fluid O2  mass fractions
  pause;cl
  use patgeo
  con ox5 z 1 fil;.001
  con ox5 x 3 fil;.001
  con ox5 z m fil;.001
  con ox5 z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(5th fluid O2  mass fractions
  pause;cl
  use patgeo
  con YN2 z 1 fil;.001
  con YN2 x 3 fil;.001
  con YN2 z m fil;.001
  con YN2 z 10 fil;.001
  vec x 3
  msg(Gas composition:
  msg(Single-fluid N2  mass fractions

  pause;cl
  use patgeo
  con COL2 z 1 fil;.001
  con COL2 x 3 fil;.001
  con COL2 z m fil;.001
  con COL2 z 10 fil;.001
  vec x 3
  msg(Particle composition:
  msg(Raw COAL particle content
  pause;cl
  use patgeo
  con CHA2 z 1 fil;.001
  con CHA2 x 3 fil;.001
  con CHA2 z m fil;.001
  con CHA2 z 10 fil;.001
  vec x 3
  msg(Particle composition:
  msg(CHAR particle content
  pause;cl
  use patgeo
  con ASH2 z 1 fil;.001
  con ASH2 x 3 fil;.001
  con ASH2 z m fil;.001
  con ASH2 z 10 fil;.001
  vec x 3
  msg(Particle composition:
  msg(ASH particle content
  pause;cl
  use patgeo
  con R2 z 1 fil;.001
  con R2 x 3 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 3 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 3 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 3 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 3 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 3 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 3 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 3 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 3 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
  use patgeo
  con avef z 1 fil;.001
  con avef x 3 fil;.001
  con avef z m fil;.001
  con avef z 10 fil;.001
  vec x 3
  msg(Fluid-averaged mixture fraction
  pause;cl
  use patgeo
  con fmix z 1 fil;.001
  con fmix x 3 fil;.001
  con fmix z m fil;.001
  con fmix z 10 fil;.001
  vec x 3
  msg(Fluid-averaged mixture fraction
  pause;cl
  use patgeo
  con mnsq z 1 fil;.001
  con mnsq x 3 fil;.001
  con mnsq z m fil;.001
  con mnsq z 10 fil;.001
  vec x 3
  msg(RMS of  mixture fraction fluctuations
  pause;cl
  use patgeo
  con f1 z 1 fil;.001
  con f1 x 3 fil;.001
  con f1 z m fil;.001
  con f1 z 10 fil;.001
  vec x 3
  msg(1st Fluid Mass Fraction
  pause;cl
  use patgeo
  con f2 z 1 fil;.001
  con f2 x 3 fil;.001
  con f2 z m fil;.001
  con f2 z 10 fil;.001
  vec x 3
  msg(2nd Fluid Mass Fraction
  pause;cl
  use patgeo
  con f3 z 1 fil;.001
  con f3 x 3 fil;.001
  con f3 z m fil;.001
  con f3 z 10 fil;.001
  vec x 3
  msg(3rd Fluid Mass Fraction
  pause;cl
  use patgeo
  con f4 z 1 fil;.001
  con f4 x 3 fil;.001
  con f4 z m fil;.001
  con f4 z 10 fil;.001
  vec x 3
  msg(4th Fluid Mass Fraction
  pause;cl
  use patgeo
  con f5 z 1 fil;.001
  con f5 x 3 fil;.001
  con f5 z m fil;.001
  con f5 z 10 fil;.001
  vec x 3
  msg(5th Fluid Mass Fraction
  pause;cl
  ENDUSE
  DISPLAY

  IPSA and 5 fluids for enthalpy CVA

  Dry coal, no water in solid phase .

  1 Kinetically released volatiles.
  2 Combustion of

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

  ENDDIS
 ************************************************************
 IRUNN   =       1 ;LIBREF =       0
 ************************************************************
  
  Group 1. Run Title
 TEXT(Coal combustor and 5 fluids for enthalpy)
 ************************************************************
  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( 34) =VCM5 ; NAME( 35) =VCM4
 NAME( 36) =VCM3 ; NAME( 37) =VCM2
 NAME( 38) =VCM1 ; NAME( 39) =CCO5
 NAME( 40) =CCO4 ; NAME( 41) =CCO3
 NAME( 42) =CCO2 ; NAME( 43) =CCO1
 NAME( 44) =COAV ; NAME( 45) =FUAV
 NAME( 46) =VCO4 ; NAME( 47) =VCO3
 NAME( 48) =VCO2 ; NAME( 49) =VFU4
 NAME( 50) =VFU3 ; NAME( 51) =VFU2
 NAME( 52) =CO5  ; NAME( 53) =CO4
 NAME( 54) =CO3  ; NAME( 55) =CO2
 NAME( 56) =CO1  ; NAME( 57) =FU5
 NAME( 58) =FU4  ; NAME( 59) =FU3
 NAME( 60) =FU2  ; NAME( 61) =FU1
 NAME( 62) =RNO5 ; NAME( 63) =RNO4
 NAME( 64) =RNO3 ; NAME( 65) =RNO2
 NAME( 66) =RNO1 ; NAME( 67) =VNO4
 NAME( 68) =VNO3 ; NAME( 69) =VNO2
 NAME( 70) =NOAV ; NAME( 71) =NO5
 NAME( 72) =NO4  ; NAME( 73) =NO3
 NAME( 74) =NO2  ; NAME( 75) =NO1
 NAME( 76) =VOX4 ; NAME( 77) =VOX3
 NAME( 78) =VOX2 ; NAME( 79) =OXAV
 NAME( 80) =OX5  ; NAME( 81) =OX4
 NAME( 82) =OX3  ; NAME( 83) =OX2
 NAME( 84) =MNSQ ; NAME( 85) =AVEF
 NAME( 86) =FMIX ; NAME( 87) =FSUM
 NAME( 88) =OX1  ; NAME( 89) =TFAV
 NAME( 90) =HFAV ; NAME( 91) =VHF4
 NAME( 92) =VHF3 ; NAME( 93) =VHF2
 NAME( 94) =COF4 ; NAME( 95) =COF3
 NAME( 96) =COF2 ; NAME( 97) =TF5
 NAME( 98) =TF4  ; NAME( 99) =TF3
 NAME(100) =TF2  ; NAME(101) =TF1
 NAME(102) =HF5  ; NAME(103) =HF4
 NAME(104) =F5   ; NAME(105) =F4
 NAME(106) =F3   ; NAME(107) =F2
 NAME(108) =F1   ; NAME(109) =HF3
 NAME(110) =HF2  ; NAME(111) =RNOX
 NAME(112) =MOX  ; NAME(113) =FK1
 NAME(114) =YNOX ; NAME(115) =HF1
 NAME(116) =CCOK ; NAME(117) =CCOD
 NAME(118) =VCOD ; NAME(119) =BURN
 NAME(120) =WMIX ; NAME(121) =CCOM
 NAME(122) =VCOM ; NAME(123) =YN2
 NAME(124) =YH2O ; NAME(125) =YCO2
 NAME(126) =YCO  ; NAME(127) =YO2
 NAME(128) =YCH4 ; NAME(129) =ASH2
 NAME(130) =VCOK ; NAME(131) =HCOF
 NAME(132) =KC   ; NAME(133) =KD
 NAME(134) =SIZE ; NAME(135) =COL2
 NAME(136) =CHA2 ; NAME(137) =DEN1
 NAME(138) =APRJ ; NAME(139) =REYN
 NAME(140) =CD   ; NAME(141) =VREL
 NAME(142) =NUSS ; NAME(143) =TMP2
 NAME(144) =TMP1 ; NAME(145) =CFIP
 NAME(146) =VRAT ; NAME(147) =CMDO
 NAME(148) =EL1  ; NAME(149) =ENUT
 NAME(150) =EPKE
    * Solved variables list
 SOLVE(P1  ,U1  ,U2  ,V1  ,V2  ,W1  ,W2  ,R1  )
 SOLVE(R2  ,RS  ,H1  ,H2  ,CO5 ,CO4 ,CO3 ,CO2 )
 SOLVE(CO1 ,FU5 ,FU4 ,FU3 ,FU2 ,FU1 ,NO5 ,NO4 )
 SOLVE(NO3 ,NO2 ,NO1 ,OX5 ,OX4 ,OX3 ,OX2 ,OX1 )
 SOLVE(HF5 ,HF4 ,F5  ,F4  ,F3  ,F2  ,F1  ,HF3 )
 SOLVE(HF2 ,YNOX,HF1 ,YH2O,YCO2,YCO ,YO2 ,YCH4)
 SOLVE(COL2,CHA2)
    * 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,TF1 )
 STORE(TF2 ,TF3 ,TF4 ,TF5 ,COF2,COF3,COF4,VHF2)
 STORE(VHF3,VHF4,HFAV,TFAV,FSUM,FMIX,AVEF,MNSQ)
 STORE(OXAV,VOX2,VOX3,VOX4,NOAV,VNO2,VNO3,VNO4)
 STORE(RNO1,RNO2,RNO3,RNO4,RNO5,VFU2,VFU3,VFU4)
 STORE(VCO2,VCO3,VCO4,FUAV,COAV,CCO1,CCO2,CCO3)
 STORE(CCO4,CCO5)
 STORE(VCM1,VCM2,VCM3,VCM4,VCM5)
 STORE(BRN1,BRN2,BRN3,BRN4,BRN5)
    * 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 (CO5 ,N,Y,Y,Y,Y,N)
 TERMS (CO4 ,N,Y,Y,Y,Y,N)
 TERMS (CO3 ,N,Y,Y,Y,Y,N)
 TERMS (CO2 ,N,Y,Y,Y,Y,N)
 TERMS (CO1 ,N,Y,Y,Y,Y,N)
 TERMS (FU5 ,N,Y,Y,Y,Y,N)
 TERMS (FU4 ,N,Y,Y,Y,Y,N)
 TERMS (FU3 ,N,Y,Y,Y,Y,N)
 TERMS (FU2 ,N,Y,Y,Y,Y,N)
 TERMS (FU1 ,N,Y,Y,Y,Y,N)
 TERMS (NO5 ,N,Y,Y,Y,Y,N)
 TERMS (NO4 ,N,Y,Y,Y,Y,N)
 TERMS (NO3 ,N,Y,Y,Y,Y,N)
 TERMS (NO2 ,N,Y,Y,Y,Y,N)
 TERMS (NO1 ,N,Y,Y,Y,Y,N)
 TERMS (OX5 ,N,Y,Y,Y,Y,N)
 TERMS (OX4 ,N,Y,Y,Y,Y,N)
 TERMS (OX3 ,N,Y,Y,Y,Y,N)
 TERMS (OX2 ,N,Y,Y,Y,Y,N)
 TERMS (OX1 ,N,Y,Y,Y,Y,N)
 TERMS (HF5 ,N,Y,Y,Y,Y,N)
 TERMS (HF4 ,N,Y,Y,Y,Y,N)
 TERMS (F5  ,N,Y,Y,Y,Y,N)
 TERMS (F4  ,N,Y,Y,Y,Y,N)
 TERMS (F3  ,N,Y,Y,Y,Y,N)
 TERMS (F2  ,N,Y,Y,Y,Y,N)
 TERMS (F1  ,N,Y,Y,Y,Y,N)
 TERMS (HF3 ,N,Y,Y,Y,Y,N)
 TERMS (HF2 ,N,Y,Y,Y,Y,N)
 TERMS (YNOX,N,Y,Y,Y,Y,N)
 TERMS (HF1 ,N,Y,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)
 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
 CP1     = 1.050000E+03
 RHO2    = 1.000000E+03
 RHO2A   = 3.000000E-01 ;RHO2B  = 0.000000E+00
 RHO2C   = 0.000000E+00
 TMP2    = GRND
 CP2     = 1.800000E+03
 ENUL    = 1.000000E-05
 DRH1DP  = GRND5
 PRNDTL(R1  ) =  1.000000E+09 ;PRNDTL(R2  ) =  1.000000E+09
 PRT   (R1  ) =  1.000000E+09 ;PRT   (R2  ) =  1.000000E+09
 PRT   (EP  ) =  1.314000E+00
 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  = 1.000000E-04
 CFIPC   = 1.000000E+02 ;CFIPD  = 0.000000E+00
 CINT  (H1  ) =  1.000000E-20 ;CINT  (H2  ) =  1.000000E-20
 ************************************************************
  Echo PLANT settings
  
  PLANTBEGIN
NAMSAT=MOSG
    * Interphase mass transfer
       VRAT-volatile release; BURN-char burning
   INTMDT=VRAT+BURN
    * Phase temperatures
   TEMP1=H1/1300.
   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)
   VAL=VRAT*1.0
COVAL(PLSR5,F5,FIXFLU,GRND)

   Reaction:   2CH4  + 3O2 ->  2CO +4H2O
PATCH(PLSR6,PHASEM,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)  + O2 ->  CO2
    * 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.*BURN
COVAL(PLSR9,YCO2,FIXFLU,GRND)
   VAL=-32./12.*BURN
COVAL(PLSR9,YO2,FIXFLU,GRND)
   VAL=BURN*1.0
COVAL(PLSR9,F1,FIXFLU,GRND)
     ===========================================
         Carbon monoxide oxidation
   Reaction:   2CO +  O2 -> 2CO2
PATCH(PLSR7,PHASEM,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,PHASEM,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*3.28e7
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/1300.
   VAL=H2*1300./1800.
COVAL(PLSR4,H1,GRND,GRND)
   CO=HCOF/1800.
   VAL=1800./1300.*H1
COVAL(PLSR4,H2,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)
   VAL=CMDO*YO2
COVAL(PLSR90,YO2,FIXFLU,GRND)
   VAL=CMDO*YCO
COVAL(PLSR90,YCO,FIXFLU,GRND)
   VAL=CMDO*YCO2
COVAL(PLSR90,YCO2,FIXFLU,GRND)
   VAL=CMDO*YH2O
COVAL(PLSR90,YH2O,FIXFLU,GRND)
     ===========================================
PATCH(NOXSOR,VOLUME,1,NX,1,NY,1,NZ,1,1)
     VAL=RNOX
COVAL(NOXSOR,YNOX,FIXFLU,GRND)
     ===========================================
       * Micromixing contributions to CVA transport
         ** Contributions to enthalpies
PATCH(MICRHF,PHASEM,1,NX,1,NY,1,NZ,1,1)
   CO = 5.*EPKE*COF2
   VAL= VHF2
COVAL(MICRHF,HF2,GRND,GRND)
   CO = 5.*EPKE*COF3
   VAL= VHF3
COVAL(MICRHF,HF3,GRND,GRND)
   CO = 5.*EPKE*COF4
   VAL= VHF4
COVAL(MICRHF,HF4,GRND,GRND)
         ** Contributions to O2 mass fractions
PATCH(MICROX,PHASEM,1,NX,1,NY,1,NZ,1,1)
   CO = 5.*EPKE*COF2
   VAL= VOX2
COVAL(MICROX,OX2,GRND,GRND)
   CO = 5.*EPKE*COF3
   VAL= VOX3
COVAL(MICROX,OX3,GRND,GRND)
   CO = 5.*EPKE*COF4
   VAL= VOX4
COVAL(MICROX,OX4,GRND,GRND)
         ** Contributions to NO mass fractions
PATCH(MICRNO,PHASEM,1,NX,1,NY,1,NZ,1,1)
   CO = 5.*EPKE*COF2
   VAL= VNO2
COVAL(MICRNO,NO2,GRND,GRND)
   CO = 5.*EPKE*COF3
   VAL= VNO3
COVAL(MICRNO,NO3,GRND,GRND)
   CO = 5.*EPKE*COF4
   VAL= VNO4
COVAL(MICRNO,NO4,GRND,GRND)
         ** Contributions to CH4 mass fractions
PATCH(MICRFU,PHASEM,1,NX,1,NY,1,NZ,1,1)
   CO = 5.*EPKE*COF2
   VAL= VFU2
COVAL(MICRFU,FU2,GRND,GRND)
   CO = 5.*EPKE*COF3
   VAL= VFU3
COVAL(MICRFU,FU3,GRND,GRND)
   CO = 5.*EPKE*COF4
   VAL= VFU4
COVAL(MICRFU,FU4,GRND,GRND)
         ** Contributions to CO mass fractions
PATCH(MICCO,PHASEM,1,NX,1,NY,1,NZ,1,1)
   CO = 5.*EPKE*COF2
   VAL= VCO2
COVAL(MICCO,CO2,GRND,GRND)
   CO = 5.*EPKE*COF3
   VAL= VCO3
COVAL(MICCO,CO3,GRND,GRND)
   CO = 5.*EPKE*COF4
   VAL= VCO4
COVAL(MICCO,CO4,GRND,GRND)
       * Within-fluid enthalpy sources due combustion heats
PATCH(COMHEATF,CELL,1,NX,1,NY,1,NZ,1,1)
  VAL=((VCM1*5.5e7+CCO1*7.45e6)*VOL*DEN1*R1+BRN1*3.28e7)*F1
COVAL(COMHEATF,HF1,FIXFLU,GRND)
  VAL=((VCM2*5.5e7+CCO2*7.45e6)*VOL*DEN1*R1+BRN2*3.28e7)*F2
COVAL(COMHEATF,HF2,FIXFLU,GRND)
  VAL=((VCM3*5.5e7+CCO3*7.45e6)*VOL*DEN1*R1+BRN3*3.28e7)*F3
COVAL(COMHEATF,HF3,FIXFLU,GRND)
  VAL=((VCM4*5.5e7+CCO4*7.45e6)*VOL*DEN1*R1+BRN4*3.28e7)*F4
COVAL(COMHEATF,HF4,FIXFLU,GRND)
  VAL=((VCM5*5.5e7+CCO5*7.45e6)*VOL*DEN1*R1+BRN5*3.28e7)*F5
COVAL(COMHEATF,HF5,FIXFLU,GRND)
    * Within-fluid enthalpy sources due interphase heat transfer
PATCH(INTHEATF,VOLUME,1,NX,1,NY,1,NZ,1,1)
   CO=HCOF/1300.*F1
   VAL=H2*1300./1800.
COVAL(INTHEATF,HF1,GRND,GRND)
   CO=HCOF/1300.*F2
   VAL=H2*1300./1800.
COVAL(INTHEATF,HF2,GRND,GRND)
   CO=HCOF/1300.*F3
   VAL=H2*1300./1800.
COVAL(INTHEATF,HF3,GRND,GRND)
   CO=HCOF/1300.*F4
   VAL=H2*1300./1800.
COVAL(INTHEATF,HF4,GRND,GRND)
   CO=HCOF/1300.*F5
   VAL=H2*1300./1800.
COVAL(INTHEATF,HF5,GRND,GRND)
       * Within-fluid O2 sources
PATCH(OXYGSOUR,CELL,1,NX,1,NY,1,NZ,1,1)
  CO=((.57*CCO1+3.*VCM1)*VOL*DEN1*R1+2.67*BRN1)/           $
   (OX1+tiny)*F1
COVAL(OXYGSOUR,OX1,GRND,0.0)
  CO=((.57*CCO2+3.*VCM2)*VOL*DEN1*R1+2.67*BRN2)/           $
   (OX2+tiny)*F2
COVAL(OXYGSOUR,OX2,GRND,0.0)
  CO=((.57*CCO3+3.*VCM3)*VOL*DEN1*R1+2.67*BRN3)/           $
   (OX3+tiny)*F3
COVAL(OXYGSOUR,OX3,GRND,0.0)
  CO=((.57*CCO4+3.*VCM4)*VOL*DEN1*R1+2.67*BRN4)/           $
   (OX4+tiny)*F4
COVAL(OXYGSOUR,OX4,GRND,0.0)
  CO=((.57*CCO5+3.*VCM5)*VOL*DEN1*R1+2.67*BRN5)/           $
   (OX5+tiny)*F5
COVAL(OXYGSOUR,OX5,GRND,0.0)
       * Within-fluid NO sources
PATCH(NOSOUR,VOLUME,1,NX,1,NY,1,NZ,1,1)
  VAL=RNO1*F1
COVAL(NOSOUR,NO1,FIXFLU,GRND)
  VAL=RNO2*F2
COVAL(NOSOUR,NO2,FIXFLU,GRND)
  VAL=RNO3*F3
COVAL(NOSOUR,NO3,FIXFLU,GRND)
  VAL=RNO4*F4
COVAL(NOSOUR,NO4,FIXFLU,GRND)
  VAL=RNO5*F5
COVAL(NOSOUR,NO5,FIXFLU,GRND)
       * Within-fluid CH4 sources
PATCH(CH4SOUR,CELL,1,NX,1,NY,1,NZ,1,1)
  VAL=VRAT*F1
COVAL(CH4SOUR,FU1,FIXFLU,GRND)
  VAL=VRAT*F2
COVAL(CH4SOUR,FU2,FIXFLU,GRND)
  VAL=VRAT*F3
COVAL(CH4SOUR,FU3,FIXFLU,GRND)
  VAL=VRAT*F4
COVAL(CH4SOUR,FU4,FIXFLU,GRND)
  VAL=VRAT*F5
COVAL(CH4SOUR,FU5,FIXFLU,GRND)

PATCH(CH4SINK,CELL,1,NX,1,NY,1,NZ,1,1)
  CO=VCM1/(FU1+tiny)*F1
COVAL(CH4SINK,FU1,GRND,0.0)
  CO=VCM2/(FU2+tiny)*F2
COVAL(CH4SINK,FU2,GRND,0.0)
  CO=VCM3/(FU3+tiny)*F3
COVAL(CH4SINK,FU3,GRND,0.0)
  CO=VCM4/(FU4+tiny)*F4
COVAL(CH4SINK,FU4,GRND,0.0)
  CO=VCM5/(FU5+tiny)*F5
COVAL(CH4SINK,FU5,GRND,0.0)
       * Within-fluid CO sources
PATCH(COSOUR,PHASEM,1,NX,1,NY,1,NZ,1,1)
  VAL=56./32.*VCM1*F1
COVAL(COSOUR,CO1,FIXFLU,GRND)
  VAL=56./32.*VCM2*F2
COVAL(COSOUR,CO2,FIXFLU,GRND)
  VAL=56./32.*VCM3*F3
COVAL(COSOUR,CO3,FIXFLU,GRND)
  VAL=56./32.*VCM4*F4
COVAL(COSOUR,CO4,FIXFLU,GRND)
  VAL=56./32.*VCM5*F5
COVAL(COSOUR,CO5,FIXFLU,GRND)

PATCH(COSINK,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)
   CO=CCO1/(CO1+tiny)*F1
COVAL(COSINK,CO1,GRND,0.0)
   CO=CCO2/(CO2+tiny)*F2
COVAL(COSINK,CO2,GRND,0.0)
   CO=CCO3/(CO3+tiny)*F3
COVAL(COSINK,CO3,GRND,0.0)
   CO=CCO4/(CO4+tiny)*F4
COVAL(COSINK,CO4,GRND,0.0)
   CO=CCO5/(CO5+tiny)*F5
COVAL(COSINK,CO5,GRND,0.0)

     ===========================================
store(HCOF,VRAT,VCOM,VCOD,VCOK,CCOM,CCOD,CCOK,WMIX,KD,KC,KT,BURN)

        # Volatile release, VRAT, kg/s
   VRAT=2.e3*RHO2A*COL2/EXP(2829./AMAX1(300.,TMP2))*       $
         RHO2*R2*VOL*RHO2A
   ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
        # Volatile combustion, VCOM, 1/s
         * Single fluid
           ** 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+0.*1./VCOK)
         * 1st fluid
   VCOD=4.*EPKE*AMIN1(FU1,32.*OX1/96.)
   VCOK=1.15e9*DEN1*OX1**1.3/EXP(24444./AMAX1(300.,TF1))/  $
               (FU1**0.3+tiny)
   VCM1=1./(1./VCOD+0.*1./VCOK)
         * 2nd fluid
   VCOD=4.*EPKE*AMIN1(FU2,32.*OX2/96.)
   VCOK=1.15e9*DEN1*OX2**1.3/EXP(24444./AMAX1(300.,TF2))/  $
               (FU2**0.3+tiny)
   VCM2=1./(1./VCOD+0.*1./VCOK)
         * 3rd fluid
   VCOD=4.*EPKE*AMIN1(FU3,32.*OX3/96.)
   VCOK=1.15e9*DEN1*OX3**1.3/EXP(24444./AMAX1(300.,TF3))/  $
               (FU3**0.3+tiny)
   VCM3=1./(1./VCOD+0.*1./VCOK+tiny)
         * 4th fluid
   VCOD=4.*EPKE*AMIN1(FU4,32.*OX4/96.)
   VCOK=1.15e9*DEN1*OX4**1.3/EXP(24444./AMAX1(300.,TF4))/  $
               (FU4**0.3+tiny)
   VCM4=1./(1./VCOD+0.*1./VCOK+tiny)
         * 5th fluid
   VCOD=4.*EPKE*AMIN1(FU5,32.*OX5/96.)
   VCOK=1.15e9*DEN1*OX5**1.3/EXP(24444./AMAX1(300.,TF5))/  $
               (FU5**0.3+tiny)
   VCM5=1./(1./VCOD+0.*1./VCOK+tiny)

    ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
        # Carbon monoxide (CO) combustion, CCOM, 1/s
         * Single fluid
           ** EBU
   CCOD=4.*EPKE*AMIN1(YCO,32.*YO2/56.)
           ** Kinetic
  CCOK=5.4e9*DEN1*YCO*YO2**0.25/EXP(15152./AMAX1(300.,TMP1))
           ** Harmonic blending
   CCOM=1./(1./CCOD+1./CCOK)
         * 1st fluid
   CCOD=4.*EPKE*AMIN1(CO1,32.*OX1/56.)
   CCOK=5.4e9*DEN1*CO1*OX1**0.25/EXP(15152./AMAX1(300.,TF1))
   CCO1=1./(1./CCOD+1./CCOK)
         * 2nd fluid
   CCOD=4.*EPKE*AMIN1(CO2,32.*OX2/56.)
   CCOK=5.4e9*DEN1*CO2*OX2**0.25/EXP(15152./AMAX1(300.,TF2))
   CCO2=1./(1./CCOD+1./CCOK)
         * 3rd fluid
   CCOD=4.*EPKE*AMIN1(CO3,32.*OX3/56.)
   CCOK=5.4e9*DEN1*CO3*OX3**0.25/EXP(15152./AMAX1(300.,TF3))
   CCO3=1./(1./CCOD+1./CCOK)
         * 4th fluid
   CCOD=4.*EPKE*AMIN1(CO4,32.*OX4/56.)
   CCOK=5.4e9*DEN1*CO4*OX4**0.25/EXP(15152./AMAX1(300.,TF4))
   CCO4=1./(1./CCOD+1./CCOK)
         * 5th fluid
   CCOD=4.*EPKE*AMIN1(CO5,32.*OX5/56.)
   CCOK=5.4e9*DEN1*CO5*OX5**0.25/EXP(15152./AMAX1(300.,TF5))
   CCO5=1./(1./CCOD+1./CCOK)

   ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

    # Char burning rate constants:
           ** Mixture moelecular mass
   WMIX=28.8
           ** Diffusion
   KD=2.*1.78e-5*12./8314./AMAX1(300.,TMP2)/SIZE
           ** Kinetic
   KC=0.1309/EXP(26850./AMAX1(300.,TMP2))
           ** Harmonic blending
   KT=1./(1./KD+1./KC)
    # Char (C) burning out, BURN, kg/s
         * Single fluid
   BURN=6./SIZE*KT*(P1+PRESS0)*YO2*WMIX/32.*R2*VOL
         * Each of 5 fluids
   BRN1=6./SIZE*KT*(P1+PRESS0)*OX1*WMIX/32.*R2*VOL
   BRN2=6./SIZE*KT*(P1+PRESS0)*OX2*WMIX/32.*R2*VOL
   BRN3=6./SIZE*KT*(P1+PRESS0)*OX3*WMIX/32.*R2*VOL
   BRN4=6./SIZE*KT*(P1+PRESS0)*OX4*WMIX/32.*R2*VOL
   BRN5=6./SIZE*KT*(P1+PRESS0)*OX5*WMIX/32.*R2*VOL
       =================================================
    # 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)
         * Single-fluid 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.
         * 1st fluid
   FK1=1.8e8*EXP(-38370./AMAX1(TF1,300.))
   MOX=3.97e5*(OX1*DEN1/32./AMAX1(TF1,300.))**.5           $
    *EXP(-31090./AMAX1(TF1,300.))
   RNO1=30.*2.*FK1*MOX*YN2*DEN1/28.
         * 2nd fluid
   FK1=1.8e8*EXP(-38370./AMAX1(TF2,300.))
   MOX=3.97e5*(OX2*DEN1/32./AMAX1(TF2,300.))**.5           $
    *EXP(-31090./AMAX1(TF2,300.))
   RNO2=30.*2.*FK1*MOX*YN2*DEN1/28.
         * 3rd fluid
   FK1=1.8e8*EXP(-38370./AMAX1(TF3,300.))
   MOX=3.97e5*(OX3*DEN1/32./AMAX1(TF3,300.))**.5           $
    *EXP(-31090./AMAX1(TF3,300.))
   RNO3=30.*2.*FK1*MOX*YN2*DEN1/28.
         * 4th fluid
   FK1=1.8e8*EXP(-38370./AMAX1(TF4,300.))
   MOX=3.97e5*(OX4*DEN1/32./AMAX1(TF4,300.))**.5           $
    *EXP(-31090./AMAX1(TF4,300.))
   RNO4=30.*2.*FK1*MOX*YN2*DEN1/28.
         * 5th fluid
   FK1=1.8e8*EXP(-38370./AMAX1(TF5,300.))
   MOX=3.97e5*(OX5*DEN1/32./AMAX1(TF5,300.))**.5           $
    *EXP(-31090./AMAX1(TF5,300.))
   RNO5=30.*2.*FK1*MOX*YN2*DEN1/28.

    ===================================================
       * COefficients for child fluids
  COF2=F1*F3+F1*F4/2.+F1*F5/3.+tiny
  COF3=F2*F4+F1*F4/2.+F2*F5/2.+F1*F5/3.+tiny
  COF4=F3*F5+F2*F5/2.+F1*F5/3.+tiny
       * VALues for child enthalpies
  VHF2=HF1/2.+(F1*F3/2.*HF3+F1*F4/4.*HF4+F1*F5/6.*HF5)/COF2
  VHF3=(F1*F4/4.+F1*F5/6.)*HF1+(F2*F4/2.+F2*F5/4.)*HF2
  VHF3=VHF3+(F1*F4/4.+F2*F4/2.)*HF4+(F1*F5/6.+F2*F5/4.)*HF5
  VHF3=VHF3/COF3
  VHF4=HF5/2.+(F3*F5/2.*HF3+F2*F5/4.*HF2+F1*F5/6.*HF1)/COF4
       * VALues for child O2
  VOX2=OX1/2.+(F1*F3/2.*OX3+F1*F4/4.*OX4+F1*F5/6.*OX5)/COF2
  VOX3=(F1*F4/4.+F1*F5/6.)*OX1+(F2*F4/2.+F2*F5/4.)*OX2
  VOX3=VOX3+(F1*F4/4.+F2*F4/2.)*OX4+(F1*F5/6.+F2*F5/4.)*OX5
  VOX3=VOX3/COF3
  VOX4=OX5/2.+(F3*F5/2.*OX3+F2*F5/4.*OX2+F1*F5/6.*OX1)/COF4
       * VALues for child NO
  VNO2=NO1/2.+(F1*F3/2.*NO3+F1*F4/4.*NO4+F1*F5/6.*NO5)/COF2
  VNO3=(F1*F4/4.+F1*F5/6.)*NO1+(F2*F4/2.+F2*F5/4.)*NO2
  VNO3=VNO3+(F1*F4/4.+F2*F4/2.)*NO4+(F1*F5/6.+F2*F5/4.)*NO5
  VNO3=VNO3/COF3
  VNO4=NO5/2.+(F3*F5/2.*NO3+F2*F5/4.*NO2+F1*F5/6.*NO1)/COF4
       * VALues for child CH4
  VFU2=FU1/2.+(F1*F3/2.*FU3+F1*F4/4.*FU4+F1*F5/6.*FU5)/COF2
  VFU3=(F1*F4/4.+F1*F5/6.)*FU1+(F2*F4/2.+F2*F5/4.)*FU2
  VFU3=VFU3+(F1*F4/4.+F2*F4/2.)*FU4+(F1*F5/6.+F2*F5/4.)*FU5
  VFU3=VFU3/COF3
  VFU4=FU5/2.+(F3*F5/2.*FU3+F2*F5/4.*FU2+F1*F5/6.*FU1)/COF4
       * VALues for child CO
  VCO2=CO1/2.+(F1*F3/2.*CO3+F1*F4/4.*CO4+F1*F5/6.*CO5)/COF2
  VCO3=(F1*F4/4.+F1*F5/6.)*CO1+(F2*F4/2.+F2*F5/4.)*CO2
  VCO3=VCO3+(F1*F4/4.+F2*F4/2.)*CO4+(F1*F5/6.+F2*F5/4.)*CO5
  VCO3=VCO3/COF3
  VCO4=CO5/2.+(F3*F5/2.*CO3+F2*F5/4.*CO2+F1*F5/6.*CO1)/COF4
       * Temperatures of the fluids
  TF1=HF1/1300.
  TF2=HF2/1300.
  TF3=HF3/1300.
  TF4=HF4/1300.
  TF5=HF5/1300.
       * Fluid-averaged properties
  HFAV=F1*HF1+F2*HF2+F3*HF3+F4*HF4+F5*HF5
  TFAV=F1*TF1+F2*TF2+F3*TF3+F4*TF4+F5*TF5
  OXAV=F1*OX1+F2*OX2+F3*OX3+F4*OX4+F5*OX5
  NOAV=F1*NO1+F2*NO2+F3*NO3+F4*NO4+F5*NO5
  FUAV=F1*FU1+F2*FU2+F3*FU3+F4*FU4+F5*FU5
  COAV=F1*CO1+F2*CO2+F3*CO3+F4*CO4+F5*CO5
store(FSUM,FMIX,AVEF,MNSQ)
   FSUM=F1+F2+F3+F4+F5
   FMIX=F1*0.0+F2*0.25+F3*0.5+F4*0.75+F5*1.

  PLANTEND
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(W1  ) =  2.000000E+01 ;FIINIT(W2  ) =  2.000000E+01
 FIINIT(R1  ) =  9.990000E-01 ;FIINIT(R2  ) =  1.000000E-03
 FIINIT(RS  ) =  5.000000E-01 ;FIINIT(YO2 ) =  2.320000E-01
 FIINIT(COL2) =  1.000000E+00 ;FIINIT(EPKE) =  1.000000E+00
   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

 PATCH (AGRAV   ,PHASEM,9,0,0,0,0,0,1,1)
 COVAL (AGRAV   ,V2  , FIXFLU      ,-9.810000E+00)

 PATCH (MICROMIX,PHASEM,0,0,0,0,0,0,1,1)
 COVAL (MICROMIX,F5  , GRND        , SAME        )
 COVAL (MICROMIX,F4  , GRND        , SAME        )
 COVAL (MICROMIX,F3  , GRND        , SAME        )
 COVAL (MICROMIX,F2  , GRND        , SAME        )
 COVAL (MICROMIX,F1  , GRND        , SAME        )

 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   ,F5  , FIXFLU      , GRND        )
 COVAL (PLSR5   ,YCH4, FIXFLU      , GRND        )

 PATCH (PLSR6   ,PHASEM,1,5,1,5,1,20,1,1)
 COVAL (PLSR6   ,YH2O, FIXFLU      , GRND        )
 COVAL (PLSR6   ,YCO , FIXFLU      , GRND        )
 COVAL (PLSR6   ,YO2 , GRND        , 0.000000E+00)
 COVAL (PLSR6   ,YCH4, GRND        , 0.000000E+00)

 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   ,F1  , FIXFLU      , GRND        )
 COVAL (PLSR9   ,YCO2, FIXFLU      , GRND        )
 COVAL (PLSR9   ,YO2 , FIXFLU      , GRND        )

 PATCH (PLSR7   ,PHASEM,1,5,1,5,1,20,1,1)
 COVAL (PLSR7   ,YCO2, FIXFLU      , GRND        )
 COVAL (PLSR7   ,YCO , GRND        , 0.000000E+00)
 COVAL (PLSR7   ,YO2 , GRND        , 0.000000E+00)

 PATCH (PLSR3   ,PHASEM,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 (PLSR90  ,CELL  ,1,5,1,5,1,20,1,1)
 COVAL (PLSR90  ,YH2O, FIXFLU      , GRND        )
 COVAL (PLSR90  ,YCO2, FIXFLU      , GRND        )
 COVAL (PLSR90  ,YCO , FIXFLU      , GRND        )
 COVAL (PLSR90  ,YO2 , FIXFLU      , GRND        )
 COVAL (PLSR90  ,COL2, FIXFLU      , GRND        )
 COVAL (PLSR90  ,CHA2, FIXFLU      , GRND        )

 PATCH (NOXSOR  ,VOLUME,1,5,1,5,1,20,1,1)
 COVAL (NOXSOR  ,YNOX, FIXFLU      , GRND        )

 PATCH (MICRHF  ,PHASEM,1,5,1,5,1,20,1,1)
 COVAL (MICRHF  ,HF4 , GRND        , GRND        )
 COVAL (MICRHF  ,HF3 , GRND        , GRND        )
 COVAL (MICRHF  ,HF2 , GRND        , GRND        )

 PATCH (MICROX  ,PHASEM,1,5,1,5,1,20,1,1)
 COVAL (MICROX  ,OX4 , GRND        , GRND        )
 COVAL (MICROX  ,OX3 , GRND        , GRND        )
 COVAL (MICROX  ,OX2 , GRND        , GRND        )

 PATCH (MICRNO  ,PHASEM,1,5,1,5,1,20,1,1)
 COVAL (MICRNO  ,NO4 , GRND        , GRND        )
 COVAL (MICRNO  ,NO3 , GRND        , GRND        )
 COVAL (MICRNO  ,NO2 , GRND        , GRND        )

 PATCH (MICRFU  ,PHASEM,1,5,1,5,1,20,1,1)
 COVAL (MICRFU  ,FU4 , GRND        , GRND        )
 COVAL (MICRFU  ,FU3 , GRND        , GRND        )
 COVAL (MICRFU  ,FU2 , GRND        , GRND        )

 PATCH (MICCO   ,PHASEM,1,5,1,5,1,20,1,1)
 COVAL (MICCO   ,CO4 , GRND        , GRND        )
 COVAL (MICCO   ,CO3 , GRND        , GRND        )
 COVAL (MICCO   ,CO2 , GRND        , GRND        )

 PATCH (COMHEATF,CELL  ,1,5,1,5,1,20,1,1)
 COVAL (COMHEATF,HF5 , FIXFLU      , GRND        )
 COVAL (COMHEATF,HF4 , FIXFLU      , GRND        )
 COVAL (COMHEATF,HF3 , FIXFLU      , GRND        )
 COVAL (COMHEATF,HF2 , FIXFLU      , GRND        )
 COVAL (COMHEATF,HF1 , FIXFLU      , GRND        )

 PATCH (INTHEATF,VOLUME,1,5,1,5,1,20,1,1)
 COVAL (INTHEATF,HF5 , GRND        , GRND        )
 COVAL (INTHEATF,HF4 , GRND        , GRND        )
 COVAL (INTHEATF,HF3 , GRND        , GRND        )
 COVAL (INTHEATF,HF2 , GRND        , GRND        )
 COVAL (INTHEATF,HF1 , GRND        , GRND        )

 PATCH (OXYGSOUR,CELL  ,1,5,1,5,1,20,1,1)
 COVAL (OXYGSOUR,OX5 , GRND        , 0.000000E+00)
 COVAL (OXYGSOUR,OX4 , GRND        , 0.000000E+00)
 COVAL (OXYGSOUR,OX3 , GRND        , 0.000000E+00)
 COVAL (OXYGSOUR,OX2 , GRND        , 0.000000E+00)
 COVAL (OXYGSOUR,OX1 , GRND        , 0.000000E+00)

 PATCH (NOSOUR  ,VOLUME,1,5,1,5,1,20,1,1)
 COVAL (NOSOUR  ,NO5 , FIXFLU      , GRND        )
 COVAL (NOSOUR  ,NO4 , FIXFLU      , GRND        )
 COVAL (NOSOUR  ,NO3 , FIXFLU      , GRND        )
 COVAL (NOSOUR  ,NO2 , FIXFLU      , GRND        )
 COVAL (NOSOUR  ,NO1 , FIXFLU      , GRND        )

 PATCH (CH4SOUR ,CELL  ,1,5,1,5,1,20,1,1)
 COVAL (CH4SOUR ,FU5 , FIXFLU      , GRND        )
 COVAL (CH4SOUR ,FU4 , FIXFLU      , GRND        )
 COVAL (CH4SOUR ,FU3 , FIXFLU      , GRND        )
 COVAL (CH4SOUR ,FU2 , FIXFLU      , GRND        )
 COVAL (CH4SOUR ,FU1 , FIXFLU      , GRND        )

 PATCH (CH4SINK ,CELL  ,1,5,1,5,1,20,1,1)
 COVAL (CH4SINK ,FU5 , GRND        , 0.000000E+00)
 COVAL (CH4SINK ,FU4 , GRND        , 0.000000E+00)
 COVAL (CH4SINK ,FU3 , GRND        , 0.000000E+00)
 COVAL (CH4SINK ,FU2 , GRND        , 0.000000E+00)
 COVAL (CH4SINK ,FU1 , GRND        , 0.000000E+00)

 PATCH (COSOUR  ,PHASEM,1,5,1,5,1,20,1,1)
 COVAL (COSOUR  ,CO5 , FIXFLU      , GRND        )
 COVAL (COSOUR  ,CO4 , FIXFLU      , GRND        )
 COVAL (COSOUR  ,CO3 , FIXFLU      , GRND        )
 COVAL (COSOUR  ,CO2 , FIXFLU      , GRND        )
 COVAL (COSOUR  ,CO1 , FIXFLU      , GRND        )

 PATCH (COSINK  ,PHASEM,1,5,1,5,1,20,1,1)
 COVAL (COSINK  ,CO5 , GRND        , 0.000000E+00)
 COVAL (COSINK  ,CO4 , GRND        , 0.000000E+00)
 COVAL (COSINK  ,CO3 , GRND        , 0.000000E+00)
 COVAL (COSINK  ,CO2 , GRND        , 0.000000E+00)
 COVAL (COSINK  ,CO1 , GRND        , 0.000000E+00)

 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP  =    5000
 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  ,FALSDT, 1.000000E-02)
 RELAX(EP  ,FALSDT, 1.000000E-02)
 RELAX(H1  ,FALSDT, 1.000000E-02)
 RELAX(H2  ,FALSDT, 1.000000E-02)
 RELAX(CO5 ,FALSDT, 3.000000E-03)
 RELAX(CO4 ,FALSDT, 3.000000E-03)
 RELAX(CO3 ,FALSDT, 3.000000E-03)
 RELAX(CO2 ,FALSDT, 3.000000E-03)
 RELAX(CO1 ,FALSDT, 3.000000E-03)
 RELAX(FU5 ,FALSDT, 3.000000E-03)
 RELAX(FU4 ,FALSDT, 3.000000E-03)
 RELAX(FU3 ,FALSDT, 3.000000E-03)
 RELAX(FU2 ,FALSDT, 3.000000E-03)
 RELAX(FU1 ,FALSDT, 3.000000E-03)
 RELAX(NO5 ,FALSDT, 1.000000E-03)
 RELAX(NO4 ,FALSDT, 1.000000E-03)
 RELAX(NO3 ,FALSDT, 1.000000E-03)
 RELAX(NO2 ,FALSDT, 1.000000E-03)
 RELAX(NO1 ,FALSDT, 1.000000E-03)
 RELAX(OX5 ,FALSDT, 3.000000E-03)
 RELAX(OX4 ,FALSDT, 3.000000E-03)
 RELAX(OX3 ,FALSDT, 3.000000E-03)
 RELAX(OX2 ,FALSDT, 3.000000E-03)
 RELAX(OX1 ,FALSDT, 3.000000E-03)
 RELAX(HF5 ,FALSDT, 3.000000E-03)
 RELAX(HF4 ,FALSDT, 3.000000E-03)
 RELAX(F5  ,FALSDT, 3.000000E-03)
 RELAX(F4  ,FALSDT, 3.000000E-03)
 RELAX(F3  ,FALSDT, 3.000000E-03)
 RELAX(F2  ,FALSDT, 3.000000E-03)
 RELAX(F1  ,FALSDT, 3.000000E-03)
 RELAX(HF3 ,FALSDT, 3.000000E-03)
 RELAX(HF2 ,FALSDT, 3.000000E-03)
 RELAX(YNOX,FALSDT, 1.000000E-03)
 RELAX(HF1 ,FALSDT, 3.000000E-03)
 RELAX(YH2O,FALSDT, 3.000000E-03)
 RELAX(YCO2,FALSDT, 3.000000E-03)
 RELAX(YCO ,FALSDT, 3.000000E-03)
 RELAX(YO2 ,FALSDT, 3.000000E-03)
 RELAX(YCH4,FALSDT, 3.000000E-03)
 RELAX(COL2,FALSDT, 1.000000E-02)
 RELAX(CHA2,FALSDT, 1.000000E-02)
 RELAX(DEN1,LINRLX, 1.000000E-01)
 RELAX(CMDO,LINRLX, 3.000000E-01)
 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(CO5 ) = 1.000000E+00 ;VARMIN(CO5 ) =.0
 VARMAX(CO4 ) = 1.000000E+00 ;VARMIN(CO4 ) =.0
 VARMAX(CO3 ) = 1.000000E+00 ;VARMIN(CO3 ) =.0
 VARMAX(CO2 ) = 1.000000E+00 ;VARMIN(CO2 ) =.0
 VARMAX(CO1 ) = 1.000000E+00 ;VARMIN(CO1 ) =.0
 VARMAX(FU5 ) = 1.000000E+00 ;VARMIN(FU5 ) =.0
 VARMAX(FU4 ) = 1.000000E+00 ;VARMIN(FU4 ) =.0
 VARMAX(FU3 ) = 1.000000E+00 ;VARMIN(FU3 ) =.0
 VARMAX(FU2 ) = 1.000000E+00 ;VARMIN(FU2 ) =.0
 VARMAX(FU1 ) = 1.000000E+00 ;VARMIN(FU1 ) =.0
 VARMAX(NO5 ) = 1.000000E+00 ;VARMIN(NO5 ) = 0.000000E+00
 VARMAX(NO4 ) = 1.000000E+00 ;VARMIN(NO4 ) = 0.000000E+00
 VARMAX(NO3 ) = 1.000000E+00 ;VARMIN(NO3 ) = 0.000000E+00
 VARMAX(NO2 ) = 1.000000E+00 ;VARMIN(NO2 ) = 0.000000E+00
 VARMAX(NO1 ) = 1.000000E+00 ;VARMIN(NO1 ) = 0.000000E+00
 VARMAX(OX5 ) = 1.000000E+00 ;VARMIN(OX5 ) = 0.000000E+00
 VARMAX(OX4 ) = 1.000000E+00 ;VARMIN(OX4 ) = 0.000000E+00
 VARMAX(OX3 ) = 1.000000E+00 ;VARMIN(OX3 ) = 0.000000E+00
 VARMAX(OX2 ) = 1.000000E+00 ;VARMIN(OX2 ) = 0.000000E+00
 VARMAX(OX1 ) = 1.000000E+00 ;VARMIN(OX1 ) = 0.000000E+00
 VARMAX(F5  ) = 1.000000E+00 ;VARMIN(F5  ) = 0.000000E+00
 VARMAX(F4  ) = 1.000000E+00 ;VARMIN(F4  ) = 0.000000E+00
 VARMAX(F3  ) = 1.000000E+00 ;VARMIN(F3  ) = 0.000000E+00
 VARMAX(F2  ) = 1.000000E+00 ;VARMIN(F2  ) = 0.000000E+00
 VARMAX(F1  ) = 1.000000E+00 ;VARMIN(F1  ) = 0.000000E+00
 VARMAX(YNOX) = 1.000000E+00 ;VARMIN(YNOX) = 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(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(ASH2) = 1.000000E+00 ;VARMIN(ASH2) = 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
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 USEGRD  =    T  ;USEGRX =    T
 NAMSAT  =MOSG
 GENK    =    T
 ASAP    =    T
 SPEDAT(SET,MFM,MFMMOD,C,MFM)
 SPEDAT(SET,MFM,NFLUIDS,I,5)
 SPEDAT(SET,MFM,NFLR,I,1)
 SPEDAT(SET,MFM,NFLF,I,5)
 SPEDAT(SET,MFM,CONREA,R,0.00000E+00)
 SPEDAT(SET,MFM,CONMIX,R,5.00000E+00)
 SPEDAT(SET,MFM,VISCON,R,0.00000E+00)
 ************************************************************
  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   =       5 ;IYMON  =       3 ;IZMON  =      15
 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,-5.320389E-01,4.837199E-01,-6.949458E-01)
 GVIEW(UP,3.406668E-01,8.736793E-01,3.473192E-01)

> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+01
> DOM,    MONIT,       9.000000E-01, 5.000000E-01, 7.250000E+00
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    SNAPSIZE,    1.000000E-02
> GRID,   RSET_Y_1,      5, 1.000000E+00
> DOM,    RELAX,       3.000000E-01

> OBJ,    NAME,        E-WALL
> OBJ,    POSITION,    1.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.000000E+00, 1.000000E+01
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    ADIABATI1,     0.000000E+00, 0.000000E+00
> OBJ,    ADIABATI2,     0.000000E+00, 0.000000E+00

> OBJ,    NAME,        W-WALL
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, 1.000000E+00, 1.000000E+01
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    VISIBLE,          -1
> OBJ,    TYPE,        PLATE
> OBJ,    ADIABATI1,     0.000000E+00, 0.000000E+00
> OBJ,    ADIABATI2,     0.000000E+00, 0.000000E+00

> OBJ,    NAME,        H-WALL
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 1.000000E+01
> OBJ,    SIZE,        1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    ADIABATI1,     0.000000E+00, 0.000000E+00
> OBJ,    ADIABATI2,     0.000000E+00, 0.000000E+00

> OBJ,    NAME,        L-WALL
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    VISIBLE,          -1
> OBJ,    TYPE,        PLATE
> OBJ,    ADIABATI1,     0.000000E+00, 0.000000E+00
> OBJ,    ADIABATI2,     0.000000E+00, 0.000000E+00

> OBJ,    NAME,        N-WALL
> OBJ,    POSITION,    0.000000E+00, 1.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 0.000000E+00, 9.000000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    ADIABATI1,     0.000000E+00, 0.000000E+00
> OBJ,    ADIABATI2,     0.000000E+00, 0.000000E+00

> OBJ,    NAME,        S-WALL
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 1.000000E+00
> OBJ,    SIZE,        1.000000E+00, 0.000000E+00, 9.000000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    ADIABATI1,     0.000000E+00, 0.000000E+00
> OBJ,    ADIABATI2,     0.000000E+00, 0.000000E+00

> OBJ,    NAME,        OUTLET
> OBJ,    POSITION,    0.000000E+00, 1.000000E+00, 9.000000E+00
> OBJ,    SIZE,        6.000000E-01, 0.000000E+00, 1.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,        INLET
> OBJ,    POSITION,    4.000000E-01, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        6.000000E-01, 0.000000E+00, 1.000000E+00
> OBJ,    CLIPART,     cube3t
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        INLET
> OBJ,    DENSITY,       4.620000E-01
> OBJ,    DENSITY-2,     1.000000E-01
> 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_OX5,     2.320000E-01
> OBJ,    INLET_OX4,     2.320000E-01
> OBJ,    INLET_OX3,     2.320000E-01
> OBJ,    INLET_OX2,     2.320000E-01
> OBJ,    INLET_OX1,     2.320000E-01
> OBJ,    INLET_HF5,     5.250000E+05
> OBJ,    INLET_HF4,     5.250000E+05
> OBJ,    INLET_F1,      1.000000E+00
> OBJ,    INLET_HF3,     5.250000E+05
> OBJ,    INLET_HF2,     5.250000E+05
> OBJ,    INLET_HF1,     5.250000E+05
> OBJ,    INLET_YO2,     2.320000E-01
> OBJ,    INLET_COL2,    8.000000E-01
> OBJ,    TURB-INTENS,   5.000000E+00

> OBJ,    NAME,        GRAVITY
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+01
> OBJ,    CLIPART,     cube14
> OBJ,    ROTATION24,        1
> OBJ,    VISIBLE,          -1
> OBJ,    TYPE,        USER_DEFINED

> OBJ,    NAME,        N-WALL2
> OBJ,    POSITION,    6.000000E-01, 1.000000E+00, 9.000000E+00
> OBJ,    SIZE,        4.000000E-01, 0.000000E+00, 1.000000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    ADIABATI1,     0.000000E+00, 0.000000E+00
> OBJ,    ADIABATI2,     0.000000E+00, 0.000000E+00

> OBJ,    NAME,        S-WALL2
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        4.000000E-01, 0.000000E+00, 1.000000E+00
> OBJ,    CLIPART,     cube11
> OBJ,    ROTATION24,        1
> OBJ,    TYPE,        PLATE
> OBJ,    ADIABATI1,     0.000000E+00, 0.000000E+00
> OBJ,    ADIABATI2,     0.000000E+00, 0.000000E+00
STOP