TALK=T;RUN(1,1)
  DISPLAY
  The problem considered is an axisymmetric, supersonic, rocket
  exhaust plume discharging into a subsonic moving air stream.

  The exhaust gas has a specific-heat ratio of 1.21, a Mach number
  of 2.39, a total pressure of 5.786MPa, and a total temperature
  of 3462 K.

  The free-stream has a Mach number of 0.029, a static temperature
  of 288K, and a static pressure of 1.01325 bar.

  The plume is underexpanded with a nozzle-to-ambient static
  pressure ratio of 3.8.

  In addition to the continuity and momentum equations,
  conservation equations are solved for the total enthalpy equation
  and the mass fraction of the exhaust gas.

  The turbulence is represented by means of the k-e turbulence mode
  The results of this elliptic-solver case may be compared with
  those produced for library case 161 by the much more economical
  parabolic solver.

  ENDDIS
   PHOTON USE
     P


       0.20443E+04 0.15633E+04 CR
     CON MACH X 1 FI;.5
     msg Mach number contours
     msg Press  to continue
     pau;cl
     vec x 1 sh;
     pau;cl
     con tmp1 x 1 fi;.5
     pau;cl
     con rho1 x 1 fi;.5
     pau;cl
     con c1 x 1 fi;.5
     pause
   ENDUSE
   AUTOPLOT USE
     FILE
     PHI 5

     D 1 MACH Y 1;PLOT;REDR;LEVEL Y 1
     msg Mach number distribution along flow axis
     msg Press  to continue
     pause
     cl
     D 1 P1 Y 1;PLOT
     shift y
     1.01325e5 1
     scale;redr
     msg static pressure distribution along flow axis
     msg Press  to continue
     pause
     cl
     D 1 RHO1 Y 1;PLOT;REDR
     msg static density distribution along flow axis
     msg Press  to continue
     pause
     cl

   ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(2D Supersonic Rocket Exhaust Plume      )
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 14
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = F
 NX = 1
 XULAST =0.1
 XFRAC(1)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 40
 YVLAST =0.181625
 YFRAC(1)=1.0E-02 ;YFRAC(2)=0.02
 YFRAC(3)=0.03 ;YFRAC(4)=0.04
 YFRAC(5)=0.05 ;YFRAC(6)=0.06
 YFRAC(7)=0.07 ;YFRAC(8)=0.08
 YFRAC(9)=0.09 ;YFRAC(10)=0.1
 YFRAC(11)=0.11 ;YFRAC(12)=0.12
 YFRAC(13)=0.13 ;YFRAC(14)=0.14
 YFRAC(15)=0.15 ;YFRAC(16)=0.16
 YFRAC(17)=0.17 ;YFRAC(18)=0.18
 YFRAC(19)=0.19 ;YFRAC(20)=0.2
 YFRAC(21)=0.212068 ;YFRAC(22)=0.231849
 YFRAC(23)=0.256185 ;YFRAC(24)=0.284049
 YFRAC(25)=0.31487 ;YFRAC(26)=0.348272
 YFRAC(27)=0.383986 ;YFRAC(28)=0.421806
 YFRAC(29)=0.46157 ;YFRAC(30)=0.503143
 YFRAC(31)=0.546416 ;YFRAC(32)=0.591293
 YFRAC(33)=0.637692 ;YFRAC(34)=0.685542
 YFRAC(35)=0.734781 ;YFRAC(36)=0.78535
 YFRAC(37)=0.837201 ;YFRAC(38)=0.890287
 YFRAC(39)=0.944566 ;YFRAC(40)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 80
 ZWLAST =0.5812
 ZFRAC(1)=0.0125 ;ZFRAC(2)=0.025
 ZFRAC(3)=0.0375 ;ZFRAC(4)=0.05
 ZFRAC(5)=0.0625 ;ZFRAC(6)=0.075
 ZFRAC(7)=0.0875 ;ZFRAC(8)=0.1
 ZFRAC(9)=0.1125 ;ZFRAC(10)=0.125
 ZFRAC(11)=0.1375 ;ZFRAC(12)=0.15
 ZFRAC(13)=0.1625 ;ZFRAC(14)=0.175
 ZFRAC(15)=0.1875 ;ZFRAC(16)=0.2
 ZFRAC(17)=0.2125 ;ZFRAC(18)=0.225
 ZFRAC(19)=0.2375 ;ZFRAC(20)=0.25
 ZFRAC(21)=0.2625 ;ZFRAC(22)=0.275
 ZFRAC(23)=0.2875 ;ZFRAC(24)=0.3
 ZFRAC(25)=0.3125 ;ZFRAC(26)=0.325
 ZFRAC(27)=0.3375 ;ZFRAC(28)=0.35
 ZFRAC(29)=0.3625 ;ZFRAC(30)=0.375
 ZFRAC(31)=0.3875 ;ZFRAC(32)=0.4
 ZFRAC(33)=0.4125 ;ZFRAC(34)=0.425
 ZFRAC(35)=0.4375 ;ZFRAC(36)=0.45
 ZFRAC(37)=0.4625 ;ZFRAC(38)=0.475
 ZFRAC(39)=0.4875 ;ZFRAC(40)=0.5
 ZFRAC(41)=0.5125 ;ZFRAC(42)=0.525
 ZFRAC(43)=0.5375 ;ZFRAC(44)=0.55
 ZFRAC(45)=0.5625 ;ZFRAC(46)=0.575
 ZFRAC(47)=0.5875 ;ZFRAC(48)=0.6
 ZFRAC(49)=0.6125 ;ZFRAC(50)=0.625
 ZFRAC(51)=0.6375 ;ZFRAC(52)=0.65
 ZFRAC(53)=0.6625 ;ZFRAC(54)=0.675
 ZFRAC(55)=0.6875 ;ZFRAC(56)=0.7
 ZFRAC(57)=0.7125 ;ZFRAC(58)=0.725
 ZFRAC(59)=0.7375 ;ZFRAC(60)=0.75
 ZFRAC(61)=0.7625 ;ZFRAC(62)=0.775
 ZFRAC(63)=0.7875 ;ZFRAC(64)=0.8
 ZFRAC(65)=0.8125 ;ZFRAC(66)=0.825
 ZFRAC(67)=0.8375 ;ZFRAC(68)=0.85
 ZFRAC(69)=0.8625 ;ZFRAC(70)=0.875
 ZFRAC(71)=0.8875 ;ZFRAC(72)=0.9
 ZFRAC(73)=0.9125 ;ZFRAC(74)=0.925
 ZFRAC(75)=0.9375 ;ZFRAC(76)=0.95
 ZFRAC(77)=0.9625 ;ZFRAC(78)=0.975
 ZFRAC(79)=0.9875 ;ZFRAC(80)=1.
 ************************************************************
  Group 6. Body-Fitted Coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
 ONEPHS = T
 NAME(1)=P1 ;NAME(5)=V1
 NAME(7)=W1 ;NAME(12)=KE
 NAME(13)=EP ;NAME(14)=H1
 NAME(16)=C1 ;NAME(144)=EPKE
 NAME(145)=CP1 ;NAME(146)=TMP1
 NAME(147)=LEN1 ;NAME(148)=ENUT
 NAME(149)=MACH ;NAME(150)=RHO1
    * 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(V1,Y,Y,N,N,N,N)
 SOLUTN(W1,Y,Y,N,N,N,N)
 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(C1,Y,Y,N,N,N,Y)
 SOLUTN(EPKE,Y,N,N,N,N,N)
 SOLUTN(CP1,Y,N,N,N,N,Y)
 SOLUTN(TMP1,Y,N,N,N,N,Y)
 SOLUTN(LEN1,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 SOLUTN(MACH,Y,N,N,N,N,Y)
 SOLUTN(RHO1,Y,N,N,N,N,Y)
 DEN1 = 150
 VIST = 148
 TEMP1 = 146
 LEN1 = 147
 ************************************************************
  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(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,Y,Y,Y)
 TERMS(C1,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWRH1 = T
 NEWENT = T
 UCONV = 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 = GRND5 ;TMP1 = GRND6
 EL1 = GRND4
 TSURR =0. ;TEMP0 =0.
 PRESS0 =1.01325E+05
 DVO1DT =0. ;DRH1DP = GRND5
 RHO1A =289.096985 ;RHO1B =337.162598
 RHO1C =0.714286
 TMP1A =0. ;TMP1B =0.
 TMP1C =0.
 TMP2A =0. ;TMP2B =0.
 TMP2C =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0. ;EL1B =0.
 EL1C =0.
 ENUL =1.0E-05 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0.
 ENUTC =0.
 IENUTA = 0
 PRNDTL(V1)=1. ;PRNDTL(W1)=1.
 PRNDTL(KE)=1. ;PRNDTL(EP)=1.
 PRNDTL(H1)=1. ;PRNDTL(C1)=1.
 PRT(V1)=1. ;PRT(W1)=1.
 PRT(KE)=1. ;PRT(EP)=1.314
 PRT(H1)=1. ;PRT(C1)=1.
 CP1 = GRND7 ;CP2 =1.
 CP1A =1011.839478 ;CP1B =1942.698486
 CP1C =0. ;CP1D =0.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=2.83775E+05 ;FIINIT(V1)=0.
 FIINIT(W1)=2246.340088 ;FIINIT(KE)=5.046044E+04
 FIINIT(EP)=5.126943E+08 ;FIINIT(H1)=6.725079E+06
 FIINIT(C1)=1. ;FIINIT(EPKE)=1.0E-10
 FIINIT(CP1)=1.0E-10 ;FIINIT(TMP1)=2163.
 FIINIT(LEN1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
 FIINIT(MACH)=1.0E-10 ;FIINIT(RHO1)=0.528053
 
 PATCH(INITA ,INIVAL, 1, 1, 21, 40, 1, 80, 1, 1)
 INIT(INITA ,P1 ,0. ,0. )
 INIT(INITA ,W1 ,0. ,10. )
 INIT(INITA ,KE ,0. ,1.0E-10 )
 INIT(INITA ,EP ,0. ,9.000001E-15 )
 INIT(INITA ,H1 ,0. ,2.914598E+05 )
 INIT(INITA ,C1 ,0. ,0. )
 INIT(INITA ,TMP1,0. ,288. )
 INIT(INITA ,RHO1,0. ,1.216972 )
 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(IN ,LOW , 1, 1, 1, 20, 1, 1, 1, 1)
 COVAL(IN ,P1 , FIXFLU ,1186.18689 )
 COVAL(IN ,V1 ,0. ,0. )
 COVAL(IN ,W1 ,0. ,2246.340088 )
 COVAL(IN ,KE ,0. ,5.046044E+04 )
 COVAL(IN ,EP ,0. ,5.126943E+08 )
 COVAL(IN ,H1 ,0. ,6.725079E+06 )
 COVAL(IN ,C1 ,0. ,1. )
 
 PATCH(LB ,LOW , 1, 1, 21, 40, 1, 1, 1, 1)
 COVAL(LB ,P1 , FIXFLU ,12.16972 )
 COVAL(LB ,W1 ,0. ,10. )
 COVAL(LB ,KE ,0. ,1.0E-10 )
 COVAL(LB ,EP ,0. ,9.000001E-15 )
 COVAL(LB ,H1 ,0. ,2.914598E+05 )
 
 PATCH(NB ,NORTH , 1, 1, 40, 40, 1, 80, 1, 1)
 COVAL(NB ,P1 ,1000. ,0. )
 COVAL(NB ,W1 ,0. ,10. )
 COVAL(NB ,KE ,0. ,1.0E-10 )
 COVAL(NB ,EP ,0. ,9.000001E-15 )
 COVAL(NB ,H1 ,0. ,2.914598E+05 )
 COVAL(NB ,C1 ,0. ,0. )
 
 PATCH(OUT ,HIGH , 1, 1, 1, 40, 80, 80, 1, 1)
 COVAL(OUT ,P1 ,1000. ,0. )
 COVAL(OUT ,KE ,0. , SAME )
 COVAL(OUT ,EP ,0. , SAME )
 COVAL(OUT ,H1 ,0. , SAME )
 COVAL(OUT ,C1 ,0. , SAME )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 700 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = F
 RESREF(P1)=7.825901E-14 ;RESREF(V1)=1.757964E-10
 RESREF(W1)=1.757964E-10 ;RESREF(KE)=3.948984E-09
 RESREF(EP)=4.012295E-05 ;RESREF(H1)=5.26298E-07
 RESREF(C1)=7.825901E-14
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=20 ;LITER(V1)=10
 LITER(W1)=10 ;LITER(KE)=20
 LITER(EP)=20 ;LITER(H1)=20
 LITER(C1)=20
 ENDIT(P1)=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(C1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(V1,FALSDT,4.56098E-06)
 RELAX(W1,FALSDT,4.56098E-06)
 RELAX(KE,FALSDT,4.56098E-06)
 RELAX(EP,FALSDT,4.56098E-06)
 RELAX(H1,FALSDT,4.56098E-06)
 RELAX(C1,LINRLX,1.)
 RELAX(EPKE,LINRLX,1.)
 RELAX(CP1,LINRLX,1.)
 RELAX(TMP1,LINRLX,1.)
 RELAX(LEN1,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 RELAX(MACH,LINRLX,1.)
 RELAX(RHO1,LINRLX,1.)
 KELIN = 3
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 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(C1)=1. ;VARMIN(C1)=1.0E-10
 VARMAX(EPKE)=1.0E+10 ;VARMIN(EPKE)=-1.0E+10
 VARMAX(CP1)=1.0E+10 ;VARMIN(CP1)=-1.0E+10
 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=200.
 VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 VARMAX(MACH)=1.0E+10 ;VARMIN(MACH)=-1.0E+10
 VARMAX(RHO1)=1.0E+10 ;VARMIN(RHO1)=1.0E-05
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMGRD =CONV
 GENK = T
 PARSOL = F
 ISG62 = 1
 ************************************************************
  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(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(C1,Y,N,Y,Y,Y,Y)
 OUTPUT(EPKE,Y,N,Y,N,N,N)
 OUTPUT(CP1,Y,N,Y,N,N,N)
 OUTPUT(TMP1,Y,N,Y,N,N,N)
 OUTPUT(LEN1,Y,N,Y,N,N,N)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 OUTPUT(MACH,Y,N,Y,N,Y,N)
 OUTPUT(RHO1,Y,N,Y,N,Y,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 20 ;IZMON = 40
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 700 ;NUMCLS = 5
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 IPLTF = 1 ;IPLTL = -1 ;NPLT = 20
 ISWPRF = 1 ;ISWPRL = 100000
 ITABL = 3 ;IPROF = 1
 ABSIZ =0.5 ;ORSIZ =0.4
 NTZPRF = 1 ;NCOLPF = 50
 ICHR = 2 ;NCOLCO = 45 ;NROWCO = 20
   No PATCHes yet used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP