************************************************************

   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