************************************************************
   Q1 created by VDI menu, Version 2017, Date 20/01/17
 CPVNAM=VDI; SPPNAM=Core
 ************************************************************
  Echo DISPLAY / USE settings
  DISPLAY
  The case considered is the simulation of the constant-stress-
  layer region of a 2d turbulent neutral atmospheric boundary
  layer using the Realisable k-e model and fully-rough logarithmic
  wall functions. In this region the velocity and turbulence
  variables are given by:
 
      u/u*=ln(z/z0)/K       ; k=u**2/sqrt(cmucd)
 
      e=u**3/(K*z)          ; enut=K*z*u*
 
  where u*=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 u*=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 a uniform pressure field.
  ENDDIS
 ************************************************************
 IRUNN = 1 ;LIBREF = 314
 ************************************************************
  Group 1. Run Title
 TEXT(NEUTRAL ABL WIND FLOW-REAL.K-E MOD:T314 )
 ************************************************************
  Group 2. Transience
 STEADY = T
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,30,5,30)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
    * Non-default variable names
 NAME(135)=DWDZ ;NAME(136)=DWDY
 NAME(137)=DWDX ;NAME(138)=DVDZ
 NAME(139)=DVDY ;NAME(140)=DVDX
 NAME(141)=DUDZ ;NAME(142)=DUDY
 NAME(143)=DUDX ;NAME(144)=EPKE
 NAME(145)=VABS ;NAME(146)=STRS
 NAME(147)=CMU ;NAME(148)=DEN1
 NAME(149)=EL1 ;NAME(150)=ENUT
    * Solved variables list
 SOLVE(P1,U1,V1,W1)
    * Stored variables list
 STORE(ENUT,EL1,DEN1,CMU,STRS,VABS,EPKE,DUDX)
 STORE(DUDY,DUDZ,DVDX,DVDY,DVDZ,DWDX,DWDY,DWDZ)
    * Additional solver options
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 TURMOD(KEREAL)
 
 ************************************************************
  Group 8. Terms & Devices
 ************************************************************
  Group 9. Properties
 PRESS0 =1.01325E+05 ;TEMP0 =273.
    * Domain material index is   0 signifying:
    * Air at 20 deg C, 1 atm, treated as incompressible
 SETPRPS(1, 0)
 DVO1DT =3.41E-03
 PRT(EP)=1.2
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(P1)=0. ;FIINIT(CMU)=0.09
   No PATCHes used for this Group
 
 
 INIADD = F
 ************************************************************
  Group 12. Convection and diffusion adjustments
   No PATCHes used for this Group
 ************************************************************
  Group 13. Boundary & Special Sources
   No PATCHes used for this Group
 
 BUOYD =1.189
 BUOYE =20.
 EGWF = T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 200
 RESFAC =9.999999E-08
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=200
 ************************************************************
  Group 17. Relaxation
 RELAX(P1 ,LINRLX,1. )
 RELAX(KE ,LINRLX,0.5 )
 RELAX(EP ,LINRLX,0.5 )
 KELIN = 3
 ************************************************************
  Group 18. Limits
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 GENK = T
 PARSOL = F
 IENUTA = 14
 ISG62 = 1
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,KECONST,C2E,R,1.9)
 SPEDAT(SET,MAXINC,KE,R,0.025)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1  )=3.056E-04;  EX(U1  )=5.768E+00
 EX(V1  )=2.603E-04;  EX(W1  )=1.418E-04
 EX(KE  )=4.160E-01;  EX(EP  )=1.158E-02
 EX(DWDZ)=2.363E-05;  EX(DWDY)=1.129E-05
 EX(DWDX)=3.513E-06;  EX(DVDZ)=2.045E-05
 EX(DVDY)=2.650E-05;  EX(DVDX)=4.369E-06
 EX(DUDZ)=8.167E-02;  EX(DUDY)=3.206E-05
 EX(DUDX)=6.347E-05;  EX(EPKE)=2.803E-02
 EX(VABS)=5.768E+00;  EX(STRS)=3.826E-03
 EX(CMU )=9.044E-02;  EX(DEN1)=1.189E+00
 EX(EL1 )=1.561E+01;  EX(ENUT)=5.514E+00
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(DWDZ,Y,N,N,N,N,N)
 OUTPUT(DWDY,Y,N,N,N,N,N)
 OUTPUT(DWDX,Y,N,N,N,N,N)
 OUTPUT(DVDZ,Y,N,N,N,N,N)
 OUTPUT(DVDY,Y,N,N,N,N,N)
 OUTPUT(DVDX,Y,N,N,N,N,N)
 OUTPUT(DUDZ,Y,N,N,N,N,N)
 OUTPUT(DUDY,Y,N,N,N,N,N)
 OUTPUT(DUDX,Y,N,N,N,N,N)
 OUTPUT(EPKE,Y,N,N,N,N,N)
 OUTPUT(CMU ,Y,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 23 ;IYMON = 1 ;IZMON = 5
 NPRMON = 100000
 NPRMNT = 1
 TSTSWP = -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000
 ISWPRF = 1 ;ISWPRL = 100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 
 GVIEW(P,-4.999917E-03,-0.999975,4.999979E-03)
 GVIEW(UP,2.499979E-05,4.999917E-03,0.999987)
 GVIEW(DEPTH,3.0E+04)
 GVIEW(VDIS,109.981827)
 GVIEW(CENTRE,100.,25.,50.)
 
> DOM,    SIZE,        2.000000E+02, 5.000000E+01, 1.000000E+02
> DOM,    MONIT,       1.500000E+02, 5.000000E+00, 7.402118E+00
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM,    INCREMENT,   1.000000E-02, 1.000000E-02, 1.000000E-02
> GRID,   RSET_X_1,     30, 1.000000E+00,G
> GRID,   RSET_Y_1,      5, 1.000000E+00
> GRID,   RSET_Z_1,     30, 1.050000E+00,G
> DOM,    INI_AMB,    YES
> DOM,    INI_BUOY,   YES
 
> OBJ,    NAME,        WIND1
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        TO_END,       TO_END,       TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    windgr
> OBJ,    ROTATION24,  7
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        WIND
> OBJ,    WIREFRAME,   YES
> OBJ,    PRESSURE,    1.01325E+05
> OBJ,    COEFFICIENT, 1.0E+04
> OBJ,    VELOCITY,    5.
> OBJ,    WIND_DIR,    270.
> OBJ,    AXIS_DIR,    0.
> OBJ,    PROFILE,      Logarithmic
> OBJ,    REF_HEIGHT,  10.
> OBJ,    RGHNS_HEIGHT,0.03
> OBJ,    UP-DIR,       Z
> OBJ,    SKY,          YES
> OBJ,    GROUND,       YES
> OBJ,    GROUND-TEMP,  ADIABATIC