&HEAD CHID='2d_VegaHigh_HeadAndBack_U1ms',TITLE='Boundary Fuel model' /
&TIME TWFIN=300 /

- 2D
&GRID IBAR=200,JBAR=1,KBAR=100 /
&PDIM XBAR0=0,XBAR=300,YBAR0=0,YBAR=3,ZBAR0=0.0,ZBAR=120.0 /
&TRNZ IDERIV=1,CC=0.,PC=0.5 /
&TRNZ IDERIV=2,CC=0.,PC=0.0 /

- Boundary condtions
&VENT XB=0,0,0,3,0,120,SURF_ID='INFLOW',VENT_COLOR='INVISIBLE' /
cVENT XB=0,0,0,3,0,120,SURF_ID='OPEN' /
&VENT XB=300,300,0,3,0,120,SURF_ID='OPEN' /
&VENT XB=0,300,0,3,120,120,SURF_ID='OPEN' /

&MISC NSTRATA=1,REACTION='WOOD',RADIATION=.TRUE.,BAROCLINIC=.TRUE.,
      DT0DZ=-0.01,U0=0.,TMPA=20.0,
      UZ_PROFILE=.TRUE.,UZ_ZS=1.,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'
       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. /
 
&PART ID='from grass',MASSLESS=.TRUE.,QUANTITY='GREEN' /
 
- For fuel break
&SURF ID='NO BURN SHRUB VEGA',
      NOBURN_SURF_VEG = .TRUE.,
      VEGETATION_ELEMENT_DENSITY = 512.,
      VEGETATION_SVRATIO = 5000.,
      VEGETATION_HEIGHT = 0.55,
      VEGETATION_LOAD = 1.89,
      VEGETATION_CDRAG = 1,
      VBC=-1,
      RGB=0.141,0.644,0.523 /
&VENT XB=  0, 99,0,3,0,0,SURF_ID='NO BURN SHRUB VEGA' /
&VENT XB=198,300,0,3,0,0,SURF_ID='NO BURN SHRUB VEGA' /

- Themally thin homogeneous surface fuel Cheney F19
&SURF ID='SHRUB VEGA',
      BURN_SURF_VEG = .TRUE.,
      VEGETATION_CHAR_FRACTION = 0.2,
      VEGETATION_ELEMENT_DENSITY = 512.,
      VEGETATION_SVRATIO = 5000.,
      VEGETATION_MOISTURE = 0.08,
      VEGETATION_HEIGHT = 0.55,
      VEGETATION_LOAD = 1.89,
      VEGETATION_CDRAG = 1,
      HEAT_OF_COMBUSTION=15585.,
      HEAT_OF_VAPORIZATION=200.,
      TMPIGN=320.,
      EMISSIVITY = 0.99,
      VBC=-1, 
      RGB = 0.48,0.46,0.19,
      PART_ID='from grass'/
&VENT XB=99,198,0,3,0,0,SURF_ID='SHRUB VEGA' /
 
- Ignitor fire
&SURF ID='IGN FIRE',HRRPUA=1460.,RAMP_Q='RAMPFIRE' 
                    NOBURN_SURF_VEG=.TRUE.,
                    VEGETATION_ELEMENT_DENSITY = 512.,
                    VEGETATION_SVRATIO = 12240.,
                    VEGETATION_HEIGHT = 0.51,
                    VEGETATION_LOAD = 0.313,
                    VEGETATION_CDRAG = 1,
                    VBC=-1, /
&RAMP ID='RAMPFIRE',T= 0.0,F=0.0 /
&RAMP ID='RAMPFIRE',T=10.0,F=0.0 /
&RAMP ID='RAMPFIRE',T=11.0,F=1.0 /
&RAMP ID='RAMPFIRE',T=19.0,F=1.0 /
&RAMP ID='RAMPFIRE',T=20.0,F=0.0 /
&VENT XB= 99,102,0,3,0.0,0.0,SURF_ID='IGN FIRE' / head fire
&VENT XB=195,198,0,3,0.0,0.0,SURF_ID='IGN FIRE' / back fire
 
- Inflow
&PART ID='inflow particles',SAMPLING_FACTOR=2,MASSLESS=.TRUE.,QUANTITY='BLUE' /
&SURF ID='INFLOW',VEL=-1,PROFILE='ATMOSPHERIC',RAMP_V='RAMPVEL',Z0=2.,PLE=0.143
                  PART_ID='inflow particles' /
cSURF ID='INFLOW',VEL=-3,RAMP_V='RAMPVEL'  /
&RAMP ID='RAMPVEL',T=0.0,F=1.0 /
&RAMP ID='RAMPVEL',T=0.5,F=1 /
 
- Outputs
&SLCF PBY= 3, QUANTITY='TEMPERATURE',VECTOR=.TRUE.,DTSAM=0.1 / 
&SLCF PBY= 3, QUANTITY='HRRPUV',DTSAM=0.1/ 
&SLCF PBY= 3, QUANTITY='MIXTURE_FRACTION',DTSAM=0.1 / 
&SLCF PBY= 3, QUANTITY='oxygen mass fraction',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 /

     - temperature and heat fluxes at back fire ignition
&OBST XB= 196.5,198,0,3,1.8,2.4 /
&THCP XYZ=197,2,2,QUANTITY='TEMPERATURE',LABEL='Gas_T1' /
&THCP XYZ=197,2,2,QUANTITY='WALL_TEMPERATURE',IOR=3,LABEL='Wall_T1_Z+' /
&THCP XYZ=197,2,2.4,QUANTITY='RADIATIVE_FLUX',IOR=3,LABEL='Qrad_1_Z+' /
&THCP XYZ=197,2,2.4,QUANTITY='CONVECTIVE_FLUX',IOR=3,LABEL='Qc_1_Z+' /
&THCP XYZ=196.5,2,2.0,QUANTITY='RADIATIVE_FLUX',IOR=-1,LABEL='Qrad_1_X-' /
&THCP XYZ=196.5,2,2.0,QUANTITY='CONVECTIVE_FLUX',IOR=-1,LABEL='Qc_1_X-' /

     - temperature and heat fluxes about 5 m downwind of back fire ignition
&OBST XB= 199.5,201,0,3,1.8,2.4 /
&THCP XYZ=200,2,2,QUANTITY='TEMPERATURE',LABEL='Gas_T2' /
&THCP XYZ=200,2,2,QUANTITY='WALL_TEMPERATURE',IOR=3,LABEL='Wall_T2_Z+' /
&THCP XYZ=200,2,2.4,QUANTITY='RADIATIVE_FLUX',IOR=3,LABEL='Qrad_2_Z+' /
&THCP XYZ=200,2,2.4,QUANTITY='CONVECTIVE_FLUX',IOR=3,LABEL='Qc_2_Z+' /
&THCP XYZ=199.5,2,2.0,QUANTITY='RADIATIVE_FLUX',IOR=-1,LABEL='Qrad_2_X-' /
&THCP XYZ=199.5,2,2.0,QUANTITY='CONVECTIVE_FLUX',IOR=-1,LABEL='Qc_2_X-' /

&PLT3D DTSAM=25. /