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