Added crown area by size class and PFT to the output.

mpaiao · Dec 6, 2023 · bf99977 · bf99977
1 parent 1874511
commit bf99977
Showing 1 changed file with 24 additions and 0 deletions.
diff --git a/main/FatesHistoryInterfaceMod.F90 b/main/FatesHistoryInterfaceMod.F90
@@ -277,6 +277,9 @@ module FatesHistoryInterfaceMod
   !  Size-class x PFT LAI states
   integer :: ih_lai_canopy_si_scpf
   integer :: ih_lai_understory_si_scpf
+  !  Size-class x PFT LAI states
+  integer :: ih_crownarea_canopy_si_scpf
+  integer :: ih_crownarea_understory_si_scpf
 
   integer :: ih_totvegc_scpf
   integer :: ih_leafc_scpf
@@ -2318,6 +2321,8 @@ subroutine update_history_dyn(this,nc,nsites,sites,bc_in)
                hio_bleaf_understory_si_scpf  => this%hvars(ih_bleaf_understory_si_scpf)%r82d, &
                hio_lai_canopy_si_scpf        => this%hvars(ih_lai_canopy_si_scpf)%r82d, &
                hio_lai_understory_si_scpf    => this%hvars(ih_lai_understory_si_scpf)%r82d, &
+               hio_crownarea_canopy_si_scpf     => this%hvars(ih_crownarea_canopy_si_scpf)%r82d, &
+               hio_crownarea_understory_si_scpf => this%hvars(ih_crownarea_understory_si_scpf)%r82d, &
                hio_mortality_canopy_si_scpf         => this%hvars(ih_mortality_canopy_si_scpf)%r82d, &
                hio_mortality_canopy_secondary_si_scls      => this%hvars(ih_mortality_canopy_secondary_si_scls)%r82d, &
                hio_mortality_understory_si_scpf     => this%hvars(ih_mortality_understory_si_scpf)%r82d, &
@@ -3359,6 +3364,9 @@ subroutine update_history_dyn(this,nc,nsites,sites,bc_in)
                      leaf_m * ccohort%n / m2_per_ha
                   hio_lai_canopy_si_scpf(io_si,scpf) = hio_lai_canopy_si_scpf(io_si,scpf) + &
                       ccohort%treelai*ccohort%c_area * AREA_INV
+                  hio_crownarea_canopy_si_scpf(io_si,scpf) = hio_crownarea_canopy_si_scpf(io_si,scpf) + &
+                      ccohort%c_area * AREA_INV
+
 
                   hio_canopy_biomass_si(io_si) = hio_canopy_biomass_si(io_si) + n_perm2 * total_m
 
@@ -3502,6 +3510,8 @@ subroutine update_history_dyn(this,nc,nsites,sites,bc_in)
                      n_perm2 * total_m
                   hio_lai_understory_si_scpf(io_si,scpf) = hio_lai_understory_si_scpf(io_si,scpf) + &
                       ccohort%treelai*ccohort%c_area  * AREA_INV
+                  hio_crownarea_understory_si_scpf(io_si,scpf) = hio_crownarea_understory_si_scpf(io_si,scpf) + &
+                      ccohort%c_area * AREA_INV
 
                   !hio_mortality_understory_si_scpf(io_si,scpf) = hio_mortality_understory_si_scpf(io_si,scpf)+ &
                   !    (ccohort%bmort + ccohort%hmort + ccohort%cmort +
@@ -7227,6 +7237,13 @@ subroutine define_history_vars(this, initialize_variables)
           ivar=ivar, initialize=initialize_variables,                          &
           index = ih_lai_canopy_si_scpf )          
 
+    call this%set_history_var(vname='FATES_CROWNAREA_CANOPY_SZPF',             &
+          units = 'm2 m-2',                                                    &
+          long='Total crown area of canopy plants by pft/size',                &
+          use_default='inactive', avgflag='A', vtype=site_size_pft_r8,         &
+          hlms='CLM:ALM', upfreq=1, ivar=ivar,                                 &
+          initialize=initialize_variables, index = ih_crownarea_canopy_si_scpf )
+
     call this%set_history_var(vname='FATES_NPLANT_CANOPY_SZPF', units = 'm-2', &
           long='number of canopy plants by size/pft per m2',                   &
           use_default='inactive', avgflag='A', vtype=site_size_pft_r8,         &
@@ -7262,6 +7279,13 @@ subroutine define_history_vars(this, initialize_variables)
           hlms='CLM:ALM', upfreq=1, ivar=ivar,                                 &
           initialize=initialize_variables, index = ih_lai_understory_si_scpf )
 
+    call this%set_history_var(vname='FATES_CROWNAREA_USTORY_SZPF',             &
+          units = 'm2 m-2',                                                    &
+          long='Total crown area of understory plants by pft/size',            &
+          use_default='inactive', avgflag='A', vtype=site_size_pft_r8,         &
+          hlms='CLM:ALM', upfreq=1, ivar=ivar,                                 &
+          initialize=initialize_variables, index = ih_crownarea_understory_si_scpf )
+
     call this%set_history_var(vname='FATES_NPLANT_USTORY_SZPF',            &
           units = 'm-2',                                                       &
           long='density of understory plants by pft/size in number of plants per m2', &