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(Case5a-3d-Yup-Zflo-PasqF3-DensD ) ************************************************************ Group 2. Transience STEADY = T ************************************************************ Groups 3, 4, 5 Grid Information * Overall number of cells, RSET(M,NX,NY,NZ,tolerance) RSET(M,4,8,10) ************************************************************ Group 6. Body-Fitted coordinates ************************************************************ Group 7. Variables: STOREd,SOLVEd,NAMEd * Non-default variable names NAME(134)=RHIN ;NAME(135)=PIN NAME(136)=TREF ;NAME(137)=ZETA NAME(138)=PSIF ;NAME(139)=PSIU NAME(140)=EPSA ;NAME(141)=TKEA NAME(142)=UANA ;NAME(143)=WIN 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,WIN) STORE(PSIF,TREF,PIN,RHIN) * 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 NEWRH1 = T NEWENL = T ************************************************************ Group 9. Properties PRESS0 =1.01325E+05 ;TEMP0 =273.149994 * Domain material index is 2 signifying: * Air using Ideal Gas Law, STP SETPRPS(1, 2) DRH1DP = GRND5 DVO1DT =3.33E-03 PRNDTL(TEM1)=-0.0263 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, GRND2) COVAL(BUOYANCY, V1, FIXFLU, GRND2) COVAL(BUOYANCY, W1, FIXFLU, GRND2) BUOYA =0. ; BUOYB =-9.81 ; BUOYC =0. BUOYD =1.225219 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 LITER(TEM1)=50 ************************************************************ 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(P1)=1.0E+10 ;VARMIN(P1)=-7.599375E+04 VARMAX(TEM1)=3000. ;VARMIN(TEM1)=-30. ************************************************************ 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 = 1 ;IZMON = 2 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.792077,0.462073,0.398877) GVIEW(UP,0.33566,0.875484,-0.347649) GVIEW(DEPTH,3.0E+04) GVIEW(VDIS,408.677094) GVIEW(CENTRE,300.,35.,300.) > DOM, SIZE, 6.000000E+02, 7.000000E+01, 6.000000E+02 > DOM, MONIT, 7.500000E+01, 4.375000E+00, 9.000000E+01 > 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, 4, 1.000000E+00,G > GRID, RSET_Y_1, 8, 1.000000E+00,G > GRID, RSET_Z_1, 10, 1.000000E+00 > 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, 3 > OBJ, SELECTABLE, NO > 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, 0. > OBJ, AXIS_DIR, 90. > OBJ, TEMPERATURE, 15. > OBJ, PROFILE, Logarithmic > OBJ, REF_HEIGHT, 10. > OBJ, RGHNS_HEIGHT,0.1 > OBJ, UP-DIR, Y > OBJ, GROUND-TEMP, 15 > OBJ, PASQUILL, Pasquill F - Moderately stable > OBJ, MONIN-OBUKHOV, Golder STOP DISTIL=T EX(RHIN)=1.221E+00 EX(PIN )=2.097E+02 EX(PSIF)=1.000E-10 EX(EPKE)=3.090E-02 EX(DEN1)=1.222E+00 EX(WIN )=3.946E+00 EX(P1 )=1.440E-01 EX(U1 )=4.758E-02 EX(V1 )=1.165E-01 EX(W1 )=6.629E+00 EX(KE )=4.117E-01 EX(EP )=1.317E-02 EX(TREF)=1.888E+01 EX(TIN )=1.242E+01 EX(TEM1)=1.911E+01 EX(EL1 )=3.076E+00 EX(ENUT)=1.189E+00