TALK=T;RUN(1,1)
  DISPLAY
  The problem considered is an axisymmetric sonic jet discharging
  into stagnant surroundings from a nozzle at a pressure 3.56
  times higher than the ambient pressure. The stagnation enthalpy of
  the nozzle fluid is equal to that of the free stream, so that with
  the assumption of unit Prandtl numbers the energy equation need
  not be solved. The turbulence is represented by means of the k-e
  turbulence model. The calculation may be made with the van-Leer
  higher-order scheme or the hybrid differencing scheme. A
  relatively coarse mesh of 30 radial by 90 axial cells is used in
  the computations.
  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
     pause
   ENDUSE
   AUTOPLOT USE
     FILE
     PHI 5

     D 1 MACH Y 1;PLOT;REDR;LEVEL Y 1;level x 1.54
     msg Mach number distribution along flow axis
     msg Press  to continue
     pause
   ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(2D TRANSONIC UNDEREXPANDED JET: N113    )
 
 ************************************************************
 ************************************************************
 
 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 = 30
 YVLAST =1.
 YFRAC(1)=0.025 ;YFRAC(2)=0.05
 YFRAC(3)=0.075 ;YFRAC(4)=0.1
 YFRAC(5)=0.125 ;YFRAC(6)=0.15
 YFRAC(7)=0.175 ;YFRAC(8)=0.2
 YFRAC(9)=0.225 ;YFRAC(10)=0.25
 YFRAC(11)=0.275 ;YFRAC(12)=0.3
 YFRAC(13)=0.325 ;YFRAC(14)=0.35
 YFRAC(15)=0.375 ;YFRAC(16)=0.4
 YFRAC(17)=0.425 ;YFRAC(18)=0.45
 YFRAC(19)=0.475 ;YFRAC(20)=0.5
 YFRAC(21)=0.525059 ;YFRAC(22)=0.561703
 YFRAC(23)=0.604527 ;YFRAC(24)=0.651932
 YFRAC(25)=0.703063 ;YFRAC(26)=0.757375
 YFRAC(27)=0.814483 ;YFRAC(28)=0.874099
 YFRAC(29)=0.935999 ;YFRAC(30)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 90
 ZWLAST =5.
 ZFRAC(1)=7.692308E-03 ;ZFRAC(2)=0.015385
 ZFRAC(3)=0.023077 ;ZFRAC(4)=0.030769
 ZFRAC(5)=0.038462 ;ZFRAC(6)=0.046154
 ZFRAC(7)=0.053846 ;ZFRAC(8)=0.061538
 ZFRAC(9)=0.069231 ;ZFRAC(10)=0.076923
 ZFRAC(11)=0.084615 ;ZFRAC(12)=0.092308
 ZFRAC(13)=0.1 ;ZFRAC(14)=0.107692
 ZFRAC(15)=0.115385 ;ZFRAC(16)=0.123077
 ZFRAC(17)=0.130769 ;ZFRAC(18)=0.138462
 ZFRAC(19)=0.146154 ;ZFRAC(20)=0.153846
 ZFRAC(21)=0.161538 ;ZFRAC(22)=0.169231
 ZFRAC(23)=0.176923 ;ZFRAC(24)=0.184615
 ZFRAC(25)=0.192308 ;ZFRAC(26)=0.2
 ZFRAC(27)=0.207692 ;ZFRAC(28)=0.215385
 ZFRAC(29)=0.223077 ;ZFRAC(30)=0.230769
 ZFRAC(31)=0.238462 ;ZFRAC(32)=0.246154
 ZFRAC(33)=0.253846 ;ZFRAC(34)=0.261538
 ZFRAC(35)=0.269231 ;ZFRAC(36)=0.276923
 ZFRAC(37)=0.284615 ;ZFRAC(38)=0.292308
 ZFRAC(39)=0.3 ;ZFRAC(40)=0.307692
 ZFRAC(41)=0.315385 ;ZFRAC(42)=0.323077
 ZFRAC(43)=0.330769 ;ZFRAC(44)=0.338462
 ZFRAC(45)=0.346154 ;ZFRAC(46)=0.353846
 ZFRAC(47)=0.361538 ;ZFRAC(48)=0.369231
 ZFRAC(49)=0.376923 ;ZFRAC(50)=0.384615
 ZFRAC(51)=0.392308 ;ZFRAC(52)=0.4
 ZFRAC(53)=0.407692 ;ZFRAC(54)=0.415385
 ZFRAC(55)=0.423077 ;ZFRAC(56)=0.430769
 ZFRAC(57)=0.438462 ;ZFRAC(58)=0.446154
 ZFRAC(59)=0.453846 ;ZFRAC(60)=0.461538
 ZFRAC(61)=0.469231 ;ZFRAC(62)=0.476923
 ZFRAC(63)=0.484615 ;ZFRAC(64)=0.492308
 ZFRAC(65)=0.5 ;ZFRAC(66)=0.511636
 ZFRAC(67)=0.526733 ;ZFRAC(68)=0.543487
 ZFRAC(69)=0.561416 ;ZFRAC(70)=0.580274
 ZFRAC(71)=0.599906 ;ZFRAC(72)=0.620207
 ZFRAC(73)=0.641098 ;ZFRAC(74)=0.662519
 ZFRAC(75)=0.684422 ;ZFRAC(76)=0.706768
 ZFRAC(77)=0.729525 ;ZFRAC(78)=0.752664
 ZFRAC(79)=0.776163 ;ZFRAC(80)=0.8
 ZFRAC(81)=0.812619 ;ZFRAC(82)=0.828991
 ZFRAC(83)=0.84716 ;ZFRAC(84)=0.866604
 ZFRAC(85)=0.887055 ;ZFRAC(86)=0.908346
 ZFRAC(87)=0.930361 ;ZFRAC(88)=0.953016
 ZFRAC(89)=0.976247 ;ZFRAC(90)=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(146)=TMP1 ;NAME(147)=VABS
 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,N,N,N,N,Y)
 SOLUTN(TMP1,Y,N,N,N,N,Y)
 SOLUTN(VABS,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
 ************************************************************
  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)
 DIFCUT =0. ;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 =0.
 DVO1DT =0. ;DRH1DP = GRND5
 RHO1A =0. ;RHO1B =1.
 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 =2.0E-05 ;ENUT = GRND3
 ENUTA =0. ;ENUTB =0.
 ENUTC =0.
 IENUTA = 0
 PRNDTL(V1)=1. ;PRNDTL(W1)=1.
 PRNDTL(KE)=1. ;PRNDTL(EP)=1.
 PRT(V1)=1. ;PRT(W1)=1.
 PRT(KE)=1. ;PRT(EP)=1.314
 CP1 =3.5 ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=0.338275 ;FIINIT(V1)=0.
 FIINIT(W1)=1.080123 ;FIINIT(KE)=2.916666E-03
 FIINIT(EP)=5.176041E-04 ;FIINIT(H1)=3.5
 FIINIT(TMP1)=1.0E-10 ;FIINIT(VABS)=1.0E-10
 FIINIT(ENUT)=1.0E-10 ;FIINIT(MACH)=1.0E-10
 FIINIT(RHO1)=0.316969
   No PATCHes yet used for this Group
 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(HOCS ,CELL , 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(HOCS ,V1 , FIXFLU , GRND2 )
 COVAL(HOCS ,W1 , FIXFLU , GRND2 )
 COVAL(HOCS ,KE , FIXFLU , GRND2 )
 COVAL(HOCS ,EP , FIXFLU , GRND2 )
 
 PATCH(IN ,LOW , 1, 1, 1, 15, 1, 1, 1, 1)
 COVAL(IN ,P1 , FIXFLU ,0.684731 )
 COVAL(IN ,V1 ,0. ,0. )
 COVAL(IN ,W1 ,0. ,1.080123 )
 COVAL(IN ,KE ,0. ,2.916666E-03 )
 COVAL(IN ,EP ,0. ,5.176041E-04 )
 
 PATCH(LB ,LOW , 1, 1, 16, 30, 1, 1, 1, 1)
 COVAL(LB ,P1 ,1000. ,0.148269 )
 COVAL(LB ,KE ,0. ,1.0E-10 )
 COVAL(LB ,EP ,0. ,4.5E-15 )
 
 PATCH(NB ,NORTH , 1, 1, 30, 30, 1, 90, 1, 1)
 COVAL(NB ,P1 ,1000. ,0.148269 )
 COVAL(NB ,KE ,0. ,1.0E-10 )
 COVAL(NB ,EP ,0. ,4.5E-15 )
 
 PATCH(OUT ,HIGH , 1, 1, 1, 30, 90, 90, 1, 1)
 COVAL(OUT ,P1 ,1000. ,0.148269 )
 COVAL(OUT ,KE ,0. , SAME )
 COVAL(OUT ,EP ,0. , SAME )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 1600 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-04
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=10 ;LITER(V1)=10
 LITER(W1)=10 ;LITER(KE)=20
 LITER(EP)=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
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,0.7)
 RELAX(V1,FALSDT,0.012859)
 RELAX(W1,FALSDT,0.012859)
 RELAX(KE,FALSDT,0.012859)
 RELAX(EP,FALSDT,0.012859)
 RELAX(H1,LINRLX,1.)
 RELAX(TMP1,LINRLX,1.)
 RELAX(VABS,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 RELAX(MACH,LINRLX,1.)
 RELAX(RHO1,LINRLX,1.)
 KELIN = 1
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=1.0E-10
 VARMAX(V1)=100. ;VARMIN(V1)=-100.
 VARMAX(W1)=100. ;VARMIN(W1)=-100.
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(H1)=3.5 ;VARMIN(H1)=1.0E-10
 VARMAX(TMP1)=1.0E+10 ;VARMIN(TMP1)=-1.0E+10
 VARMAX(VABS)=1.0E+10 ;VARMIN(VABS)=-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-10
 ************************************************************
  Group 19. Data transmitted to GROUND
 NAMGRD =CONV
 GENK = T
 PARSOL = F
 ISG62 = 1
 SPEDAT(SET,SCHEME,INLCS005,I,8)
 SPEDAT(SET,SCHEME,INLCS007,I,8)
 SPEDAT(SET,SCHEME,INLCS012,I,8)
 SPEDAT(SET,SCHEME,INLCS013,I,8)
 ************************************************************
  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,N,N,N)
 OUTPUT(TMP1,Y,N,Y,N,N,N)
 OUTPUT(VABS,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 = 1 ;IZMON = 78
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 1600 ;NUMCLS = 5
 NYPRIN = 1 ;IYPRF = 1 ;IYPRL = 5
 NZPRIN = 1 ;IZPRF = 43 ;IZPRL = 55
 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
 
 PATCH(PLOT4 ,PROFIL, 1, 1, 1, 1, 1, 90, 1, 1)
 PLOT(PLOT4 ,W1 ,0. ,0. )
 
 PATCH(PLOT5 ,PROFIL, 1, 1, 1, 1, 1, 90, 1, 1)
 PLOT(PLOT5 ,MACH,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP