TALK=T;RUN(1,1)
  DISPLAY
  The problem considered is the calculation of incompressible
  turbulent flow and heat transfer in a sudden pipe expansion
  with a diameter ratio Do/Di=2.5. The Reynolds number of the
  larger pipe is 4.075E4, and the laminar Prandtl number is 0.7.
  Downstream of the expansion, the fluid is heated by a uniform
  heat flux through the outer wall. This case was studied
  experimentally by Baugh et al [1984] to determine the local
  Nusselt number distributions for various Reynolds numbers.

  Following the expansion, the flow enters the larger pipe in the
  form of a circular jet, and then separates from the expansion
  corner to generate to a primary recirculation zone, whose reverse
  flow is driven by the adverse pressure gradient associated with
  the expansion geometry. At the point of reattachment, the heat
  transfer coefficients (htc) are known to be several times greater
  than the corresponding fully-developed values. For this case, the
  measurements show that the peak htc downstream of the expansion
  is roughly 4 times the fully-developed htc for the same Reynolds
  number. These high htcs occur in the region of the reattachment
  of the shear layer to the tube wall.

  The calculations are started at z= -4h and terminated at z=41h,
  where a fixed-pressure boundary condition is applied. Here, h is
  the step height given by h=0.5*(Do-Di).The calculation may be
  performed with one of the following low-Re turbulence models: the
  Lam-Bremhorst, Chen-Kim k-e & 2-Layer k-e models; the Wilcox
  1988 & 2008 k-w models; and the Menter & SST k-w models. The
  turbulent Prandtl number is taken as 0.86. The numerical
  integration is taken right down to the wall, and a non-uniform
  mesh of NY=110 by NX=200 is used which concentrates grid cells
  close to the wall. Typically, about 5000 sweeps are required
  for a converged solution with this mesh density.
  ENDDIS
   AUTOPLOT USE
   file
   phida 3


   da 1 w1 y m
   shift x -1.2 1
   divide x .3 1
   col9 1
   level y 0;level x 0
   scale x 0 40
   msg Velocity (W1) profile
   pause
   cl
   da 1 ypls y m
   shift x -1.2 1
   divide x .3 1
   col3 1
   level y 0;level x 0
   scale x 0 40
   pause
   cl
   da 1 ypls y 50
   divide x .3 1
   col3 1
   level y 0;level x 0
   scale x 0 5
   pause
   cl
   da 1 nuss y m
   shift x -1.2 1
   divide x .3 1
   plot 1
   level y 0;level x 0
   scale x 0 40
   pause
   cl
   da 1 nusc y m
   shift x -1.2 1
   divide x .3 1
   plot 1
   level y 0;level x 0
   scale x 0 40
   msg Press e to END
   ENDUSE
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(SST k-w__2D ABRUPT PIPE EXPANS  :T213   )
 
 ************************************************************
 ************************************************************
 
 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 = 110
 YVLAST =0.5
 YFRAC(1)=0.030275 ;YFRAC(2)=0.058308
 YFRAC(3)=0.084264 ;YFRAC(4)=0.108297
 YFRAC(5)=0.13055 ;YFRAC(6)=0.151155
 YFRAC(7)=0.170233 ;YFRAC(8)=0.187899
 YFRAC(9)=0.204255 ;YFRAC(10)=0.219401
 YFRAC(11)=0.233424 ;YFRAC(12)=0.246408
 YFRAC(13)=0.258431 ;YFRAC(14)=0.269563
 YFRAC(15)=0.279871 ;YFRAC(16)=0.289415
 YFRAC(17)=0.298252 ;YFRAC(18)=0.306434
 YFRAC(19)=0.31401 ;YFRAC(20)=0.321025
 YFRAC(21)=0.327521 ;YFRAC(22)=0.333535
 YFRAC(23)=0.339104 ;YFRAC(24)=0.34426
 YFRAC(25)=0.349035 ;YFRAC(26)=0.353455
 YFRAC(27)=0.357549 ;YFRAC(28)=0.361339
 YFRAC(29)=0.364848 ;YFRAC(30)=0.368097
 YFRAC(31)=0.371106 ;YFRAC(32)=0.373892
 YFRAC(33)=0.376471 ;YFRAC(34)=0.37886
 YFRAC(35)=0.381071 ;YFRAC(36)=0.383119
 YFRAC(37)=0.385015 ;YFRAC(38)=0.38677
 YFRAC(39)=0.388396 ;YFRAC(40)=0.389901
 YFRAC(41)=0.391294 ;YFRAC(42)=0.392585
 YFRAC(43)=0.39378 ;YFRAC(44)=0.394886
 YFRAC(45)=0.39591 ;YFRAC(46)=0.396859
 YFRAC(47)=0.397737 ;YFRAC(48)=0.39855
 YFRAC(49)=0.399303 ;YFRAC(50)=0.4
 YFRAC(51)=0.401824 ;YFRAC(52)=0.40383
 YFRAC(53)=0.406037 ;YFRAC(54)=0.408464
 YFRAC(55)=0.411134 ;YFRAC(56)=0.414072
 YFRAC(57)=0.417302 ;YFRAC(58)=0.420856
 YFRAC(59)=0.424766 ;YFRAC(60)=0.429066
 YFRAC(61)=0.433797 ;YFRAC(62)=0.439
 YFRAC(63)=0.444724 ;YFRAC(64)=0.45102
 YFRAC(65)=0.457946 ;YFRAC(66)=0.465564
 YFRAC(67)=0.473944 ;YFRAC(68)=0.483163
 YFRAC(69)=0.493303 ;YFRAC(70)=0.504457
 YFRAC(71)=0.516726 ;YFRAC(72)=0.530223
 YFRAC(73)=0.545069 ;YFRAC(74)=0.561399
 YFRAC(75)=0.579363 ;YFRAC(76)=0.599123
 YFRAC(77)=0.620859 ;YFRAC(78)=0.644769
 YFRAC(79)=0.671069 ;YFRAC(80)=0.7
 YFRAC(81)=0.728931 ;YFRAC(82)=0.755231
 YFRAC(83)=0.779141 ;YFRAC(84)=0.800877
 YFRAC(85)=0.820637 ;YFRAC(86)=0.838601
 YFRAC(87)=0.854931 ;YFRAC(88)=0.869778
 YFRAC(89)=0.883274 ;YFRAC(90)=0.895543
 YFRAC(91)=0.906697 ;YFRAC(92)=0.916837
 YFRAC(93)=0.926056 ;YFRAC(94)=0.934436
 YFRAC(95)=0.942054 ;YFRAC(96)=0.94898
 YFRAC(97)=0.955276 ;YFRAC(98)=0.961
 YFRAC(99)=0.966203 ;YFRAC(100)=0.970934
 YFRAC(101)=0.975234 ;YFRAC(102)=0.979144
 YFRAC(103)=0.982698 ;YFRAC(104)=0.985929
 YFRAC(105)=0.988866 ;YFRAC(106)=0.991536
 YFRAC(107)=0.993963 ;YFRAC(108)=0.99617
 YFRAC(109)=0.998176 ;YFRAC(110)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 200
 ZWLAST =13.5
 ZFRAC(1)=4.762894E-03 ;ZFRAC(2)=9.298985E-03
 ZFRAC(3)=0.013619 ;ZFRAC(4)=0.017733
 ZFRAC(5)=0.021652 ;ZFRAC(6)=0.025384
 ZFRAC(7)=0.028938 ;ZFRAC(8)=0.032323
 ZFRAC(9)=0.035547 ;ZFRAC(10)=0.038617
 ZFRAC(11)=0.041541 ;ZFRAC(12)=0.044325
 ZFRAC(13)=0.046978 ;ZFRAC(14)=0.049504
 ZFRAC(15)=0.051909 ;ZFRAC(16)=0.0542
 ZFRAC(17)=0.056382 ;ZFRAC(18)=0.05846
 ZFRAC(19)=0.060439 ;ZFRAC(20)=0.062324
 ZFRAC(21)=0.064119 ;ZFRAC(22)=0.065829
 ZFRAC(23)=0.067457 ;ZFRAC(24)=0.069008
 ZFRAC(25)=0.070484 ;ZFRAC(26)=0.071891
 ZFRAC(27)=0.07323 ;ZFRAC(28)=0.074506
 ZFRAC(29)=0.075721 ;ZFRAC(30)=0.076878
 ZFRAC(31)=0.07798 ;ZFRAC(32)=0.07903
 ZFRAC(33)=0.080029 ;ZFRAC(34)=0.080981
 ZFRAC(35)=0.081888 ;ZFRAC(36)=0.082752
 ZFRAC(37)=0.083574 ;ZFRAC(38)=0.084357
 ZFRAC(39)=0.085103 ;ZFRAC(40)=0.085813
 ZFRAC(41)=0.08649 ;ZFRAC(42)=0.087134
 ZFRAC(43)=0.087748 ;ZFRAC(44)=0.088332
 ZFRAC(45)=0.088889 ;ZFRAC(46)=0.089088
 ZFRAC(47)=0.0893 ;ZFRAC(48)=0.089524
 ZFRAC(49)=0.089761 ;ZFRAC(50)=0.090013
 ZFRAC(51)=0.09028 ;ZFRAC(52)=0.090563
 ZFRAC(53)=0.090863 ;ZFRAC(54)=0.09118
 ZFRAC(55)=0.091517 ;ZFRAC(56)=0.091875
 ZFRAC(57)=0.092253 ;ZFRAC(58)=0.092654
 ZFRAC(59)=0.09308 ;ZFRAC(60)=0.093531
 ZFRAC(61)=0.094008 ;ZFRAC(62)=0.094515
 ZFRAC(63)=0.095052 ;ZFRAC(64)=0.095621
 ZFRAC(65)=0.096225 ;ZFRAC(66)=0.096864
 ZFRAC(67)=0.097542 ;ZFRAC(68)=0.098261
 ZFRAC(69)=0.099022 ;ZFRAC(70)=0.09983
 ZFRAC(71)=0.100686 ;ZFRAC(72)=0.101593
 ZFRAC(73)=0.102555 ;ZFRAC(74)=0.103574
 ZFRAC(75)=0.104655 ;ZFRAC(76)=0.1058
 ZFRAC(77)=0.107014 ;ZFRAC(78)=0.108301
 ZFRAC(79)=0.109665 ;ZFRAC(80)=0.111111
 ZFRAC(81)=0.113294 ;ZFRAC(82)=0.115509
 ZFRAC(83)=0.117758 ;ZFRAC(84)=0.120041
 ZFRAC(85)=0.122357 ;ZFRAC(86)=0.124709
 ZFRAC(87)=0.127096 ;ZFRAC(88)=0.129518
 ZFRAC(89)=0.131977 ;ZFRAC(90)=0.134473
 ZFRAC(91)=0.137006 ;ZFRAC(92)=0.139577
 ZFRAC(93)=0.142187 ;ZFRAC(94)=0.144836
 ZFRAC(95)=0.147524 ;ZFRAC(96)=0.150253
 ZFRAC(97)=0.153023 ;ZFRAC(98)=0.155835
 ZFRAC(99)=0.158688 ;ZFRAC(100)=0.161585
 ZFRAC(101)=0.164525 ;ZFRAC(102)=0.167509
 ZFRAC(103)=0.170537 ;ZFRAC(104)=0.173612
 ZFRAC(105)=0.176732 ;ZFRAC(106)=0.179899
 ZFRAC(107)=0.183114 ;ZFRAC(108)=0.186376
 ZFRAC(109)=0.189688 ;ZFRAC(110)=0.19305
 ZFRAC(111)=0.196461 ;ZFRAC(112)=0.199924
 ZFRAC(113)=0.203439 ;ZFRAC(114)=0.207007
 ZFRAC(115)=0.210628 ;ZFRAC(116)=0.214304
 ZFRAC(117)=0.218035 ;ZFRAC(118)=0.221821
 ZFRAC(119)=0.225665 ;ZFRAC(120)=0.229566
 ZFRAC(121)=0.233525 ;ZFRAC(122)=0.237544
 ZFRAC(123)=0.241624 ;ZFRAC(124)=0.245764
 ZFRAC(125)=0.249967 ;ZFRAC(126)=0.254232
 ZFRAC(127)=0.258562 ;ZFRAC(128)=0.262956
 ZFRAC(129)=0.267417 ;ZFRAC(130)=0.271944
 ZFRAC(131)=0.276539 ;ZFRAC(132)=0.281204
 ZFRAC(133)=0.285938 ;ZFRAC(134)=0.290743
 ZFRAC(135)=0.29562 ;ZFRAC(136)=0.300571
 ZFRAC(137)=0.305595 ;ZFRAC(138)=0.310695
 ZFRAC(139)=0.315872 ;ZFRAC(140)=0.321126
 ZFRAC(141)=0.326459 ;ZFRAC(142)=0.331872
 ZFRAC(143)=0.337366 ;ZFRAC(144)=0.342943
 ZFRAC(145)=0.348603 ;ZFRAC(146)=0.354348
 ZFRAC(147)=0.36018 ;ZFRAC(148)=0.366098
 ZFRAC(149)=0.372106 ;ZFRAC(150)=0.378204
 ZFRAC(151)=0.384393 ;ZFRAC(152)=0.390675
 ZFRAC(153)=0.397051 ;ZFRAC(154)=0.403523
 ZFRAC(155)=0.410092 ;ZFRAC(156)=0.416759
 ZFRAC(157)=0.423527 ;ZFRAC(158)=0.430396
 ZFRAC(159)=0.437368 ;ZFRAC(160)=0.444444
 ZFRAC(161)=0.453642 ;ZFRAC(162)=0.463024
 ZFRAC(163)=0.472593 ;ZFRAC(164)=0.482354
 ZFRAC(165)=0.492309 ;ZFRAC(166)=0.502464
 ZFRAC(167)=0.512822 ;ZFRAC(168)=0.523388
 ZFRAC(169)=0.534164 ;ZFRAC(170)=0.545156
 ZFRAC(171)=0.556368 ;ZFRAC(172)=0.567804
 ZFRAC(173)=0.579469 ;ZFRAC(174)=0.591367
 ZFRAC(175)=0.603503 ;ZFRAC(176)=0.615882
 ZFRAC(177)=0.628508 ;ZFRAC(178)=0.641387
 ZFRAC(179)=0.654524 ;ZFRAC(180)=0.667923
 ZFRAC(181)=0.68159 ;ZFRAC(182)=0.695531
 ZFRAC(183)=0.70975 ;ZFRAC(184)=0.724254
 ZFRAC(185)=0.739048 ;ZFRAC(186)=0.754137
 ZFRAC(187)=0.769529 ;ZFRAC(188)=0.785228
 ZFRAC(189)=0.801241 ;ZFRAC(190)=0.817575
 ZFRAC(191)=0.834235 ;ZFRAC(192)=0.851229
 ZFRAC(193)=0.868562 ;ZFRAC(194)=0.886242
 ZFRAC(195)=0.904276 ;ZFRAC(196)=0.92267
 ZFRAC(197)=0.941432 ;ZFRAC(198)=0.960569
 ZFRAC(199)=0.98009 ;ZFRAC(200)=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(125)=OMRS
 NAME(126)=OMCR ;NAME(127)=PRPS
 NAME(128)=NUSC ;NAME(129)=NUSS
 NAME(130)=HTCB ;NAME(131)=TWAL
 NAME(132)=TAVE ;NAME(133)=CNH1
 NAME(134)=LEN1 ;NAME(135)=SKIN
 NAME(136)=SPH1 ;NAME(137)=ENUL
 NAME(138)=KOND ;NAME(139)=STRS
 NAME(140)=YPLS ;NAME(141)=LTLS
 NAME(142)=WDIS ;NAME(143)=GEN1
 NAME(144)=BF2 ;NAME(145)=BF1
 NAME(146)=OMEG ;NAME(147)=ENUT
 NAME(148)=TEM1 ;NAME(149)=TDIF
 NAME(150)=CP
    * 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,Y,N,N,N)
 SOLUTN(W1,Y,Y,Y,N,N,N)
 SOLUTN(KE,Y,Y,N,N,N,N)
 SOLUTN(EP,Y,N,N,N,N,N)
 SOLUTN(OMRS,Y,N,N,N,N,N)
 SOLUTN(OMCR,Y,N,N,N,N,N)
 SOLUTN(PRPS,Y,N,N,N,N,Y)
 SOLUTN(NUSC,Y,N,N,N,N,N)
 SOLUTN(NUSS,Y,N,N,N,N,N)
 SOLUTN(HTCB,Y,N,N,N,N,N)
 SOLUTN(TWAL,Y,N,N,N,N,N)
 SOLUTN(TAVE,Y,N,N,N,N,N)
 SOLUTN(CNH1,Y,N,N,N,N,N)
 SOLUTN(LEN1,Y,N,N,N,N,Y)
 SOLUTN(SKIN,Y,N,N,N,N,Y)
 SOLUTN(SPH1,Y,N,N,N,N,Y)
 SOLUTN(ENUL,Y,N,N,N,N,Y)
 SOLUTN(KOND,Y,N,N,N,N,Y)
 SOLUTN(STRS,Y,N,N,N,N,Y)
 SOLUTN(YPLS,Y,N,N,N,N,Y)
 SOLUTN(LTLS,Y,Y,Y,N,N,Y)
 SOLUTN(WDIS,Y,N,N,N,N,N)
 SOLUTN(GEN1,Y,N,N,N,N,N)
 SOLUTN(BF2,Y,N,N,N,N,N)
 SOLUTN(BF1,Y,N,N,N,N,N)
 SOLUTN(OMEG,Y,Y,N,N,N,N)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 SOLUTN(TEM1,Y,Y,Y,N,N,N)
 SOLUTN(TDIF,Y,N,N,N,N,N)
 SOLUTN(CP,Y,N,N,N,N,N)
 VISL = 137
 VIST = 147
 LEN1 = 134
 PRPS = 127
 ************************************************************
  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(LTLS,N,N,Y,N,Y,Y)
 TERMS(OMEG,N,Y,Y,Y,Y,N)
 TERMS(TEM1,N,Y,Y,Y,Y,Y)
 DIFCUT =0.5 ;ZDIFAC =1.
 GALA = F ;ADDDIF = F
 NEWENT = T
 ISOLX = -1 ;ISOLY = -1 ;ISOLZ = -1
 ************************************************************
  Group 9. Properties used if PRPS is not
  stored, and where PRPS = -1.0 if it is!
 RHO1 =1. ;TMP1 =0. ;EL1 = GRND5
 TSURR =0. ;TEMP0 =0. ;PRESS0 =0.
 DVO1DT =0. ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0.6 ;EL1B =0. ;EL1C =0.
 ENUL =3.926381E-05 ;ENUT = GRND7
 ENUTA =0. ;ENUTB =0. ;ENUTC =0.
 IENUTA = 20
 PRNDTL(V1)=1. ;PRNDTL(W1)=1.
 PRNDTL(KE)=GRND6 ;PRNDTL(LTLS)=1.
 PRNDTL(OMEG)=GRND6 ;PRNDTL(TEM1)=0.7
 PRT(V1)=1. ;PRT(W1)=1.
 PRT(KE)=1.0E+10 ;PRT(LTLS)=1.
 PRT(OMEG)=1.0E+10 ;PRT(TEM1)=0.86
 CP1 =1000. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=1.3E-04 ;FIINIT(V1)=0.
 FIINIT(W1)=10. ;FIINIT(KE)=0.44746
 FIINIT(EP)=1.0E-10 ;FIINIT(OMRS)=1.0E-10
 FIINIT(OMCR)=1.0E-10 ;FIINIT(PRPS)=-1.
 FIINIT(NUSC)=1.0E-10 ;FIINIT(NUSS)=1.0E-10
 FIINIT(HTCB)=1.0E-10 ;FIINIT(TWAL)=1.0E-10
 FIINIT(TAVE)=1.0E-10 ;FIINIT(CNH1)=1.0E-10
 FIINIT(LEN1)=1.0E-10 ;FIINIT(SKIN)=1.0E-10
 FIINIT(SPH1)=1.0E-10 ;FIINIT(ENUL)=1.0E-10
 FIINIT(KOND)=1.0E-10 ;FIINIT(STRS)=1.0E-10
 FIINIT(YPLS)=1.0E-10 ;FIINIT(LTLS)=1.0E-10
 FIINIT(WDIS)=1.0E-02 ;FIINIT(GEN1)=1.0E-10
 FIINIT(BF2)=1. ;FIINIT(BF1)=1.
 FIINIT(OMEG)=67.842171 ;FIINIT(ENUT)=1.0E-10
 FIINIT(TEM1)=20. ;FIINIT(TDIF)=1.0E-10
 FIINIT(CP)=1.0E-10
 
 PATCH(STEP ,INIVAL, 1, 1, 51, 110, 1, 45, 1, 1)
 INIT(STEP ,PRPS,0. ,198. )
 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(KWSOURCE,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(KWSOURCE,KE , GRND4 , GRND4 )
 COVAL(KWSOURCE,OMEG, GRND4 , GRND4 )
 
 PATCH(KWSOGD ,PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
 COVAL(KWSOGD ,OMEG, GRND4 , GRND4 )
 
 PATCH(PATCH1 ,CELL , 1, 1, 1, 110, 1, 200, 1, 1)
 
 PATCH(INLET ,LOW , 1, 1, 1, 50, 1, 1, 1, 1)
 COVAL(INLET ,P1 , FIXFLU ,10. )
 COVAL(INLET ,V1 ,0. ,0. )
 COVAL(INLET ,W1 ,0. ,10. )
 COVAL(INLET ,KE ,0. ,0.44746 )
 COVAL(INLET ,OMEG,0. ,67.842171 )
 COVAL(INLET ,TEM1,0. ,20. )
 
 PATCH(OUTLET ,HIGH , 1, 1, 1, 110, 200, 200, 1, 1)
 COVAL(OUTLET ,P1 ,1.0E+05 ,0. )
 COVAL(OUTLET ,TEM1,0. , SAME )
 
 PATCH(WFUNNORT,NWALL , 1, 1, 110, 110, 46, 200, 1, 1)
 COVAL(WFUNNORT,W1 , GRND2 ,0. )
 COVAL(WFUNNORT,KE ,1. ,0. )
 COVAL(WFUNNORT,OMEG, GRND2 , GRND2 )
 
 PATCH(HEATIN ,NORTH , 1, 1, 110, 110, 46, 200, 1, 1)
 COVAL(HEATIN ,TEM1, FIXFLU ,32.520325 )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 5000 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-05
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=50 ;LITER(V1)=10
 LITER(W1)=10 ;LITER(KE)=20
 LITER(LTLS)=20 ;LITER(OMEG)=20
 LITER(TEM1)=50
 ENDIT(P1)=1.0E-03 ;ENDIT(V1)=1.0E-03
 ENDIT(W1)=1.0E-03 ;ENDIT(KE)=1.0E-03
 ENDIT(LTLS)=1.0E-03 ;ENDIT(OMEG)=1.0E-03
 ENDIT(TEM1)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(V1,FALSDT,3.375E-03)
 RELAX(W1,FALSDT,3.375E-03)
 RELAX(KE,FALSDT,8.4375E-04)
 RELAX(EP,LINRLX,1.)
 RELAX(OMRS,LINRLX,1.)
 RELAX(OMCR,LINRLX,1.)
 RELAX(PRPS,LINRLX,1.)
 RELAX(NUSC,LINRLX,1.)
 RELAX(NUSS,LINRLX,1.)
 RELAX(HTCB,LINRLX,1.)
 RELAX(TWAL,LINRLX,1.)
 RELAX(TAVE,LINRLX,1.)
 RELAX(CNH1,LINRLX,1.)
 RELAX(LEN1,LINRLX,1.)
 RELAX(SKIN,LINRLX,1.)
 RELAX(SPH1,LINRLX,1.)
 RELAX(ENUL,LINRLX,1.)
 RELAX(KOND,LINRLX,1.)
 RELAX(STRS,LINRLX,1.)
 RELAX(YPLS,LINRLX,1.)
 RELAX(LTLS,FALSDT,1.0E+09)
 RELAX(WDIS,LINRLX,1.)
 RELAX(GEN1,LINRLX,1.)
 RELAX(BF2,LINRLX,1.)
 RELAX(BF1,LINRLX,0.05)
 RELAX(OMEG,FALSDT,8.4375E-04)
 RELAX(ENUT,LINRLX,1.)
 RELAX(TEM1,FALSDT,1.)
 RELAX(TDIF,LINRLX,1.)
 RELAX(CP,LINRLX,1.)
 OVRRLX =0.
 EXPERT = F ;NNORSL = F
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-1.0E+10
 VARMAX(V1)=50. ;VARMIN(V1)=-50.
 VARMAX(W1)=50. ;VARMIN(W1)=-50.
 VARMAX(KE)=1.0E+10 ;VARMIN(KE)=1.0E-10
 VARMAX(EP)=1.0E+10 ;VARMIN(EP)=1.0E-10
 VARMAX(OMRS)=1.0E+10 ;VARMIN(OMRS)=-1.0E+10
 VARMAX(OMCR)=1.0E+10 ;VARMIN(OMCR)=-1.0E+10
 VARMAX(PRPS)=1.0E+10 ;VARMIN(PRPS)=-1.0E+10
 VARMAX(NUSC)=1.0E+10 ;VARMIN(NUSC)=-1.0E+10
 VARMAX(NUSS)=1.0E+10 ;VARMIN(NUSS)=-1.0E+10
 VARMAX(HTCB)=1.0E+10 ;VARMIN(HTCB)=-1.0E+10
 VARMAX(TWAL)=1.0E+10 ;VARMIN(TWAL)=-1.0E+10
 VARMAX(TAVE)=1.0E+10 ;VARMIN(TAVE)=-1.0E+10
 VARMAX(CNH1)=1.0E+10 ;VARMIN(CNH1)=-1.0E+10
 VARMAX(LEN1)=1.0E+10 ;VARMIN(LEN1)=-1.0E+10
 VARMAX(SKIN)=1.0E+10 ;VARMIN(SKIN)=-1.0E+10
 VARMAX(SPH1)=1.0E+10 ;VARMIN(SPH1)=-1.0E+10
 VARMAX(ENUL)=1.0E+10 ;VARMIN(ENUL)=-1.0E+10
 VARMAX(KOND)=1.0E+10 ;VARMIN(KOND)=-1.0E+10
 VARMAX(STRS)=1.0E+10 ;VARMIN(STRS)=-1.0E+10
 VARMAX(YPLS)=1.0E+10 ;VARMIN(YPLS)=-1.0E+10
 VARMAX(LTLS)=1.0E+10 ;VARMIN(LTLS)=-1.0E+10
 VARMAX(WDIS)=1.0E+10 ;VARMIN(WDIS)=-1.0E+10
 VARMAX(GEN1)=1.0E+10 ;VARMIN(GEN1)=-1.0E+10
 VARMAX(BF2)=1.0E+10 ;VARMIN(BF2)=-1.0E+10
 VARMAX(BF1)=1.0E+10 ;VARMIN(BF1)=-1.0E+10
 VARMAX(OMEG)=1.0E+10 ;VARMIN(OMEG)=1.0E-10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=1.0E-10
 VARMAX(TEM1)=1.0E+09 ;VARMIN(TEM1)=19.
 VARMAX(TDIF)=1.0E+10 ;VARMIN(TDIF)=-1.0E+10
 VARMAX(CP)=1.0E+10 ;VARMIN(CP)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 IENUTA = 20
 ISG62 = 1
 SPEDAT(SET,STORED,CP,C,=(P1-P1[1&1&1])/50.!IMAT<100!ZSLFIN)
 SPEDAT(SET,STORED,TDIF,C,=(TEM1-20.)!IMAT<100!ZSLFIN)
 SPEDAT(SET,STORED,CNH1,C,=CNH1[&&-1]!IF(IZ.EQ.200)!ZSLFIN)
 SPEDAT(SET,MAKE,TSUM,C,=0.)
 SPEDAT(SET,MAKE,ASUM,C,=0.)
 SPEDAT(SET,STORED,ASUM!PATCH1,C,=SSUM(CNH1*CP1)!IMAT<100!ZSLFIN)
 SPEDAT(SET,STORED,TSUM!PATCH1,C,=SSUM(CNH1*CP1*TEM1)!IMAT<100!ZSL$)
 SPEDAT(SET,STORED,TSUM!PATCH1,C,FIN)
 SPEDAT(SET,STORED,TAVE!PATCH1,C,=TSUM/ASUM!ZSLFIN)
 SPEDAT(SET,STORED,TWAL!HEATIN,C,=TEM1+32.520325*0.5*DYV/KOND!ZSLFIN)
 SPEDAT(SET,STORED,HTCB!HEATIN,C,=32.520325/(TWAL-TAVE+TINY)!ZSLFIN)
 SPEDAT(SET,STORED,NUSS!HEATIN,C,=(HTCB*1./KOND)!ZSLFIN)
 SPEDAT(SET,STORED,NUSC!HEATIN,C,=NUSS/97.24543!ZSLFIN)
 SPEDAT(SET,STORED,OMCR,C,=CORR(OMEG)!IMAT<100!CORREC)
 SPEDAT(SET,STORED,OMRS,C,=RESI(OMEG)!IMAT<100!RESIDU)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,MATERIAL,198,L,T)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=6.9 ;EX(V1)=0.05786
 EX(W1)=3.583 ;EX(KE)=1.087
 EX(EP)=20.15 ;EX(OMRS)=1.067E-06
 EX(OMCR)=1.415E-04 ;EX(PRPS)=0.8773
 EX(NUSC)=0.01408 ;EX(NUSS)=1.369
 EX(HTCB)=0.07681 ;EX(TWAL)=0.1778
 EX(TAVE)=17.75 ;EX(CNH1)=2.053E-04
 EX(LEN1)=0.03346 ;EX(SKIN)=1.659
 EX(SPH1)=877.299988 ;EX(ENUL)=3.445E-05
 EX(KOND)=0.04921 ;EX(STRS)=6.137E-04
 EX(YPLS)=0.01447 ;EX(LTLS)=0.02731
 EX(WDIS)=0.1313 ;EX(GEN1)=5.891E+07
 EX(BF2)=0.7881 ;EX(BF1)=0.6633
 EX(OMEG)=568.599976 ;EX(ENUT)=0.016
 EX(TEM1)=18.18 ;EX(TDIF)=0.6343
 EX(CP)=0.09523
 ************************************************************
  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(OMRS,Y,N,Y,N,N,N)
 OUTPUT(OMCR,Y,N,Y,N,N,N)
 OUTPUT(PRPS,Y,N,Y,N,N,N)
 OUTPUT(NUSC,Y,N,Y,N,N,N)
 OUTPUT(NUSS,Y,N,Y,N,N,N)
 OUTPUT(HTCB,Y,N,Y,N,N,N)
 OUTPUT(TWAL,Y,N,Y,N,N,N)
 OUTPUT(TAVE,Y,N,Y,N,N,N)
 OUTPUT(CNH1,Y,N,Y,N,N,N)
 OUTPUT(LEN1,Y,N,Y,N,N,N)
 OUTPUT(SKIN,Y,N,Y,N,N,N)
 OUTPUT(SPH1,Y,N,Y,N,N,N)
 OUTPUT(ENUL,Y,N,Y,N,N,N)
 OUTPUT(KOND,Y,N,Y,N,N,N)
 OUTPUT(STRS,Y,N,Y,N,N,N)
 OUTPUT(YPLS,Y,N,Y,N,N,N)
 OUTPUT(LTLS,Y,N,Y,Y,Y,Y)
 OUTPUT(WDIS,Y,N,Y,N,N,N)
 OUTPUT(GEN1,Y,N,Y,Y,Y,Y)
 OUTPUT(BF2,Y,N,Y,Y,Y,Y)
 OUTPUT(BF1,Y,N,Y,Y,Y,Y)
 OUTPUT(OMEG,Y,N,Y,Y,Y,Y)
 OUTPUT(ENUT,Y,N,Y,N,Y,Y)
 OUTPUT(TEM1,Y,N,Y,Y,Y,Y)
 OUTPUT(TDIF,Y,N,Y,N,N,N)
 OUTPUT(CP,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 1 ;IYMON = 52 ;IZMON = 47
 NPRMON = 10000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 IPLTF = 1 ;IPLTL = 2000 ;NPLT = 50
 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(NUWALL ,PROFIL, 1, 1, 110, 110, 46, 200, 1, 1)
 PLOT(NUWALL ,NUSC,0. ,0. )
 PLOT(NUWALL ,NUSS,0. ,0. )
 PLOT(NUWALL ,YPLS,0. ,0. )
 ************************************************************
  Group 24. Dumps For Restarts
 SAVE = T ;NOWIPE = F
 NSAVE =CHAM
STOP