Encyclopaedia Index


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

 

 

   ----------------------------------------------------------  

   PHOENICS May 2012 - EARTH                                   

   (C) Copyright 2011                                          

   Concentration Heat and Momentum Limited                     

   All rights reserved.                                        

   Address:  Bakery House, 40 High Street                      

   Wimbledon, London, SW19 5AU, England                        

   tel:       +44 (0) 20-8947-7651                             

   fax:       +44 (0) 20-8879-3497                             

   e-mail:    phoenics@cham.co.uk                              

   web:       www.cham.co.uk                                   

   ----------------------------------------------------------  

   Code expiry date is the end  : oct 2012                     

   ----------------------------------------------------------  

   Running with 64-bit Single Precision executable             

 

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

 

 Initial estimated storage requirement is                  10000000

 Information about material properties

 

 Total number of SPEDATs is                    13

 number of materials specified by SPEDATs is                      1

 solprp = 100 porprp = 198 vacprp = 199

 !!!! The properties file is PROPS

 Properties being read from  PROPS

 Properties have been read from  PROPS

 Formula used for setting RHO1

 Formula used for setting TMP1

 

 PRPS is stored with initial value=

 =-1.

 Material properties used for phase 1 are

 density

 laminar viscosity

 turbulent viscosity

 mixing length

 specific heat

 >>>   End of property-related data   <<<

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

 

 

 >>> Data read in from the Q1 file    <<<

 

   The following text is provided as an example of what, by use of   

   the readq1 command, can be transmitted to EARTH for writing near  

   the top of the RESULT file. Such text should leave columns 1 and 2

   blank; and it should not extend beyond the 68th column.           

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

     * This library case dates from the earliest days of PHOENICS,  *

     * when Professor WU Chung-Hua ('turbomachinery Wu', who had    *

     * returned to China from the USA) visited CHAM in 1982.        *

     *                                                              *

     * The configuration of the combustion chamber, and its being   *

     * supplied with premixed fuel vapour and air, was proposed by  *

     * Professor Wu's accompanying assistant.                       *

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

 

 >>> End of data read in from Q1      <<<

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

 

 Number of F-array locations available is                  10000000

 Number used before BFC allowance is                          68628

 Number used after BFC allowance is                           68628

 

 >> Current turbulence model constants <<

 They may be changed by inserting in Q1

 SPEDAT(KECONST,name of constant,R,value)

 

 CMU =0.5478 CD =0.1643

 CMUCD =0.090004

 C1E =1.44 C2E =1.92

 AK =0.41 EWAL =8.6

 

 biggest cell volume divided by average is     2.469999      at:

 ix = 1 iy = 10 iz = 7

 xg =0.05236 yg =0.06175

 zg =0.08575

 smallest cell volume divided by average is   2.1224484E-02  at:

 ix = 1 iy = 1 iz = 3

 xg =0.05236 yg =3.25E-03

 zg =0.027

 ratio of smallest to biggest is   8.5929111E-03

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

  Group 1. Run Title and Number

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

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

 

 TEXT(Idealised Gas-Turbine Combustion Chamber)

 

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

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

 

 IRUNN = 1 ;LIBREF = 492

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

  Group 2. Time dependence

 STEADY = T

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

  Group 3. X-Direction Grid Spacing

 CARTES = F

 NX = 6

 XULAST =0.628318

   Method of pairs used for grid setting.

 XFRAC(1)=-6. ;XFRAC(2)=0.166666

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

  Group 4. Y-Direction Grid Spacing

 NY = 10

 YVLAST =1.

   Method of pairs used for grid setting.

 YFRAC(1)=-5. ;YFRAC(3)=1.

 YFRAC(5)=1. ;YFRAC(7)=3.

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

  Group 5. Z-Direction Grid Spacing

 PARAB = F

 NZ = 13

 ZWLAST =1.

   Method of pairs used for grid setting.

 ZFRAC(1)=-2. ;ZFRAC(3)=1.

 ZFRAC(5)=1. ;ZFRAC(7)=1.

 ZFRAC(9)=1. ;ZFRAC(11)=7.

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

  Group 6. Body-Fitted Coordinates

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

  Group 7. Variables: STOREd,SOLVEd,NAMEd

 ONEPHS = T

 NAME(1)=P1 ;NAME(3)=U1                                                         

 NAME(5)=V1 ;NAME(7)=W1                                                         

 NAME(12)=KE ;NAME(13)=EP                                                       

 NAME(14)=H1 ;NAME(136)=PRPS                                                    

 NAME(137)=DEGF ;NAME(138)=RRAT                                                 

 NAME(139)=EPKE ;NAME(140)=HPOR                                                 

 NAME(141)=NPOR ;NAME(142)=EPOR                                                 

 NAME(143)=VPOR ;NAME(144)=ENUT                                                 

 NAME(145)=RHO1 ;NAME(146)=TMP1                                                 

 NAME(147)=PROD ;NAME(148)=OXID                                                 

 NAME(149)=FUEL ;NAME(150)=MIXF                                                 

    * Y in SOLUTN argument list denotes:

    * 1-stored 2-solved 3-whole-field

    * 4-point-by-point 5-explicit 6-harmonic averaging 

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

 SOLUTN(U1,Y,Y,N,N,N,Y)

 SOLUTN(V1,Y,Y,N,N,N,Y)

 SOLUTN(W1,Y,Y,N,N,N,Y)

 SOLUTN(KE,Y,Y,N,N,N,N)

 SOLUTN(EP,Y,Y,N,N,N,N)

 SOLUTN(H1,Y,Y,N,N,N,Y)

 SOLUTN(PRPS,Y,N,N,N,N,N)

 SOLUTN(DEGF,Y,N,N,N,N,N)

 SOLUTN(RRAT,Y,N,N,N,N,N)

 SOLUTN(EPKE,Y,N,N,N,N,N)

 SOLUTN(HPOR,Y,N,N,N,N,N)

 SOLUTN(NPOR,Y,N,N,N,N,N)

 SOLUTN(EPOR,Y,N,N,N,N,N)

 SOLUTN(VPOR,Y,N,N,N,N,N)

 SOLUTN(ENUT,Y,N,N,N,N,Y)

 SOLUTN(RHO1,Y,N,N,N,N,Y)

 SOLUTN(TMP1,Y,N,N,N,N,Y)

 SOLUTN(PROD,Y,N,N,N,N,Y)

 SOLUTN(OXID,Y,N,N,N,N,Y)

 SOLUTN(FUEL,Y,Y,N,N,N,Y)

 SOLUTN(MIXF,Y,Y,N,N,N,Y)

 DEN1 = 145

 VIST = 144

 EPOR = 142 ;HPOR = 140 ;NPOR = 141 ;VPOR = 143

 TEMP1 = 146

 PRPS = 136

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

  Group 8. Terms & Devices

    * Y in TERMS argument list denotes:

    * 1-built-in source 2-convection 3-diffusion 4-transient

    * 5-first phase variable 6-interphase transport         

 TERMS(P1,Y,Y,Y,N,Y,N)

 TERMS(U1,Y,Y,Y,N,Y,N)

 TERMS(V1,Y,Y,Y,N,Y,N)

 TERMS(W1,Y,Y,Y,N,Y,N)

 TERMS(KE,N,Y,Y,N,Y,N)

 TERMS(EP,N,Y,Y,N,Y,N)

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

 TERMS(FUEL,N,Y,Y,N,Y,N)

 TERMS(MIXF,N,Y,Y,N,Y,N)

 DIFCUT =0.5 ;ZDIFAC =1.

 GALA = F ;ADDDIF = F

 NEWRH1 = T

 NEWENT = T

 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1

 DENPCO = T

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

  Group 9. Properties used if PRPS is not

  stored, and where PRPS = -1.0 if it is!

 RHO1 = GRND6 ;TMP1 = GRND8

 EL1 = GRND4

 TSURR =0. ;TEMP0 =0.

 PRESS0 =8.0E+05

 DVO1DT =0. ;DRH1DP =0.

 RHO1A =16. ;RHO1B =29.

 RHO1C =28.

 TMP1A =0. ;TMP1B =0.

 TMP1C =0.

 TMP2A =0.054825 ;TMP2B =4.9E+07

 TMP2C =0.

 EMISS =0. ;SCATT =0.

 RADIA =0. ;RADIB =0.

 EL1A =0. ;EL1B =0.

 EL1C =0.

 ENUL =9.983361E-06 ;ENUT = GRND3

 ENUTA =0. ;ENUTB =0.

 ENUTC =0.

 IENUTA = 0

 PRNDTL(U1)=1. ;PRNDTL(V1)=1.

 PRNDTL(W1)=1. ;PRNDTL(KE)=1.

 PRNDTL(EP)=1. ;PRNDTL(H1)=1.

 PRNDTL(FUEL)=1. ;PRNDTL(MIXF)=1.

 PRT(U1)=1. ;PRT(V1)=1.

 PRT(W1)=1. ;PRT(KE)=1.

 PRT(EP)=1.314 ;PRT(H1)=1.

 PRT(FUEL)=1. ;PRT(MIXF)=1.

 CP1 = GRND10 ;CP2 =1.

 CP1A =1500. ;CP1B =1500.

 CP1C =1500. ;CP1D =0.

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

  Group 10.Inter-Phase Transfer Processes

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

  Group 11.Initial field variables (PHIs)

 FIINIT(P1)=1.0E-10 ;FIINIT(U1)=1.0E-10

 FIINIT(V1)=1.0E-10 ;FIINIT(W1)=10.

 FIINIT(KE)=11.25 ;FIINIT(EP)=353.255707

 FIINIT(H1)=1.1595E+06 ;FIINIT(PRPS)=-1.

 FIINIT(DEGF)=1.0E-10 ;FIINIT(RRAT)=1.0E-10

 FIINIT(EPKE)=1.0E-10 ;FIINIT(HPOR)=1.

 FIINIT(NPOR)=1. ;FIINIT(EPOR)=1.

 FIINIT(VPOR)=1. ;FIINIT(ENUT)=1.0E-10

 FIINIT(RHO1)=3.606 ;FIINIT(TMP1)=773.

 FIINIT(PROD)=1.0E-10 ;FIINIT(OXID)=1.0E-10

 FIINIT(FUEL)=0.054825 ;FIINIT(MIXF)=0.054825

 

 PATCH(CMP1 ,INIVAL, 1, 6, 1, 2, 1, 3, 1, 1)

 INIT(CMP1 ,PRPS,0. ,199. )

 

 PATCH(CMP2 ,INIVAL, 1, 5, 7, 7, 1, 1, 1, 1)

 INIT(CMP2 ,EPOR,0. ,0.709 )

 

 PATCH(CMP3 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1)

 INIT(CMP3 ,VPOR,0. ,0.709 )

 

 PATCH(CMP4 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1)

 INIT(CMP4 ,NPOR,0. ,0.418 )

 

 PATCH(CMP5 ,INIVAL, 1, 6, 7, 7, 1, 1, 1, 1)

 INIT(CMP5 ,HPOR,0. ,1. )

 

 PATCH(CMP6 ,INIVAL, 1, 5, 8, 8, 1, 1, 1, 1)

 INIT(CMP6 ,EPOR,0. ,0.185 )

 

 PATCH(CMP7 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1)

 INIT(CMP7 ,VPOR,0. ,0.185 )

 

 PATCH(CMP8 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1)

 INIT(CMP8 ,V1 ,0. ,0. )

 INIT(CMP8 ,NPOR,0. ,0. )

 

 PATCH(CMP9 ,INIVAL, 1, 6, 8, 8, 1, 1, 1, 1)

 INIT(CMP9 ,HPOR,0. ,0.877 )

 

 PATCH(CMP10 ,INIVAL, 1, 5, 9, 10, 1, 1, 1, 1)

 INIT(CMP10 ,U1 ,0. ,0. )

 INIT(CMP10 ,EPOR,0. ,0. )

 

 PATCH(CMP11 ,INIVAL, 1, 6, 9, 10, 1, 1, 1, 1)

 INIT(CMP11 ,PRPS,0. ,199. )

 

 PATCH(CMP12 ,INIVAL, 1, 6, 9, 9, 1, 1, 1, 1)

 INIT(CMP12 ,V1 ,0. ,0. )

 INIT(CMP12 ,NPOR,0. ,0. )

 

 PATCH(CMP13 ,INIVAL, 1, 6, 9, 10, 1, 1, 1, 1)

 INIT(CMP13 ,W1 ,0. ,0. )

 INIT(CMP13 ,HPOR,0. ,0. )

 

 PATCH(CMP14 ,INIVAL, 1, 5, 8, 8, 2, 2, 1, 1)

 INIT(CMP14 ,EPOR,0. ,0.997 )

 

 PATCH(CMP15 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1)

 INIT(CMP15 ,VPOR,0. ,0.997 )

 

 PATCH(CMP16 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1)

 INIT(CMP16 ,NPOR,0. ,0.946 )

 

 PATCH(CMP17 ,INIVAL, 1, 6, 8, 8, 2, 2, 1, 1)

 INIT(CMP17 ,HPOR,0. ,1. )

 

 PATCH(CMP18 ,INIVAL, 1, 5, 9, 9, 2, 2, 1, 1)

 INIT(CMP18 ,EPOR,0. ,0.706 )

 

 PATCH(CMP19 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1)

 INIT(CMP19 ,VPOR,0. ,0.706 )

 

 PATCH(CMP20 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1)

 INIT(CMP20 ,NPOR,0. ,0.473 )

 

 PATCH(CMP21 ,INIVAL, 1, 6, 9, 9, 2, 2, 1, 1)

 INIT(CMP21 ,HPOR,0. ,1. )

 

 PATCH(CMP22 ,INIVAL, 1, 5, 10, 10, 2, 2, 1, 1)

 INIT(CMP22 ,EPOR,0. ,0.236 )

 

 PATCH(CMP23 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1)

 INIT(CMP23 ,VPOR,0. ,0.236 )

 

 PATCH(CMP24 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1)

 

 PATCH(CMP25 ,INIVAL, 1, 6, 10, 10, 2, 2, 1, 1)

 INIT(CMP25 ,HPOR,0. ,1. )

 

 PATCH(WALL4 ,INIVAL, 1, 6, 2, 2, 1, 3, 1, 1)

 INIT(WALL4 ,V1 ,0. ,0. )

 INIT(WALL4 ,NPOR,0. ,0. )

 

 PATCH(WALL6 ,INIVAL, 1, 6, 1, 2, 3, 3, 1, 1)

 INIT(WALL6 ,W1 ,0. ,0. )

 INIT(WALL6 ,HPOR,0. ,0. )

 INIADD = F

 FSWEEP = 1

 NAMFI =CHAM

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

  Group 12. Patchwise adjustment of terms

  Patches for this group are printed with those

  for Group 13.

  Their names begin either with GP12 or &

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

  Group 13. Boundary & Special Sources

 

 PATCH(FOIN ,SOUTH , 1, 6, 3, 3, 3, 3, 1, 1)

 COVAL(FOIN ,P1 , FIXFLU ,360.600006 )

 COVAL(FOIN ,U1 ,0. ,0. )

 COVAL(FOIN ,V1 ,0. ,100. )

 COVAL(FOIN ,W1 ,0. ,0. )

 COVAL(FOIN ,KE ,0. ,0.25 )

 COVAL(FOIN ,EP ,0. ,58.752609 )

 COVAL(FOIN ,H1 ,0. ,3.845904E+06 )

 COVAL(FOIN ,FUEL,0. ,0.054825 )

 COVAL(FOIN ,MIXF,0. ,0.054825 )

 

 PATCH(SOIN1 ,CELL , 1, 1, 10, 10, 7, 7, 1, 1)

 COVAL(SOIN1 ,P1 , FIXFLU ,8.157E-03 )

 COVAL(SOIN1 ,U1 ,0. ,0. )

 COVAL(SOIN1 ,V1 ,0. ,-40. )

 COVAL(SOIN1 ,W1 ,0. ,0. )

 COVAL(SOIN1 ,KE ,0. ,0.123 )

 COVAL(SOIN1 ,EP ,0. ,8.429416 )

 COVAL(SOIN1 ,H1 ,0. ,1.1595E+06 )

 COVAL(SOIN1 ,FUEL,0. ,0. )

 COVAL(SOIN1 ,MIXF,0. ,0. )

 

 PATCH(SOIN2 ,CELL , 4, 4, 10, 10, 7, 7, 1, 1)

 COVAL(SOIN2 ,P1 , FIXFLU ,8.157E-03 )

 COVAL(SOIN2 ,U1 ,0. ,0. )

 COVAL(SOIN2 ,V1 ,0. ,-40. )

 COVAL(SOIN2 ,W1 ,0. ,0. )

 COVAL(SOIN2 ,KE ,0. ,0.123 )

 COVAL(SOIN2 ,EP ,0. ,8.429416 )

 COVAL(SOIN2 ,H1 ,0. ,1.1595E+06 )

 COVAL(SOIN2 ,FUEL,0. ,0. )

 COVAL(SOIN2 ,MIXF,0. ,0. )

 

 PATCH(DILUIN ,CELL , 1, 1, 10, 10, 10, 10, 1, 1)

 COVAL(DILUIN ,P1 , FIXFLU ,0.0145 )

 COVAL(DILUIN ,U1 ,0. ,0. )

 COVAL(DILUIN ,V1 ,0. ,-40. )

 COVAL(DILUIN ,W1 ,0. ,0. )

 COVAL(DILUIN ,KE ,0. ,0.123 )

 COVAL(DILUIN ,EP ,0. ,6.299674 )

 COVAL(DILUIN ,H1 ,0. ,1.1595E+06 )

 COVAL(DILUIN ,FUEL,0. ,0. )

 COVAL(DILUIN ,MIXF,0. ,0. )

 

 PATCH(OUTLET ,HIGH , 1, 6, 1, 10, 13, 13, 1, 1)

 COVAL(OUTLET ,P1 ,1000. ,0. )

 COVAL(OUTLET ,U1 ,0. ,0. )

 COVAL(OUTLET ,V1 ,0. ,0. )

 COVAL(OUTLET ,W1 ,0. ,0. )

 COVAL(OUTLET ,KE ,0. , SAME )

 COVAL(OUTLET ,EP ,0. , SAME )

 COVAL(OUTLET ,H1 ,0. , SAME )

 COVAL(OUTLET ,FUEL,0. , SAME )

 COVAL(OUTLET ,MIXF,0. , SAME )

 

 PATCH(WALL1 ,NWALL , 1, 6, 10, 10, 3, 6, 1, 1)

 COVAL(WALL1 ,U1 , GRND2 ,0. )

 COVAL(WALL1 ,W1 , GRND2 ,0. )

 COVAL(WALL1 ,KE , GRND2 , GRND2 )

 COVAL(WALL1 ,EP , GRND2 , GRND2 )

 

 PATCH(WALL2 ,NWALL , 1, 6, 10, 10, 8, 9, 1, 1)

 COVAL(WALL2 ,U1 , GRND2 ,0. )

 COVAL(WALL2 ,W1 , GRND2 ,0. )

 COVAL(WALL2 ,KE , GRND2 , GRND2 )

 COVAL(WALL2 ,EP , GRND2 , GRND2 )

 

 PATCH(WALL3 ,NWALL , 1, 6, 10, 10, 11, 13, 1, 1)

 COVAL(WALL3 ,U1 , GRND2 ,0. )

 COVAL(WALL3 ,W1 , GRND2 ,0. )

 COVAL(WALL3 ,KE , GRND2 , GRND2 )

 COVAL(WALL3 ,EP , GRND2 , GRND2 )

 

 PATCH(WALL4-NW,NWALL , 1, 6, 2, 2, 1, 3, 1, 1)

 COVAL(WALL4-NW,U1 , GRND2 ,0. )

 COVAL(WALL4-NW,W1 , GRND2 ,0. )

 COVAL(WALL4-NW,KE , GRND2 , GRND2 )

 COVAL(WALL4-NW,EP , GRND2 , GRND2 )

 

 PATCH(WALL4-SW,SWALL , 1, 6, 3, 3, 1, 3, 1, 1)

 COVAL(WALL4-SW,U1 , GRND2 ,0. )

 COVAL(WALL4-SW,W1 , GRND2 ,0. )

 COVAL(WALL4-SW,KE , GRND2 , GRND2 )

 COVAL(WALL4-SW,EP , GRND2 , GRND2 )

 

 PATCH(WALL5 ,LWALL , 1, 6, 3, 6, 1, 1, 1, 1)

 COVAL(WALL5 ,U1 , GRND2 ,0. )

 COVAL(WALL5 ,V1 , GRND2 ,0. )

 COVAL(WALL5 ,KE , GRND2 , GRND2 )

 COVAL(WALL5 ,EP , GRND2 , GRND2 )

 

 PATCH(WALL6-HW,HWALL , 1, 6, 1, 2, 3, 3, 1, 1)

 COVAL(WALL6-HW,U1 , GRND2 ,0. )

 COVAL(WALL6-HW,V1 , GRND2 ,0. )

 COVAL(WALL6-HW,KE , GRND2 , GRND2 )

 COVAL(WALL6-HW,EP , GRND2 , GRND2 )

 

 PATCH(WALL6-LW,LWALL , 1, 6, 1, 2, 4, 4, 1, 1)

 COVAL(WALL6-LW,U1 , GRND2 ,0. )

 COVAL(WALL6-LW,V1 , GRND2 ,0. )

 COVAL(WALL6-LW,KE , GRND2 , GRND2 )

 COVAL(WALL6-LW,EP , GRND2 , GRND2 )

 

 PATCH(KESOURCE,PHASEM, 1, 6, 1, 10, 1, 13, 1, 1)

 COVAL(KESOURCE,KE , GRND4 , GRND4 )

 COVAL(KESOURCE,EP , GRND4 , GRND4 )

 

 PATCH(WHOLE ,PHASEM, 1, 6, 1, 10, 1, 13, 1, 1)

 COVAL(WHOLE ,FUEL,In-Form:COVAL(1.*EPKE&0.0))

 XCYCLE = F

 EGWF = T

 WALLCO = GRND2

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

  Group 14. Downstream Pressure For PARAB

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

  Group 15. Terminate Sweeps

 LSWEEP = 100 ;ISWC1 = 1

 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1

 SARAH =1.

 SELREF = T

 RESFAC =1.0E-03

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

  Group 16. Terminate Iterations

 LITER(P1)=20 ;LITER(U1)=20

 LITER(V1)=20 ;LITER(W1)=20

 LITER(KE)=20 ;LITER(EP)=20

 LITER(H1)=20 ;LITER(FUEL)=20

 LITER(MIXF)=20

 ENDIT(P1)=1.0E-03 ;ENDIT(U1)=1.0E-03

 ENDIT(V1)=1.0E-03 ;ENDIT(W1)=1.0E-03

 ENDIT(KE)=1.0E-03 ;ENDIT(EP)=1.0E-03

 ENDIT(H1)=1.0E-03 ;ENDIT(FUEL)=1.0E-03

 ENDIT(MIXF)=1.0E-03

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

  Group 17. Relaxation

 RELAX(P1,LINRLX,1.)                                                            

 RELAX(U1,FALSDT,0.02)                                                          

 RELAX(V1,FALSDT,0.02)                                                          

 RELAX(W1,FALSDT,0.02)                                                          

 RELAX(KE,LINRLX,0.5)                                                           

 RELAX(EP,LINRLX,0.5)                                                           

 RELAX(H1,FALSDT,10.)                                                           

 RELAX(PRPS,LINRLX,1.)                                                          

 RELAX(DEGF,LINRLX,1.)                                                          

 RELAX(RRAT,LINRLX,1.)                                                          

 RELAX(EPKE,LINRLX,1.)                                                          

 RELAX(ENUT,LINRLX,1.)                                                          

 RELAX(RHO1,LINRLX,1.)                                                          

 RELAX(TMP1,LINRLX,1.)                                                          

 RELAX(PROD,LINRLX,1.)                                                          

 RELAX(OXID,LINRLX,1.)                                                          

 RELAX(FUEL,FALSDT,1.0E+09)                                                     

 RELAX(MIXF,FALSDT,1.0E+09)                                                     

 KELIN = 3

 OVRRLX =0.

 EXPERT = F ;NNORSL = F

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

  Group 18. Limits

 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10                                        

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

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

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

 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10                                         

 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10                                         

 VARMAX(H1)=1.0E+10 ;VARMIN(H1)=-1.0E+10                                        

 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10                                    

 VARMAX(DEGF)=1.0E+10 ;VARMIN(DEGF)=-1.0E+10                                    

 VARMAX(RRAT)=1.0E+10 ;VARMIN(RRAT)=-1.0E+10                                    

 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=0.                                          

 VARMAX(HPOR)=1.0E+10 ;VARMIN(HPOR)=-1.0E+10                                    

 VARMAX(NPOR)=1.0E+10 ;VARMIN(NPOR)=-1.0E+10                                    

 VARMAX(EPOR)=1.0E+10 ;VARMIN(EPOR)=-1.0E+10                                    

 VARMAX(VPOR)=1.0E+10 ;VARMIN(VPOR)=-1.0E+10                                    

 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10                                    

 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=-1.0E+10                                    

 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=773.                                        

 VARMAX(PROD)=1.0E+10 ;VARMIN(PROD)=0.                                          

 VARMAX(OXID)=1.0E+10 ;VARMIN(OXID)=0.                                          

 VARMAX(FUEL)=1.0E+10 ;VARMIN(FUEL)=0.                                          

 VARMAX(MIXF)=1.0E+10 ;VARMIN(MIXF)=-1.0E+10                                    

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

  Group 19. Data transmitted to GROUND

 GENK = T

 READQ1 = T

 PARSOL = F

 GEN1 = 681

 CHSOA =0.054825 ;CHSOB =1.

 CSG3    ='CNGR'

 SPEDAT(SET,PROPERTY,TMP1,C,=MAX(773.&(H1-4.9E+07*FUEL)/(1500.*FUEL+1500.*OXID+15$)

 SPEDAT(SET,PROPERTY,TMP1,C,00.*PROD))!IMAT<100)

 SPEDAT(SET,PROPERTY,RHO1,C,=(P1+8.0E+05)*1.0/(FUEL/16.+OXID/29.+PROD/28.)/(TMP1*$)

 SPEDAT(SET,PROPERTY,RHO1,C,8313.4))

 SPEDAT(SET,STORED,PROD,C,=(MIXF-FUEL)*(1+17.24)!ZSLFIN)

 SPEDAT(SET,STORED,OXID,C,=1-FUEL-PROD!ZSLFIN)

 SPEDAT(SET,SOURCE,FUEL!WHOLE,C,=COVAL(1.*EPKE&0.0))

 SPEDAT(SET,STORED,RRAT,C,=1.*EPKE*FUEL!ZSLFIN)

 SPEDAT(SET,LONGNAME,TMP1,C,absolute_temperature_of_the_gas_Kelvin)

 SPEDAT(SET,STORED,DEGF,C,=(9./5.)*(TMP1-273)!ZSLFIN)

 SPEDAT(SET,LONGNAME,DEGF,C,degrees_Fahrenheit)

 SPEDAT(SET,LONGNAME,RRAT,C,rate_of_consumption_of_fuel_kg/m^3sec)

 SPEDAT(SET,MATERIAL,199,L,T)

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

  Group 20. Preliminary Printout

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

  Group 21. Print-out of Variables

 INIFLD = F ;SUBWGR = F

    * Y in OUTPUT argument list denotes:

    * 1-field 2-correction-eq. monitor 3-selective dumping      

    * 4-whole-field residual 5-spot-value table 6-residual table

 OUTPUT(P1,Y,N,Y,Y,Y,Y)

 OUTPUT(U1,Y,N,Y,Y,Y,Y)

 OUTPUT(V1,Y,N,Y,Y,Y,Y)

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

 OUTPUT(KE,Y,N,Y,Y,Y,Y)

 OUTPUT(EP,Y,N,Y,Y,Y,Y)

 OUTPUT(H1,Y,N,Y,Y,Y,Y)

 OUTPUT(PRPS,Y,N,Y,N,N,N)

 OUTPUT(DEGF,Y,N,Y,N,N,N)

 OUTPUT(RRAT,Y,N,Y,N,N,N)

 OUTPUT(EPKE,Y,N,Y,N,N,N)

 OUTPUT(HPOR,Y,N,Y,N,N,N)

 OUTPUT(NPOR,Y,N,Y,N,N,N)

 OUTPUT(EPOR,Y,N,Y,N,N,N)

 OUTPUT(VPOR,Y,N,Y,N,N,N)

 OUTPUT(ENUT,Y,N,Y,N,N,N)

 OUTPUT(RHO1,Y,N,Y,N,N,N)

 OUTPUT(TMP1,Y,N,Y,N,N,N)

 OUTPUT(PROD,Y,N,Y,N,N,N)

 OUTPUT(OXID,Y,N,Y,N,N,N)

 OUTPUT(FUEL,Y,N,Y,Y,Y,Y)

 OUTPUT(MIXF,Y,N,Y,Y,Y,Y)

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

  Group 22. Monitor Print-Out

 IXMON = 3 ;IYMON = 5 ;IZMON = 5

 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = 10005

 UWATCH = F ;USTEER = F

 HIGHLO = F

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

  Group 23.Field Print-Out & Plot Control

 NPRINT = 100 ;NUMCLS = 5

 NXPRIN = 1 ;IXPRF = 3 ;IXPRL = 3

 NYPRIN = 2 ;IYPRF = 1 ;IYPRL = 10

 NZPRIN = 2 ;IZPRF = 1 ;IZPRL = 10000

 XZPR = F ;YZPR = T

 IPLTF = 2 ;IPLTL = 100 ;NPLT = 5

 ISWPRF = 1 ;ISWPRL = 100000

 ITABL = 3 ;IPROF = 1

 ABSIZ =0.5 ;ORSIZ =0.4

 NTZPRF = 1 ;NCOLPF = 50

 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 40

 

 PATCH(IZ4 ,CONTUR, 1, 6, 1, 10, 4, 4, 1, 1)

 PLOT(IZ4 ,TMP1,1. ,10. )

 

 PATCH(XSECIN1 ,CONTUR, 1, 1, 1, 10, 1, 13, 1, 1)

 PLOT(XSECIN1 ,W1 ,0. ,10. )

 PLOT(XSECIN1 ,TMP1,0. ,10. )

 PLOT(XSECIN1 ,FUEL,0. ,10. )

 PLOT(XSECIN1 ,MIXF,0. ,10. )

 

 PATCH(XSECIN2 ,CONTUR, 4, 4, 1, 10, 1, 13, 1, 1)

 PLOT(XSECIN2 ,W1 ,0. ,10. )

 PLOT(XSECIN2 ,TMP1,0. ,10. )

 PLOT(XSECIN2 ,FUEL,0. ,10. )

 PLOT(XSECIN2 ,MIXF,0. ,10. )

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

  Group 24. Dumps For Restarts

 SAVE = T ;NOWIPE = F

 NSAVE =CHAM

 

 *** grid-geometry information ***

 X-coordinates of the cell centres

    2.618E-01

 Y-coordinates of the cell centres

    3.250E-03   1.625E-02   2.925E-02   4.150E-02   5.525E-02



 Z-coordinates of the cell centres

    6.250E-03   2.700E-02   4.500E-02   8.575E-02   1.347E-01

    1.837E-01   2.327E-01

 

 X-coordinates of the (higher) cell faces

    3.142E-01

 Y-coordinates of the (higher) cell faces

    6.500E-03   1.950E-02   3.250E-02   4.550E-02   5.850E-02



 Z-coordinates of the (higher) cell faces

    1.250E-02   2.900E-02   5.000E-02   9.800E-02   1.470E-01

    1.960E-01   2.450E-01

 

 Total number of F-array elements used is                     72255

 --- INTEGRATION OF EQUATIONS BEGINS ---

 

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

 *** Start of sweep 100

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

 

 Flow field at ITHYD=   1, ISWEEP=   100, ISTEP=   1

 YZPR IX= 3 isweep = 100 istep = 1

 

 Field Values of P1

 IY=  10    none       9.796E+02   3.570E+02  -1.081E+00   2.133E+02

 IY=   8   4.790E+02   8.482E+02   4.196E+01   3.236E+02   2.253E+02

 IY=   6   4.395E+01   2.847E+02  -1.122E+02   4.080E+02   2.852E+02

 IY=   4   4.858E+01   1.952E+02   1.183E+01   4.599E+02   3.657E+02

 IY=   2    none        none       1.909E+02   5.193E+02   4.919E+02

 IZ=          1           3           5           7           9

 IY=  10  -4.674E+01   3.576E-02

 IY=   8  -8.783E+01   1.937E-02

 IY=   6  -1.352E+01   1.986E-02

 IY=   4   1.880E+02   3.985E-02

 IY=   2   3.445E+02   2.588E-03

 IZ=         11          13

 

 Field Values of U1

 IY=  10    none       2.883E-02   6.426E-02   2.576E+00  -2.158E+00

 IY=   8  -1.359E-02  -1.897E-02  -1.380E-01  -2.574E+00   6.440E-01

 IY=   6  -1.124E-02  -3.861E-02  -3.383E-01   1.045E+00   8.897E-01

 IY=   4   3.789E-03  -2.515E-02  -2.512E-01   2.164E+00  -1.591E-02

 IY=   2    none        none      -1.170E-02   3.564E-01  -5.555E-02

 IZ=          1           3           5           7           9

 IY=  10  -7.426E+00  -4.195E+00

 IY=   8  -2.859E+00  -8.512E-01

 IY=   6   2.922E+00   1.959E+00

 IY=   4   3.441E+00   2.477E+00

 IY=   2   5.226E-01   4.560E+00

 IZ=         11          13

 

 Field Values of V1

 IY=   9    none       1.675E+01   2.216E+00   1.526E+00   5.192E-01

 IY=   7  -1.105E+01   5.495E+01   1.575E+00  -8.835E+00   2.755E+00

 IY=   5  -1.264E+01   8.061E+01   3.091E+00  -1.247E+01   4.359E+00

 IY=   3  -5.461E+00   9.071E+01   3.670E+00  -9.780E+00   4.472E+00

 IY=   1    none        none       1.475E+00  -3.379E+00   2.002E+00

 IZ=          1           3           5           7           9

 IY=   9   7.010E-02   9.604E-01

 IY=   7  -6.000E-01   5.482E-02

 IY=   5  -2.385E+00  -9.551E-01

 IY=   3  -1.684E+00   5.573E-01

 IY=   1  -3.233E-01   5.726E-02

 IZ=         11          13

 

 Field Values of W1

 IY=  10         none       1.359E+01   3.082E+01   2.373E+01   2.195E+01

 IY=   8       -9.175E+00   7.894E+00   1.214E+01   1.214E+01   1.512E+01

 IY=   6       -5.642E-01   9.216E-01  -7.048E+00   1.014E+01   1.153E+01

 IY=   4        1.371E+01  -2.930E+00  -1.373E+01   9.003E+00   6.591E+00

 IY=   2         none        none      -1.644E+01   9.173E+00  -1.837E+00

 IZ=               1           3           5           7           9

 IY=  10        1.861E+01

 IY=   8        1.468E+01

 IY=   6        1.224E+01

 IY=   4        1.405E+01

 IY=   2        1.833E+01

 IZ=              11

 

 Field Values of KE

 IY=  10    none       8.041E-01   6.824E+00   9.949E+01   3.717E+00

 IY=   8   7.849E+02   4.485E+02   2.648E+02   1.206E+02   4.274E+01

 IY=   6   1.173E+00   1.226E+02   1.122E+02   8.509E+01   4.912E+01

 IY=   4   5.841E-01   4.290E+01   8.091E+01   8.837E+01   6.460E+01

 IY=   2    none        none       1.016E+02   1.693E+02   9.427E+01

 IZ=          1           3           5           7           9

 IY=  10   3.999E+00   3.114E+00

 IY=   8   3.892E+01   2.582E+01

 IY=   6   4.773E+01   2.793E+01

 IY=   4   5.387E+01   4.317E+01

 IY=   2   7.239E+01   1.068E+02

 IZ=         11          13

 

 Field Values of EP

 IY=  10    none       8.891E+01   2.198E+03   6.335E+04   8.837E+02

 IY=   8   6.408E+05   3.924E+05   1.556E+05   6.444E+04   1.169E+04

 IY=   6   8.147E+01   1.043E+05   3.466E+04   3.655E+04   1.119E+04

 IY=   4   2.863E+01   4.210E+04   1.952E+04   3.355E+04   1.477E+04

 IY=   2    none        none       2.860E+04   6.474E+04   2.623E+04

 IZ=          1           3           5           7           9

 IY=  10   9.861E+02   6.775E+02

 IY=   8   1.073E+04   4.999E+03

 IY=   6   1.840E+04   6.436E+03

 IY=   4   2.150E+04   1.023E+04

 IY=   2   2.528E+04   3.394E+04

 IZ=         11          13

 

 Field Values of H1

 IY=  10    none       3.471E+06   3.296E+06   2.990E+06   2.799E+06

 IY=   8   3.496E+06   3.482E+06   2.977E+06   2.212E+06   2.150E+06

 IY=   6   3.502E+06   3.568E+06   2.669E+06   2.070E+06   2.018E+06

 IY=   4   3.502E+06   3.695E+06   2.287E+06   1.995E+06   1.974E+06

 IY=   2    none        none       2.087E+06   1.988E+06   1.916E+06

 IZ=          1           3           5           7           9

 IY=  10   2.437E+06   2.157E+06

 IY=   8   1.966E+06   1.910E+06

 IY=   6   1.771E+06   1.859E+06

 IY=   4   1.645E+06   1.699E+06

 IY=   2   1.555E+06   1.602E+06

 IZ=         11          13

 

 Field Values of PRPS

 IY=  10    empty       pil prop    pil prop    pil prop    pil prop

 IY=   8    pil prop    pil prop    pil prop    pil prop    pil prop

 IY=   6    pil prop    pil prop    pil prop    pil prop    pil prop

 IY=   4    pil prop    pil prop    pil prop    pil prop    pil prop

 IY=   2    empty       empty       pil prop    pil prop    pil prop

 IZ=          1           3           5           7           9

 IY=  10    pil prop    pil prop

 IY=   8    pil prop    pil prop

 IY=   6    pil prop    pil prop

 IY=   4    pil prop    pil prop

 IY=   2    pil prop    pil prop

 IZ=         11          13

 

 Field Values of DEGF: degrees_Fahrenheit

 IY=  10    none       2.923E+03   2.950E+03   2.864E+03   2.749E+03

 IY=   8   3.238E+03   2.424E+03   2.753E+03   2.078E+03   2.060E+03

 IY=   6   3.312E+03   1.884E+03   2.615E+03   1.937E+03   1.911E+03

 IY=   4   3.334E+03   1.478E+03   2.217E+03   1.865E+03   1.864E+03

 IY=   2    none        none       1.993E+03   1.870E+03   1.799E+03

 IZ=          1           3           5           7           9

 IY=  10   2.386E+03   2.079E+03

 IY=   8   1.856E+03   1.794E+03

 IY=   6   1.627E+03   1.735E+03

 IY=   4   1.479E+03   1.545E+03

 IY=   2   1.372E+03   1.430E+03

 IZ=         11          13

 

 Field Values of RRAT: rate_of_consumption_of_fuel_kg/m^3sec

 IY=  10    none       1.411E+00   2.814E+00   2.514E+00   4.799E-01

 IY=   8   6.465E+00   1.878E+01   3.276E+00   7.714E-01   1.318E-01

 IY=   6   4.711E-01   2.758E+01   5.099E-01   4.045E-01   7.326E-02

 IY=   4   3.136E-01   4.115E+01   1.469E-01   2.442E-01   5.323E-02

 IY=   2    none        none       9.347E-02   1.555E-01   4.510E-02

 IZ=          1           3           5           7           9

 IY=  10   1.989E-01   6.573E-02

 IY=   8   5.833E-02   2.025E-02

 IY=   6   4.608E-02   1.755E-02

 IY=   4   3.059E-02   1.044E-02

 IY=   2   1.471E-02   7.240E-03

 IZ=         11          13

 

 Field Values of EPKE

 IY=  10    none       1.105E+02   3.220E+02   6.367E+02   2.377E+02

 IY=   8   8.163E+02   8.749E+02   5.877E+02   5.343E+02   2.735E+02

 IY=   6   6.942E+01   8.507E+02   3.089E+02   4.296E+02   2.279E+02

 IY=   4   4.896E+01   9.815E+02   2.413E+02   3.796E+02   2.286E+02

 IY=   2    none        none       2.813E+02   3.824E+02   2.782E+02

 IZ=          1           3           5           7           9

 IY=  10   2.466E+02   2.175E+02

 IY=   8   2.757E+02   1.936E+02

 IY=   6   3.855E+02   2.305E+02

 IY=   4   3.991E+02   2.370E+02

 IY=   2   3.492E+02   3.170E+02

 IZ=         11          13

 

 Field Values of HPOR

 IY=  10   1.000E-11   1.0         1.0         1.0         1.0

 IY=   8   8.770E-01   1.0         1.0         1.0         1.0

 IY=   6   1.0         1.0         1.0         1.0         1.0

 IY=   4   1.0         1.0         1.0         1.0         1.0

 IY=   2   1.0         1.000E-11   1.0         1.0         1.0

 IZ=          1           3           5           7           9

 IY=  10   1.0         1.0

 IY=   8   1.0         1.0

 IY=   6   1.0         1.0

 IY=   4   1.0         1.0

 IY=   2   1.0         1.0

 IZ=         11          13

 

 Field Values of NPOR

 IY=  10   1.0         1.0         1.0         1.0         1.0

 IY=   8   1.000E-11   1.0         1.0         1.0         1.0

 IY=   6   1.0         1.0         1.0         1.0         1.0

 IY=   4   1.0         1.0         1.0         1.0         1.0

 IY=   2   1.000E-11   1.000E-11   1.0         1.0         1.0

 IZ=          1           3           5           7           9

 IY=  10   1.0         1.0

 IY=   8   1.0         1.0

 IY=   6   1.0         1.0

 IY=   4   1.0         1.0

 IY=   2   1.0         1.0

 IZ=         11          13

 

 Field Values of EPOR

 IY=  10   1.000E-11   1.0         1.0         1.0         1.0

 IY=   8   1.850E-01   1.0         1.0         1.0         1.0

 IY=   6   1.0         1.0         1.0         1.0         1.0

 IY=   4   1.0         1.0         1.0         1.0         1.0

 IY=   2   1.0         1.0         1.0         1.0         1.0

 IZ=          1           3           5           7           9

 IY=  10   1.0         1.0

 IY=   8   1.0         1.0

 IY=   6   1.0         1.0

 IY=   4   1.0         1.0

 IY=   2   1.0         1.0

 IZ=         11          13

 

 Field Values of VPOR

 IY=  10   1.000E-11   1.0         1.0         1.0         1.0

 IY=   8   1.850E-01   1.0         1.0         1.0         1.0

 IY=   6   1.0         1.0         1.0         1.0         1.0

 IY=   4   1.0         1.0         1.0         1.0         1.0

 IY=   2   1.000E-11   1.000E-11   1.0         1.0         1.0

 IZ=          1           3           5           7           9

 IY=  10   1.0         1.0

 IY=   8   1.0         1.0

 IY=   6   1.0         1.0

 IY=   4   1.0         1.0

 IY=   2   1.0         1.0

 IZ=         11          13

 

 Field Values of ENUT

 IY=  10    none       6.546E-04   1.907E-03   1.406E-02   1.407E-03

 IY=   8   8.653E-02   4.614E-02   4.055E-02   2.032E-02   1.407E-02

 IY=   6   1.520E-03   1.297E-02   3.270E-02   1.783E-02   1.940E-02

 IY=   4   1.073E-03   3.934E-03   3.018E-02   2.095E-02   2.544E-02

 IY=   2    none        none       3.252E-02   3.984E-02   3.049E-02

 IZ=          1           3           5           7           9

 IY=  10   1.460E-03   1.288E-03

 IY=   8   1.271E-02   1.201E-02

 IY=   6   1.114E-02   1.091E-02

 IY=   4   1.215E-02   1.639E-02

 IY=   2   1.866E-02   3.024E-02

 IZ=         11          13

 

 Field Values of RHO1

 IY=  10    none       1.426E+00   1.418E+00   1.460E+00   1.517E+00

 IY=   8   1.305E+00   1.667E+00   1.510E+00   1.929E+00   1.944E+00

 IY=   6   1.280E+00   2.039E+00   1.585E+00   2.044E+00   2.068E+00

 IY=   4   1.272E+00   2.453E+00   1.827E+00   2.110E+00   2.111E+00

 IY=   2    none        none       1.997E+00   2.105E+00   2.173E+00

 IZ=          1           3           5           7           9

 IY=  10   1.716E+00   1.929E+00

 IY=   8   2.117E+00   2.176E+00

 IY=   6   2.352E+00   2.236E+00

 IY=   4   2.534E+00   2.449E+00

 IY=   2   2.683E+00   2.600E+00

 IZ=         11          13

 

 Field Values of TMP1: absolute_temperature_of_the_gas_Kelvin

 IY=  10    none       1.897E+03   1.912E+03   1.864E+03   1.800E+03

 IY=   8   2.072E+03   1.620E+03   1.803E+03   1.427E+03   1.417E+03

 IY=   6   2.113E+03   1.320E+03   1.726E+03   1.349E+03   1.335E+03

 IY=   4   2.125E+03   1.094E+03   1.505E+03   1.309E+03   1.309E+03

 IY=   2    none        none       1.380E+03   1.312E+03   1.272E+03

 IZ=          1           3           5           7           9

 IY=  10   1.598E+03   1.428E+03

 IY=   8   1.304E+03   1.270E+03

 IY=   6   1.177E+03   1.237E+03

 IY=   4   1.094E+03   1.131E+03

 IY=   2   1.035E+03   1.067E+03

 IZ=         11          13

 

 Field Values of PROD

 IY=  10    none       6.274E-01   6.358E-01   6.094E-01   5.737E-01

 IY=   8   7.253E-01   4.728E-01   5.750E-01   3.654E-01   3.598E-01

 IY=   6   7.482E-01   3.054E-01   5.319E-01   3.217E-01   3.136E-01

 IY=   4   7.551E-01   1.792E-01   4.085E-01   2.993E-01   2.991E-01

 IY=   2    none        none       3.392E-01   3.009E-01   2.788E-01

 IZ=          1           3           5           7           9

 IY=  10   4.609E-01   3.656E-01

 IY=   8   2.965E-01   2.774E-01

 IY=   6   2.255E-01   2.589E-01

 IY=   4   1.795E-01   2.001E-01

 IY=   2   1.463E-01   1.643E-01

 IZ=         11          13

 

 Field Values of OXID

 IY=  10    none       3.598E-01   3.554E-01   3.867E-01   4.243E-01

 IY=   8   2.667E-01   5.057E-01   4.195E-01   6.332E-01   6.397E-01

 IY=   6   2.450E-01   6.622E-01   4.664E-01   6.773E-01   6.861E-01

 IY=   4   2.385E-01   7.789E-01   5.909E-01   7.001E-01   7.007E-01

 IY=   2    none        none       6.605E-01   6.987E-01   7.211E-01

 IZ=          1           3           5           7           9

 IY=  10   5.383E-01   6.341E-01

 IY=   8   7.033E-01   7.225E-01

 IY=   6   7.744E-01   7.410E-01

 IY=   4   8.204E-01   7.999E-01

 IY=   2   8.536E-01   8.357E-01

 IZ=         11          13

 

 Field Values of FUEL

 IY=  10    none       1.277E-02   8.737E-03   3.948E-03   2.019E-03

 IY=   8   7.919E-03   2.147E-02   5.574E-03   1.444E-03   4.819E-04

 IY=   6   6.786E-03   3.242E-02   1.651E-03   9.417E-04   3.215E-04

 IY=   4   6.406E-03   4.193E-02   6.086E-04   6.433E-04   2.329E-04

 IY=   2    none        none       3.323E-04   4.067E-04   1.621E-04

 IZ=          1           3           5           7           9

 IY=  10   8.065E-04   3.022E-04

 IY=   8   2.115E-04   1.046E-04

 IY=   6   1.195E-04   7.615E-05

 IY=   4   7.663E-05   4.403E-05

 IY=   2   4.214E-05   2.284E-05

 IZ=         11          13

 

 Field Values of MIXF

 IY=  10    none       4.717E-02   4.360E-02   3.736E-02   3.347E-02

 IY=   8   4.769E-02   4.739E-02   3.710E-02   2.147E-02   2.021E-02

 IY=   6   4.780E-02   4.916E-02   3.081E-02   1.858E-02   1.752E-02

 IY=   4   4.780E-02   5.175E-02   2.301E-02   1.705E-02   1.663E-02

 IY=   2    none        none       1.893E-02   1.690E-02   1.545E-02

 IZ=          1           3           5           7           9

 IY=  10   2.607E-02   2.035E-02

 IY=   8   1.646E-02   1.531E-02

 IY=   6   1.248E-02   1.427E-02

 IY=   4   9.918E-03   1.101E-02

 IY=   2   8.065E-03   9.030E-03

 IZ=         11          13

 

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

 Whole-field residuals before solution at sweep 100

 with resref values determined by EARTH

 & resfac=1.0E-03

  variable   resref  (res sum)/resref  (res sum)

     P1     1.841E-01   2.083E-03      3.833E-04

     U1     1.151E+01   1.041E-03      1.198E-02

     V1     1.151E+01   2.281E-04      2.625E-03

     W1     1.151E+01   1.638E-04      1.885E-03

     KE     5.023E+01   8.412E-04      4.225E-02

     EP     6.144E+04   6.554E-04      4.027E+01

     H1     3.501E+05   7.901E-05      2.766E+01

     FUEL   1.146E-03   1.613E-04      1.849E-07

     MIXF   3.519E-03   1.603E-04      5.642E-07

 

 SARAH = Self-Adjusting Relxtn. algoritHm

 DTFALS values set in EARTH by SARAH with

 relaxfac=1.

 

 U1 dtf=3.301611E-03

 V1 dtf=3.301611E-03

 W1 dtf=3.301611E-03

 

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

 Sources and sinks

 !! Zero nett sources are not printed !!!

 

    Nett Sources have units of mass_per_unit_time * variable

    Average values have units of the variable

 

    Typically the units of the sources are:

        U1,V1,W1 - Force  - Newtons

        R1       - Mass   - kg/s

        TEM1     - Energy - Watts

 

 Nett source of U1   at patch named: OUTLET   =-1.782140E-02

 Nett source of U1   at patch named: WALL1    =-6.723997E-05

 Nett source of U1   at patch named: WALL2    = 3.873256E-04

 Nett source of U1   at patch named: WALL3    = 1.171623E-03

 Nett source of U1   at patch named: WALL4-SW =-6.572206E-08

 Nett source of U1   at patch named: WALL5    = 1.989909E-08

 Nett source of U1   at patch named: WALL6-LW = 3.912061E-09

 pos. sum=1.558973E-03 neg. sum=-0.017889

 nett sum=-0.01633

 

 Nett source of V1   at patch named: FOIN     = 1.177630E+00

 Nett source of V1   at patch named: SOIN1    =-3.262800E-01

 Nett source of V1   at patch named: SOIN2    =-3.262800E-01

 Nett source of V1   at patch named: DILUIN   =-5.800000E-01

 Nett source of V1   at patch named: OUTLET   =-7.067809E-02

 Nett source of V1   at patch named: WALL5    = 1.571128E-04

 Nett source of V1   at patch named: WALL6-LW =-1.599715E-05

 pos. sum=1.177787 neg. sum=-1.303254

 nett sum=-0.125467

 

 Nett source of W1   at patch named: OUTLET   =-6.941262E-01

 Nett source of W1   at patch named: WALL1    =-5.542239E-03

 Nett source of W1   at patch named: WALL2    =-2.597433E-03

 Nett source of W1   at patch named: WALL3    =-3.171513E-03

 Nett source of W1   at patch named: WALL4-SW =-7.088535E-05

 pos. sum=0. neg. sum=-0.705508

 nett sum=-0.705508

 

 Nett source of R1   at patch named: FOIN     = 1.177630E-02

 Nett source of R1   at patch named: SOIN1    = 8.157000E-03

 Nett source of R1   at patch named: SOIN2    = 8.157000E-03

 Nett source of R1   at patch named: DILUIN   = 1.450000E-02

 Nett source of R1   at patch named: OUTLET   =-4.257892E-02 (Mass Out -4.821704E-02 In 5.638121E-03)

 pos. sum=0.04259 neg. sum=-0.042579

 nett sum=1.136959E-05

 

 Nett source of KE   at patch named: FOIN     = 2.944074E-03 (Average  2.500000E-01)

 Nett source of KE   at patch named: SOIN1    = 1.003311E-03 (Average  1.230000E-01)

 Nett source of KE   at patch named: SOIN2    = 1.003311E-03 (Average  1.230000E-01)

 Nett source of KE   at patch named: DILUIN   = 1.783500E-03 (Average  1.230000E-01)

 Nett source of KE   at patch named: OUTLET   =-1.146152E+00 (Ave Out  4.018373E+01 In  1.402622E+02)

 Nett source of KE   at patch named: KESOURCE =-3.491978E+01

 pos. sum=6.734196E-03 neg. sum=-36.065937

 nett sum=-36.059204

 

 Nett source of EP   at patch named: FOIN     = 6.918881E-01 (Average  5.875261E+01)

 Nett source of EP   at patch named: SOIN1    = 6.875874E-02 (Average  8.429416E+00)

 Nett source of EP   at patch named: SOIN2    = 6.875874E-02 (Average  8.429416E+00)

 Nett source of EP   at patch named: DILUIN   = 9.134528E-02 (Average  6.299675E+00)

 Nett source of EP   at patch named: OUTLET   =-2.699967E+02 (Ave Out  1.183352E+04 In  5.325760E+04)

 Nett source of EP   at patch named: KESOURCE =-5.170584E+04

 pos. sum=0.920751 neg. sum=-5.197583E+04

 nett sum=-5.197491E+04

 

 Nett source of H1   at patch named: FOIN     = 4.529050E+04 (Average  3.845904E+06)

 Nett source of H1   at patch named: SOIN1    = 9.458041E+03 (Average  1.159500E+06)

 Nett source of H1   at patch named: SOIN2    = 9.458041E+03 (Average  1.159500E+06)

 Nett source of H1   at patch named: DILUIN   = 1.681275E+04 (Average  1.159500E+06)

 Nett source of H1   at patch named: OUTLET   =-8.100595E+04 (Ave Out  1.867542E+06 In  1.604206E+06)

 pos. sum=8.101934E+04 neg. sum=-8.100595E+04

 nett sum=13.382812

 

 Nett source of FUEL at patch named: FOIN     = 6.456302E-04 (Average  5.482456E-02)

 Nett source of FUEL at patch named: OUTLET   =-5.677206E-06 (Ave Out  1.202788E-04 In  2.173634E-05)

 Nett source of FUEL at patch named: WHOLE    =-6.400137E-04

 pos. sum=6.456302E-04 neg. sum=-6.456909E-04

 nett sum=-6.071059E-08

 

 Nett source of MIXF at patch named: FOIN     = 6.456302E-04 (Average  5.482456E-02)

 Nett source of MIXF at patch named: OUTLET   =-6.456262E-04 (Ave Out  1.444982E-02 In  9.075633E-03)

 pos. sum=6.456302E-04 neg. sum=-6.456262E-04

 nett sum=3.958121E-09

 

 

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

 spot values vs sweep or iteration number

 IXMON= 3 IYMON= 5 IZMON= 5 TIMESTEP= 1

 

 Tabulation of abscissa and ordinates...

   ISWP       P1           U1           V1           W1           KE  

      2   1.0000E-10   1.1373E-02   1.5868E+00   4.8148E+00   8.2805E+01

      7   2.7085E+02  -1.9238E-02   3.4977E+00   3.6640E-01   3.9605E+02

     12   2.7328E+02  -4.6909E-03   2.3285E+00   2.7930E-01   4.9987E+02

     17   2.7439E+02   4.1345E-03   2.1552E+00   3.3566E-01   5.7223E+02

     22   3.2096E+02   1.2746E-03   2.6231E+00   6.0980E-01   5.8231E+02

     27   2.8791E+02  -6.6709E-02   1.9868E+00  -1.3522E+00   4.7680E+02

     32   5.2028E+01  -1.1800E-01  -4.2655E-01  -7.1366E+00   3.1569E+02

     37  -4.3144E+01  -1.8577E-01   3.5410E-01  -1.0739E+01   2.0624E+02

     42  -8.7473E+01  -2.5942E-01   1.5747E+00  -1.1368E+01   1.4783E+02

     47  -8.0000E+01  -3.1146E-01   2.2079E+00  -1.1137E+01   1.2698E+02

     52  -7.7364E+01  -3.2648E-01   2.5590E+00  -1.1044E+01   1.1381E+02

     57  -7.2674E+01  -3.3998E-01   2.7698E+00  -1.1083E+01   1.0341E+02

     62  -6.7921E+01  -3.5029E-01   2.9428E+00  -1.1106E+01   9.6034E+01

     67  -6.5745E+01  -3.5614E-01   2.9980E+00  -1.1094E+01   9.2266E+01

     72  -6.5868E+01  -3.5798E-01   3.0121E+00  -1.1093E+01   9.0532E+01

     77  -6.5626E+01  -3.5873E-01   3.0331E+00  -1.1086E+01   8.9694E+01

     82  -6.3653E+01  -3.5957E-01   3.0570E+00  -1.1084E+01   8.9036E+01

     87  -6.2357E+01  -3.6060E-01   3.0744E+00  -1.1087E+01   8.8440E+01

     92  -6.1785E+01  -3.6147E-01   3.0841E+00  -1.1088E+01   8.7991E+01

     97  -6.1500E+01  -3.6171E-01   3.0896E+00  -1.1088E+01   8.7720E+01

   ISWP       EP           H1           FUEL         MIXF

      2   4.0633E+03   2.6063E+06   2.6416E-02   5.4825E-02

      7   4.0407E+04   2.4547E+06   2.2418E-02   5.0959E-02

     12   8.1255E+04   2.4342E+06   1.8959E-02   4.9531E-02

     17   1.6280E+05   2.5712E+06   1.6428E-02   4.8260E-02

     22   2.5804E+05   2.7958E+06   1.2351E-02   4.7560E-02

     27   2.5206E+05   3.0475E+06   8.5119E-03   4.6662E-02

     32   1.4422E+05   3.0904E+06   4.6993E-03   4.3260E-02

     37   6.6135E+04   2.9826E+06   2.6903E-03   3.7934E-02

     42   4.1846E+04   2.7910E+06   1.7385E-03   3.3174E-02

     47   3.3091E+04   2.6707E+06   1.4597E-03   3.0805E-02

     52   2.9094E+04   2.6084E+06   1.2901E-03   2.9566E-02

     57   2.6212E+04   2.5667E+06   1.1583E-03   2.8717E-02

     62   2.4192E+04   2.5353E+06   1.0619E-03   2.8078E-02

     67   2.3067E+04   2.5196E+06   1.0231E-03   2.7757E-02

     72   2.2549E+04   2.5122E+06   1.0097E-03   2.7607E-02

     77   2.2288E+04   2.5070E+06   1.0022E-03   2.7499E-02

     82   2.2085E+04   2.5027E+06   9.9502E-04   2.7412E-02

     87   2.1905E+04   2.4995E+06   9.8851E-04   2.7348E-02

     92   2.1772E+04   2.4977E+06   9.8442E-04   2.7310E-02

     97   2.1690E+04   2.4968E+06   9.8267E-04   2.7291E-02

   Variable    1 = P1     2 = U1     3 = V1     4 = W1     5 = KE  

     Minval= -8.747E+01 -3.617E-01 -4.266E-01 -1.137E+01  8.281E+01

     Maxval=  3.210E+02  1.137E-02  3.498E+00  4.815E+00  5.823E+02

     Cellav=  3.022E+01 -2.241E-01  2.375E+00 -7.309E+00  2.173E+02

   Variable    6 = EP     7 = H1     8 = FUEL   9 = MIXF

     Minval=  4.063E+03  2.434E+06  9.827E-04  2.729E-02

     Maxval=  2.580E+05  3.090E+06  2.642E-02  5.482E-02

     Cellav=  6.594E+04  2.635E+06  6.308E-03  3.608E-02

 1.00 9..3.+..5.+6.6.+7...+....+....+....+....+....+....+

      .    2     2 7                                    .

 0.90 +  2 1  1    1                       3 3  3  3 3  3

      .  9 5            7             3 3               .

 0.80 +    9     3 5               3                    +

      .  4    9  9              3                       .

 0.70 +    8  4    9                                    +

      .  5    6       2       3                         .

 0.60 +       8    4  9                                 +

      .          7    6    7                            .

 0.50 3                    3                            +

      .          8    5 2                               .

 0.40 +                 9                               +

      .               1       7                         .

 0.30 +    6       8                                    +

      7               4 6  2    7                       .

 0.20 1       7         3  9       7                    +

      .  6            8    6  9       7 7               .

 0.10 +                 1     6 9  6  6    7 7  7  7 7  7

      .  7              8  8    5  9  9 6  6 6  6  6 6  6

 0.00 6....7....+....+3...+4..8+8..8+.8.9+.9.9+.9..9.9..9

      0   .1   .2   .3   .4   .5   .6   .7   .8   .9  1.0

 the abscissa is      ISWP.  min= 2.00E+00 max= 9.70E+01

 

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

 

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

 residuals vs sweep or iteration number

 

 Tabulation of abscissa and ordinates...

   ISWP       P1           U1           V1           W1           KE  

      2   2.7498E+00   1.1559E+00   1.8774E+00   8.6199E-01   3.5539E+01

      7   5.0283E-01   1.0854E-01   1.6545E-01   2.0491E-01   1.0534E+00

     12   3.3892E-01   3.1227E-02   7.8816E-02   7.4995E-02   6.1569E-01

     17   3.7875E-01   3.3065E-02   8.3938E-02   6.8393E-02   3.8218E-01

     22   2.8214E-01   2.5133E-02   8.7674E-02   5.9639E-02   3.4040E-01

     27   2.2481E-01   1.8958E-02   7.9198E-02   4.6124E-02   3.1365E-01

     32   1.5733E-01   1.4108E-02   5.0053E-02   3.1947E-02   2.2951E-01

     37   9.1134E-02   1.0925E-02   2.6319E-02   1.9183E-02   1.7407E-01

     42   5.5244E-02   8.2981E-03   1.4997E-02   1.3579E-02   9.0523E-02

     47   4.0905E-02   6.8946E-03   1.1460E-02   7.5317E-03   6.0835E-02

     52   3.4229E-02   5.7004E-03   8.9586E-03   5.0034E-03   5.1946E-02

     57   2.4859E-02   4.6175E-03   6.3617E-03   3.2589E-03   3.7047E-02

     62   1.7111E-02   3.6640E-03   3.9697E-03   2.2672E-03   2.2080E-02

     67   1.1924E-02   2.8376E-03   2.0358E-03   1.6583E-03   9.2269E-03

     72   7.4577E-03   2.2588E-03   1.0126E-03   9.8703E-04   6.5301E-03

     77   5.8739E-03   1.8962E-03   9.8518E-04   6.9460E-04   5.8846E-03

     82   5.1959E-03   1.6428E-03   8.7275E-04   5.5791E-04   4.6176E-03

     87   4.1243E-03   1.4488E-03   6.1094E-04   4.0175E-04   3.3887E-03

     92   3.2642E-03   1.2689E-03   4.1488E-04   3.1452E-04   2.2406E-03

     97   2.3718E-03   1.0972E-03   2.6551E-04   2.0456E-04   1.0568E-03

   ISWP       EP           H1           FUEL         MIXF

      2   3.0966E+01   1.1457E+00   2.9159E+00   2.4634E+00

      7   2.4883E+00   2.8752E-01   5.9149E-01   1.1038E-01

     12   1.2131E+00   1.2309E-01   2.6755E-01   9.6982E-02

     17   9.2877E-01   9.9313E-02   2.7266E-01   1.0082E-01

     22   5.1616E-01   1.0261E-01   2.2904E-01   7.6989E-02

     27   4.1424E-01   8.5540E-02   1.2966E-01   6.1008E-02

     32   2.4160E-01   3.9599E-02   8.2510E-02   4.0441E-02

     37   1.6738E-01   1.7976E-02   3.5614E-02   3.6160E-02

     42   7.8407E-02   1.3660E-02   1.8975E-02   2.9526E-02

     47   4.9735E-02   9.9836E-03   1.2850E-02   2.0713E-02

     52   4.9733E-02   5.8387E-03   1.0435E-02   1.1490E-02

     57   4.1307E-02   3.4927E-03   7.2182E-03   6.9200E-03

     62   2.7321E-02   1.9735E-03   4.5726E-03   3.9475E-03

     67   1.2626E-02   1.2521E-03   2.5912E-03   2.5155E-03

     72   3.9474E-03   8.1504E-04   1.0638E-03   1.6424E-03

     77   4.5277E-03   5.0547E-04   1.4786E-03   1.0209E-03

     82   4.5468E-03   2.9277E-04   1.2852E-03   5.9256E-04

     87   3.3259E-03   1.8885E-04   7.2980E-04   3.8275E-04

     92   1.8609E-03   1.1749E-04   3.2105E-04   2.3830E-04

     97   8.6473E-04   8.9906E-05   1.4423E-04   1.8251E-04

   Variable    1 = P1     2 = U1     3 = V1     4 = W1     5 = KE  

     Minval= -6.044E+00 -6.815E+00 -8.234E+00 -8.495E+00 -6.853E+00

     Maxval=  1.012E+00  1.449E-01  6.299E-01 -1.485E-01  3.571E+00

   Variable    6 = EP     7 = H1     8 = FUEL   9 = MIXF

     Minval= -7.053E+00 -9.317E+00 -8.844E+00 -8.609E+00

     Maxval=  3.433E+00  1.360E-01  1.070E+00  9.015E-01

 1.00 9....+....+....+....+....+....+....+....+....+....+

      .                                                 .

 0.90 +                                                 +

      .  8                                              .

 0.80 +                                                 +

      .  6 8  8  8 7                                    .

 0.70 +    6  4  4 8                                    +

      .  9 9  9  9 4  8                                 .

 0.60 +    5     6 9  4                                 +

      .       5  5 5  9 9  9                            .

 0.50 +    2  2       5 6  8  9                         +

      .          2         6  8 9                       .

 0.40 +            2          6 6  9                    +

      .               2 2       5  6  8                 .

 0.30 +                    2          9 9               +

      .                       2 2       6  9 8          .

 0.20 +                            2    5  8 9  8       +

      .                               2 2  6 6  6  8    .

 0.10 +                                    2 2  9  9 8  +

      .                                         2  2 9  .

 0.00 +....+....+....+....+....+....+....+....+....+.2..9

      0   .1   .2   .3   .4   .5   .6   .7   .8   .9  1.0

 the abscissa is      ISWP.  min= 2.00E+00 max= 9.70E+01

 

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

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

 PATCH(IZ4     ,CONTUR,   1,   6,   1,  10,   4,   4,     1,     1)

 PLOT(IZ4     ,TMP1, 1.000E+00, 1.000E+01)

 Contour plot of TMP1

 Total=0.119788 Cell ave=1574.848022

 Mass ave=1662.279541

 Min valu=1366.047607 Max valu=1885.914673

 PHOENICS polar plot: Angle X,radius Y

 Zero angle is at 3 o-clock

 

 30 I:::::::::::::::::::::::::::::::::::   .......I

 31 I:::::::::::::::::::::::::::::::: 888    .....I

 32 I::::::::::::::::::::::::::::666   888    ....I

 33 I:::::::::::::::::::::::::    6666  888    ...I

 34 I:::::::::::::::::::::44444    6666  888   ...I

 35 I::::::::::::::::::     4444    6666  888   ..I

 36 I:::::::::::::: 2222     4444    666   88    .I

 37 I:::::::::::     2222    4444    6666  888   .I

 38 I:::::::00000     222     4444    666  888   .I

 39 I.:::000000000    2222    4444    666   88   .I

 40 I...::::::::::::::::::::::::::::::::::::::::::I

 40   2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4  

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

 PATCH(XSECIN1 ,CONTUR,   1,   1,   1,  10,   1,  13,     1,     1)

 PLOT(XSECIN1 ,W1  , 0.000E+00, 1.000E+01)

 Contour plot of W1  

   X-plane; Y horizontal, Z vertical

 Total=567.842896 Cell ave=4.368022

 Min valu=-16.920452 Max valu=30.489183

 N,S,E,W... denote: North, South, East ..

 

  0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH

  1 S                                             N

  2 S     4444444444444444444444444444444444444444N

  3 S 4444444444444444444444444444444444444       N

  4 S444444                                     66N

  5 S4444       6666666              6666666666666N

  6 S4      666666666666             6666666666666N

  7 S   666666666666666                  666666666N

  8 S66666666666666666                      666666N

  9 S        66666666                          666N

 10 S         666666                              N

 11 S       66666666                              N

 12 S      66666666               4444444444444444N

 13 S    666666666             4444444444444444444N

 14 S66666666666             44444444444444444444 N

 15 S                          4444444444444444444N

 16 S44444444                   444444444444444444N

 17 S     44444444               44444444444444444N

 18 S22      444444                444444444444444N

 19 S      444444                   44444444444444N

 20 S444444444                       44444444444  N

 21 S                                 444444444   N

 22 S                                  4444444    N

 23 S                                  444444444  N

 24 S                                   4444444444N

 25 S                                   4444444444N

 26 S                                    444444444N

 27 S4444444444444444444                    444444N

 28 S               4444444444                  44N

 29 S22222222222222      44444444                 N

 30 S000000       22222    4444444          666666N

 31 S00000000000     2222     44444     6666666666N

 32 S0000000000000      22222    444    6666    88N

 33 S0000000000000000      2222    44  666   888  N

 34 S00000000000000000       222   44  666   888  N

 35 S0000000000000000         222  44  666   888  N

 36 S222222222222222222222222     4444    66666   N

 37 S                                4444444      N

 38 S       4444444444444          222222      444N

 39 S      444          444    222222222222       N

 40  LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 

 

 40   2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4  

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

 PATCH(XSECIN1 ,CONTUR,   1,   1,   1,  10,   1,  13,     1,     1)

 PLOT(XSECIN1 ,TMP1, 0.000E+00, 1.000E+01)

 Contour plot of TMP1

   X-plane; Y horizontal, Z vertical

 Total=1.821506E+05 Cell ave=1401.158569

 Min valu=773. Max valu=2212.91748

 N,S,E,W... denote: North, South, East ..

 

  0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH

  1 S222222222222222222            444444444444   N

  2 S     2222222222222            444444444444   N

  3 S       22222222222            444444444444   N

  4 S        22222222222           44444444444    N

  5 S         2222222222           44444444444    N

  6 S          2222222222           4444444444    N

  7 S          22222222222            44444444    N

  8 S          222222222222             4444444   N

  9 S          22222222222222            444444   N

 10 S            222222222222222222          44444N

 11 S                  2222222222222222222222     N

 12 S                                    222222222N

 13 S                                             N

 14 S                                             N

 15 S2222222222222222222                2222222222N

 16 S222222222222222222222222222222222222         N

 17 S                                    4444444  N

 18 S                             44444444     666N

 19 S                              44444444     66N

 20 S                              444444444    66N

 21 S                               444444444    6N

 22 S                                 4444444    6N

 23 S                                      444444 N

 24 S                                            4N

 25 S                               22222222222222N

 26 S                         22222222222222222   N

 27 S                                       222222N

 28 S            444444444444444444444444444444444N

 29 S      444444444444         66666666666       N

 30 S 44444444444444     66666                    N

 31 S44444444444444     6666                888888N

 32 S44444444444444     66666                   88N

 33 S44444444444444      66666666                 N

 34 S444444444444444       66666666666666         N

 35 S44444444444444444           6666666666       N

 36 S         222222      4444444      66666      N

 37 S000000  2 444             666666666      8888N

 38 S0000002468              888888888888888888888N

 39 S000000246                   88888888 6 444444N

 40  LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 

 

 40   2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4  

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

 PATCH(XSECIN1 ,CONTUR,   1,   1,   1,  10,   1,  13,     1,     1)

 PLOT(XSECIN1 ,FUEL, 0.000E+00, 1.000E+01)

 Contour plot of FUEL

   X-plane; Y horizontal, Z vertical

 Total=0.455206 Cell ave=3.501582E-03

 Min valu=0. Max valu=0.04833

 N,S,E,W... denote: North, South, East ..

 

  0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH

  1 S000000000000000000000000000000000000000000000N

  2 S000000000000000000000000000000000000000000000N

  3 S000000000000000000000000000000000000000000000N

  4 S000000000000000000000000000000000000000000000N

  5 S000000000000000000000000000000000000000000000N

  6 S000000000000000000000000000000000000000000000N

  7 S000000000000000000000000000000000000000000000N

  8 S000000000000000000000000000000000000000000000N

  9 S000000000000000000000000000000000000000000000N

 10 S000000000000000000000000000000000000000000000N

 11 S000000000000000000000000000000000000000000000N

 12 S000000000000000000000000000000000000000000000N

 13 S000000000000000000000000000000000000000000000N

 14 S000000000000000000000000000000000000000000000N

 15 S000000000000000000000000000000000000000000000N

 16 S000000000000000000000000000000000000000000000N

 17 S000000000000000000000000000000000000000000000N

 18 S000000000000000000000000000000000000000000000N

 19 S000000000000000000000000000000000000000000000N

 20 S000000000000000000000000000000000000000000000N

 21 S000000000000000000000000000000000000000000000N

 22 S000000000000000000000000000000000000000000000N

 23 S000000000000000000000000000000000000000000000N

 24 S000000000000000000000000000000000000000000000N

 25 S000000000000000000000000000000000000000000000N

 26 S000000000000000000000000000000000000000000000N

 27 S000000000000000000000000000000000000000000000N

 28 S000000000000000000000000000000000000000000000N

 29 S00000000000000000000000000000000000000000000 N

 30 S00000000000000000000000000000000000000000    N

 31 S000000000000000000000000000000000000000      N

 32 S00000000000000000000000000000000000          N

 33 S000000000000000000000000000000000            N

 34 S000000000000000000000000000000               N

 35 S00000000000000000000000000                 22N

 36 S000000 24 6666        444444444       2222222N

 37 S000000 24 666       4444444444       22222222N

 38 S0000000000   0000           22222            N

 39 S000000000                              000000N

 40  LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 

 

 40   2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4  

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

 PATCH(XSECIN1 ,CONTUR,   1,   1,   1,  10,   1,  13,     1,     1)

 PLOT(XSECIN1 ,MIXF, 0.000E+00, 1.000E+01)

 Contour plot of MIXF

   X-plane; Y horizontal, Z vertical

 Total=2.955019 Cell ave=0.022731

 Min valu=0. Max valu=0.0533

 N,S,E,W... denote: North, South, East ..

 

  0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH

  1 S             2222222222                  4444N

  2 S             2222222222                 44444N

  3 S              222222222                 44444N

  4 S              222222222                444444N

  5 S               222222222               444444N

  6 S               2222222222               44444N

  7 S               2222222222222            44444N

  8 S                222222222222222          4444N

  9 S                2222222222222222         4444N

 10 S                    22222222222222222      44N

 11 S                              2222222222222  N

 12 S                                          222N

 13 S                                             N

 14 S                            00000000000000000N

 15 S                                          222N

 16 S22222222222222222222222222222222222222222    N

 17 S22222222222222222222222222222222         4444N

 18 S222222222                           44444    N

 19 S222222                              444444   N

 20 S                                     44444   N

 21 S                                      4444   N

 22 S                                      44444  N

 23 S                                         4444N

 24 S                  22222222222222222222       N

 25 S               222222222222222222222222222222N

 26 S             2222222222222222222222222       N

 27 S                                  22222222222N

 28 S                      44444444444444444444444N

 29 S                4444444                      N

 30 S             444444     6666666666666666     N

 31 S            444444     66666666            88N

 32 S            444444      666666666          88N

 33 S            4444444      6666666666        88N

 34 S             4444444       66666666        88N

 35 S              4444444       6666666       888N

 36 S      24 6                8888888888888888888N

 37 S0000002468                   8888888888888888N

 38 S0000002   888888888         88888888888888888N

 39 S0000002   88888888888888 88888888888 6 444444N

 40  LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 

 

 40   2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4  

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

 PATCH(XSECIN2 ,CONTUR,   4,   4,   1,  10,   1,  13,     1,     1)

 PLOT(XSECIN2 ,W1  , 0.000E+00, 1.000E+01)

 Contour plot of W1  

   X-plane; Y horizontal, Z vertical

 Total=748.194763 Cell ave=5.755344

 Min valu=-16.921331 Max valu=30.721508

 N,S,E,W... denote: North, South, East ..

 

  0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH

  1 S                                             N

  2 S    44444444444444444444444444444444444444444N

  3 S 444444444         4444444444444444          N

  4 S4444                                      666N

  5 S444    6666666666666666666666666666666666666 N

  6 S     6666666666666666666666666666666666666   N

  7 S  666666666666666666666666666666666666666    N

  8 S66666666666666666666666666666666666666666    N

  9 S            66666666666666666666666666666    N

 10 S           66666666666666666666666666666     N

 11 S          666666666666666666666666666666     N

 12 S       666666666666666666666666666666        N

 13 S     66666666666666666666666666666           N

 14 S6666666666666666666666666666666666           N

 15 S                 6666666666666666666         N

 16 S444444444            66666666666666666       N

 17 S      444444444         6666666666666666     N

 18 S222        444444        6666666666666666    N

 19 S       44444444          666666666666666666  N

 20 S44444444444               6666666666666666666N

 21 S                                         6666N

 22 S               6666666                     66N

 23 S                                             N

 24 S                                             N

 25 S                                          444N

 26 S                                        44444N

 27 S4444444444444444444                      4444N

 28 S                 444444444                   N

 29 S2222222222222222      444444                 N

 30 S0000000         22222   444444       666666  N

 31 S0000000000000      2222    4444    6666666666N

 32 S000000000000000      2222    444   6666     8N

 33 S000000000000000000     2222   44   666   88  N

 34 S000000000000000000       222  44   66   888  N

 35 S00000000000000000        222  44  666   8888 N

 36 S222222222222222222222222     4444    66666   N

 37 S                                4444444      N

 38 S       4444444444444         2222222       44N

 39 S      444          444    222222222222       N

 40  LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 

 

 40   2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4  

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

 PATCH(XSECIN2 ,CONTUR,   4,   4,   1,  10,   1,  13,     1,     1)

 PLOT(XSECIN2 ,TMP1, 0.000E+00, 1.000E+01)

 Contour plot of TMP1

   X-plane; Y horizontal, Z vertical

 Total=1.83871E+05 Cell ave=1414.39209

 Min valu=773. Max valu=2212.382568

 N,S,E,W... denote: North, South, East ..

 

  0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH

  1 S222222222222222222222222222                44N

  2 S      2222222222222222222222               44N

  3 S        22222222222222222222              444N

  4 S         22222222222222222222             444N

  5 S           222222222222222222             444N

  6 S           222222222222222222            4444N

  7 S           222222222222222222            4444N

  8 S           222222222222222222           4444 N

  9 S           222222222222222222           4444 N

 10 S            2222222222222222           4444  N

 11 S            2222222222222222          4444   N

 12 S             22222222222222          44444   N

 13 S             22222222222222         44444    N

 14 S           222222222222222         444444   6N

 15 S2222222222222222222222            444444    6N

 16 S222222222222222                 44444444   66N

 17 S                              444444444    66N

 18 S                          4444444444444   666N

 19 S                         4444444444444    666N

 20 S                        44444444444444    666N

 21 S                      4444444444444444    666N

 22 S                     4444444444444444     666N

 23 S                      4444444444444444444    N

 24 S                                       444444N

 25 S                                             N

 26 S                                  22222222222N

 27 S              44444444444444444             2N

 28 S         4444444444444            44444444444N

 29 S     444444444444       666666666666666666666N

 30 S 44444444444444     66666                    N

 31 S44444444444444     66666               888888N

 32 S44444444444444      66666                  88N

 33 S444444444444444     66666666                 N

 34 S444444444444444        6666666666666         N

 35 S44444444444444444           6666666666       N

 36 S         222222      4444444      66666      N

 37 S000000  2 444             666666666      8888N

 38 S0000002468              888888888888888888888N

 39 S000000246                   88888888 6 444444N

 40  LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 

 

 40   2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4  

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

 PATCH(XSECIN2 ,CONTUR,   4,   4,   1,  10,   1,  13,     1,     1)

 PLOT(XSECIN2 ,FUEL, 0.000E+00, 1.000E+01)

 Contour plot of FUEL

   X-plane; Y horizontal, Z vertical

 Total=0.456564 Cell ave=3.512027E-03

 Min valu=0. Max valu=0.048329

 N,S,E,W... denote: North, South, East ..

 

  0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH

  1 S000000000000000000000000000000000000000000000N

  2 S000000000000000000000000000000000000000000000N

  3 S000000000000000000000000000000000000000000000N

  4 S000000000000000000000000000000000000000000000N

  5 S000000000000000000000000000000000000000000000N

  6 S000000000000000000000000000000000000000000000N

  7 S000000000000000000000000000000000000000000000N

  8 S000000000000000000000000000000000000000000000N

  9 S000000000000000000000000000000000000000000000N

 10 S000000000000000000000000000000000000000000000N

 11 S000000000000000000000000000000000000000000000N

 12 S000000000000000000000000000000000000000000000N

 13 S000000000000000000000000000000000000000000000N

 14 S000000000000000000000000000000000000000000000N

 15 S000000000000000000000000000000000000000000000N

 16 S000000000000000000000000000000000000000000000N

 17 S000000000000000000000000000000000000000000000N

 18 S000000000000000000000000000000000000000000000N

 19 S000000000000000000000000000000000000000000000N

 20 S000000000000000000000000000000000000000000000N

 21 S000000000000000000000000000000000000000000000N

 22 S000000000000000000000000000000000000000000000N

 23 S000000000000000000000000000000000000000000000N

 24 S000000000000000000000000000000000000000000000N

 25 S000000000000000000000000000000000000000000000N

 26 S000000000000000000000000000000000000000000000N

 27 S000000000000000000000000000000000000000000000N

 28 S000000000000000000000000000000000000000000000N

 29 S00000000000000000000000000000000000000000000 N

 30 S00000000000000000000000000000000000000000    N

 31 S000000000000000000000000000000000000000      N

 32 S000000000000000000000000000000000000         N

 33 S000000000000000000000000000000000            N

 34 S000000000000000000000000000000               N

 35 S00000000000000000000000000                 22N

 36 S000000 24 6666        444444444       2222222N

 37 S000000 24 666       4444444444       22222222N

 38 S0000000000   0000           22222            N

 39 S000000000                              000000N

 40  LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 

 

 40   2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4  

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

 PATCH(XSECIN2 ,CONTUR,   4,   4,   1,  10,   1,  13,     1,     1)

 PLOT(XSECIN2 ,MIXF, 0.000E+00, 1.000E+01)

 Contour plot of MIXF

   X-plane; Y horizontal, Z vertical

 Total=3.009044 Cell ave=0.023146

 Min valu=0. Max valu=0.053297

 N,S,E,W... denote: North, South, East ..

 

  0 HHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHHH

  1 S               22222222222222222222222222    N

  2 S                2222222222222222222222222    N

  3 S                 22222222222222222222222     N

  4 S                 22222222222222222222222     N

  5 S                  2222222222222222222222     N

  6 S                  222222222222222222222     4N

  7 S                  22222222222222222222      4N

  8 S                   222222222222222222      44N

  9 S                   22222222222222222       44N

 10 S                   2222222222222222       444N

 11 S                   222222222222222       4444N

 12 S                   22222222222222        4444N

 13 S                   2222222222222        4444 N

 14 S                 22222222222222         4444 N

 15 S          2222222222222222222          4444  N

 16 S22222222222222222222222222            4444   N

 17 S222222222222222222222                44444   N

 18 S222222222                           444444   N

 19 S222222                              44444    N

 20 S                                   444444    N

 21 S                                   444444    N

 22 S                                   444444    N

 23 S                                      444444 N

 24 S                                            4N

 25 S                                     2222    N

 26 S                               22222222222222N

 27 S                                         2222N

 28 S                  444444444444444444444444444N

 29 S               4444444                       N

 30 S             444444     6666666666666666     N

 31 S            444444     66666666            88N

 32 S            4444444     666666666          88N

 33 S            4444444       666666666        88N

 34 S             4444444       66666666        88N

 35 S              4444444       6666666       888N

 36 S      24 6                8888888888888888888N

 37 S0000002468                   8888888888888888N

 38 S0000002   888888888         88888888888888888N

 39 S0000002   88888888888888888888888888 6 444444N

 40  LLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL 

 

 40   2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4 6 8 0 2 4  

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

 

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

 SATLIT RUN NUMBER =   1 ; LIBRARY REF.=  492

 Run started at 15:01:01 on Friday, 25 May 2012                                

 Run completed at 15:01:02 on Friday, 25 May 2012                              

 CPU time of run 2 s

 This includes 1 seconds of user time and 0 seconds of system time.

 TIME/(VARIABLES*CELLS*TSTEPS*SWEEPS*ITS) = 2.849E-06

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