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

   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

  rot z ang -10

  con tfav x 1 fil;.001

  con tfav z 10 fil;.001

  con tfav x 5 fil;.001

  vec x m; vec z m; vec y m

  msg( Fluid-average temperature

  pause;con cl;vec cl;red

  con tmp1 x 1 fil;.001

  con tmp1 z 10 fil;.001

  con tmp1 x 5 fil;.001

  msg(  Single-fluid temperature

  pause;con cl;red

  con crav x 1 fil;.001

  con crav z 10 fil;.001

  con crav x 5 fil;.001

  msg(  Fluid-average radiation flux

  pause;con cl;red

  con crad x 1 fil;.001

  con crad z 10 fil;.001

  con crad x 5 fil;.001

  msg(  Single-fluid radiation flux

  pause;con cl;red

  con cr1 x 1 fil;.001

  con cr1 z 10 fil;.001

  con cr1 x 5 fil;.001

  msg(  1st-fluid radiation

  pause;con cl;red

  con cr2 x 1 fil;.001

  con cr2 z 10 fil;.001

  con cr2 x 5 fil;.001

  msg(  2nd-fluid radiation

  pause;con cl;red

  con cr3 x 1 fil;.001

  con cr3 z 10 fil;.001

  con cr3 x 5 fil;.001

  msg(  3rd-fluid radiation

  pause;con cl;red

  con cr4 x 1 fil;.001

  con cr4 z 10 fil;.001

  con cr4 x 5 fil;.001

  msg(  4th-fluid radiation

  pause;con cl;red

  con cr5 x 1 fil;.001

  con cr5 z 10 fil;.001

  con cr5 x 5 fil;.001

  msg(  5th-fluid radiation

  pause;con cl;red

  con f1 x 1 fil;.001

  con f1 z 10 fil;.001

  con f1 x 5 fil;.001

  msg(  1st-fluid 'presence probability'

  pause;con cl;red

  con f2 x 1 fil;.001

  con f2 z 10 fil;.001

  con f2 x 5 fil;.001

  msg(  2nd-fluid 'presence probability'

  pause;con cl;red

  con f3 x 1 fil;.001

  con f3 z 10 fil;.001

  con f3 x 5 fil;.001

  msg(  3th-fluid 'presence probability'

  pause;con cl;red

  con f4 x 1 fil;.001

  con f4 z 10 fil;.001

  con f4 x 5 fil;.001

  msg(  4th-fluid 'presence probability'

  pause;con cl;red

  con f5 x 1 fil;.001

  con f5 z 10 fil;.001

  con f5 x 5 fil;.001

  msg(  5th-fluid 'presence probability'

  ENDUSE



  DISPLAY



    MFM: 5 fluid population



    Reaction: CO   + 0.5 (O2+0.86N2) -> CO2 + 0.43N2

         i.e. FUEL + 0.5 OXID        -> PROD



    Stoichiometric ratio: S=1.



    Rate: Rfu = 4*rho*EP/KE*min(FUEL,OXID)



    Radiation: P-1 model



  ENDDIS



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

 IRUNN   =       1 ;LIBREF =     302

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

  

  Group 1. Run Title

 TEXT(P-1 radiation model for 5-fluid MFM      )

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

  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  =    T

    * Non-default variable names

 NAME( 70) =RFU5 ; NAME( 71) =RFU4

 NAME( 72) =RFU3 ; NAME( 73) =RFU2

 NAME( 74) =RFU1 ; NAME( 75) =SUMA

 NAME( 76) =SUM5 ; NAME( 77) =SUM4

 NAME( 78) =SUM3 ; NAME( 79) =SUM2

 NAME( 80) =SUM1 ; NAME( 81) =SUMS

 NAME( 82) =PR5  ; NAME( 83) =PR4

 NAME( 84) =PR3  ; NAME( 85) =PR2

 NAME( 86) =PR1  ; NAME( 87) =PRAV

 NAME( 88) =VPR4 ; NAME( 89) =VPR3

 NAME( 90) =VPR2 ; NAME( 91) =OXAV

 NAME( 92) =VOX4 ; NAME( 93) =VOX3

 NAME( 94) =VOX2 ; NAME( 95) =OX5

 NAME( 96) =OX4  ; NAME( 97) =OX3

 NAME( 98) =OX2  ; NAME( 99) =OX1

 NAME(100) =VFU5 ; NAME(101) =VFU4

 NAME(102) =VFU3 ; NAME(103) =VFU2

 NAME(104) =FUAV ; NAME(105) =FU5

 NAME(106) =FU4  ; NAME(107) =FU3

 NAME(108) =FU2  ; NAME(109) =FU1

 NAME(110) =VHF4 ; NAME(111) =VHF3

 NAME(112) =VHF2 ; NAME(113) =COF4

 NAME(114) =COF3 ; NAME(115) =COF2

 NAME(116) =TF5  ; NAME(117) =TF4

 NAME(118) =TF3  ; NAME(119) =TF2

 NAME(120) =TF1  ; NAME(121) =TFAV

 NAME(122) =HFAV ; NAME(123) =HF5

 NAME(124) =HF4  ; NAME(125) =HF3

 NAME(126) =HF2  ; NAME(127) =HF1

 NAME(128) =CRAV ; NAME(129) =CR5

 NAME(130) =CR4  ; NAME(131) =CR3

 NAME(132) =CR2  ; NAME(133) =CR1

 NAME(134) =FAV  ; NAME(135) =F5

 NAME(136) =F4   ; NAME(137) =F3

 NAME(138) =F2   ; NAME(139) =F1

 NAME(140) =ENUT ; NAME(141) =EL1

 NAME(142) =EPKE ; NAME(143) =DEN1

 NAME(144) =FSUM ; NAME(145) =FUEL

 NAME(146) =TMP1 ; NAME(147) =PROD

 NAME(148) =OXID ; NAME(149) =CRAD

 NAME(150) =MIXF

    * Solved variables list

 SOLVE(P1  ,U1  ,V1  ,W1  ,H1  ,PR5 ,PR4 ,PR3 )

 SOLVE(PR2 ,PR1 ,OX5 ,OX4 ,OX3 ,OX2 ,OX1 ,FU5 )

 SOLVE(FU4 ,FU3 ,FU2 ,FU1 ,HF5 ,HF4 ,HF3 ,HF2 )

 SOLVE(HF1 ,CR5 ,CR4 ,CR3 ,CR2 ,CR1 ,F5  ,F4  )

 SOLVE(F3  ,F2  ,F1  ,FUEL,CRAD,MIXF)

    * Stored variables list

 STORE(OXID,PROD,TMP1,FSUM,DEN1,EPKE,EL1 ,ENUT)

 STORE(FAV ,CRAV,HFAV,TFAV,TF1 ,TF2 ,TF3 ,TF4 )

 STORE(TF5 ,COF2,COF3,COF4,VHF2,VHF3,VHF4,FUAV)

 STORE(VFU2,VFU3,VFU4,VFU5,VOX2,VOX3,VOX4,OXAV)

 STORE(VPR2,VPR3,VPR4,PRAV,SUMS,SUM1,SUM2,SUM3)

 STORE(SUM4,SUM5,SUMA,RFU1,RFU2,RFU3,RFU4,RFU5)

    * Additional solver options

 SOLUTN(P1  ,Y,Y,Y,N,N,Y)

 TURMOD(KEMODL)



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

  Group 8. Terms & Devices

 TERMS (H1  ,N,Y,Y,Y,Y,Y)

 TERMS (CR5 ,N,N,Y,Y,N,Y)

 TERMS (CR4 ,N,N,Y,Y,Y,Y)

 TERMS (CR3 ,N,N,Y,Y,N,Y)

 TERMS (CR2 ,N,N,Y,Y,Y,Y)

 TERMS (CR1 ,N,N,Y,Y,N,Y)

 TERMS (CRAD,N,N,Y,Y,Y,N)

 NEWRH1  =    T

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

  Group 9. Properties

 RHO1    = GRND6

 RHO1A   = 2.800000E+01 ;RHO1B  = 5.600000E+01

 RHO1C   = 5.600000E+01

 PRESS0  = 1.000000E+05

 TMP1    = GRND8

 TMP1A   = 1.500000E+03 ;TMP1B  = 1.500000E+03

 TMP1C   = 1.500000E+03

 TMP2A   = 5.000000E-01 ;TMP2B  = 7.450000E+06

 TMP2C   = 0.000000E+00

 CP1     = 1.500000E+03

 ENUL    = 1.000000E-05

 DVO1DT  = 3.330000E-03

 PRNDTL(CR5 ) =  GRND         ;PRNDTL(CR4 ) =  GRND

 PRNDTL(CR3 ) =  GRND         ;PRNDTL(CR2 ) =  GRND

 PRNDTL(CR1 ) =  GRND         ;PRNDTL(CRAD) =  GRND

 PRT   (EP  ) =  1.314000E+00 ;PRT   (CR5 ) =  1.000000E+09

 PRT   (CR4 ) =  1.000000E+09 ;PRT   (CR3 ) =  1.000000E+09

 PRT   (CR2 ) =  1.000000E+09 ;PRT   (CR1 ) =  1.000000E+09

 PRT   (CRAD) =  1.000000E+09

 EMISS   = 4.000000E-01 ;SCATT  = 1.000000E-01

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

  Group 10.Inter-Phase Transfer Processes

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

  Echo PLANT settings

  

    PLANTBEGIN



    * Exchange coefficients for incident radiations

   LAMPR(CRAD)=:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)

   LAMPR(CR1) =:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)

   LAMPR(CR2) =:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)

   LAMPR(CR3) =:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)

   LAMPR(CR4) =:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)

   LAMPR(CR5) =:ENUL:*DEN1*3.*(:EMISS:+:SCATT:)



    * Single-fluid sources due to radiation

PATCH(RADSOR,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)

   VAL=:EMISS:*(CRAD-4.*5.68e-8*TMP1**4)

COVAL(RADSOR,H1,FIXFLU,GRND)

   VAL=4.*5.68e-8*TEMP1**4

COVAL(RADSOR,CRAD,EMISS,GRND)



    * Within-fluid sources due to radiation

PATCH(RADFLUID,VOLUME,1,NX,1,NY,1,NZ,1,LSTEP)

      ** 1st fluid

   VAL=:EMISS:*(CR1-4.*5.68e-8*TF1**4)*F1

COVAL(RADFLUID,HF1,FIXFLU,GRND)

   VAL=4.*5.68e-8*TF1**4*F1

COVAL(RADFLUID,CR1,EMISS,GRND)

      ** 2nd fluid

   VAL=:EMISS:*(CR2-4.*5.68e-8*TF2**4)*F2

COVAL(RADFLUID,HF2,FIXFLU,GRND)

   VAL=4.*5.68e-8*TF2**4*F2

COVAL(RADFLUID,CR2,EMISS,GRND)

      ** 3rd fluid

   VAL=:EMISS:*(CR3-4.*5.68e-8*TF3**4)*F3

COVAL(RADFLUID,HF3,FIXFLU,GRND)

   VAL=4.*5.68e-8*TF3**4*F3

COVAL(RADFLUID,CR3,EMISS,GRND)

      ** 4th fluid

   VAL=:EMISS:*(CR4-4.*5.68e-8*TF4**4)*F4

COVAL(RADFLUID,HF4,FIXFLU,GRND)

   VAL=4.*5.68e-8*TF4**4*F4

COVAL(RADFLUID,CR4,EMISS,GRND)

      ** 5th fluid

   VAL=:EMISS:*(CR5-4.*5.68e-8*TF5**4)*F5

COVAL(RADFLUID,HF5,FIXFLU,GRND)

   VAL=4.*5.68e-8*TF5**4*F5

COVAL(RADFLUID,CR5,EMISS,GRND)



    * Within-fluid fuel sources

PATCH(FUELFLUI,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)

   VAL=-Rfu1*F1

COVAL(FUELFLUI,FU1,FIXFLU,GRND)

   VAL=-Rfu2*F2

COVAL(FUELFLUI,FU2,FIXFLU,GRND)

   VAL=-Rfu3*F3

COVAL(FUELFLUI,FU3,FIXFLU,GRND)

   VAL=-Rfu4*F4

COVAL(FUELFLUI,FU4,FIXFLU,GRND)

   VAL=-Rfu5*F5

COVAL(FUELFLUI,FU5,FIXFLU,GRND)



   * Within-fluid oxidant sources

PATCH(OXIDFLUI,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)

   VAL=-Rfu1*F1

COVAL(OXIDFLUI,OX1,FIXFLU,GRND)

   VAL=-Rfu2*F2

COVAL(OXIDFLUI,OX2,FIXFLU,GRND)

   VAL=-Rfu3*F3

COVAL(OXIDFLUI,OX3,FIXFLU,GRND)

   VAL=-Rfu4*F4

COVAL(OXIDFLUI,OX4,FIXFLU,GRND)

   VAL=-Rfu5*F5

COVAL(OXIDFLUI,OX5,FIXFLU,GRND)



  * Within-fluid product sources

PATCH(PRODFLUI,PHASEM,1,NX,1,NY,1,NZ,1,LSTEP)

   VAL=2.*Rfu1*F1

COVAL(PRODFLUI,PR1,FIXFLU,GRND)

   VAL=2.*Rfu2*F2

COVAL(PRODFLUI,PR2,FIXFLU,GRND)

   VAL=2.*Rfu3*F3

COVAL(PRODFLUI,PR3,FIXFLU,GRND)

   VAL=2.*Rfu4*F4

COVAL(PRODFLUI,PR4,FIXFLU,GRND)

   VAL=2.*Rfu5*F5

COVAL(PRODFLUI,PR5,FIXFLU,GRND)



        * Micromixing contributions to CVA transport*

PATCH(MICRCVA,PHASEM,1,NX,1,NY,1,NZ,1,1)

     ** Within-fluid enthalpy

   CO = 5.*EPKE*COF2

   VAL= VHF2

COVAL(MICRCVA,HF2,GRND,GRND)

   CO = 5.*EPKE*COF3

   VAL= VHF3

COVAL(MICRCVA,HF3,GRND,GRND)

   CO = 5.*EPKE*COF4

   VAL= VHF4

COVAL(MICRCVA,HF4,GRND,GRND)

     ** Within-fluid fuel

   CO = 5.*EPKE*COF2

   VAL= VFU2

COVAL(MICRCVA,FU2,GRND,GRND)

   CO = 5.*EPKE*COF3

   VAL= VFU3

COVAL(MICRCVA,FU3,GRND,GRND)

   CO = 5.*EPKE*COF4

   VAL= VFU4

COVAL(MICRCVA,FU4,GRND,GRND)

     ** Within-fluid oxidant

   CO = 5.*EPKE*COF2

   VAL= VOX2

COVAL(MICRCVA,OX2,GRND,GRND)

   CO = 5.*EPKE*COF3

   VAL= VOX3

COVAL(MICRCVA,OX3,GRND,GRND)

   CO = 5.*EPKE*COF4

   VAL= VOX4

COVAL(MICRCVA,OX4,GRND,GRND)

     ** Within-fluid products

   CO = 5.*EPKE*COF2

   VAL= VPR2

COVAL(MICRCVA,PR2,GRND,GRND)

   CO = 5.*EPKE*COF3

   VAL= VPR3

COVAL(MICRCVA,PR3,GRND,GRND)

   CO = 5.*EPKE*COF4

   VAL= VPR4

COVAL(MICRCVA,PR4,GRND,GRND)



        * 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 fluids

          ** Enthalpy

   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

          ** Fuel

   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

          ** Oxidant

   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

          ** Products

   VPR2=PR1/2.+(F1*F3/2.*PR3+F1*F4/4.*PR4+F1*F5/6.*PR5)/COF2

   VPR3=(F1*F4/4.+F1*F5/6.)*PR1+(F2*F4/2.+F2*F5/4.)*PR2

   VPR3=VPR3+(F1*F4/4.+F2*F4/2.)*PR4+(F1*F5/6.+F2*F5/4.)*PR5

   VPR3=VPR3/COF3

   VPR4=PR5/2.+(F3*F5/2.*PR3+F2*F5/4.*PR2+F1*F5/6.*PR1)/COF4

          ** Within-fluid temperatures

              * TMP2B is the neat of combustion

              * TMP1A is specific heat

   TF1=(HF1-:TMP2B:*FU1)/:TMP1A:

   TF2=(HF2-:TMP2B:*FU2)/:TMP1A:

   TF3=(HF3-:TMP2B:*FU3)/:TMP1A:

   TF4=(HF4-:TMP2B:*FU4)/:TMP1A:

   TF5=(HF5-:TMP2B:*FU5)/:TMP1A:

          ** Within-fluid reaction rates

    Rfu1=4.*EPKE*AMIN1(FU1,OX1)

    Rfu2=4.*EPKE*AMIN1(FU2,OX2)

    Rfu3=4.*EPKE*AMIN1(FU3,OX3)

    Rfu4=4.*EPKE*AMIN1(FU4,OX4)

    Rfu5=4.*EPKE*AMIN1(FU5,OX5)



          ** Population averaged properties

             * Mixture fraction

   FAV =F1*0.0+F2*0.25+F3*0.5+F4*0.75+F5*1.

             * Radiation

   CRAV=F1*CR1+F2*CR2+F3*CR3+F4*CR4+F5*CR5

             * Enthalpy

   HFAV=F1*HF1+F2*HF2+F3*HF3+F4*HF4+F5*HF5

             * Temperature

   TFAV=F1*TF1+F2*TF2+F3*TF3+F4*TF4+F5*TF5

             * Fuel mass fraction

   FUAV=F1*FU1+F2*FU2+F3*FU3+F4*FU4+F5*FU5

             * Oxygen mass fraction

   OXAV=F1*OX1+F2*OX2+F3*OX3+F4*OX4+F5*OX5

             * Product mass fraction

   PRAV=F1*PR1+F2*PR2+F3*PR3+F4*PR4+F5*PR5

          ** Checks for continuity

             * Sum of fluid mass fractions

   FSUM=F1+F2+F3+F4+F5

             * Sum of single-fluid components

   SUMS=OXID+FUEL+PROD

             * Sum of 1st-fluid components

   SUM1=OX1+FU1+PR1

             * Sum of 2nd-fluid components

   SUM2=OX2+FU2+PR2

             * Sum of 3rd-fluid components

   SUM3=OX3+FU3+PR3

             * Sum of 4th-fluid components

   SUM4=OX4+FU4+PR4

             * Sum of 5th-fluid components

   SUM5=OX5+FU5+PR5

             * Sum of fluid-averged components

   SUMA=OXAV+FUAV+PRAV



    PLANTEND

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

  Group 11.Initialise Var/Porosity Fields

   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 (MICROMIX,PHASEM,1,5,1,5,1,20,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 (CHSO    ,PHASEM,0,0,0,0,0,0,1,1)

 COVAL (CHSO    ,FUEL, GRND9       , GRND9       )



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

  Group 14. Downstream Pressure For PARAB

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

  Group 15. Terminate Sweeps

 LSWEEP  =    2000

 RESFAC  = 1.000000E-03

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

  Group 16. Terminate Iterations

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

  Group 17. Relaxation

 RELAX(P1  ,LINRLX, 3.000000E-01)

 RELAX(U1  ,FALSDT, 1.000000E-01)

 RELAX(V1  ,FALSDT, 1.000000E-01)

 RELAX(W1  ,FALSDT, 1.000000E-01)

 RELAX(KE  ,FALSDT, 1.000000E-01)

 RELAX(EP  ,FALSDT, 1.000000E-01)

 RELAX(H1  ,FALSDT, 1.000000E-01)

 RELAX(PR5 ,FALSDT, 1.000000E-01)

 RELAX(PR4 ,FALSDT, 1.000000E-01)

 RELAX(PR3 ,FALSDT, 1.000000E-01)

 RELAX(PR2 ,FALSDT, 1.000000E-01)

 RELAX(PR1 ,FALSDT, 1.000000E-01)

 RELAX(OX5 ,FALSDT, 1.000000E-01)

 RELAX(OX4 ,FALSDT, 1.000000E-01)

 RELAX(OX3 ,FALSDT, 1.000000E-01)

 RELAX(OX2 ,FALSDT, 1.000000E-01)

 RELAX(OX1 ,FALSDT, 1.000000E-01)

 RELAX(FU5 ,FALSDT, 1.000000E-01)

 RELAX(FU4 ,FALSDT, 1.000000E-01)

 RELAX(FU3 ,FALSDT, 1.000000E-01)

 RELAX(FU2 ,FALSDT, 1.000000E-01)

 RELAX(FU1 ,FALSDT, 1.000000E-01)

 RELAX(HF5 ,FALSDT, 1.000000E-01)

 RELAX(HF4 ,FALSDT, 1.000000E-01)

 RELAX(HF3 ,FALSDT, 1.000000E-01)

 RELAX(HF2 ,FALSDT, 1.000000E-01)

 RELAX(HF1 ,FALSDT, 1.000000E-01)

 RELAX(CR5 ,FALSDT, 1.000000E-01)

 RELAX(CR4 ,FALSDT, 1.000000E-01)

 RELAX(CR3 ,FALSDT, 1.000000E-01)

 RELAX(CR2 ,FALSDT, 1.000000E-01)

 RELAX(CR1 ,FALSDT, 1.000000E-01)

 RELAX(F5  ,FALSDT, 1.000000E-01)

 RELAX(F4  ,FALSDT, 1.000000E-01)

 RELAX(F3  ,FALSDT, 1.000000E-01)

 RELAX(F2  ,FALSDT, 1.000000E-01)

 RELAX(F1  ,FALSDT, 1.000000E-01)

 RELAX(DEN1,LINRLX, 1.000000E-01)

 RELAX(FUEL,FALSDT, 1.000000E-01)

 RELAX(CRAD,FALSDT, 1.000000E+01)

 RELAX(MIXF,FALSDT, 1.000000E-01)

 KELIN   =       3

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

  Group 18. Limits

 VARMAX(U1  ) = 1.000000E+06 ;VARMIN(U1  ) =-1.000000E+06

 VARMAX(V1  ) = 1.000000E+06 ;VARMIN(V1  ) =-1.000000E+06

 VARMAX(W1  ) = 1.000000E+06 ;VARMIN(W1  ) =-1.000000E+06

 VARMAX(H1  ) = 1.000000E+10 ;VARMIN(H1  ) = 3.000000E+05

 VARMAX(PR5 ) = 1.000000E+00 ;VARMIN(PR5 ) = 0.000000E+00

 VARMAX(PR4 ) = 1.000000E+00 ;VARMIN(PR4 ) = 0.000000E+00

 VARMAX(PR3 ) = 1.000000E+00 ;VARMIN(PR3 ) = 0.000000E+00

 VARMAX(PR2 ) = 1.000000E+00 ;VARMIN(PR2 ) = 0.000000E+00

 VARMAX(PR1 ) = 1.000000E+00 ;VARMIN(PR1 ) = 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(FU4 ) = 1.000000E+00 ;VARMIN(FU4 ) = 0.000000E+00

 VARMAX(FU3 ) = 1.000000E+00 ;VARMIN(FU3 ) = 0.000000E+00

 VARMAX(FU2 ) = 1.000000E+00 ;VARMIN(FU2 ) = 0.000000E+00

 VARMAX(FU1 ) = 1.000000E+00 ;VARMIN(FU1 ) = 0.000000E+00

 VARMAX(HF5 ) = 1.000000E+10 ;VARMIN(HF5 ) = 3.000000E+05

 VARMAX(HF4 ) = 1.000000E+10 ;VARMIN(HF4 ) = 3.000000E+05

 VARMAX(HF3 ) = 1.000000E+10 ;VARMIN(HF3 ) = 3.000000E+05

 VARMAX(HF2 ) = 1.000000E+10 ;VARMIN(HF2 ) = 3.000000E+05

 VARMAX(HF1 ) = 1.000000E+10 ;VARMIN(HF1 ) = 3.000000E+05

 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(DEN1) = 2.000000E+00 ;VARMIN(DEN1) = 1.000000E-03

 VARMAX(FUEL) = 1.000000E+00 ;VARMIN(FUEL) = 0.000000E+00

 VARMAX(TMP1) = 1.000000E+04 ;VARMIN(TMP1) = 3.000000E+02

 VARMAX(PROD) = 1.000000E+00 ;VARMIN(PROD) = 0.000000E+00

 VARMAX(OXID) = 1.000000E+00 ;VARMIN(OXID) = 0.000000E+00

 VARMAX(MIXF) = 1.000000E+00 ;VARMIN(MIXF) = 0.000000E+00

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

  Group 19. EARTH Calls To GROUND Station

 USEGRD  =    T  ;USEGRX =    T

 NAMSAT  =MOSG

 GENK    =    T

 ASAP    =    T

 CHSOA   = 5.000000E-01 ;CHSOB  = 4.000000E+00

 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(V1  ,Y,Y,Y,N,Y,Y)

 OUTPUT(W1  ,Y,Y,Y,Y,Y,Y)

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

  Group 22. Monitor Print-Out

 IXMON   =       3 ;IYMON  =       5 ;IZMON  =      20

 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.953568E-01,6.788709E-01,-4.297495E-01)

 GVIEW(UP,6.290849E-01,7.265853E-01,2.762717E-01)



> DOM,    SIZE,        1.000000E+00, 1.000000E+00, 1.000000E+01

> DOM,    MONIT,       5.000000E-01, 9.000000E-01, 9.750000E+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,        FUELIN

> OBJ,    POSITION,    6.000000E-01, 0.000000E+00, 0.000000E+00

> OBJ,    SIZE,        4.000000E-01, 0.000000E+00, 5.000000E-01

> OBJ,    CLIPART,     cube3t

> OBJ,    ROTATION24,        1

> OBJ,    TYPE,        INLET

> OBJ,    PRESSURE,      0.000000E+00

> OBJ,    VELOCITY,      0.000000E+00, 2.000000E+00, 0.000000E+00

> OBJ,    TEMPERATURE,   3.000000E+02

> OBJ,    INLET_FU5,     1.000000E+00

> OBJ,    INLET_FU4,     1.000000E+00

> OBJ,    INLET_FU3,     1.000000E+00

> OBJ,    INLET_FU2,     1.000000E+00

> OBJ,    INLET_FU1,     1.000000E+00

> OBJ,    INLET_HF5,     7.900000E+06

> OBJ,    INLET_HF4,     7.900000E+06

> OBJ,    INLET_HF3,     7.900000E+06

> OBJ,    INLET_HF2,     7.900000E+06

> OBJ,    INLET_HF1,     7.900000E+06

> OBJ,    INLET_F5,      1.000000E+00

> OBJ,    INLET_FUEL,    1.000000E+00

> OBJ,    INLET_MIXF,    1.000000E+00

> OBJ,    TURB-INTENS,   5.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+03

> OBJ,    TURBULENCE,   -1.026000E+04,-1.026000E+04



> OBJ,    NAME,        AIRIN1

> OBJ,    POSITION,    6.000000E-01, 0.000000E+00, 5.000000E-01

> OBJ,    SIZE,        4.000000E-01, 0.000000E+00, 5.000000E-01

> OBJ,    CLIPART,     cube10

> OBJ,    ROTATION24,        1

> OBJ,    TYPE,        INLET

> OBJ,    PRESSURE,      0.000000E+00

> OBJ,    VELOCITY,      0.000000E+00, 2.000000E+00, 0.000000E+00

> OBJ,    TEMPERATURE,   3.000000E+02

> OBJ,    INLET_OX5,     1.000000E+00

> OBJ,    INLET_OX4,     1.000000E+00

> OBJ,    INLET_OX3,     1.000000E+00

> OBJ,    INLET_OX2,     1.000000E+00

> OBJ,    INLET_OX1,     1.000000E+00

> OBJ,    INLET_HF5,     4.500000E+05

> OBJ,    INLET_HF4,     4.500000E+05

> OBJ,    INLET_HF3,     4.500000E+05

> OBJ,    INLET_HF2,     4.500000E+05

> OBJ,    INLET_HF1,     4.500000E+05

> OBJ,    INLET_F1,      1.000000E+00

> OBJ,    TURB-INTENS,   5.000000E+00



> OBJ,    NAME,        AIRIN2

> OBJ,    POSITION,    4.000000E-01, 0.000000E+00, 0.000000E+00

> OBJ,    SIZE,        2.000000E-01, 0.000000E+00, 5.000000E-01

> OBJ,    CLIPART,     cube10

> OBJ,    ROTATION24,        1

> OBJ,    TYPE,        INLET

> OBJ,    PRESSURE,      0.000000E+00

> OBJ,    VELOCITY,      0.000000E+00, 2.000000E+00, 0.000000E+00

> OBJ,    TEMPERATURE,   3.000000E+02

> OBJ,    INLET_OX5,     1.000000E+00

> OBJ,    INLET_OX4,     1.000000E+00

> OBJ,    INLET_OX3,     1.000000E+00

> OBJ,    INLET_OX2,     1.000000E+00

> OBJ,    INLET_OX1,     1.000000E+00

> OBJ,    INLET_HF5,     4.500000E+05

> OBJ,    INLET_HF4,     4.500000E+05

> OBJ,    INLET_HF3,     4.500000E+05

> OBJ,    INLET_HF2,     4.500000E+05

> OBJ,    INLET_HF1,     4.500000E+05

> OBJ,    INLET_F1,      1.000000E+00

> OBJ,    TURB-INTENS,   5.000000E+00

STOP