&HEAD CHID='firec064_1grid',TITLE='Cheney/Gould grass fire case C064' /

&GRID IBAR=180,JBAR=180,KBAR=72 /
&PDIM XBAR0=0,XBAR=300.,YBAR0=-150,YBAR= 150,ZBAR0=0.0,ZBAR=80.0 /
&TRNZ IDERIV=1,CC=0.,PC=0.5 /
&TRNZ IDERIV=2,CC=0.,PC=0.0 /

&TIME TWFIN=125. /

&MISC NSTRATA=1,REACTION='WOOD',RADIATION=.TRUE.,BAROCLINIC=.TRUE.,
      DT0DZ=-0.01,U0=0.,TMPA=31.5,
      UZ_PROFILE=.TRUE.,UZ_ZS=7.,ZS_UZ=2.,EXP_UZ=0.143 /

&RADI NUMBER_RADIATION_ANGLES=200 /

- Specify parameters for combustion of fuel gases
  from pyrolysis of the solid fuel(s)
&REAC ID='WOOD'
      FYI='Ritchie, et al., 5th IAFSS, C_3.4 H_6.2 O_2.5, dHc = 15MW/kg'
       RADIATIVE_FRACTION = 0.35
       HRRPUA_SHEET = 200.
       SOOT_YIELD = 0.02
       NU_O2      = 3.7
       NU_CO2     = 3.4
       NU_H2O     = 3.1
       MW_FUEL    = 87.
       EPUMO2     = 11962. /
 
- Themally thin homogeneous surface fuel
&SURF ID='NO BURN GRASS FUEL SURFACE',
      NOBURN_SURF_VEG = .TRUE.,
      VEGETATION_ELEMENT_DENSITY = 512.,
      VEGETATION_SVRATIO = 9770.,
      VEGETATION_HEIGHT = 0.21,
      VEGETATION_LOAD = 0.283,
      VEGETATION_CDRAG = 0.41,
      VBC = -1,
      RGB = 0.141,0.180,0.086 / 

&SURF ID='BURN GRASS FUEL SURFACE',
      BURN_SURF_VEG = .TRUE.,
      VEGETATION_CHAR_FRACTION = 0.2,
      VEGETATION_ELEMENT_DENSITY = 512.,
      VEGETATION_SVRATIO = 9770.,
      VEGETATION_MOISTURE = 0.063,
      VEGETATION_HEIGHT = 0.21,
      VEGETATION_LOAD = 0.283,
      VEGETATION_CDRAG = 0.41,
      HEAT_OF_COMBUSTION=15585.,
      TMPIGN=320.,
      EMISSIVITY = 0.99,
      VBC = -1,
      RGB = 0.48,0.46,0.19 /

&VENT XB= 100,202, -54,  54,0,0,SURF_ID='BURN GRASS FUEL SURFACE'/
&VENT XB=  73,100, -79,  79,0,0,SURF_ID='NO BURN GRASS FUEL SURFACE'/
&VENT XB= 100,202,  54,  79,0,0,SURF_ID='NO BURN GRASS FUEL SURFACE'/
&VENT XB= 202,227, -79,  79,0,0,SURF_ID='NO BURN GRASS FUEL SURFACE'/
&VENT XB= 100,202, -79, -54,0,0,SURF_ID='NO BURN GRASS FUEL SURFACE'/
 
- Time spatially dependent ignitor fire
- Positive y locations of ignitor
&SURF ID='IGN FIRE P1',HRRPUA=1460.,RAMP_Q='RAMPFIRE P1' 
                    NOBURN_SURF_VEG=.TRUE.,
                    VEGETATION_ELEMENT_DENSITY = 512.,
                    VEGETATION_SVRATIO = 9770.,
                    VEGETATION_HEIGHT = 0.21,
                    VEGETATION_LOAD = 0.283,
                    VBC = -1,
                    VEGETATION_CDRAG = 0.41 /
&RAMP ID='RAMPFIRE P1',T=0.0,F=0.0 /
&RAMP ID='RAMPFIRE P1',T=1.0,F=1.0 /
&RAMP ID='RAMPFIRE P1',T=6.0,F=1.0 /
&RAMP ID='RAMPFIRE P1',T=7.0,F=0.0 /
&VENT XB=100,101, 0,5,0.0,0.0,SURF_ID='IGN FIRE P1' /

&SURF ID='IGN FIRE P2',HRRPUA=1460.,RAMP_Q='RAMPFIRE P2' 
                    NOBURN_SURF_VEG=.TRUE.,
                    VEGETATION_ELEMENT_DENSITY = 512.,
                    VEGETATION_SVRATIO = 9770.,
                    VEGETATION_HEIGHT = 0.21,
                    VEGETATION_LOAD = 0.283,
                    VBC = -1,
                    VEGETATION_CDRAG = 0.41 /
&RAMP ID='RAMPFIRE P2',T=5.0,F=0.0 /
&RAMP ID='RAMPFIRE P2',T=6.0,F=1.0 /
&RAMP ID='RAMPFIRE P2',T=11.0,F=1.0 /
&RAMP ID='RAMPFIRE P2',T=12.0,F=0.0 /
&VENT XB=100,101,5,10,0.0,0.0,SURF_ID='IGN FIRE P2' /

&SURF ID='IGN FIRE P3',HRRPUA=1460.,RAMP_Q='RAMPFIRE P3' 
                    NOBURN_SURF_VEG=.TRUE.,
                    VEGETATION_ELEMENT_DENSITY = 512.,
                    VEGETATION_SVRATIO = 9770.,
                    VEGETATION_HEIGHT = 0.21,
                    VEGETATION_LOAD = 0.283,
                    VBC = -1,
                    VEGETATION_CDRAG = 0.41 /
&RAMP ID='RAMPFIRE P3',T=10.0,F=0.0 /
&RAMP ID='RAMPFIRE P3',T=11.0,F=1.0 /
&RAMP ID='RAMPFIRE P3',T=16.0,F=1.0 /
&RAMP ID='RAMPFIRE P3',T=17.0,F=0.0 /
&VENT XB=100,101,10,15,0.0,0.0,SURF_ID='IGN FIRE P3' /

&SURF ID='IGN FIRE P4',HRRPUA=1460.,RAMP_Q='RAMPFIRE P4' 
                    NOBURN_SURF_VEG=.TRUE.,
                    VEGETATION_ELEMENT_DENSITY = 512.,
                    VEGETATION_SVRATIO = 9770.,
                    VEGETATION_HEIGHT = 0.21,
                    VEGETATION_LOAD = 0.283,
                    VBC = -1,
                    VEGETATION_CDRAG = 0.41 /
&RAMP ID='RAMPFIRE P4',T=15.0,F=0.0 /
&RAMP ID='RAMPFIRE P4',T=16.0,F=1.0 /
&RAMP ID='RAMPFIRE P4',T=21.0,F=1.0 /
&RAMP ID='RAMPFIRE P4',T=22.0,F=0.0 /
&VENT XB=100,101,15,20,0.0,0.0,SURF_ID='IGN FIRE P4' /

&SURF ID='IGN FIRE P5',HRRPUA=1460.,RAMP_Q='RAMPFIRE P5' 
                    NOBURN_SURF_VEG=.TRUE.,
                    VEGETATION_ELEMENT_DENSITY = 512.,
                    VEGETATION_SVRATIO = 9770.,
                    VEGETATION_HEIGHT = 0.21,
                    VEGETATION_LOAD = 0.283,
                    VBC = -1,
                    VEGETATION_CDRAG = 0.41 /
&RAMP ID='RAMPFIRE P5',T=20.0,F=0.0 /
&RAMP ID='RAMPFIRE P5',T=21.0,F=1.0 /
&RAMP ID='RAMPFIRE P5',T=26.0,F=1.0 /
&RAMP ID='RAMPFIRE P5',T=27.0,F=0.0 /
&VENT XB=100,101,20,25,0.0,0.0,SURF_ID='IGN FIRE P5' /
 
- Negative y locations of ignitor
&SURF ID='IGN FIRE M1',HRRPUA=1460.,RAMP_Q='RAMPFIRE M1' 
                    NOBURN_SURF_VEG=.TRUE.,
                    VEGETATION_ELEMENT_DENSITY = 512.,
                    VEGETATION_SVRATIO = 9770.,
                    VEGETATION_HEIGHT = 0.21,
                    VEGETATION_LOAD = 0.283,
                    VBC = -1,
                    VEGETATION_CDRAG = 0.41 /
&RAMP ID='RAMPFIRE M1',T=0.0,F=0.0 /
&RAMP ID='RAMPFIRE M1',T=1.0,F=1.0 /
&RAMP ID='RAMPFIRE M1',T=6.0,F=1.0 /
&RAMP ID='RAMPFIRE M1',T=7.0,F=0.0 /
&VENT XB=100,101,-5,0,0.0,0.0,SURF_ID='IGN FIRE M1' /

&SURF ID='IGN FIRE M2',HRRPUA=1460.,RAMP_Q='RAMPFIRE M2' 
                    NOBURN_SURF_VEG=.TRUE.,
                    VEGETATION_ELEMENT_DENSITY = 512.,
                    VEGETATION_SVRATIO = 9770.,
                    VEGETATION_HEIGHT = 0.21,
                    VEGETATION_LOAD = 0.283,
                    VBC = -1,
                    VEGETATION_CDRAG = 0.41 /
&RAMP ID='RAMPFIRE M2',T=5.0,F=0.0 /
&RAMP ID='RAMPFIRE M2',T=6.0,F=1.0 /
&RAMP ID='RAMPFIRE M2',T=11.0,F=1.0 /
&RAMP ID='RAMPFIRE M2',T=12.0,F=0.0 /
&VENT XB=100,101,-10,-5,0.0,0.0,SURF_ID='IGN FIRE M2' /

&SURF ID='IGN FIRE M3',HRRPUA=1460.,RAMP_Q='RAMPFIRE M3' 
                    NOBURN_SURF_VEG=.TRUE.,
                    VEGETATION_ELEMENT_DENSITY = 512.,
                    VEGETATION_SVRATIO = 9770.,
                    VEGETATION_HEIGHT = 0.21,
                    VEGETATION_LOAD = 0.283,
                    VBC = -1,
                    VEGETATION_CDRAG = 0.41 /
&RAMP ID='RAMPFIRE M3',T=10.0,F=0.0 /
&RAMP ID='RAMPFIRE M3',T=11.0,F=1.0 /
&RAMP ID='RAMPFIRE M3',T=16.0,F=1.0 /
&RAMP ID='RAMPFIRE M3',T=17.0,F=0.0 /
&VENT XB=100,101,-15,-10,0.0,0.0,SURF_ID='IGN FIRE M3' /

&SURF ID='IGN FIRE M4',HRRPUA=1460.,RAMP_Q='RAMPFIRE M4' 
                    NOBURN_SURF_VEG=.TRUE.,
                    VEGETATION_ELEMENT_DENSITY = 512.,
                    VEGETATION_SVRATIO = 9770.,
                    VEGETATION_HEIGHT = 0.21,
                    VEGETATION_LOAD = 0.283,
                    VBC = -1,
                    VEGETATION_CDRAG = 0.41 /
&RAMP ID='RAMPFIRE M4',T=15.0,F=0.0 /
&RAMP ID='RAMPFIRE M4',T=16.0,F=1.0 /
&RAMP ID='RAMPFIRE M4',T=21.0,F=1.0 /
&RAMP ID='RAMPFIRE M4',T=22.0,F=0.0 /
&VENT XB=100,101,-20,-15,0.0,0.0,SURF_ID='IGN FIRE M4' /

&SURF ID='IGN FIRE M5',HRRPUA=1460.,RAMP_Q='RAMPFIRE M5' 
                    NOBURN_SURF_VEG=.TRUE.,
                    VEGETATION_ELEMENT_DENSITY = 512.,
                    VEGETATION_SVRATIO = 9770.,
                    VEGETATION_HEIGHT = 0.21,
                    VEGETATION_LOAD = 0.283,
                    VBC = -1,
                    VEGETATION_CDRAG = 0.41 /
&RAMP ID='RAMPFIRE M5',T=20.0,F=0.0 /
&RAMP ID='RAMPFIRE M5',T=21.0,F=1.0 /
&RAMP ID='RAMPFIRE M5',T=26.0,F=1.0 /
&RAMP ID='RAMPFIRE M5',T=27.0,F=0.0 /
&VENT XB=100,101,-25,-20,0.0,0.0,SURF_ID='IGN FIRE M5' /
 
- Inflow
&SURF ID='INFLOW',VEL=-7,RAMP_V='RAMPVEL',PROFILE='ATMOSPHERIC',Z0=2.,PLE=0.143 /
&RAMP ID='RAMPVEL',T=0.0,F=1.0 /
&RAMP ID='RAMPVEL',T=0.5,F=1.0 /
 
- Boundary conditions
&VENT XB=  0,  0,-150, 150,  0,80,SURF_ID='INFLOW',VENT_COLOR='INVISIBLE' /
&VENT XB=300,300,-150, 150,  0,80,SURF_ID='OPEN' /
&VENT XB=  0,300,-150,-150,  0,80,SURF_ID='OPEN' /
&VENT XB=  0,300, 150, 150,  0,80,SURF_ID='OPEN' /
&VENT XB=  0,300,-150, 150, 80,80,SURF_ID='OPEN' /
 
- Outputs
&SLCF PBY= 0, QUANTITY='TEMPERATURE',VECTOR=.TRUE.,DTSAM=0.1 / 
&SLCF PBY= 0, QUANTITY='MIXTURE_FRACTION',DTSAM=0.1 / 
&SLCF PBZ= 1.11, QUANTITY='U-VELOCITY',DTSAM=0.1/ 
&SLCF PBZ= 2.22, QUANTITY='U-VELOCITY',DTSAM=0.1/ 
&SLCF PBZ= 4.49, QUANTITY='U-VELOCITY',DTSAM=0.1/ 
&SLCF PBZ= 2, QUANTITY='TEMPERATURE',DTSAM=0.1/ 
&SLCF PBZ= 6, QUANTITY='TEMPERATURE',DTSAM=0.1/ 

&BNDF QUANTITY='INCIDENT_HEAT_FLUX',DTSAM=0.1/
&BNDF QUANTITY='HEAT_FLUX',DTSAM=0.1/
&BNDF QUANTITY='RADIATIVE_FLUX',DTSAM=0.1/
&BNDF QUANTITY='CONVECTIVE_FLUX',DTSAM=0.1/
&BNDF QUANTITY='WALL_TEMPERATURE',DTSAM=0.1/
&BNDF QUANTITY='BURNING_RATE',DTSAM=0.1/
&BNDF QUANTITY='MASS_LOSS',DTSAM=0.1/
&PL3D DTSAM=50 /