TALK=T;RUN( 1, 1)
 
 ************************************************************
   Q1 created by VDI menu, Version 2020, Date 16/06/20
 CPVNAM=VDI; SPPNAM=FLAIR
 ************************************************************
  Echo DISPLAY / USE settings
 ************************************************************
 IRUNN = 1 ;LIBREF = 2
 ************************************************************
  Group 1. Run Title
 TEXT(Fire-Spray in a compartment             )
 ************************************************************
  Group 2. Transience
 STEADY=F
    * Set overall time and no. of steps
 RSET(U,0.,120.,12)
    * Modify regions
 ************************************************************
  Groups 3, 4, 5  Grid Information
    * Overall number of cells, RSET(M,NX,NY,NZ,tolerance)
 RSET(M,31,21,17)
 ************************************************************
  Group 6. Body-Fitted coordinates
 ************************************************************
  Group 7. Variables: STOREd,SOLVEd,NAMEd
    * Non-default variable names
 NAME(128)=SLEN ;NAME(130)=ENUL
 NAME(137)=MASS ;NAME(138)=HEAT
 NAME(139)=MOMZ ;NAME(140)=MOMY
 NAME(141)=MOMX ;NAME(142)=EPKE
 NAME(143)=DEN1 ;NAME(144)=EL1
 NAME(145)=ENUT ;NAME(146)=SPH1
 NAME(147)=TEM1 ;NAME(149)=SMOK
 NAME(150) =MH2O
    * Solved variables list
 SOLVE(P1,U1,V1,W1,TEM1,SMOK,MH2O)
    * Stored variables list
 STORE(SPH1,ENUT,EL1,DEN1,EPKE,MOMX,MOMY,MOMZ)
 STORE(HEAT,MASS,ENUL,SLEN)
    * Additional solver options
 SOLUTN(P1,Y,Y,Y,N,N,Y)
 SOLUTN(TEM1,Y,Y,Y,N,N,Y)
 SOLUTN(SMOK,Y,Y,Y,N,N,Y)
 SOLUTN(MH2O,Y,Y,Y,N,N,N)
 TURMOD(KECHEN)
 
 ************************************************************
  Group 8. Terms & Devices
 NEWRH1 = T
 NEWENL = T
 ************************************************************
  Group 9. Properties
 PRESS0 =1.01325E+05 ;TEMP0 =273.
    * 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(KE)=0.75 ;PRT(EP)=1.15
 ************************************************************
  Group 10.Inter-Phase Transfer Processes
 ************************************************************
  Group 11.Initialise Var/Porosity Fields
 FIINIT(P1)=0. ;FIINIT(U1)=0.
 FIINIT(V1)=0. ;FIINIT(W1)=0.
 FIINIT(KE)=1.0E-03 ;FIINIT(EP)=9.0E-03
 FIINIT(TEM1)=20. ;FIINIT(SMOK)=0.
   No PATCHes used for this Group
 
 
 INIADD = F
 NAMFI   ='stdy'
 ************************************************************
  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, 12)
 COVAL(BUOYANCY, U1, FIXFLU, GRND2)
 COVAL(BUOYANCY, V1, FIXFLU, GRND2)
 COVAL(BUOYANCY, W1, FIXFLU, GRND2)
 
 PATCH(GENPAT, CELL, 0, 0, 0, 0, 0, 0, 1, 12)
 COVAL(GENPAT, U1, FIXFLU, GRND)
 COVAL(GENPAT, V1, FIXFLU, GRND)
 COVAL(GENPAT, W1, FIXFLU, GRND)
 COVAL(GENPAT, TEM1, FIXFLU, GRND)
 COVAL(GENPAT, MH2O, FIXFLU, GRND)
 
 PATCH(GENMAS, CELL, 0, 0, 0, 0, 0, 0, 1, 12)
 COVAL(GENMAS, P1, FIXFLU, GRND)
 
 BUOYA =0. ; BUOYB =0. ; BUOYC =-9.81
 BUOYD =1.204938
 EGWF = T
 ************************************************************
  Echo save-block settings for Group 13
  SAVE13BEGIN
  SAVE13END
 ************************************************************
  Group 14. Downstream Pressure For PARAB
 ************************************************************
  Group 15. Terminate Sweeps
 LSWEEP = 100
 RESFAC =1.0E-03
 ************************************************************
  Group 16. Terminate Iterations
 LITER(P1)=200 ;LITER(TEM1)=50
 ************************************************************
  Group 17. Relaxation
 RELAX(P1 ,LINRLX,1. )
 RELAX(KE ,LINRLX,0.5 )
 RELAX(EP ,LINRLX,0.5 )
 RELAX(MASS,LINRLX,0.3 )
 RELAX(HEAT,LINRLX,0.3 )
 RELAX(MOMZ,LINRLX,0.3 )
 RELAX(MOMY,LINRLX,0.3 )
 RELAX(MOMX,LINRLX,0.3 )
 RELAX(DEN1,LINRLX,0.3 )
 RELAX(TEM1,LINRLX,0.25 )
 RELAX(SMOK,LINRLX,0.5 )
 RELAX(MH2O,LINRLX,0.3 )
 KELIN = 3
 ************************************************************
  Group 18. Limits
 VARMAX(P1)=1.0E+10 ;VARMIN(P1)=-7.5E+04
 VARMAX(TEM1)=3000. ;VARMIN(TEM1)=8.
 VARMAX(SMOK)=1. ;VARMIN(SMOK)=0.
 VARMAX(MH2O)=1. ;VARMIN(MH2O)=0.
 ************************************************************
  Group 19. EARTH Calls To GROUND Station
 NAMGRD =FLAR
 GENK = T
 PARSOL = F
 CONWIZ = T
 IENUTA = 2
 ISG62 = 1
 SPEDAT(SET,FACETDAT,NUMOBJ,I,14)
 SPEDAT(SET,SLEN,ACON,R,3.)
 SPEDAT(SET,SLEN,BCON,R,150.)
 SPEDAT(SET,SLEN,CCON,R,30.)
 SPEDAT(SET,GXMONI,PLOTALL,L,T)
 SPEDAT(SET,IMMERSOL,EMCON,R,5.)
 SPEDAT(SET,FLAIR,HCOMB,R,1.3E+07)
 SPEDAT(SET,FLAIR,RADFAC,R,0.)
 SPEDAT(SET,FLAIR,SYIELD,R,0.02)
 SPEDAT(SET,FLAIR,STOIC,R,0.992366)
 SPEDAT(SET,FLAIR,KEXT,R,7600.)
 SPEDAT(SET,FLAIR,COYIELD,R,0.06)
 SPEDAT(SET,MAXINC,TEM1,R,10.)
 SPEDAT(SET,MAXINC,MH2O,R,0.5)
L($G001
 GENTR = T
  *------------------------------------------------------
  * GENTRA GROUP 1: Particle physics
  *------------------------------------------------------
  * Particle type -  60
 GPTYPE = 60
  * Gravity components in GENTRA Cartesian system
 GGRAX =0. ;GGRAY =0. ;GGRAZ =-9.81
  * Buoyancy forces
 GBUOYA = T ;GSURPR = F
  * Stochastic model of turbulence
 GSTOCH = F
  * Data for vaporising particles
  * Data for melting/solidifying particles
 GDRAG = GRND1 ;GKONC =0.0263 ;GNUSS = GRND1
 GKONV = GRND1 ;GCPCON =1007. ;GCPVAP = GRND1
 GMWCON =28.9 ;GMWVAP =18. ;GLIQST =273.149994
 GCPLIQ =4131.799805 ;GHLIQD = GRND1 ;GVAPST = GRND1
 GLATVP = GRND1 ;GSTPRE = GRND1 ;GDTRCT =0.
  *-------------------------------------------------------
  * GENTRA GROUP 2: Boundary conditions for particles
  *-------------------------------------------------------
  * Inlet-data file-name
 GINFIL  ='spin'
  * Wall treatment, and rest coefficient if appropriate
 GWALLC = 1
  * Porosity threshold
 GPOROS =0.
 GTSTRT =0. ;GTEND =0. ;GFLENG =0.
  *-------------------------------------------------------
  * GENTRA GROUP 3: Numerical controls
  *-------------------------------------------------------
  * 1st GENTRA sweep; frequency of calls
 GSWEP1 = 1 ;GSWEPF = 1
  * Maximum Lagrangian time-step; time step size multplier
 GDTMAX =1. ;GRTFRC =1.
  * Min # of t-steps per cell; max # of t-steps; timeout
 GLAGTS = 5 ;GSTEMX = 100
 GTIMMX =100.
  *-------------------------------------------------------
  * GENTRA GROUP 4: Output controls
  *-------------------------------------------------------
  * Restart-file, history-file and frequency for output
 GRSFIL  ='NONE'
 GHFILE  ='GHIS'
 GOUTFR = 1
  * The identifier of the individual history and
    trajectory files
 GH1STC  ='NONE'
 GT1STC  ='NONE'
  * The first, last particles and the interval for
    writing history and trajectory files
 NGWSTR = 0 ;NGWEND = 0 ;NGWINT = 1
 GSWOUT = 1
  * Particle number for residence-time calculation
 GRESTI = 0
L(GENSET
 ************************************************************
  Group 20. Preliminary Printout
 ************************************************************
  Group 21. Print-out of Variables
 OUTPUT(ENUL,Y,N,N,N,N,N)
 ************************************************************
  Group 22. Monitor Print-Out
 IXMON = 14 ;IYMON = 11 ;IZMON = 8
 NPRMON = 100000
 NPRMNT = 1
 TSTSWP = -1
 ************************************************************
  Group 23.Field Print-Out & Plot Control
 NPRINT = 100000
 NTPRIN = 1 ;ISTPRF = 1 ;ISTPRL = 10000
 ISWPRF = 1 ;ISWPRL = 100000
   No PATCHes used for this Group
 ************************************************************
  Group 24. Dumps For Restarts
 IDISPA = 1 ;IDISPB = 0 ;IDISPC = 0
 CSG1    ='M'
 
 GVIEW(P,0.15279,-0.927805,0.340344)
 GVIEW(UP,-0.053857,0.336058,0.9403)
 GVIEW(VDIS,2.946434)
 GVIEW(CENTRE,2.5,1.4,1.1)
 
> DOM,    SIZE,        5.000000E+00, 2.800000E+00, 2.200000E+00
> DOM,    MONIT,       1.592500E+00, 1.400000E+00, 9.259092E-01
> DOM,    SCALE,       1.000000E+00, 1.000000E+00, 1.000000E+00
> GRID,   MINCELL,     2.5E-03 ,5.0E-03 ,5.0E-03
> GRID,   BOUNDS,       T F T T T T
> GRID,   RSET_X_1,    -10, 1.200000E+00
> GRID,   RSET_X_2,     -2, 1.000000E+00
> GRID,   RSET_X_3,     -1, 1.000000E+00
> GRID,   RSET_X_4,     -2, 1.000000E+00
> GRID,   RSET_X_5,    -10, 1.200000E+00
> GRID,   RSET_X_6,     -1, 1.000000E+00
> GRID,   RSET_X_7,      5, 1.500000E+00
> GRID,   RSET_Y_1,     -7, 1.000000E+00
> GRID,   RSET_Y_2,     -1, 1.000000E+00
> GRID,   RSET_Y_3,     -2, 1.000000E+00
> GRID,   RSET_Y_4,     -1, 1.000000E+00
> GRID,   RSET_Y_5,     -2, 1.000000E+00
> GRID,   RSET_Y_6,     -1, 1.000000E+00
> GRID,   RSET_Y_7,     -7, 1.200000E+00,G
> GRID,   RSET_Z_1,      3, 1.000000E+00
> GRID,   RSET_Z_2,     11, 1.000000E+00
> GRID,   RSET_Z_3,      2, 1.000000E+00
> GRID,   RSET_Z_4,      1, 1.000000E+00
> DOM,    INI_AMB,    YES
> DOM,    INI_BUOY,   YES
 
> OBJ,    NAME,        BLK1
> OBJ,    POSITION,    2.800000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        1.000000E-01, 1.030000E+00, TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    box
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    199,Solid allowing fluid-slip at walls
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
 
> OBJ,    NAME,        BLK2
> OBJ,    POSITION,    2.800000E+00, AT_END,       0.000000E+00
> OBJ,    SIZE,        1.000000E-01, 1.050000E+00, TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    box
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    199,Solid allowing fluid-slip at walls
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
 
> OBJ,    NAME,        BLK3
> OBJ,    POSITION,    2.800000E+00, 1.030000E+00, AT_END
> OBJ,    SIZE,        1.000000E-01, 7.200000E-01, 3.700000E-01
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    box
> OBJ,    TYPE,        BLOCKAGE
> OBJ,    MATERIAL,    199,Solid allowing fluid-slip at walls
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
 
> OBJ,    NAME,        OPN1
> OBJ,    POSITION,    2.900000E+00, AT_END,       0.000000E+00
> OBJ,    SIZE,        TO_END,       0.000000E+00, TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    box
> OBJ,    TYPE,        OPENING
> OBJ,    PRESSURE,     P_AMBIENT
> OBJ,    TEMPERATURE,  T_AMBIENT
> OBJ,    COEFFICIENT, 1.
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    TURBULENCE,  1.0E-03 ,9.0E-03
 
> OBJ,    NAME,        OPN2
> OBJ,    POSITION,    2.900000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        TO_END,       0.000000E+00, TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    box
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        OPENING
> OBJ,    PRESSURE,     P_AMBIENT
> OBJ,    TEMPERATURE,  T_AMBIENT
> OBJ,    COEFFICIENT, 1.
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    TURBULENCE,  1.0E-03 ,9.0E-03
 
> OBJ,    NAME,        OPN3
> OBJ,    POSITION,    AT_END,       0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, TO_END,       TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    box
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        OPENING
> OBJ,    PRESSURE,     P_AMBIENT
> OBJ,    TEMPERATURE,  T_AMBIENT
> OBJ,    COEFFICIENT, 1.
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    TURBULENCE,  1.0E-03 ,9.0E-03
 
> OBJ,    NAME,        OPN4
> OBJ,    POSITION,    2.900000E+00, 0.000000E+00, AT_END
> OBJ,    SIZE,        TO_END,       TO_END,       0.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    box
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        OPENING
> OBJ,    PRESSURE,     P_AMBIENT
> OBJ,    TEMPERATURE,  T_AMBIENT
> OBJ,    COEFFICIENT, 1.
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    TURBULENCE,  1.0E-03 ,9.0E-03
 
> OBJ,    NAME,        WALL1
> OBJ,    POSITION,    0.000000E+00, AT_END,       0.000000E+00
> OBJ,    SIZE,        2.800000E+00, 0.000000E+00, TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    box
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        PLATE
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    SURF_TEMP,   0. ,20.
 
> OBJ,    NAME,        WALL2
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        2.800000E+00, 0.000000E+00, 2.200000E+00
> OBJ,    GEOMETRY,    box
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        PLATE
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    SURF_TEMP,   0. ,20.
 
> OBJ,    NAME,        WALL3
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        0.000000E+00, TO_END,       TO_END
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    box
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        PLATE
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    SURF_TEMP,   0. ,20.
 
> OBJ,    NAME,        FLOOR
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, 0.000000E+00
> OBJ,    SIZE,        TO_END,       TO_END,       0.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    box
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        PLATE
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    SURF_TEMP,   0. ,20.
 
> OBJ,    NAME,        CEILING
> OBJ,    POSITION,    0.000000E+00, 0.000000E+00, AT_END
> OBJ,    SIZE,        2.800000E+00, TO_END,       0.000000E+00
> OBJ,    DOMCLIP,     NO
> OBJ,    GEOMETRY,    box
> OBJ,    VISIBLE,     NO
> OBJ,    TYPE,        PLATE
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    SURF_TEMP,   0. ,20.
 
> OBJ,    NAME,        FIRE
> OBJ,    POSITION,    1.250000E+00, 1.150000E+00, 0.000000E+00
> OBJ,    SIZE,        5.000000E-01, 5.000000E-01, 3.000000E-01
> OBJ,    GEOMETRY,    box
> OBJ,    TYPE,        FIRE
> OBJ,    INFSRC_P1,   (exp(-500.*MH2O)*46.9*tim)*tim*(1+rox)/hcom$
> OBJ,    INFSRC_P1,   b with whol
> OBJ,    INFSRC_TEM1, (exp(-500.*MH2O)*46.9*tim)*tim with whol
> OBJ,    INFSRC_SMOK, 1 with onlyms
> OBJ,    TIME_LIMITS, ALWAYS_ACTIVE
> OBJ,    PRE-TEMP,     T_AMBIENT
> OBJ,    MASS-SOURCE,   No Source
> OBJ,    HEAT-SOURCE,   Fixed Power
> OBJ,    FXD-FLU,     0.
> OBJ,    SCALAR-SOURCE, No Source
> OBJ,    INLET_SMOK,  1.
 
> OBJ,    NAME,        SPRAY
> OBJ,    POSITION,    1.440000E+00, 1.350000E+00, 2.100000E+00
> OBJ,    SIZE,        1.000000E-01, 1.000000E-01, 0.000000E+00
> OBJ,    GEOMETRY,    spray
> OBJ,    TYPE,        SPRAY-HEAD
> OBJ,    TIME_SPRAY,  120.
> OBJ,    LINK_T,       YES
> OBJ,    T_ACTIVATE,  50.
> OBJ,    RTI,         160.
> OBJ,    AXIS,          Z
> OBJ,    ORIGIN,      1.49 ,1.4 ,2.1
> OBJ,    RADIUS,      0.05
> OBJ,    PORTS,       8.
> OBJ,    TOTAL-VOL,   60.
> OBJ,    TOTAL-VEL,   1.
> OBJ,    ANGLE,       135.
> OBJ,    TEMPERATURE, 8.
> OBJ,    DMEAN,       0.8
> OBJ,    NSIZE,       1.
> OBJ,    DMIN,        0.
> OBJ,    DMAX,        2.5
> OBJ,    SPREAD_EXP,  1.7
STOP