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" ;