Crown damage module

biogeochem/EDCanopyStructureMod.F90

@@ -381,7 +381,8 @@ subroutine DemoteFromLayer(currentSite,currentPatch,i_lyr,bc_in)
        currentCohort => currentPatch%shortest
        do while (associated(currentCohort))
           call carea_allom(currentCohort%dbh,currentCohort%n, &
-               currentSite%spread,currentCohort%pft,currentCohort%c_area)
+               currentSite%spread,currentCohort%pft,&
+               currentCohort%crowndamage, currentCohort%c_area)
 
           if(debug) then
              if(currentCohort%c_area<0._r8)then
@@ -686,9 +687,10 @@ subroutine DemoteFromLayer(currentSite,currentPatch,i_lyr,bc_in)
                 currentSite%demotion_carbonflux = currentSite%demotion_carbonflux + &
                      (leaf_c + store_c + fnrt_c + sapw_c + struct_c) * currentCohort%n
 
-                call carea_allom(copyc%dbh,copyc%n,currentSite%spread,copyc%pft,copyc%c_area)
+                call carea_allom(copyc%dbh,copyc%n,currentSite%spread,copyc%pft, &
+                     copyc%crowndamage, copyc%c_area)
                 call carea_allom(currentCohort%dbh,currentCohort%n,currentSite%spread, &
-                     currentCohort%pft,currentCohort%c_area)
+                     currentCohort%pft,currentCohort%crowndamage, currentCohort%c_area)
 
                 !----------- Insert copy into linked list ------------------------!
                 copyc%shorter => currentCohort
@@ -727,7 +729,8 @@ subroutine DemoteFromLayer(currentSite,currentPatch,i_lyr,bc_in)
              end if
 
              call carea_allom(currentCohort%dbh,currentCohort%n, &
-                  currentSite%spread,currentCohort%pft,currentCohort%c_area)
+                  currentSite%spread,currentCohort%pft,&
+                  currentCohort%crowndamage, currentCohort%c_area)
 
           endif !canopy layer = i_ly
 
@@ -833,7 +836,7 @@ subroutine PromoteIntoLayer(currentSite,currentPatch,i_lyr)
 
                 currentCohort%canopy_layer = i_lyr
                 call carea_allom(currentCohort%dbh,currentCohort%n,currentSite%spread, &
-                     currentCohort%pft,currentCohort%c_area)
+                     currentCohort%pft,currentCohort%crowndamage, currentCohort%c_area)
                 ! keep track of number and biomass of promoted cohort
                 currentSite%promotion_rate(currentCohort%size_class) = &
                      currentSite%promotion_rate(currentCohort%size_class) + currentCohort%n
@@ -859,7 +862,7 @@ subroutine PromoteIntoLayer(currentSite,currentPatch,i_lyr)
           currentCohort => currentPatch%tallest
           do while (associated(currentCohort))
              call carea_allom(currentCohort%dbh,currentCohort%n,currentSite%spread, &
-                  currentCohort%pft,currentCohort%c_area)
+                  currentCohort%pft,currentCohort%crowndamage, currentCohort%c_area)
              if(currentCohort%canopy_layer == i_lyr+1)then !look at the cohorts in the canopy layer below...
 
                 if (ED_val_comp_excln .ge. 0.0_r8 ) then
@@ -1126,7 +1129,7 @@ subroutine PromoteIntoLayer(currentSite,currentPatch,i_lyr)
                    newarea = currentCohort%c_area - cc_gain !new area of existing cohort
 
                    call carea_allom(currentCohort%dbh,currentCohort%n,currentSite%spread, &
-                        currentCohort%pft,currentCohort%c_area)
+                        currentCohort%pft,currentCohort%crowndamage, currentCohort%c_area)
 
                    ! number of individuals in promoted cohort.
                    copyc%n = currentCohort%n*cc_gain/currentCohort%c_area
@@ -1145,8 +1148,9 @@ subroutine PromoteIntoLayer(currentSite,currentPatch,i_lyr)
                         (leaf_c + fnrt_c + store_c + sapw_c + struct_c) * copyc%n
 
                    call carea_allom(currentCohort%dbh,currentCohort%n,currentSite%spread, &
-                        currentCohort%pft,currentCohort%c_area)
-                   call carea_allom(copyc%dbh,copyc%n,currentSite%spread,copyc%pft,copyc%c_area)
+                        currentCohort%pft,currentCohort%crowndamage, currentCohort%c_area)
+                   call carea_allom(copyc%dbh,copyc%n,currentSite%spread,copyc%pft,&
+                        copyc%crowndamage, copyc%c_area)
 
                    !----------- Insert copy into linked list ------------------------!
                    copyc%shorter => currentCohort
@@ -1223,7 +1227,9 @@ subroutine canopy_spread( currentSite )
        currentCohort => currentPatch%tallest
        do while (associated(currentCohort))
           call carea_allom(currentCohort%dbh,currentCohort%n, &
-               currentSite%spread,currentCohort%pft,currentCohort%c_area)
+               currentSite%spread,currentCohort%pft, &
+               currentCohort%crowndamage, currentCohort%c_area)
+
           if( ( int(prt_params%woody(currentCohort%pft)) .eq. itrue ) .and. &
                (currentCohort%canopy_layer .eq. 1 ) ) then
              sitelevel_canopyarea = sitelevel_canopyarea + currentCohort%c_area
@@ -1323,7 +1329,7 @@ subroutine canopy_summarization( nsites, sites, bc_in )
              ! Update the cohort's index within the size bin classes
              ! Update the cohort's index within the SCPF classification system
              call sizetype_class_index(currentCohort%dbh,currentCohort%pft, &
-                  currentCohort%size_class,currentCohort%size_by_pft_class)
+                                       currentCohort%size_class,currentCohort%size_by_pft_class)
 
              if (hlm_use_cohort_age_tracking .eq. itrue) then
                 call coagetype_class_index(currentCohort%coage,currentCohort%pft, &
@@ -1332,7 +1338,7 @@ subroutine canopy_summarization( nsites, sites, bc_in )
 
              if(hlm_use_sp.eq.ifalse)then
                 call carea_allom(currentCohort%dbh,currentCohort%n,sites(s)%spread,&
-                     currentCohort%pft,currentCohort%c_area)
+                     currentCohort%pft,currentCohort%crowndamage,currentCohort%c_area)
              endif
              currentCohort%treelai = tree_lai(leaf_c,             &
                   currentCohort%pft, currentCohort%c_area, currentCohort%n, &
@@ -1445,7 +1451,7 @@ subroutine leaf_area_profile( currentSite )
     ! currentCohort%treesai    ! SAI per unit crown area  (m2/m2)
     ! currentCohort%lai        ! LAI per unit canopy area (m2/m2)
     ! currentCohort%sai        ! SAI per unit canopy area (m2/m2)
-    ! currentCohort%NV         ! The number of discrete vegetation
+    ! currentCohort%nv         ! The number of discrete vegetation
     !                          ! layers needed to describe this crown
     !
     ! The following patch level diagnostics are updated here:
@@ -1471,7 +1477,6 @@ subroutine leaf_area_profile( currentSite )
     ! !USES:
 
     use EDtypesMod           , only : area, dinc_ed, hitemax, n_hite_bins
-
     !
     ! !ARGUMENTS
     type(ed_site_type)     , intent(inout) :: currentSite
@@ -1498,8 +1503,9 @@ subroutine leaf_area_profile( currentSite )
     real(r8) :: min_chite                ! bottom of cohort canopy  (m)
     real(r8) :: max_chite                ! top of cohort canopy      (m)
     real(r8) :: lai                      ! summed lai for checking m2 m-2
-    real(r8) :: leaf_c                   ! leaf carbon [kg]
-
+    real(r8) :: snow_depth_avg           ! avg snow over whole site
+    real(r8) :: leaf_c                   ! leaf carbon [kgC]
+    real(r8) :: target_c_area            ! crown area of undamaged cohort given dbh
     !----------------------------------------------------------------------
 
 
@@ -1554,13 +1560,14 @@ subroutine leaf_area_profile( currentSite )
              currentCohort%treelai = tree_lai(leaf_c, currentCohort%pft, currentCohort%c_area, &
                   currentCohort%n, currentCohort%canopy_layer,               &
                   currentPatch%canopy_layer_tlai,currentCohort%vcmax25top )
-
-            if (hlm_use_sp .eq. ifalse) then
-             currentCohort%treesai = tree_sai(currentCohort%pft, currentCohort%dbh, currentCohort%canopy_trim, &
-                  currentCohort%c_area, currentCohort%n, currentCohort%canopy_layer, &
-                  currentPatch%canopy_layer_tlai, currentCohort%treelai , &
-                  currentCohort%vcmax25top,4)
-            end if
+
+             if (hlm_use_sp .eq. ifalse) then
+                currentCohort%treesai = tree_sai(currentCohort%pft, currentCohort%dbh,&
+                     currentSite%spread, currentCohort%canopy_trim, &
+                     currentCohort%c_area, currentCohort%n, currentCohort%canopy_layer, &
+                     currentPatch%canopy_layer_tlai, currentCohort%treelai , &
+                     currentCohort%vcmax25top,4)
+             end if
 
              currentCohort%lai =  currentCohort%treelai *currentCohort%c_area/currentPatch%total_canopy_area
              currentCohort%sai =  currentCohort%treesai *currentCohort%c_area/currentPatch%total_canopy_area
@@ -2164,7 +2171,7 @@ subroutine CanopyLayerArea(currentPatch,site_spread,layer_index,layer_area)
     currentCohort => currentPatch%tallest
     do while (associated(currentCohort))
        call carea_allom(currentCohort%dbh,currentCohort%n,site_spread, &
-            currentCohort%pft,currentCohort%c_area)
+            currentCohort%pft,currentCohort%crowndamage, currentCohort%c_area)
        if (currentCohort%canopy_layer .eq. layer_index) then
           layer_area = layer_area + currentCohort%c_area
        end if
@@ -2208,7 +2215,8 @@ function NumPotentialCanopyLayers(currentPatch,site_spread,include_substory) res
        currentCohort => currentPatch%tallest
        do while (associated(currentCohort))
           if(currentCohort%canopy_layer == z) then
-             call carea_allom(currentCohort%dbh,currentCohort%n,site_spread,currentCohort%pft,c_area)
+             call carea_allom(currentCohort%dbh,currentCohort%n,site_spread,&
+                  currentCohort%pft,currentCohort%crowndamage, c_area)
              arealayer = arealayer + c_area
           end if
           currentCohort => currentCohort%shorter