&HEAD CHID='pinendls',TITLE='Laboratory experiment with pine needle bed' / 
&GRID IBAR=80,JBAR=25,KBAR=45 /
&PDIM XBAR=4.0,YBAR0=-0.05, YBAR=1.2, ZBAR=2.25 / 
&TRNZ IDERIV=1,CC=0.,PC=0.5 /
&TRNZ IDERIV=2,CC=0.,PC=0.0 /

&MISC SURF_DEFAULT='INERT',
	REACTION='WOOD',
        TMPA = 20. / 
	DATABASE='database95.data',

&REAC ID='WOOD'
      FYI='Dupuy et al., Comb Flame, 135:65-76, 2003 dHc = 16MW/kg'
      RADIATIVE_FRACTION = 0.35
      HRRPUA_SHEET = 100.
      SOOT_YIELD = 0.02
      NU_O2      = 3.7
      NU_CO2     = 3.4
      NU_H2O     = 3.1
      MW_FUEL    = 87.
      EPUMO2     = 12300. /

&TIME TWFIN= 300 / Time when finished (length of simulation)

&SURF ID='BLOW',VEL=-0.00, TMPWAL=20.0,   RGB=0.8,0.6,0.4, PART_ID='cold air',/
&PART ID='cold air', QUANTITY='BLUE', MASSLESS=.TRUE., DTPAR=0.1 / 

--- BCs
&VENT XB= 0.0, 4.0, -0.05, 1.2, 2.25, 2.25, SURF_ID='OPEN' / top
&VENT XB= 0.0, 0.0, -0.05, 1.2, 0.0, 2.25, SURF_ID='OPEN' / left side
&VENT XB= 4.0, 4.0, -0.05, 1.2, 0.0, 2.25,  SURF_ID='OPEN'/ right side 
&VENT XB= 3.4, 4.0, -0.05, -0.05, 0.0, 0.95, SURF_ID='OPEN' / front
&VENT XB= 0.0, 4.0, -0.05, -0.05, 0.95, 2.25, SURF_ID='OPEN' / 
&VENT XB= 3.4, 4.0,  1.20,  1.20, 0.0, 0.95, SURF_ID='OPEN' / back
&VENT XB= 0.0, 4.0,  1.20,  1.20, 0.95, 2.25, SURF_ID='OPEN' / 

-- Nuno M. G. Santos parameters (personal communication) except 
   for vegetation_element_density which is from Dupuy et al. 
   Comb. Flame, 135:65-76, 2003
&SURF ID='NEEDLES',
      BURN_SURF_VEG = .TRUE.,
      VEGETATION_CHAR_FRACTION = 0.25,
      VEGETATION_ELEMENT_DENSITY = 640.,
      VEGETATION_SVRATIO = 3500.,
      VEGETATION_MOISTURE = 0.11,
      VEGETATION_HEIGHT = 0.03,
      VEGETATION_LOAD = 0.586,
      VEGETATION_CDRAG = 1.
      TMPIGN=320.,
      EMISSIVITY = 0.99,
      VBC = -1.,
      RGB = 0.48,0.46,0.19 /
&VENT XB=1.4,3.4,0.25,0.95,0.0,0.0,SURF_ID="NEEDLES"/

-- Burner ignitor

&SURF ID='BURNER',HRRPUA=2400.,RAMP_Q='RAMPBURNER'/ 
&RAMP ID='RAMPBURNER',T=0.0,F=0.0 /
&RAMP ID='RAMPBURNER',T=0.0,F=1.0 /
&RAMP ID='RAMPBURNER',T=8.0,F=1.0 /
&RAMP ID='RAMPBURNER',T=9.0,F=0.0 /
&VENT XB= 1.35, 1.40, 0.25, 0.95, 0.00, 0.00, SURF_ID='BURNER'/ Burner on floor.

-- Ouputs

&SLCF PBY=0.6,QUANTITY='TEMPERATURE',VECTOR=.TRUE. /
&SLCF PBZ=0.30,QUANTITY='TEMPERATURE',VECTOR=.TRUE. /

&BNDF QUANTITY='HEAT_FLUX' /
&BNDF QUANTITY='BURNING_RATE' /
&BNDF QUANTITY='WALL_TEMPERATURE' /