TALK=T;RUN(1,1)
  DISPLAY
  The case considered is the simulation of the constant-stress-
  layer region of a 2d turbulent neutral atmospheric boundary
  layer using the Wilcox 1988 k-w model and fully-rough logarithmic
  wall functions. In this region the velocity and turbulence
  variables are given by:

      u/utau=ln(z/z0)/K       ; k=utau^2/sqrt(cmucd)

      w=utau/(sqrt(cmucd)*K*z); enut=K*z*utau

  where utau=friction velocity, K=0.41, cmucd=0.09 and z0 is the
  roughness height. A calculation is made using the following
  wind inlet-profile data set by the wind object: Wind reference
  speed: Ur=5m/s at a reference height Zr=10m, wind direction
  from the west, and terrain roughness height z0=0.03 m. The
  friction velocity is given by utau=Ur*K/ln(zr/z0)=0.353m/s
  and so the shear_stress/density=0.125 m^2/s^2. The wind object
  is used which sets the wind inlet profiles, a fixed-pressure
  condition at the downstream boundary, and a constant shear-
  stress condition at the sky boundary via an entrainment and
  diffusive boundary condition. The height(z) of the solution
  domain is 100m, and the domain length is 200m(x). The wind enters
  from the west boundary, and leaves via the east boundary. The
  flow is 2d, although the simulation is 3d. The calculation
  demonstrates that PHOENICS preserves the wind inlet profiles
  and returns what is an essentially uniform pressure field.
  ENDDIS
 ************************************************************
  Group 1. Run Title and Number
 ************************************************************
 ************************************************************
 
 TEXT(T703:2DZY ABL WIND FLOW- Wilcox 1988 k-w)
 
 ************************************************************
 ************************************************************
 
 IRUNN = 1 ;LIBREF = 0
 ************************************************************
  Group 2. Time dependence
 STEADY = T
 ************************************************************
  Group 3. X-Direction Grid Spacing
 CARTES = T
 NX = 30
 XULAST =200.
 XFRAC(1)=0.033333 ;XFRAC(2)=0.066667
 XFRAC(3)=0.1 ;XFRAC(4)=0.133333
 XFRAC(5)=0.166667 ;XFRAC(6)=0.2
 XFRAC(7)=0.233333 ;XFRAC(8)=0.266667
 XFRAC(9)=0.3 ;XFRAC(10)=0.333333
 XFRAC(11)=0.366667 ;XFRAC(12)=0.4
 XFRAC(13)=0.433333 ;XFRAC(14)=0.466667
 XFRAC(15)=0.5 ;XFRAC(16)=0.533333
 XFRAC(17)=0.566667 ;XFRAC(18)=0.6
 XFRAC(19)=0.633333 ;XFRAC(20)=0.666667
 XFRAC(21)=0.7 ;XFRAC(22)=0.733333
 XFRAC(23)=0.766667 ;XFRAC(24)=0.8
 XFRAC(25)=0.833334 ;XFRAC(26)=0.866667
 XFRAC(27)=0.9 ;XFRAC(28)=0.933334
 XFRAC(29)=0.966667 ;XFRAC(30)=1.
 ************************************************************
  Group 4. Y-Direction Grid Spacing
 NY = 5
 YVLAST =50.
 YFRAC(1)=0.2 ;YFRAC(2)=0.4
 YFRAC(3)=0.6 ;YFRAC(4)=0.8
 YFRAC(5)=1.
 ************************************************************
  Group 5. Z-Direction Grid Spacing
 PARAB = F
 NZ = 30
 ZWLAST =100.
 ZFRAC(1)=0.015051 ;ZFRAC(2)=0.030855
 ZFRAC(3)=0.04745 ;ZFRAC(4)=0.064874
 ZFRAC(5)=0.083169 ;ZFRAC(6)=0.102379
 ZFRAC(7)=0.122549 ;ZFRAC(8)=0.143728
 ZFRAC(9)=0.165966 ;ZFRAC(10)=0.189315
 ZFRAC(11)=0.213833 ;ZFRAC(12)=0.239576
 ZFRAC(13)=0.266606 ;ZFRAC(14)=0.294988
 ZFRAC(15)=0.324788 ;ZFRAC(16)=0.356079
 ZFRAC(17)=0.388935 ;ZFRAC(18)=0.423433
 ZFRAC(19)=0.459656 ;ZFRAC(20)=0.49769
 ZFRAC(21)=0.537626 ;ZFRAC(22)=0.579559
 ZFRAC(23)=0.623588 ;ZFRAC(24)=0.669819
 ZFRAC(25)=0.718361 ;ZFRAC(26)=0.769331
 ZFRAC(27)=0.822849 ;ZFRAC(28)=0.879043
 ZFRAC(29)=0.938046 ;ZFRAC(30)=1.
 ************************************************************
  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(137)=EVDS ;NAME(138)=KDIS
 NAME(139)=UDIS ;NAME(140)=EVIN
 NAME(141)=OMIN ;NAME(142)=EPIN
 NAME(143)=KEIN ;NAME(144)=UIN
 NAME(145)=VABS ;NAME(146)=STRS
 NAME(147)=OMEG ;NAME(148)=DEN1
 NAME(149)=EL1 ;NAME(150)=ENUT
    * 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,Y)
 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,N,N,N,N,N)
 SOLUTN(EVDS,Y,N,N,N,N,Y)
 SOLUTN(KDIS,Y,N,N,N,N,Y)
 SOLUTN(UDIS,Y,N,N,N,N,Y)
 SOLUTN(EVIN,Y,N,N,N,N,Y)
 SOLUTN(OMIN,Y,N,N,N,N,Y)
 SOLUTN(EPIN,Y,N,N,N,N,Y)
 SOLUTN(KEIN,Y,N,N,N,N,Y)
 SOLUTN(UIN,Y,N,N,N,N,Y)
 SOLUTN(VABS,Y,N,N,N,N,Y)
 SOLUTN(STRS,Y,N,N,N,N,Y)
 SOLUTN(OMEG,Y,Y,N,N,N,N)
 SOLUTN(DEN1,Y,N,N,N,N,Y)
 SOLUTN(EL1,Y,N,N,N,N,Y)
 SOLUTN(ENUT,Y,N,N,N,N,Y)
 DEN1 = 148
 VIST = 150
 LEN1 = 149
 ************************************************************
  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(OMEG,N,Y,Y,Y,Y,N)
 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.189 ;TMP1 =0. ;EL1 = GRND5
 TSURR =0. ;TEMP0 =273. ;PRESS0 =1.01325E+05
 DVO1DT =3.41E-03 ;DRH1DP =0.
 EMISS =0. ;SCATT =0.
 RADIA =0. ;RADIB =0.
 EL1A =0.6 ;EL1B =0. ;EL1C =0.
 ENUL =1.544E-05 ;ENUT = GRND7
 ENUTA =0. ;ENUTB =0. ;ENUTC =0.
 IENUTA = 10
 PRNDTL(U1)=1. ;PRNDTL(V1)=1.
 PRNDTL(W1)=1. ;PRNDTL(KE)=1.
 PRNDTL(OMEG)=1.
 PRT(U1)=1. ;PRT(V1)=1.
 PRT(W1)=1. ;PRT(KE)=2.
 PRT(OMEG)=2.
 CP1 =1005. ;CP2 =1.
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initial field variables (PHIs)
 FIINIT(P1)=0. ;FIINIT(U1)=5.
 FIINIT(V1)=-5.96244E-08 ;FIINIT(W1)=1.0E-10
 FIINIT(KE)=0.415109 ;FIINIT(EP)=2.143737E-03
 FIINIT(EVDS)=1.0E-10 ;FIINIT(KDIS)=1.0E-10
 FIINIT(UDIS)=1.0E-10 ;FIINIT(EVIN)=1.0E-10
 FIINIT(OMIN)=1.0E-10 ;FIINIT(EPIN)=1.0E-10
 FIINIT(KEIN)=1.0E-10 ;FIINIT(UIN)=1.0E-10
 FIINIT(VABS)=1.0E-10 ;FIINIT(STRS)=1.0E-10
 FIINIT(OMEG)=0.05738 ;FIINIT(DEN1)=1.0E-10
 FIINIT(EL1)=1.0E-10 ;FIINIT(ENUT)=1.0E-10
   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(KWSOURCE,PHASEM, 1, 30, 1, 5, 1, 30, 1, 1)
 COVAL(KWSOURCE,KE , GRND4 , GRND4 )
 COVAL(KWSOURCE,OMEG, GRND4 , GRND4 )
 
 PATCH(BLIN1 ,NORTH , 1, 30, 5, 5, 1, 30, 1, 1)
 COVAL(BLIN1 ,P1 , GRND8 , GRND7 )
 COVAL(BLIN1 ,U1 ,0. , GRND7 )
 COVAL(BLIN1 ,V1 ,0. , GRND7 )
 COVAL(BLIN1 ,KE ,0. , GRND7 )
 COVAL(BLIN1 ,OMEG,0. , GRND7 )
 
 PATCH(BLIN2 ,SOUTH , 1, 30, 1, 1, 1, 30, 1, 1)
 COVAL(BLIN2 ,P1 , GRND8 , GRND7 )
 COVAL(BLIN2 ,U1 ,0. , GRND7 )
 COVAL(BLIN2 ,V1 ,0. , GRND7 )
 COVAL(BLIN2 ,KE ,0. , GRND7 )
 COVAL(BLIN2 ,OMEG,0. , GRND7 )
 
 PATCH(BLIN3 ,EAST , 30, 30, 1, 5, 1, 30, 1, 1)
 COVAL(BLIN3 ,P1 , GRND8 , GRND7 )
 COVAL(BLIN3 ,U1 ,0. , GRND7 )
 COVAL(BLIN3 ,V1 ,0. , GRND7 )
 COVAL(BLIN3 ,KE ,0. , GRND7 )
 COVAL(BLIN3 ,OMEG,0. , GRND7 )
 
 PATCH(BLIN4 ,WEST , 1, 1, 1, 5, 1, 30, 1, 1)
 COVAL(BLIN4 ,P1 , GRND8 , GRND7 )
 COVAL(BLIN4 ,U1 ,0. , GRND7 )
 COVAL(BLIN4 ,V1 ,0. , GRND7 )
 COVAL(BLIN4 ,KE ,0. , GRND7 )
 COVAL(BLIN4 ,OMEG,0. , GRND7 )
 
 PATCH(BLIN5 ,HIGH , 1, 30, 1, 5, 30, 30, 1, 1)
 COVAL(BLIN5 ,P1 ,1.0E+04 ,0. )
 COVAL(BLIN5 ,U1 , GRND7 , GRND7 )
 COVAL(BLIN5 ,V1 , GRND7 , GRND7 )
 COVAL(BLIN5 ,KE , GRND7 , GRND7 )
 COVAL(BLIN5 ,OMEG, GRND7 , GRND7 )
 
 PATCH(BLIN6 ,LWALL , 1, 30, 1, 5, 1, 1, 1, 1)
 COVAL(BLIN6 ,U1 , GRND5 ,0. )
 COVAL(BLIN6 ,V1 , GRND5 ,0. )
 COVAL(BLIN6 ,KE , GRND5 , GRND5 )
 COVAL(BLIN6 ,OMEG, GRND5 , GRND5 )
 XCYCLE = F
 EGWF = T
 WALLCO = GRND2
 BUOYD =1.189
 BUOYE =20.
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 200 ;ISWC1 = 1
 LITHYD = 1 ;LITFLX = 1 ;LITC = 1 ;ITHC1 = 1
 SELREF = T
 RESFAC =1.0E-05
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=200 ;LITER(U1)=10
 LITER(V1)=10 ;LITER(W1)=10
 LITER(KE)=20 ;LITER(OMEG)=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(OMEG)=1.0E-03
 ************************************************************
  Group 17. Relaxation
 RELAX(P1,LINRLX,1.)
 RELAX(U1,FALSDT,1.)
 RELAX(V1,FALSDT,1.)
 RELAX(W1,FALSDT,1.)
 RELAX(KE,FALSDT,1.)
 RELAX(EP,LINRLX,1.)
 RELAX(EVDS,LINRLX,1.)
 RELAX(KDIS,LINRLX,1.)
 RELAX(UDIS,LINRLX,1.)
 RELAX(EVIN,LINRLX,1.)
 RELAX(OMIN,LINRLX,1.)
 RELAX(EPIN,LINRLX,1.)
 RELAX(KEIN,LINRLX,1.)
 RELAX(UIN,LINRLX,1.)
 RELAX(VABS,LINRLX,1.)
 RELAX(STRS,LINRLX,1.)
 RELAX(OMEG,FALSDT,1.)
 RELAX(DEN1,LINRLX,1.)
 RELAX(EL1,LINRLX,1.)
 RELAX(ENUT,LINRLX,1.)
 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(EVDS)=1.0E+10 ;VARMIN(EVDS)=-1.0E+10
 VARMAX(KDIS)=1.0E+10 ;VARMIN(KDIS)=-1.0E+10
 VARMAX(UDIS)=1.0E+10 ;VARMIN(UDIS)=-1.0E+10
 VARMAX(EVIN)=1.0E+10 ;VARMIN(EVIN)=-1.0E+10
 VARMAX(OMIN)=1.0E+10 ;VARMIN(OMIN)=-1.0E+10
 VARMAX(EPIN)=1.0E+10 ;VARMIN(EPIN)=-1.0E+10
 VARMAX(KEIN)=1.0E+10 ;VARMIN(KEIN)=-1.0E+10
 VARMAX(UIN)=1.0E+10 ;VARMIN(UIN)=-1.0E+10
 VARMAX(VABS)=1.0E+10 ;VARMIN(VABS)=-1.0E+10
 VARMAX(STRS)=1.0E+10 ;VARMIN(STRS)=-1.0E+10
 VARMAX(OMEG)=1.0E+10 ;VARMIN(OMEG)=1.0E-10
 VARMAX(DEN1)=1.0E+10 ;VARMIN(DEN1)=-1.0E+10
 VARMAX(EL1)=1.0E+10 ;VARMIN(EL1)=-1.0E+10
 VARMAX(ENUT)=1.0E+10 ;VARMIN(ENUT)=-1.0E+10
 ************************************************************
  Group 19. Data transmitted to GROUND
 GENK = T
 PARSOL = F
 IENUTA = 10
 ISG62 = 1
 SPEDAT(SET,DOMAIN,PHASE_1_MAT,I,0)
 SPEDAT(SET,MAXINC,KE,R,0.025)
 SPEDAT(SET,STORED,UIN,C,=(0.352892*LOGE(ZG/0.03)/0.41)!ZSLSTR)
 SPEDAT(SET,STORED,KEIN,C,=0.415109!ZSLSTR)
 SPEDAT(SET,STORED,EPIN,C,=(0.107187/ZG)!ZSLSTR)
 SPEDAT(SET,STORED,EVIN,C,=0.41*0.352892*ZG!ZSLSTR)
 SPEDAT(SET,STORED,OMIN,C,=0.107187/(ZG*0.09*0.415109)!ZSLSTR)
 SPEDAT(SET,STORED,UDIS,C,=100.*(U1/UIN-1.0)!IF(IX.LT.30)!ZSLFIN)
 SPEDAT(SET,STORED,KDIS,C,=100.*(KE/KEIN-1.0)!ZSLFIN)
 SPEDAT(SET,STORED,EVDS,C,=100.*(ENUT/EVIN-1.0)!ZSLFIN)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,BLIN,INIBUOY,I,1)
 SPEDAT(SET,INITIAL,U1,C,=0.860712*LOGE(MAX(ZG/0.03&2.0))!IMAT<100)
 SPEDAT(SET,INITIAL,V1,C,=-1.026389E-08*LOGE(MAX(ZG/0.03&2.0))!IMA$)
 SPEDAT(SET,INITIAL,V1,C,T<100)
 SPEDAT(SET,INITIAL,KE,C,=0.415109!IMAT<100)
 SPEDAT(SET,INITIAL,EP,C,=0.043947/(0.41*ZG)!IMAT<100)
 SPEDAT(SET,INITIAL,OMEG,C,=0.124533/(0.123002*ZG)!IMAT<100)
 SPEDAT(SET,BLIN,WDIR,R,270.)
 SPEDAT(SET,BLIN,AXDIR,R,0.)
 SPEDAT(SET,OBJNAM,!BLIN1,C,WIND1)
 SPEDAT(SET,OBJTYP,!BLIN1,C,USER_DEFINED)
 SPEDAT(SET,BLIN1,VELX,R,5.)
 SPEDAT(SET,BLIN1,VELY,R,-5.96244E-08)
 SPEDAT(SET,BLIN1,VELZ,R,0.)
 SPEDAT(SET,BLIN1,PCOEF,R,1.0E+04)
 SPEDAT(SET,BLIN1,PEXT,R,1.01325E+05)
 SPEDAT(SET,BLIN1,VDIR,C,Z)
 SPEDAT(SET,BLIN1,RHOIN,R,1.189)
 SPEDAT(SET,BLIN1,BLTY,C,LOGL)
 SPEDAT(SET,BLIN1,ZO,R,0.03)
 SPEDAT(SET,BLIN1,REFH,R,10.)
 SPEDAT(SET,BLIN1,HO,R,0.)
 SPEDAT(SET,OBJNAM,!BLIN2,C,WIND1)
 SPEDAT(SET,OBJTYP,!BLIN2,C,USER_DEFINED)
 SPEDAT(SET,BLIN2,VELX,R,5.)
 SPEDAT(SET,BLIN2,VELY,R,-5.96244E-08)
 SPEDAT(SET,BLIN2,VELZ,R,0.)
 SPEDAT(SET,BLIN2,PCOEF,R,1.0E+04)
 SPEDAT(SET,BLIN2,PEXT,R,1.01325E+05)
 SPEDAT(SET,BLIN2,VDIR,C,Z)
 SPEDAT(SET,BLIN2,RHOIN,R,1.189)
 SPEDAT(SET,BLIN2,BLTY,C,LOGL)
 SPEDAT(SET,BLIN2,ZO,R,0.03)
 SPEDAT(SET,BLIN2,REFH,R,10.)
 SPEDAT(SET,BLIN2,HO,R,0.)
 SPEDAT(SET,OBJNAM,!BLIN3,C,WIND1)
 SPEDAT(SET,OBJTYP,!BLIN3,C,USER_DEFINED)
 SPEDAT(SET,BLIN3,VELX,R,5.)
 SPEDAT(SET,BLIN3,VELY,R,-5.96244E-08)
 SPEDAT(SET,BLIN3,VELZ,R,0.)
 SPEDAT(SET,BLIN3,PCOEF,R,1.0E+04)
 SPEDAT(SET,BLIN3,PEXT,R,1.01325E+05)
 SPEDAT(SET,BLIN3,VDIR,C,Z)
 SPEDAT(SET,BLIN3,RHOIN,R,1.189)
 SPEDAT(SET,BLIN3,BLTY,C,LOGL)
 SPEDAT(SET,BLIN3,ZO,R,0.03)
 SPEDAT(SET,BLIN3,REFH,R,10.)
 SPEDAT(SET,BLIN3,HO,R,0.)
 SPEDAT(SET,OBJNAM,!BLIN4,C,WIND1)
 SPEDAT(SET,OBJTYP,!BLIN4,C,USER_DEFINED)
 SPEDAT(SET,BLIN4,VELX,R,5.)
 SPEDAT(SET,BLIN4,VELY,R,-5.96244E-08)
 SPEDAT(SET,BLIN4,VELZ,R,0.)
 SPEDAT(SET,BLIN4,PCOEF,R,1.0E+04)
 SPEDAT(SET,BLIN4,PEXT,R,1.01325E+05)
 SPEDAT(SET,BLIN4,VDIR,C,Z)
 SPEDAT(SET,BLIN4,RHOIN,R,1.189)
 SPEDAT(SET,BLIN4,BLTY,C,LOGL)
 SPEDAT(SET,BLIN4,ZO,R,0.03)
 SPEDAT(SET,BLIN4,REFH,R,10.)
 SPEDAT(SET,BLIN4,HO,R,0.)
 SPEDAT(SET,OBJNAM,!BLIN5,C,WIND1)
 SPEDAT(SET,OBJTYP,!BLIN5,C,USER_DEFINED)
 SPEDAT(SET,BLIN5,VELX,R,5.)
 SPEDAT(SET,BLIN5,VELY,R,-5.96244E-08)
 SPEDAT(SET,BLIN5,VELZ,R,0.)
 SPEDAT(SET,BLIN5,PCOEF,R,1.0E+04)
 SPEDAT(SET,BLIN5,PEXT,R,1.01325E+05)
 SPEDAT(SET,BLIN5,VDIR,C,Z)
 SPEDAT(SET,BLIN5,RHOIN,R,1.189)
 SPEDAT(SET,BLIN5,BLTY,C,LOGL)
 SPEDAT(SET,BLIN5,ZO,R,0.03)
 SPEDAT(SET,BLIN5,REFH,R,10.)
 SPEDAT(SET,BLIN5,HO,R,0.)
 SPEDAT(SET,OBJNAM,!BLIN6,C,WIND1)
 SPEDAT(SET,OBJTYP,!BLIN6,C,USER_DEFINED)
 SPEDAT(SET,ROUGHNESS,BLIN6,R,0.03)
 SPEDAT(SET,FACETDAT,NUMOBJ,I,1)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=2.405E-04 ;EX(U1)=5.768
 EX(V1)=2.113E-04 ;EX(W1)=1.426E-04
 EX(KE)=0.4166 ;EX(EP)=0.0118
 EX(EVDS)=1.031 ;EX(KDIS)=1.247
 EX(UDIS)=0.2055 ;EX(EVIN)=5.531
 EX(OMIN)=0.3055 ;EX(EPIN)=0.01141
 EX(KEIN)=0.4151 ;EX(UIN)=5.767
 EX(VABS)=5.768 ;EX(STRS)=3.855E-03
 EX(OMEG)=0.3141 ;EX(DEN1)=1.189
 EX(EL1)=15.64 ;EX(ENUT)=5.521
 ************************************************************
  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(EVDS,Y,N,Y,N,N,N)
 OUTPUT(KDIS,Y,N,Y,N,N,N)
 OUTPUT(UDIS,Y,N,Y,N,N,N)
 OUTPUT(EVIN,Y,N,Y,N,N,N)
 OUTPUT(OMIN,Y,N,Y,N,N,N)
 OUTPUT(EPIN,Y,N,Y,N,N,N)
 OUTPUT(KEIN,Y,N,Y,N,N,N)
 OUTPUT(UIN,Y,N,Y,N,N,N)
 OUTPUT(VABS,Y,N,Y,N,N,N)
 OUTPUT(STRS,Y,N,Y,N,N,N)
 OUTPUT(OMEG,Y,N,Y,Y,Y,Y)
 OUTPUT(DEN1,Y,N,Y,N,N,N)
 OUTPUT(EL1,Y,N,Y,N,N,N)
 OUTPUT(ENUT,Y,N,Y,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 23 ;IYMON = 1 ;IZMON = 5
 NPRMON = 100000 ;NPRMNT = 1 ;TSTSWP = -1
 UWATCH = T ;USTEER = T
 HIGHLO = F
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000 ;NUMCLS = 5
 NXPRIN = -1 ;IXPRF = 1 ;IXPRL = 10000
 NYPRIN = -1 ;IYPRF = 1 ;IYPRL = 10000
 NZPRIN = -1 ;IZPRF = 1 ;IZPRL = 10000
 XZPR = F ;YZPR = F
 IPLTF = 1 ;IPLTL = -1 ;NPLT = -1
 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