TALK=T;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2017, Date 01/09/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 k-w SST 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
   save1begin  
  The inlet velocity and turbulence log-law profiles aren't 
  preserved exactly because a spike in the k profile is produced 
  at the second cell from the wall, which originates from the 
  standard discretisation practices employed for the turbulence 
  production rate Pk in the turbulence transport equations and the 
  vertical diffusion terms in the axial momentum equation (see
  library case T306 for more details).  
   save1end 
 ************************************************************
 IRUNN = 1 ;LIBREF = 704
 ************************************************************
  Group 1. Run Title
 TEXT(T704:2DZY ABL WIND FLOW- k-w SST model  )
 ************************************************************
  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(129)=EVDS ;NAME(130)=KDIS
 NAME(131)=UDIS ;NAME(132)=EVIN
 NAME(133)=OMIN ;NAME(134)=EPIN
 NAME(135)=KEIN ;NAME(136)=UIN
 NAME(137)=SIGW ;NAME(138)=DUDZ
 NAME(139)=SIGK ;NAME(141)=WDIS
 NAME(145)=VABS ;NAME(146)=STRS
 NAME(148)=DEN1 ;NAME(149)=EL1
 NAME(150) =ENUT
    * Solved variables list
 SOLVE(P1,U1,V1,W1)
    * Stored variables list
 STORE(ENUT,EL1,DEN1,STRS,VABS,WDIS,SIGK,DUDZ)
 STORE(SIGW,UIN,KEIN,EPIN,OMIN,EVIN,UDIS,KDIS)
 STORE(EVDS)
    * Additional solver options
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 TURMOD(KWSST)
 
 ************************************************************
  Echo save-block settings for Group  7
  save7begin
REAL(VKC,XLEN,REFH,zdz0)
REAL(USTAR,USTAR2,Z0,UIN)
REAL(KCNS,USTAR3,EPCON)
UIN =5.0               ! reference velocity
REFH=10.               ! reference height
Z0=0.03                !  roughness height
VKC=0.41               ! von karman's constant
ZDZ0=REFH/Z0
USTAR=UIN*VKC/LOG(zdz0)  ! friction velocity
USTAR
USTAR2=USTAR*USTAR;KCNS=USTAR2/0.3;USTAR3=USTAR2*USTAR
KCNS
EPCON=USTAR3/VKC
  *** store fully-developed boundary-layer values
(stored of UIN is (:ustar:*loge(ZG/:Z0:)/:vkc:)!zslstr)
(stored of KEIN is :kcns:!zslstr)
(stored of EPIN is (:epcon:/ZG)!zslstr)
(stored of EVIN is :vkc:*:ustar:*ZG!zslstr)
(stored of OMIN is :epcon:/(ZG*0.09*:kcns:)!zslstr)
  ** discrepancy with log-law values
(stored of UDIS is 100.*(U1/UIN-1.0) with IF(IX.LT.NX))
(stored of KDIS is 100.*(KE/KEIN-1.0))
(stored of EVDS is 100.*(ENUT/EVIN-1.0))
    (stored of S is GEN1^0.5)
    (stored of SF2 is BF2*S)
    (stored of OLIM is 0.3*OMEG)
    (stored of SANY is STRS[&&1]/(0.41*ZG))
  save7end
 ************************************************************
  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)
 DISWAL
 DVO1DT =3.41E-03
 PRT(KE)=1.0E+10 ;PRT(OMEG)=1.0E+10
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(P1)=0. ;FIINIT(WDIS)=5.
 FIINIT(BF2)=1. ;FIINIT(BF1)=1.
   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
 
 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)
 
 BUOYD =1.189
 BUOYE =20.
 EGWF = T
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 200
 RESFAC =1.0E-05
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=200
 ************************************************************
  Group 17. Relaxation
 RELAX(P1 ,LINRLX,1. )
 RELAX(KE ,FALSDT,1. )
 RELAX(OMEG,FALSDT,1. )
 ************************************************************
  Group 18. Limits
 VARMAX(OMEG)=1.0E+10 ;VARMIN(OMEG)=1.0E-10
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 GENK = T
 PARSOL = F
 IENUTA = 19
 ISG62 = 1
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,MAXINC,KE,R,0.025)
 ************************************************************
  Group 20. Preliminary Printout
 DISTIL = T ;NULLPR = F
 NDST = 0
 DSTTOL =1.0E-02
 EX(P1)=8.046E-04 ;EX(U1)=5.768
 EX(V1)=8.26E-04 ;EX(W1)=3.82E-04
 EX(KE)=0.4166 ;EX(EP)=0.01092
 EX(SIGW)=1.996 ;EX(DUDZ)=0.08889
 EX(SIGK)=1.995 ;EX(LTLS)=2692.
 EX(WDIS)=38.25 ;EX(GEN1)=0.1035
 EX(BF2)=1. ;EX(BF1)=0.995
 EX(VABS)=5.768 ;EX(STRS)=3.23E-03
 EX(OMEG)=0.3072 ;EX(DEN1)=1.189
 EX(EL1)=15.64 ;EX(ENUT)=5.521
 EX(EVDS)=3.387E+00
 EX(KDIS)=2.433E+00
 EX(UDIS)=7.054E-01
 EX(EVIN)=5.531E+00
 EX(OMIN)=3.055E-01
 EX(EPIN)=1.141E-02
 EX(KEIN)=4.151E-01
 EX(UIN )=5.767E+00

 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(EP  ,Y,N,Y,Y,Y,Y)
 OUTPUT(WDIS,Y,N,N,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)
 ************************************************************
  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,0.,-1.,0.)
 GVIEW(UP,0.,0.,1.)
 GVIEW(VDIS,116.67971)
 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
STOP