NGEET · slevis-lmwg · Sep 4, 2019 · Sep 4, 2019 · Sep 5, 2019 · Sep 5, 2019
diff --git a/biogeochem/EDPatchDynamicsMod.F90 b/biogeochem/EDPatchDynamicsMod.F90
@@ -2124,6 +2124,12 @@ subroutine create_patch(currentSite, new_patch, age, areap, label,nocomp_pft)
  new_patch%fabi_sha_z(:,:,:) = 0._r8 
  new_patch%scorch_ht(:) = 0._r8 
  new_patch%frac_burnt = 0._r8 
+ new_patch%canopy_bulk_density = 0._r8
+ new_patch%canopy_fuel_load = 0._r8
+ new_patch%passive_crown_FI = 0._r8
+ new_patch%ros_torch = 0._r8
+ new_patch%heat_per_area = 0._r8
+ new_patch%lb = 0._r8
  new_patch%litter_moisture(:) = 0._r8
  new_patch%fuel_eff_moist = 0._r8
  new_patch%livegrass = 0._r8
@@ -2235,11 +2241,18 @@ subroutine zero_patch(cp_p)
  currentPatch%fi = nan ! average fire intensity of flaming front during day. 
  ! backward ros plays no role. kj/m/s or kw/m.
  currentPatch%fire = 999 ! sr decide_fire.1=fire hot enough to proceed. 0=stop everything- no fires today
+ currentPatch%active_crown_fire_flg = 9999 ! flag to indicate active crown fire ignition
  currentPatch%fd = nan ! fire duration (mins)
  currentPatch%ros_back = nan ! backward ros (m/min)
  currentPatch%scorch_ht(:) = nan ! scorch height of flames on a given PFT
- currentPatch%frac_burnt = nan ! fraction burnt daily 
- currentPatch%burnt_frac_litter(:) = nan 
+ currentPatch%frac_burnt = nan ! fraction burnt daily
+ currentPatch%canopy_bulk_density = nan ! available canopy fuel bulk density in patch (kg biomass/m3)
+ currentPatch%canopy_fuel_load = nan ! available canopy fuel load in patch (kg biomass)
+ currentPatch%passive_crown_FI = nan ! fire intensity for ignition of passive canopy fuel (kW/m)
+ currentPatch%ros_torch = nan ! rate of spread (ROS) for crown torch initiation (m/min)
+ currentPatch%heat_per_area = nan ! heat release per unit area (kJ/m2) for surface fuel
+ currentPatch%lb = nan ! length to breadth ratio of fire ellipse (unitless)
+ currentPatch%burnt_frac_litter(:) = nan
  currentPatch%btran_ft(:) = 0.0_r8
 
  currentPatch%canopy_layer_tlai(:) = 0.0_r8
@@ -2658,6 +2671,12 @@ subroutine fuse_2_patches(csite, dp, rp)
  rp%ros_back = (dp%ros_back*dp%area + rp%ros_back*rp%area) * inv_sum_area
  rp%scorch_ht(:) = (dp%scorch_ht(:)*dp%area + rp%scorch_ht(:)*rp%area) * inv_sum_area
  rp%frac_burnt = (dp%frac_burnt*dp%area + rp%frac_burnt*rp%area) * inv_sum_area
+ rp%canopy_bulk_density = (dp%canopy_bulk_density*dp%area + rp%canopy_bulk_density*rp%area) * inv_sum_area
+ rp%canopy_fuel_load = (dp%canopy_fuel_load*dp%area + rp%canopy_fuel_load*rp%area) * inv_sum_area
+ rp%passive_crown_FI = (dp%passive_crown_FI*dp%area + rp%passive_crown_FI*rp%area) * inv_sum_area
+ rp%ros_torch = (dp%ros_torch*dp%area + rp%ros_torch*rp%area) * inv_sum_area
+ rp%heat_per_area = (dp%heat_per_area*dp%area + rp%heat_per_area*rp%area) * inv_sum_area
+ rp%lb = (dp%lb*dp%area + rp%lb*rp%area) * inv_sum_area
  rp%burnt_frac_litter(:) = (dp%burnt_frac_litter(:)*dp%area + rp%burnt_frac_litter(:)*rp%area) * inv_sum_area
  rp%btran_ft(:) = (dp%btran_ft(:)*dp%area + rp%btran_ft(:)*rp%area) * inv_sum_area
  rp%zstar = (dp%zstar*dp%area + rp%zstar*rp%area) * inv_sum_area

diff --git a/fire/SFMainMod.F90 b/fire/SFMainMod.F90
diff --git a/main/EDInitMod.F90 b/main/EDInitMod.F90
@@ -642,6 +642,12 @@ subroutine init_patches( nsites, sites, bc_in)
  currentPatch%scorch_ht(:) = 0._r8
  currentPatch%frac_burnt = 0._r8
  currentPatch%burnt_frac_litter(:) = 0._r8
+ currentPatch%canopy_bulk_density = 0._r8
+ currentPatch%canopy_fuel_load = 0._r8
+ currentPatch%passive_crown_FI = 0._r8
+ currentPatch%ros_torch = 0._r8
+ currentPatch%heat_per_area = 0._r8
+ currentPatch%lb = 0._r8
 
  currentPatch => currentPatch%older
  enddo

diff --git a/main/EDParamsMod.F90 b/main/EDParamsMod.F90
@@ -61,10 +61,9 @@ module EDParamsMod
  real(r8),protected, public :: ED_val_patch_fusion_tol
  real(r8),protected, public :: ED_val_canopy_closure_thresh ! site-level canopy closure point where trees take on forest (narrow) versus savannah (wide) crown allometry
  integer,protected, public :: stomatal_model !switch for choosing between stomatal conductance models, 1 for Ball-Berry, 2 for Medlyn
-
- logical,protected, public :: active_crown_fire ! flag, 1=active crown fire 0=no active crown fire
- character(len=param_string_length),parameter :: fates_name_active_crown_fire = "fates_fire_active_crown_fire"
 
+ logical,protected, public :: active_crown_fire_switch ! flag, 1=active crown fire 0=no active crown fire
+ character(len=param_string_length),parameter :: fates_name_active_crown_fire = "fates_fire_active_crown_fire"
 
  real(r8), protected, public :: cg_strikes ! fraction of cloud to ground lightning strikes (0-1)
  character(len=param_string_length),parameter :: fates_name_cg_strikes="fates_fire_cg_strikes"
@@ -454,7 +453,7 @@ subroutine FatesRegisterParams(fates_params)
 
  call fates_params%RegisterParameter(name=fates_name_active_crown_fire, dimension_shape=dimension_shape_scalar, &
  dimension_names=dim_names_scalar)
- 
+
  call fates_params%RegisterParameter(name=fates_name_cg_strikes, dimension_shape=dimension_shape_scalar, &
  dimension_names=dim_names_scalar)
 
@@ -623,7 +622,7 @@ subroutine FatesReceiveParams(fates_params)
 
  call fates_params%RetreiveParameter(name=fates_name_active_crown_fire, & 
  data=tmpreal)
- active_crown_fire = (abs(tmpreal-1.0_r8)<nearzero)
+ active_crown_fire_switch = (abs(tmpreal-1.0_r8)<nearzero)
 
  call fates_params%RetreiveParameter(name=fates_name_cg_strikes, &
  data=cg_strikes)
@@ -707,7 +706,7 @@ subroutine FatesReportParams(is_master)
  write(fates_log(),fmt0) 'q10_mr = ',q10_mr
  write(fates_log(),fmt0) 'q10_froz = ',q10_froz
  write(fates_log(),fmt0) 'cg_strikes = ',cg_strikes
- write(fates_log(),'(a,L2)') 'active_crown_fire = ',active_crown_fire
+ write(fates_log(),'(a,L2)') 'active_crown_fire_switch = ',active_crown_fire_switch
  write(fates_log(),*) '------------------------------------------------------'
 
  end if

diff --git a/main/EDPftvarcon.F90 b/main/EDPftvarcon.F90
@@ -48,9 +48,8 @@ module EDPftvarcon
  real(r8), allocatable :: z0mr(:) ! ratio of roughness length of vegetation to height (-)
  real(r8), allocatable :: displar(:) ! ratio of displacement height to canopy top height
  real(r8), allocatable :: bark_scaler(:) ! scaler from dbh to bark thickness. For fire model.
- real(r8), allocatable :: crown_kill(:) ! scaler on fire death. For fire model.
- real(r8), allocatable :: initd(:) ! initial seedling density
-
+ real(r8), allocatable :: crown_kill(:) ! crown resistance to fire. (1 = none) For fire model.
+ real(r8), allocatable :: initd(:) ! initial seedling density 
  real(r8), allocatable :: seed_suppl(:) ! seeds that come from outside the gridbox.
  real(r8), allocatable :: bb_slope(:) ! ball berry slope parameter
  real(r8), allocatable :: medlyn_slope(:) ! Medlyn slope parameter KPa^0.5

diff --git a/main/EDTypesMod.F90 b/main/EDTypesMod.F90
@@ -156,6 +156,9 @@ module EDTypesMod
  integer, parameter, public :: dl_sf = 5 ! array index of dead leaf pool for spitfire (dead grass and dead leaves)
  integer, parameter, public :: lg_sf = 6 ! array index of live grass pool for spitfire
 
+! real(r8), parameter, public :: crown_fire_threshold = 200.0_r8 ! threshold for passive crown fire ignition. KWm-2 (Bessie & Johnson 1995) TODO slevis: see corresponding todo in SFMainMod.F90
+
+
  ! PATCH FUSION 
  real(r8), parameter, public :: force_patchfuse_min_biomass = 0.005_r8 ! min biomass (kg / m2 patch area) below which to force-fuse patches
  integer , parameter, public :: N_DBH_BINS = 6 ! no. of dbh bins used when comparing patches
@@ -562,8 +565,14 @@ module EDTypesMod
  real(r8) :: fuel_eff_moist ! effective avearage fuel moisture content of the ground fuel 
  ! (incl. live grasses. omits 1000hr fuels)
  real(r8) :: litter_moisture(nfsc)
+ real(r8) :: canopy_bulk_density ! available canopy fuel bulk density in patch (kg biomass/m3)
+ real(r8) :: canopy_fuel_load ! available canopy fuel load in patch (kg biomass)
+ real(r8) :: passive_crown_FI ! fire intensity for ignition of passive canopy fuel (kW/m)
+ real(r8) :: heat_per_area ! heat release per unit area (kJ/m2) for surface fuel
 
  ! FIRE SPREAD
+ real(r8) :: lb ! length to breadth ratio of fire ellipse (unitless)
+ real(r8) :: ros_torch ! rate of spread for crown torch initiation: m/min
  real(r8) :: ros_front ! rate of forward spread of fire: m/min
  real(r8) :: ros_back ! rate of backward spread of fire: m/min
  real(r8) :: effect_wspeed ! windspeed modified by fraction of relative grass and tree cover: m/min
@@ -574,9 +583,10 @@ module EDTypesMod
 
  ! FIRE EFFECTS 
  real(r8) :: scorch_ht(maxpft) ! scorch height: m 
- real(r8) :: frac_burnt ! fraction burnt: frac patch/day 
+ real(r8) :: frac_burnt ! fraction burnt: frac patch/day
  real(r8) :: tfc_ros ! total intensity-relevant fuel consumed - no trunks. KgC/m2 of burned ground/day
  real(r8) :: burnt_frac_litter(nfsc) ! fraction of each litter pool burned, conditional on it being burned
+ integer :: active_crown_fire_flg ! flag for active crown fire ignition
 
 
  ! PLANT HYDRAULICS (not currently used in hydraulics RGK 03-2018) 
@@ -1018,7 +1028,8 @@ subroutine dump_patch(cpatch)
  write(fates_log(),*) 'pa%c_lblayer = ',cpatch%c_lblayer
  write(fates_log(),*) 'pa%disturbance_rate = ',cpatch%disturbance_rate
  write(fates_log(),*) 'pa%disturbance_rates = ',cpatch%disturbance_rates(:)
- write(fates_log(),*) 'pa%anthro_disturbance_label = ',cpatch%anthro_disturbance_label
+ write(fates_log(),*) 'pa%anthro_disturbance_label = ', cpatch%anthro_disturbance_label
+ write(fates_log(),*) 'pa%active_crown_fire_flg = ', cpatch%active_crown_fire_flg
  write(fates_log(),*) '----------------------------------------'
  do el = 1,num_elements
  write(fates_log(),*) 'element id: ',element_list(el)

diff --git a/main/FatesHistoryInterfaceMod.F90 b/main/FatesHistoryInterfaceMod.F90
@@ -4742,7 +4742,7 @@ subroutine define_history_vars(this, initialize_variables)
 
  call this%set_history_var(vname='FATES_FUEL_AMOUNT_APFC', units='kg m-2', &
  long='spitfire fuel quantity in each age x fuel class in kg carbon per m2 land area', &
- use_default='inactive', avgflag='A', vtype=site_agefuel_r8, &
+ use_default='active', avgflag='A', vtype=site_agefuel_r8, &
  hlms='CLM:ALM', upfreq=1, ivar=ivar, initialize=initialize_variables, &
  index = ih_fuel_amount_age_fuel)
 
@@ -5515,7 +5515,7 @@ subroutine define_history_vars(this, initialize_variables)
 
  call this%set_history_var(vname='FATES_SCORCH_HEIGHT_APPF',units = 'm', &
  long='SPITFIRE flame Scorch Height (calculated per PFT in each patch age bin)', &
- use_default='inactive', avgflag='A', vtype=site_agepft_r8, &
+ use_default='active', avgflag='A', vtype=site_agepft_r8, &
  hlms='CLM:ALM', upfreq=1, ivar=ivar, &
  initialize=initialize_variables, index = ih_scorch_height_si_agepft)
 
@@ -5706,21 +5706,21 @@ subroutine define_history_vars(this, initialize_variables)
  call this%set_history_var(vname='FATES_MORTALITY_FIRE_SZPF', &
  units = 'm-2 yr-1', &
  long='fire mortality by pft/size in number of plants per m2 per year', &
- use_default='inactive', avgflag='A', vtype=site_size_pft_r8, &
+ use_default='active', avgflag='A', vtype=site_size_pft_r8, &
  hlms='CLM:ALM', upfreq=1, ivar=ivar, &
  initialize=initialize_variables, index = ih_m5_si_scpf)
 
  call this%set_history_var(vname='FATES_MORTALITY_CROWNSCORCH_SZPF', &
  units = 'm-2 yr-1', &
  long='fire mortality from crown scorch by pft/size in number of plants per m2 per year', &
- use_default='inactive', avgflag='A', vtype=site_size_pft_r8, &
+ use_default='active', avgflag='A', vtype=site_size_pft_r8, &
  hlms='CLM:ALM', upfreq=1, ivar=ivar, &
  initialize=initialize_variables, index = ih_crownfiremort_si_scpf)
 
  call this%set_history_var(vname='FATES_MORTALITY_CAMBIALBURN_SZPF', &
  units = 'm-2 yr-1', &
  long='fire mortality from cambial burn by pft/size in number of plants per m2 per year', &
- use_default='inactive', avgflag='A', vtype=site_size_pft_r8, &
+ use_default='active', avgflag='A', vtype=site_size_pft_r8, &
  hlms='CLM:ALM', upfreq=1, ivar=ivar, &
  initialize=initialize_variables, index = ih_cambialfiremort_si_scpf)
 

diff --git a/parameter_files/fates_params_default.cdl b/parameter_files/fates_params_default.cdl
@@ -206,7 +206,7 @@ variables:
  fates_fire_crown_depth_frac:long_name = "the depth of a cohorts crown as a fraction of its height" ;
  double fates_fire_crown_kill(fates_pft) ;
  fates_fire_crown_kill:units = "NA" ;
- fates_fire_crown_kill:long_name = "fire parameter, see equation 22 in Thonicke et al 2010" ;
+ fates_fire_crown_kill:long_name = "crown resistance to fire (1 = none), EQ 22 in Thonicke et al 2010" ;
  double fates_fnrt_prof_a(fates_pft) ;
  fates_fnrt_prof_a:units = "unitless" ;
  fates_fnrt_prof_a:long_name = "Fine root profile function, parameter a" ;