Grass phenology update

biogeochem/EDCohortDynamicsMod.F90

@@ -128,9 +128,9 @@ module EDCohortDynamicsMod
 
   !-------------------------------------------------------------------------------------!
 
-  subroutine create_cohort(currentSite, patchptr, pft, nn, hite, dbh,   &
-                           prt, laimemory, status, recruitstatus,ctrim, &
-                           clayer, spread, bc_in)
+  subroutine create_cohort(currentSite, patchptr, pft, nn, hite, dbh, &
+                           prt, laimemory, sapwmemory, structmemory, &
+                           status, recruitstatus, ctrim, clayer, spread, bc_in)
 
     !
     ! !DESCRIPTION:
@@ -164,6 +164,10 @@ subroutine create_cohort(currentSite, patchptr, pft, nn, hite, dbh,   &
                                                   ! object
     real(r8), intent(in)      :: laimemory        ! target leaf biomass- set from 
                                                   ! previous year: kGC per indiv
+    real(r8), intent(in)   :: sapwmemory          ! target sapwood biomass- set from 
+                                                  ! previous year: kGC per indiv	
+    real(r8), intent(in)   :: structmemory        ! target structural biomass- set from 
+                                                  ! previous year: kGC per indiv							 
     real(r8), intent(in)      :: ctrim            ! What is the fraction of the maximum 
                                                   ! leaf biomass that we are targeting?
     real(r8), intent(in)      :: spread           ! The community assembly effects how 
@@ -212,7 +216,8 @@ subroutine create_cohort(currentSite, patchptr, pft, nn, hite, dbh,   &
     new_cohort%canopy_layer = clayer
     new_cohort%canopy_layer_yesterday = real(clayer, r8)
     new_cohort%laimemory    = laimemory
-
+    new_cohort%sapwmemory   = sapwmemory
+    new_cohort%structmemory = structmemory
 
     ! This sets things like vcmax25top, that depend on the
     ! leaf age fractions (which are defined by PARTEH)
@@ -476,6 +481,8 @@ subroutine nan_cohort(cc_p)
     currentCohort%dbh                = nan ! 'diameter at breast height' in cm                            
     currentCohort%hite               = nan ! height: meters                   
     currentCohort%laimemory          = nan ! target leaf biomass- set from previous year: kGC per indiv
+    currentCohort%sapwmemory         = nan ! target sapwood biomass- set from previous year: kGC per indiv
+    currentCohort%structmemory       = nan ! target structural biomass- set from previous year: kGC per indiv
     currentCohort%lai                = nan ! leaf area index of cohort   m2/m2      
     currentCohort%sai                = nan ! stem area index of cohort   m2/m2
     currentCohort%g_sb_laweight      = nan ! Total leaf conductance of cohort (stomata+blayer) weighted by leaf-area [m/s]*[m2]
@@ -1046,6 +1053,12 @@ subroutine fuse_cohorts(currentSite, currentPatch, bc_in)
 
                                 currentCohort%laimemory   = (currentCohort%n*currentCohort%laimemory   &
                                       + nextc%n*nextc%laimemory)/newn
+
+                                currentCohort%sapwmemory   = (currentCohort%n*currentCohort%sapwmemory   &
+                                      + nextc%n*nextc%sapwmemory)/newn
+
+                                currentCohort%structmemory   = (currentCohort%n*currentCohort%structmemory   &
+                                      + nextc%n*nextc%structmemory)/newn				      				      
 
                                 currentCohort%canopy_trim = (currentCohort%n*currentCohort%canopy_trim &
                                       + nextc%n*nextc%canopy_trim)/newn
@@ -1572,6 +1585,8 @@ subroutine copy_cohort( currentCohort,copyc )
     n%dbh             = o%dbh                                        
     n%hite            = o%hite
     n%laimemory       = o%laimemory
+    n%sapwmemory      = o%sapwmemory
+    n%structmemory    = o%structmemory
     n%lai             = o%lai                         
     n%sai             = o%sai  
     n%g_sb_laweight   = o%g_sb_laweight

biogeochem/EDMortalityFunctionsMod.F90

@@ -59,6 +59,7 @@ subroutine mortality_rates( cohort_in,bc_in,cmort,hmort,bmort,frmort )
     real(r8),intent(out) :: hmort  ! hydraulic failure mortality
     real(r8),intent(out) :: frmort ! freezing stress mortality
 
+
     real(r8) :: frac  ! relativised stored carbohydrate
     real(r8) :: leaf_c_target      ! target leaf biomass kgC
     real(r8) :: store_c
@@ -125,8 +126,7 @@ subroutine mortality_rates( cohort_in,bc_in,cmort,hmort,bmort,frmort )
        write(fates_log(),*) 'dbh problem in mortality_rates', &
             cohort_in%dbh,cohort_in%pft,cohort_in%n,cohort_in%canopy_layer
     endif
-
-
+    !-------------------------------------------------------------------------------- 
     !    Mortality due to cold and freezing stress (frmort), based on ED2 and:           
     !      Albani, M.; D. Medvigy; G. C. Hurtt; P. R. Moorcroft, 2006: The contributions 
     !           of land-use change, CO2 fertilization, and climate variability to the