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

  view z
  msg exit smoke concentration according to single-fluid model
  con smos z m fi;0.001
  pause;con off;red
  msg exit smoke concentration according to multi-fluid model
  con smom z m fi;0.001;pause;con off;red;view x

  msg x-plane smoke concentration according to single-fluid model
  con smos x 1 fi;0.001;pause;con off;red
  msg x-plane smoke concentration according to multi-fluid model
  con smom x 1 fi;0.001;pause;con off;red
  msg x-plane mixture fraction avef according to multi-fluid model
  con avef x 1 fi;0.001;pause;con off;red
  msg x-plane mixture fraction mixf according to single-fluid model
  con mixf x 1 fi;0.001;pause;con off;red
  msg x-plane fluid-1 concentration according to multi-fluid model
  con f1 x 1 fi;0.001;pause;con off;red
  msg x-plane fluid-11 concentration according to multi-fluid model
  con f11 x 1 fi;0.001;pause;con off;red
  msg x-plane fluid-6 concentration according to multi-fluid model
  con f6 x 1 fi;0.001
  ENDUSE
  DISPLAY
  This simulation shows how the multi-fluid model of turbulence
  predicts a smoke-production pattern in a gas-turbine combustor
  which is different from, and probably more realistic than, that
  predicted by a conventional single-fluid model.

  It is explained and discussed extensively in a report which is
  published on CHAM's website www.cham.co.uk (Click on short-cuts,
  then MFM then 1998b).

  The geometry is that of library case 492, which is loaded first.
  Note the $ in #$492; it ensures that the above PHOTON USE
  commands remain in the Q1 file, and therefore take priority over
  those which appear in case 493.

  The injected fuel-air composition, which is stoichiometric in
  case 492, is made twice stoichiometric here.

  To reduce computation time, it is best to answer N to the question
  about radiation.
  ENDDIS

  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
  DISPLAY
  Notes:

  (1) four macros are provided in order reduce the repetiveness of
      MFM q1s.

  (2) they are named mfm, mfm1, mfm2 and mfm3, and are normally
      called in that order at the top of the q1 file, which then
      appears  as:

      #mfm
      settings of nflr, nflf, etc
      #mfm1
      settings of
      #mfm2
      #mfm3

  (3) the character variable mfm is declared and set in the
      always-loaded core-library macro 014.htm, which, in its turn,
      declares and sets the character variables mfm1, mfm2, mfm3.

  (4) the values are:
      mfm  = $L004
      mfm1 = $L003
      mfm2 = $L002
      mfm3 = $L001
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(MFM; 3D Cmbstr; Conmix=5.               )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 600
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 1
 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(2)=6.5E-03
 YFRAC(3)=1. ;YFRAC(4)=5.0E-03
 YFRAC(5)=1. ;YFRAC(6)=8.0E-03
 YFRAC(7)=3. ;YFRAC(8)=6.5E-03
 ************************************************************
  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(115)=PRPS
 NAME(116)=PRDM ;NAME(117)=OXDM
 NAME(118)=FRFM ;NAME(119)=TMPM
 NAME(120)=SMOS ;NAME(121)=SMOM
 NAME(122)=F1 ;NAME(123)=F2
 NAME(124)=F3 ;NAME(125)=F4
 NAME(126)=F5 ;NAME(127)=F6
 NAME(128)=F7 ;NAME(129)=F8
 NAME(130)=F9 ;NAME(131)=F10
 NAME(132)=F11 ;NAME(133)=MNSQ
 NAME(134)=AVEF ;NAME(135)=MIXL
 NAME(136)=RATE ;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,Y)
 SOLUTN(PRDM,Y,N,N,N,N,Y)
 SOLUTN(OXDM,Y,N,N,N,N,Y)
 SOLUTN(FRFM,Y,N,N,N,N,Y)
 SOLUTN(TMPM,Y,N,N,N,N,Y)
 SOLUTN(SMOS,Y,Y,N,N,N,Y)
 SOLUTN(SMOM,Y,Y,N,N,N,Y)
 SOLUTN(F1,Y,Y,N,N,N,Y)
 SOLUTN(F2,Y,Y,N,N,N,Y)
 SOLUTN(F3,Y,Y,N,N,N,Y)
 SOLUTN(F4,Y,Y,N,N,N,Y)
 SOLUTN(F5,Y,Y,N,N,N,Y)
 SOLUTN(F6,Y,Y,N,N,N,Y)
 SOLUTN(F7,Y,Y,N,N,N,Y)
 SOLUTN(F8,Y,Y,N,N,N,Y)
 SOLUTN(F9,Y,Y,N,N,N,Y)
 SOLUTN(F10,Y,Y,N,N,N,Y)
 SOLUTN(F11,Y,Y,N,N,N,Y)
 SOLUTN(MNSQ,Y,N,N,N,N,Y)
 SOLUTN(AVEF,Y,N,N,N,N,Y)
 SOLUTN(MIXL,Y,N,N,N,N,N)
 SOLUTN(RATE,Y,N,N,N,N,Y)
 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,N,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 = 115
 ************************************************************
  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(SMOS,N,Y,Y,Y,Y,Y)
 TERMS(SMOM,N,Y,Y,Y,N,Y)
 TERMS(F1,N,Y,Y,Y,Y,Y)
 TERMS(F2,N,Y,Y,Y,N,Y)
 TERMS(F3,N,Y,Y,Y,Y,Y)
 TERMS(F4,N,Y,Y,Y,N,Y)
 TERMS(F5,N,Y,Y,Y,Y,Y)
 TERMS(F6,N,Y,Y,Y,N,Y)
 TERMS(F7,N,Y,Y,Y,Y,Y)
 TERMS(F8,N,Y,Y,Y,N,Y)
 TERMS(F9,N,Y,Y,Y,Y,Y)
 TERMS(F10,N,Y,Y,Y,N,Y)
 TERMS(F11,N,Y,Y,Y,Y,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 = GRND7
 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.054825
 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(SMOS)=1. ;PRNDTL(SMOM)=1.
 PRNDTL(F1)=1. ;PRNDTL(F2)=1.
 PRNDTL(F3)=1. ;PRNDTL(F4)=1.
 PRNDTL(F5)=1. ;PRNDTL(F6)=1.
 PRNDTL(F7)=1. ;PRNDTL(F8)=1.
 PRNDTL(F9)=1. ;PRNDTL(F10)=1.
 PRNDTL(F11)=1. ;PRNDTL(MIXF)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(W1)=1. ;PRT(KE)=1.
 PRT(EP)=1.314 ;PRT(H1)=1.
 PRT(SMOS)=1. ;PRT(SMOM)=1.
 PRT(F1)=1. ;PRT(F2)=1.
 PRT(F3)=1. ;PRT(F4)=1.
 PRT(F5)=1. ;PRT(F6)=1.
 PRT(F7)=1. ;PRT(F8)=1.
 PRT(F9)=1. ;PRT(F10)=1.
 PRT(F11)=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(PRDM)=1.0E-10 ;FIINIT(OXDM)=1.0E-10
 FIINIT(FRFM)=1.0E-10 ;FIINIT(TMPM)=1.0E-10
 FIINIT(SMOS)=0. ;FIINIT(SMOM)=0.
 FIINIT(F1)=0.090909 ;FIINIT(F2)=0.090909
 FIINIT(F3)=0.090909 ;FIINIT(F4)=0.090909
 FIINIT(F5)=0.090909 ;FIINIT(F6)=0.090909
 FIINIT(F7)=0.090909 ;FIINIT(F8)=0.090909
 FIINIT(F9)=0.090909 ;FIINIT(F10)=0.090909
 FIINIT(F11)=0.090909 ;FIINIT(MNSQ)=1.0E-10
 FIINIT(AVEF)=1.0E-10 ;FIINIT(MIXL)=1.0E-10
 FIINIT(RATE)=31.400507 ;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.109649
 FIINIT(MIXF)=0.109649
 
 PATCH(CMP1 ,INIVAL, 1, 6, 1, 2, 1, 3, 1, 1)
 INIT(CMP1 ,PRPS,0. ,199. )
 
 PATCH(CMP2 ,INIVAL, 1, 1, 7, 7, 1, 1, 1, 1)
 
 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, 1, 8, 8, 1, 1, 1, 1)
 
 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, 1, 9, 10, 1, 1, 1, 1)
 
 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, 1, 8, 8, 2, 2, 1, 1)
 
 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, 1, 9, 9, 2, 2, 1, 1)
 
 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, 1, 10, 10, 2, 2, 1, 1)
 
 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 ,180.300003 )
 COVAL(FOIN ,U1 ,0. ,0. )
 COVAL(FOIN ,V1 ,0. ,50. )
 COVAL(FOIN ,W1 ,0. ,0. )
 COVAL(FOIN ,KE ,0. ,0.0625 )
 COVAL(FOIN ,EP ,0. ,3.672038 )
 COVAL(FOIN ,H1 ,0. ,4.001728E+06 )
 COVAL(FOIN ,F11 ,0. ,1. )
 COVAL(FOIN ,FUEL,0. ,0.058005 )
 COVAL(FOIN ,MIXF,0. ,0.109649 )
 
 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 ,F1 ,0. ,1. )
 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 ,F1 ,0. ,1. )
 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 ,F1 ,0. ,1. )
 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 ,SMOS,0. , SAME )
 COVAL(OUTLET ,SMOM,0. , SAME )
 COVAL(OUTLET ,F1 ,0. , SAME )
 COVAL(OUTLET ,F2 ,0. , SAME )
 COVAL(OUTLET ,F3 ,0. , SAME )
 COVAL(OUTLET ,F4 ,0. , SAME )
 COVAL(OUTLET ,F5 ,0. , SAME )
 COVAL(OUTLET ,F6 ,0. , SAME )
 COVAL(OUTLET ,F7 ,0. , SAME )
 COVAL(OUTLET ,F8 ,0. , SAME )
 COVAL(OUTLET ,F9 ,0. , SAME )
 COVAL(OUTLET ,F10 ,0. , SAME )
 COVAL(OUTLET ,F11 ,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)
 
 PATCH(SMOKE ,PHASEM, 1, 6, 1, 10, 1, 13, 1, 1)
 COVAL(SMOKE ,SMOS, FIXFLU , GRND )
 COVAL(SMOKE ,SMOM, FIXFLU , GRND )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 500 ;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(SMOS)=20
 LITER(SMOM)=20 ;LITER(F1)=20
 LITER(F2)=20 ;LITER(F3)=20
 LITER(F4)=20 ;LITER(F5)=20
 LITER(F6)=20 ;LITER(F7)=20
 LITER(F8)=20 ;LITER(F9)=20
 LITER(F10)=20 ;LITER(F11)=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(SMOS)=1.0E-03
 ENDIT(SMOM)=1.0E-03 ;ENDIT(F1)=1.0E-03
 ENDIT(F2)=1.0E-03 ;ENDIT(F3)=1.0E-03
 ENDIT(F4)=1.0E-03 ;ENDIT(F5)=1.0E-03
 ENDIT(F6)=1.0E-03 ;ENDIT(F7)=1.0E-03
 ENDIT(F8)=1.0E-03 ;ENDIT(F9)=1.0E-03
 ENDIT(F10)=1.0E-03 ;ENDIT(F11)=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(PRDM,LINRLX,1.)
 RELAX(OXDM,LINRLX,1.)
 RELAX(FRFM,LINRLX,1.)
 RELAX(TMPM,LINRLX,1.)
 RELAX(SMOS,FALSDT,1.0E+09)
 RELAX(SMOM,FALSDT,1.0E+09)
 RELAX(F1,LINRLX,1.)
 RELAX(F2,LINRLX,1.)
 RELAX(F3,LINRLX,1.)
 RELAX(F4,LINRLX,1.)
 RELAX(F5,LINRLX,1.)
 RELAX(F6,LINRLX,1.)
 RELAX(F7,LINRLX,1.)
 RELAX(F8,LINRLX,1.)
 RELAX(F9,LINRLX,1.)
 RELAX(F10,LINRLX,1.)
 RELAX(F11,LINRLX,1.)
 RELAX(MNSQ,LINRLX,1.)
 RELAX(AVEF,LINRLX,1.)
 RELAX(MIXL,LINRLX,1.)
 RELAX(RATE,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,LINRLX,1.)
 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(PRDM)=1.0E+10 ;VARMIN(PRDM)=-1.0E+10
 VARMAX(OXDM)=1.0E+10 ;VARMIN(OXDM)=-1.0E+10
 VARMAX(FRFM)=1.0E+10 ;VARMIN(FRFM)=-1.0E+10
 VARMAX(TMPM)=1.0E+10 ;VARMIN(TMPM)=-1.0E+10
 VARMAX(SMOS)=1.0E+10 ;VARMIN(SMOS)=-1.0E+10
 VARMAX(SMOM)=1.0E+10 ;VARMIN(SMOM)=-1.0E+10
 VARMAX(F1)=1. ;VARMIN(F1)=0.
 VARMAX(F2)=1. ;VARMIN(F2)=0.
 VARMAX(F3)=1. ;VARMIN(F3)=0.
 VARMAX(F4)=1. ;VARMIN(F4)=0.
 VARMAX(F5)=1. ;VARMIN(F5)=0.
 VARMAX(F6)=1. ;VARMIN(F6)=0.
 VARMAX(F7)=1. ;VARMIN(F7)=0.
 VARMAX(F8)=1. ;VARMIN(F8)=0.
 VARMAX(F9)=1. ;VARMIN(F9)=0.
 VARMAX(F10)=1. ;VARMIN(F10)=0.
 VARMAX(F11)=1. ;VARMIN(F11)=0.
 VARMAX(MNSQ)=1. ;VARMIN(MNSQ)=0.
 VARMAX(AVEF)=1. ;VARMIN(AVEF)=0.
 VARMAX(MIXL)=1.0E+10 ;VARMIN(MIXL)=-1.0E+10
 VARMAX(RATE)=1.0E+10 ;VARMIN(RATE)=-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,STORED,MIXL,C,=LEN1!ZSLFIN)
 SPEDAT(SET,STORED,RATE,C,=EPKE!ZSLFIN)
 SPEDAT(SET,MFM,FFUIN,R,0.109649)
 SPEDAT(SET,MFM,TAIRIN,R,773.)
 SPEDAT(SET,MFM,TFUIN,R,773.)
 SPEDAT(SET,MFM,TSTOIC,R,2106.90918)
 SPEDAT(SET,MFM,FACST,R,0.054825)
 SPEDAT(SET,MFM,GASCON,R,8314.)
 SPEDAT(SET,MFM,SPVLMX,R,100.)
 SPEDAT(SET,MFM,SPVLMN,R,1.0E-02)
 SPEDAT(SET,MFM,SMOCON,R,1.)
 SPEDAT(SET,MFM,SMOEXP,R,7.)
 SPEDAT(SET,MFM,MFMMOD,C,SCRS)
 SPEDAT(SET,MFM,NFLUIDS,I,11)
 SPEDAT(SET,MFM,NFLR,I,1)
 SPEDAT(SET,MFM,NFLF,I,11)
 SPEDAT(SET,MFM,CONREA,R,0.)
 SPEDAT(SET,MFM,CONMIX,R,5.)
 SPEDAT(SET,MFM,VISCON,R,0.)
 SPEDAT(SET,MATERIAL,199,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=5893. ;EX(U1)=4.551E-15
 EX(V1)=6.941 ;EX(W1)=24.77
 EX(KE)=23.209999 ;EX(EP)=1.065E+04
 EX(H1)=1.765E+06 ;EX(PRPS)=0.9385
 EX(PRDM)=0.3405 ;EX(OXDM)=0.5908
 EX(FRFM)=7.154E-03 ;EX(TMPM)=1015.
 EX(SMOS)=3.597E-06 ;EX(SMOM)=7.34E-06
 EX(F1)=0.4534 ;EX(F2)=0.1159
 EX(F3)=0.04251 ;EX(F4)=0.03241
 EX(F5)=0.03351 ;EX(F6)=0.03746
 EX(F7)=0.03901 ;EX(F8)=0.05489
 EX(F9)=0.05942 ;EX(F10)=0.02541
 EX(F11)=0.04501 ;EX(MNSQ)=0.2738
 EX(AVEF)=0.0262 ;EX(MIXL)=0.
 EX(RATE)=1176. ;EX(DEGF)=0.
 EX(RRAT)=0. ;EX(EPKE)=1176.
 EX(HPOR)=0. ;EX(NPOR)=0.
 EX(EPOR)=0. ;EX(VPOR)=0.
 EX(ENUT)=5.211E-03 ;EX(RHO1)=2.618
 EX(TMP1)=983.900024 ;EX(PROD)=0.3688
 EX(OXID)=0.5638 ;EX(FUEL)=5.906E-03
 EX(MIXF)=0.02612
 ************************************************************
  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(PRDM,Y,N,Y,N,N,N)
 OUTPUT(OXDM,Y,N,Y,N,N,N)
 OUTPUT(FRFM,Y,N,Y,N,N,N)
 OUTPUT(TMPM,Y,N,Y,N,N,N)
 OUTPUT(SMOS,Y,N,Y,Y,Y,Y)
 OUTPUT(SMOM,Y,N,Y,Y,Y,Y)
 OUTPUT(F1,Y,N,Y,Y,Y,Y)
 OUTPUT(F2,Y,N,Y,Y,Y,Y)
 OUTPUT(F3,Y,N,Y,Y,Y,Y)
 OUTPUT(F4,Y,N,Y,Y,Y,Y)
 OUTPUT(F5,Y,N,Y,Y,Y,Y)
 OUTPUT(F6,Y,N,Y,Y,Y,Y)
 OUTPUT(F7,Y,N,Y,Y,Y,Y)
 OUTPUT(F8,Y,N,Y,Y,Y,Y)
 OUTPUT(F9,Y,N,Y,Y,Y,Y)
 OUTPUT(F10,Y,N,Y,Y,Y,Y)
 OUTPUT(F11,Y,N,Y,Y,Y,Y)
 OUTPUT(MNSQ,Y,N,Y,N,N,N)
 OUTPUT(AVEF,Y,N,Y,N,N,N)
 OUTPUT(MIXL,Y,N,Y,N,N,N)
 OUTPUT(RATE,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,N,N,N,N,N,N)
 OUTPUT(NPOR,N,N,N,N,N,N)
 OUTPUT(EPOR,N,N,N,N,N,N)
 OUTPUT(VPOR,N,N,N,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 = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 500 ;NUMCLS = 5
 NYPRIN = 1 ;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 = 1 ;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. )
 
 PATCH(PROF1 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF1 ,F1 ,0. ,0. )
 
 PATCH(PROF2 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF2 ,F2 ,0. ,0. )
 
 PATCH(PROF3 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF3 ,F3 ,0. ,0. )
 
 PATCH(PROF4 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF4 ,F4 ,0. ,0. )
 
 PATCH(PROF5 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF5 ,F5 ,0. ,0. )
 
 PATCH(PROF6 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF6 ,F6 ,0. ,0. )
 
 PATCH(PROF7 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF7 ,F7 ,0. ,0. )
 
 PATCH(PROF8 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF8 ,F8 ,0. ,0. )
 
 PATCH(PROF9 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF9 ,F9 ,0. ,0. )
 
 PATCH(PROF10 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF10 ,F10 ,0. ,0. )
 
 PATCH(PROF11 ,PROFIL, 1, 1, 1, 10, 1, 13, 1, 1)
 PLOT(PROF11 ,F11 ,0. ,0. )
 
 PATCH(MIDDLE ,PROFIL, 1, 1, 5, 5, 1, 13, 1, 1)
 PLOT(MIDDLE ,MNSQ,0. ,0. )
 PLOT(MIDDLE ,AVEF,0. ,0. )
 PLOT(MIDDLE ,ENUT,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP