Active crown scott 2001

biogeochem/EDPatchDynamicsMod.F90

@@ -2500,6 +2500,7 @@ 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%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

biogeochem/FatesPatchMod.F90

@@ -196,6 +196,8 @@ module FatesPatchMod
     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)   ! moisture of litter [m3/m3]
+    real(r8)              :: canopy_bulk_density 
+
 
     ! fire spread
     real(r8)              :: ros_front               ! rate of forward  spread of fire [m/min]
@@ -204,6 +206,7 @@ module FatesPatchMod
     real(r8)              :: tau_l                   ! duration of lethal heating [min]
     real(r8)              :: fi                      ! average fire intensity of flaming front [kJ/m/s] or [kW/m]
     integer               :: fire                    ! is there a fire? [1=yes; 0=no]
+    integer               ::  active_crown_fire_flg
     real(r8)              :: fd                      ! fire duration [min]
 
     ! fire effects      
@@ -382,9 +385,11 @@ subroutine NanValues(this)
       this%tau_l                        = nan
       this%fi                           = nan 
       this%fire                         = fates_unset_int
+      this%active_crown_fire_flg        = fates_unset_int
       this%fd                           = nan 
       this%scorch_ht(:)                 = nan 
       this%frac_burnt                   = nan
+      this%canopy_bulk_density          = nan 
       this%tfc_ros                      = nan    
       this%burnt_frac_litter(:)         = nan    
 
@@ -459,7 +464,8 @@ subroutine ZeroValues(this)
       this%fi                                = 0.0_r8
       this%fd                                = 0.0_r8
       this%scorch_ht(:)                      = 0.0_r8  
-      this%frac_burnt                        = 0.0_r8  
+      this%frac_burnt                        = 0.0_r8
+      this%canopy_bulk_density               = 0.0_r8  
       this%tfc_ros                           = 0.0_r8
       this%burnt_frac_litter(:)              = 0.0_r8
 
@@ -732,6 +738,7 @@ subroutine Dump(this)
       write(fates_log(),*) 'pa%c_stomata          = ',this%c_stomata
       write(fates_log(),*) 'pa%c_lblayer          = ',this%c_lblayer
       write(fates_log(),*) 'pa%disturbance_rates  = ',this%disturbance_rates(:)
+      write(fates_log(),*) 'pa%active_crown_fire_flg = ', this%active_crown_fire_flg
       write(fates_log(),*) 'pa%anthro_disturbance_label = ',this%anthro_disturbance_label
       write(fates_log(),*) '----------------------------------------'
 

main/EDParamsMod.F90

@@ -605,9 +605,6 @@ subroutine FatesRegisterParams(fates_params)
     call fates_params%RegisterParameter(name=ED_name_history_height_bin_edges, dimension_shape=dimension_shape_1d, &
          dimension_names=dim_names_height)
 
-    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)
 
@@ -815,10 +812,6 @@ subroutine FatesReceiveParams(fates_params)
     call fates_params%RetrieveParameter(name=name_dev_arbitrary, &
          data=dev_arbitrary)
 
-    call fates_params%RetrieveParameter(name=fates_name_active_crown_fire, & 
-          data=tmpreal)
-    active_crown_fire = (abs(tmpreal-1.0_r8)<nearzero)
-
     call fates_params%RetrieveParameter(name=fates_name_cg_strikes, &
           data=cg_strikes)
 

main/EDPftvarcon.F90

@@ -50,9 +50,9 @@ 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 :: active_crown_fire(:)   ! Is plant susceptible to active crown fire?
+     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
@@ -385,6 +385,10 @@ subroutine Register_PFT(this, fates_params)
     call fates_params%RegisterParameter(name=name, dimension_shape=dimension_shape_1d, &
          dimension_names=dim_names, lower_bounds=dim_lower_bound)
 
+    name = 'fates_fire_active_crown_fire'
+    call fates_params%RegisterParameter(name=name, dimension_shape=dimension_shape_1d, &
+         dimension_names=dim_names, lower_bounds=dim_lower_bound)
+
     name = 'fates_recruit_init_density'
     call fates_params%RegisterParameter(name=name, dimension_shape=dimension_shape_1d, &
          dimension_names=dim_names, lower_bounds=dim_lower_bound)
@@ -826,6 +830,10 @@ subroutine Receive_PFT(this, fates_params)
     call fates_params%RetrieveParameterAllocate(name=name, &
          data=this%crown_kill)
 
+    name = 'fates_fire_active_crown_fire'
+    call fates_params%RetrieveParameterAllocate(name=name, &
+         data=this%active_crown_fire)
+
     name = 'fates_recruit_init_density'
     call fates_params%RetrieveParameterAllocate(name=name, &
          data=this%initd)
@@ -1651,6 +1659,7 @@ subroutine FatesReportPFTParams(is_master)
         write(fates_log(),fmt0) 'displar = ',EDPftvarcon_inst%displar
         write(fates_log(),fmt0) 'bark_scaler = ',EDPftvarcon_inst%bark_scaler
         write(fates_log(),fmt0) 'crown_kill = ',EDPftvarcon_inst%crown_kill
+        write(fates_log(),fmt0) 'active_crown_fire = ',EDPftvarcon_inst%active_crown_fire
         write(fates_log(),fmt0) 'initd = ',EDPftvarcon_inst%initd
         write(fates_log(),fmt0) 'seed_suppl = ',EDPftvarcon_inst%seed_suppl
         write(fates_log(),fmt0) 'bb_slope = ',EDPftvarcon_inst%bb_slope
@@ -1918,7 +1927,7 @@ subroutine FatesCheckParams(is_master)
         ! xl must be between -0.4 and 0.6 according to Bonan (2019) doi:10.1017/9781107339217 pg. 238
         !-----------------------------------------------------------------------------------
         if (EDPftvarcon_inst%xl(ipft) < -0.4 .or. EDPftvarcon_inst%xl(ipft) > 0.6) then
-          write(fates_log(),*) 'fates_rad_leaf_xl for pft ', ipft, ' is outside the allowed range of -0.6 to 0.4'
+          write(fates_log(),*) 'fates_rad_leaf_xl for pft ', ipft, ' is outside the allowed range of -0.4 to 0.6'
           write(fates_log(),*) 'Aborting'
           call endrun(msg=errMsg(sourcefile, __LINE__))
         end if 

main/EDTypesMod.F90

@@ -101,6 +101,9 @@ module EDTypesMod
   integer, parameter, public :: phen_dstat_timeon    = 3       ! Leaves on due to time exceedance  (drought phenology)
   integer, parameter, public :: phen_dstat_pshed     = 4 ! Leaves partially abscissing       (drought phenology)
 
+  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
   real(r8), parameter, public :: patch_fusion_tolerance_relaxation_increment = 1.1_r8 ! amount by which to increment patch fusion threshold

main/FatesHistoryInterfaceMod.F90

@@ -6926,7 +6926,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)
 
@@ -7124,14 +7124,14 @@ subroutine define_history_vars(this, initialize_variables)
     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)