TALK=T;RUN( 1, 1)
************************************************************
Q1 created by VDI menu, Version 2022, Date 21/03/22
CPVNAM=VDI; SPPNAM=Core
************************************************************
IRUNN = 1 ;LIBREF = 0
************************************************************
Group 1. Run Title
TEXT(Case6b-3d-Zup-Yflo-PasqF3-Buoss )
************************************************************
Group 2. Transience
STEADY = T
************************************************************
Groups 3, 4, 5 Grid Information
* Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
RSET(M,4,10,8)
************************************************************
Group 6. Body-Fitted coordinates
************************************************************
Group 7. Variables: STOREd,SOLVEd,NAMEd
* Non-default variable names
NAME(136)=TREF ;NAME(137)=ZETA
NAME(138)=PSIF ;NAME(139)=PSIU
NAME(140)=EPSA ;NAME(141)=TKEA
NAME(142)=UANA ;NAME(143)=VIN
NAME(144)=TANA ;NAME(145)=TIN
NAME(146)=TEM1 ;NAME(147)=EPKE
NAME(148)=DEN1 ;NAME(149)=EL1
NAME(150) =ENUT
* Solved variables list
SOLVE(P1,U1,V1,W1,TEM1)
* Stored variables list
STORE(ENUT,EL1,DEN1,EPKE,TIN,VIN)
STORE(TREF)
* Additional solver options
SOLUTN(P1,Y,Y,Y,N,N,Y)
SOLUTN(TEM1,Y,Y,Y,N,N,Y)
TURMOD(KEMODL)
************************************************************
Group 8. Terms & Devices
************************************************************
Group 9. Properties
PRESS0 =1.01325E+05 ;TEMP0 =273.149994
* Domain material index is 0 signifying:
* Air at 20 deg C, 1 atm, treated as incompressible
SETPRPS(1, 0)
RHO1 =1.2265
DVO1DT =3.41E-03
PRNDTL(TEM1)=-0.0258
PRT(EP)=1.314
************************************************************
Group 10.Inter-Phase Transfer Processes
************************************************************
Group 11.Initialise Var/Porosity Fields
FIINIT(P1)=0. ;FIINIT(TEM1)=15.
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(BUOYANCY, PHASEM, 0, 0, 0, 0, 0, 0, 1, 1)
COVAL(BUOYANCY, U1, FIXFLU, GRND3)
COVAL(BUOYANCY, V1, FIXFLU, GRND3)
COVAL(BUOYANCY, W1, FIXFLU, GRND3)
BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81
BUOYE =15.
EGWF = T
************************************************************
Group 14. Downstream Pressure For PARAB
************************************************************
Group 15. Terminate Sweeps
LSWEEP = 1000
RESFAC =1.0E-05
************************************************************
Group 16. Terminate Iterations
LITER(KE)=10 ;LITER(EP)=10
************************************************************
Group 17. Relaxation
RELAX(P1 ,LINRLX,1. )
RELAX(KE ,LINRLX,0.5 )
RELAX(EP ,LINRLX,0.5 )
RELAX(TEM1,LINRLX,0.5 )
KELIN = 3
************************************************************
Group 18. Limits
VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-204.862488
************************************************************
Group 19. EARTH Calls To GROUND Station
GENK = T
PARSOL = F
ISG62 = 1
SPEDAT(SET,OUTPUT,NOFIELD,L,T)
SPEDAT(SET,GXMONI,PLOTALL,L,T)
SPEDAT(SET,RLXFAC,REFVEL,R,3.)
SPEDAT(SET,GXMONI,MINIRES,L,F)
************************************************************
Group 20. Preliminary Printout
************************************************************
Group 21. Print-out of Variables
OUTPUT(TANA,Y,N,Y,N,Y,Y)
OUTPUT(TIN ,Y,N,Y,N,Y,Y)
************************************************************
Group 22. Monitor Print-Out
IXMON = 1 ;IYMON = 2 ;IZMON = 1
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.996459,0.070419,0.045942)
GVIEW(UP,-0.072785,0.448896,0.890615)
GVIEW(DEPTH,3.0E+04)
GVIEW(VDIS,433.565491)
GVIEW(CENTRE,300.,300.,35.)
> DOM, SIZE, 6.000000E+02, 6.000000E+02, 7.000000E+01
> DOM, MONIT, 7.500000E+01, 9.000000E+01, 4.375000E+00
> DOM, SCALE, 1.000000E+00, 1.000000E+00, 1.000000E+00
> DOM, INCREMENT, 1.000000E+00, 1.000000E+00, 9.999999E-02
> GRID, RSET_X_1, -4, 1.000000E+00,G
> GRID, RSET_Y_1, 10, 1.000000E+00,G
> GRID, RSET_Z_1, 8, 1.000000E+00,G
> DOM, T_AMBIENT, 1.500000E+01
> DOM, INI_AMB, YES
> DOM, INI_BUOY, YES
> OBJ, NAME, WIND1
> OBJ, POSITION, AT_START, AT_START, AT_START
> OBJ, SIZE, TO_END, TO_END, TO_END
> OBJ, DOMCLIP, NO
> OBJ, GEOMETRY, windgr
> OBJ, ROTATION24, 7
> OBJ, TYPE, WIND
> OBJ, COLOR-MODE, FROM_FIL
> OBJ, COLOR-VAL, 144
> OBJ, WIREFRAME, YES
> OBJ, PRESSURE, 1.01325E+05
> OBJ, COEFFICIENT, 1000.
> OBJ, VELOCITY, 3.
> OBJ, WIND_DIR, 180.
> OBJ, AXIS_DIR, 0.
> OBJ, TEMPERATURE, 15.
> OBJ, PROFILE, Logarithmic
> OBJ, REF_HEIGHT, 10.
> OBJ, RGHNS_HEIGHT,0.1
> OBJ, UP-DIR, Z
> OBJ, GROUND-TEMP, 15
> OBJ, PASQUILL, Pasquill F - Moderately stable
> OBJ, MONIN-OBUKHOV, Golder
STOP
goto end
** Replicate here Wind Object settings for this case
REAL(QREF,ZO,ZREF)
ZO=0.1 ! roughness height
QREF=3.0 ! wind speed at reference height
ZREF=10.0 ! reference height
** avoid buoyancy sequences for other 2 directions
COVAL(BUOYANCY,U1,ZERO,ZERO )
COVAL(BUOYANCY,V1,ZERO,ZERO )
! STORE(RHIN,PIN,TREF,TIN) ! Compulsory
! STORE(UIN,VIN,WIN,KEIN,EPIN,OMIN) ! Optional
** Integer control variable for temperature profiles
Default value is 0 for uniform temperature
INTEGER(ITPRO); ITPRO = 0
** Integer control variable for Monin Obukhov length
MONIN = 0 user-specified value
= 1 TNO formulae (default)
= 2 PHAST formulae
INTEGER(MOLEN); MOLEN = 2
REAL(GT0,GZT0,GALR,GQWALL,AMO,BMO)
GT0=15.0 ; ! surface temperature
SPEDAT(SET,BLIN,GT0,R,GT0)
** Select Paquill Stability Class
REAL(LS,ZS,ZODZS,ZOMO,GLMO)
CHAR(PASQSC);PASQSC=D
PASQSC=F
CASE :PASQSC: OF
WHEN A,1
** Extremely unstable
+ MESG(Pasquill A - Extremely unstable
LS=33.162;ZS=1117.0 ! TNO
AMO=-11.4;BMO=0.1 ! PHAST
ITPRO=1
WHEN B,1
** Moderately unstable
+ MESG(Pasquill B - Moderately unstable
LS=32.258;ZS=11.46 ! TNO
AMO=-26.0;BMO=0.17 ! PHAST
ITPRO=2
WHEN C,1
** Slightly unstable
+ MESG(Pasquill C - Slightly unstable
LS=51.787;ZS=1.324 ! TNO
AMO=-123.0;BMO=0.3 ! PHAST
ITPRO=3
WHEN D,1
** Neutral
+ MESG(Pasquill D - Neutral
LS=1.E10;ZS=1.0
ITPRO=4 ! uniform or logarithmic
WHEN E,1
** Slightly stable
+ MESG(Pasquill E - Slightly stable
LS=-48.33;ZS=1.262 ! TNO
AMO=123.0;BMO=0.3 ! PHAST
ITPRO=5
WHEN F,1
** Moderately unstable
+ MESG(Pasquill F - Moderately stable
LS=-31.325;ZS=19.36 ! TNO
AMO=26.0;BMO=0.17 ! PHAST
ITPRO=6
ENDCASE
SPEDAT(SET,BLIN,PASQSC,C,PASQSC)
SPEDAT(SET,BLIN,ITPRO,I,ITPRO)
SPEDAT(SET,BLIN,MOLEN,I,MOLEN)
** Compute the Monin-Obukhov length
** TNO formula
IF(MOLEN.EQ.1) THEN
IF(ZO.GT.0.5) THEN
+ ZOMO=0.5 ! TNO limitation of ZO for LMO
ELSE
+ ZOMO=ZO
ENDIF
ZOMO
ZODZS=ZOMO/ZS
GLMO = LS/LOG10(ZODZS)
ENDIF
** PHAST formula
IF(MOLEN.EQ.2) THEN
+ GLMO=AMO*(ZO**BMO)
ENDIF
GLMO
** temperature profile
SPEDAT(SET,BLIN,ITPRO,I,ITPRO)
GZT0=0.0 ; SPEDAT(SET,BLIN,GZT0,R,0.)
GALR=-9.81/CP1 ! adiabatic lapse rate
GALR
** analytical profiles
REAL(QTAU,AKA,PSIUR,ZRDZO,AMO,BMO)
REAL(TSTAR,TSDK,USDK,TKCON,EPCON,ZETAR,RHOAMB)
GT0
AKA=0.41;ZRDZO=ZREF/ZO
Stability Class F - analytical profiles
ZETAR=ZREF/GLMO
ZETAR
PSIUR=-5.*ZETAR
REAL(LOGZDZ)
LOGZDZ=LOG(ZRDZO)
LOGZDZ
QTAU= AKA*QREF/(LOG(ZRDZO)-PSIUR); USDK=QTAU/AKA
QTAU
RHOAMB=PRESS0/(286.7*(GT0+TEMP0))
RHOAMB
** Compute Ground heat flux, then set
TEM1 boundary condition
GQWALL=-(QTAU**3)*(GT0+TEMP0)*RHOAMB*CP1/(AKA*9.81*GLMO)
GQWALL
COVAL(BLIN6,TEM1,FIXFLU,GQWALL)
SPEDAT(SET,BLIN,GQWALL,R,GQWALL)
TSTAR=-GQWALL/(RHOAMB*CP1*QTAU) ; TSDK=TSTAR/AKA
TSTAR
TKCON=QTAU*QTAU
EPCON=QTAU**3/AKA
(stored of ZG3D is ZG)
(stored of ZETA is ZG/GLMO)
(stored of PSIU is -5.*ZETA)
(stored of PHIF is 1.+5.*ZETA)
(stored of UANA is USDK*(LOGE(ZG/ZO)-PSIU))
(stored of TANA is GT0 + GALR*(ZG-GZT0)+TSDK*(LOGE(ZG/ZO)-PSIU))
(stored of TKEA is TKCON*(1.-ZETA/PHIF)^0.5)
(stored of EPSA is EPCON*PHIF/ZG)
EX(PHIF)=1.096E+01
EX(ZG3D)=3.500E+01
EX(ZETA)=1.991E+00
EX(EPSA)=3.757E-03
EX(TKEA)=2.482E-02
EX(UANA)=6.262E+00
EX(TANA)=1.887E+01
EX(PSIU)=9.956E+00
EX(PSIF)=1.000E-10
label end
DISTIL=T
EX(P1 )=3.896E-02
EX(U1 )=6.787E-03
EX(V1 )=6.506E+00
EX(W1 )=1.120E-01
EX(KE )=2.843E-01
EX(EP )=8.852E-03
EX(TREF)=1.887E+01
EX(VIN )=4.329E+00
EX(TIN )=1.242E+01
EX(TEM1)=1.904E+01
EX(EPKE)=3.042E-02
EX(DEN1)=1.227E+00
EX(EL1 )=2.915E+00
EX(ENUT)=8.968E-01