TALK=T;RUN(1,1)
  PHOTON USE
  ext;;;;

  gr ou z 1; use patgeo
  msg boundary condition patches. Press RETURN
  pause
  msg temperature contours. Press RETURN
  con tmp1 y m fi;0.001;con tmp1 x m fi;0.001;con tmp1 z m fi;0.001
  gr ou y m;gr ou z m
  msg press RETURN for view x
  pause;  con off;  view x;  con tmp1 x m fi;0.001
  con tmp1 y m fi;0.001;  gr ou y m;  gr ou x m;  gr ou z 1
  msg This is view x. Press RETURN for view y
  pause;  con off;  view y;  con tmp1 y m fi;0.001
  msg This is view y. Press RETURN for view z
  pause;  con off; view z;  con tmp1 z m fi;0.001;  gr ou z m
  msg This is view z. Press RETURN for view x and velocity vectors
  pause;  con off;  view x
  msg velocity vectors. Type menu for menu and further possibilities
  msg Press e to end. Otherwise enter photon-readable commands
  vec x m sh;  gr ou y m;  gr ou x m;  gr ou z 1
  ENDUSE
  DISPLAY
  READQ1_BEGIN
  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.                       *
    ****************************************************************
  READQ1_END
  The shape of the combustion chamber is as shown.
  Pre-mixed fuel vapour and air enter near the axis on left.
  Secondary and dilution air enter through holes in outer wall.
  A 36-degree sector is simulated.

                 secondary ox. inlets   dilution inlet
           ____________ 1,2 _______________     ________
  blocked /
  region /
        /
       /                                             outlet
       |
       |_____   _
  fuel-ox. inlet|            Symmetry axis
       -- - -- -| -- - -- -  -- - -- - -- - -- - -- - -- -

  The flow is turbulent; the Simple Chemical Reaction Scheme is
  used; and the reaction-rate is physically controlled by means of
  the Eddy-Breakup Model.
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(Idealised Gas-Turbine Combustion Chamber)
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 6
 XULAST =0.628319
   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,Y)
 TERMS(U1,Y,Y,Y,Y,Y,Y)
 TERMS(V1,Y,Y,Y,Y,Y,Y)
 TERMS(W1,Y,Y,Y,Y,Y,Y)
 TERMS(KE,N,Y,Y,Y,Y,N)
 TERMS(EP,N,Y,Y,Y,Y,N)
 TERMS(H1,N,Y,Y,N,Y,N)
 TERMS(FUEL,N,Y,Y,Y,N,Y)
 TERMS(MIXF,N,Y,Y,Y,Y,Y)
 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(KESOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(KESOURCE,KE , GRND4 , GRND4 )
 COVAL(KESOURCE,EP , GRND4 , GRND4 )
 
 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.752605 )
 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(WHOLE ,PHASEM, 1, 6, 1, 10, 1, 13, 1, 1)
 COVAL(WHOLE ,FUEL,In-Form:source - see Grp 19)
 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-04
 ************************************************************
  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)=-1.0E+10
 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
 ISG62 = 1
 CHSOA =0.054825 ;CHSOB =1.
 SPEDAT(SET,PROPERTY,TMP1,C,=MAX(773.&(H1-4.9E+07*FUEL)/(1500.*FUE$)
 SPEDAT(SET,PROPERTY,TMP1,C,L+1500.*OXID+1500.*PROD))!IMAT<100)
 SPEDAT(SET,PROPERTY,RHO1,C,=(P1+8.0E+05)*1.0/(FUEL/16.+OXID/29.+P$)
 SPEDAT(SET,PROPERTY,RHO1,C,ROD/28.)/(TMP1*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 = -5
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100 ;NUMCLS = 5
 NXPRIN = -1 ;IXPRF = 3 ;IXPRL = 3
 NYPRIN = 2 ;IYPRF = 1 ;IYPRL = 10000
 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
STOP