Leaf Age Response function for all Biogenic VOCs consistent with MEGAN

input/namelist.canopy

@@ -37,6 +37,7 @@
   crop_set    =  3.0
   co2_opt     =  0
   co2_set     =  400.0
+  leafage_opt =  1
   lai_thresh  =  0.1
   frt_thresh  =  0.1
   fch_thresh  =  0.5

src/canopy_bioemi_mod.F90

@@ -8,6 +8,8 @@ module canopy_bioemi_mod
     SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
         PPFD_SHADE, TLEAF_SUN, TLEAF_SHADE, TLEAF_AVE, LU_OPT, &
         VTYPE, MODRES, CCE, VERT, CO2OPT, CO2SET, &
+        LEAFAGEOPT, PLAI, CLAI, NDAYS, TABOVECANOPY, &
+        ANEW, AGRO, AMAT, AOLD, &
         EMI_IND, EMI_OUT)
 
 !-----------------------------------------------------------------------
@@ -33,10 +35,12 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
 !     Feb 2023 P.C. Campbell: Modified for multiple biogenic species
 !     Jul 2023 P.C. Campbell: Restructured to use FSUN, TLEAF, and PPFD
 !                             as inputs
+!     Sept 2023: QZ Rasool : Modifications for  LeafAge Response for multiple
+!     BVOCs
 !-----------------------------------------------------------------------
 !-----------------------------------------------------------------------
         use canopy_const_mod, ONLY: rk,rgasuniv   !constants for canopy models
-        use canopy_utils_mod,  ONLY: interp_linear1_internal, GET_GAMMA_CO2
+        use canopy_utils_mod,  ONLY: interp_linear1_internal,GET_GAMMA_CO2,GET_GAMMA_LEAFAGE
         use canopy_bioparm_mod
 
 ! Arguments:
@@ -59,6 +63,12 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
         INTEGER,     INTENT( IN )       :: VERT            ! MEGAN vertical integration option (default = 0/no integration)
         INTEGER,     INTENT( IN )       :: CO2OPT          ! Option for co2 inhibition calculation
         REAL(RK),    INTENT( IN )       :: CO2SET          ! User set atmospheric CO2 conc [ppmv]
+
+        INTEGER,    INTENT( IN )       :: LEAFAGEOPT      ! leafage_opt (0= ON, 1= off i.e. GAMMALEAFAGE =1, in canopy_readnml.F90)
+        REAL(RK),    INTENT( IN )       :: PLAI           ! Past LAI [cm2/cm2]
+        REAL(RK),    INTENT( IN )       :: CLAI           ! Current LAI [cm2/cm2]
+        REAL(RK),    INTENT( IN )       :: NDAYS           !Number of days between the past and current LAI
+
         INTEGER,     INTENT( IN )       :: EMI_IND         ! Input biogenic emissions index
         REAL(RK),    INTENT( OUT )      :: EMI_OUT(:)      ! Output canopy layer volume emissions (kg m-3 s-1)
 
@@ -91,7 +101,19 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
         REAL(RK) :: CT1                            ! Activation energy (kJ/mol)
         REAL(RK) :: CEO                            ! Empirical coefficient
         REAL(RK) :: EF                             ! Final Mapped Emission factor (EF) (ug/m2 hr)
+
+        REAL(RK) :: TABOVECANOPY(SIZE(ZK))  ! Above Canopy Temp (assigned = TLEAF_AVE for now)
+        ! Empirical coeff.'s for Leaf Age factor calculations (see
+        ! canopy_bioparm_mod or call canopy_biop)
+        REAL(RK) :: ANEW
+        REAL(RK) :: AGRO
+        REAL(RK) :: AMAT
+        REAL(RK) :: AOLD
+
         REAL(RK) :: GAMMACO2                       ! CO2 inhibition factor (isoprene only)
+
+        REAL(RK) :: GAMMALEAFAGE(SIZE(ZK))                 ! LEAF AGE factor
+
         integer i, LAYERS
 
 ! Constant Canopy Parameters
@@ -105,7 +127,7 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
         TLEAF_OPT = 313.0_rk + (0.6_rk * (TLEAF240_AVE-297.0_rk)) !Guenther et al. (2012)
 
 ! Calculate emission species/plant-dependent mapped emission factors
-        call canopy_biop(EMI_IND, LU_OPT, VTYPE, EF, CT1, CEO)
+        call canopy_biop(EMI_IND, LU_OPT, VTYPE, EF, CT1, CEO, ANEW, AGRO, AMAT, AOLD)  !Update: Added leaf age empirical coeff's
 
         E_OPT = CEO * EXP(0.05_rk * (TLEAF24_AVE-297.0_rk)) * EXP(0.05_rk * (TLEAF240_AVE-297.0_rk))
 
@@ -143,6 +165,12 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
             GAMMACO2 = 1.0_rk
         end if
 
+! Get LEAF AGE factor
+        TABOVECANOPY  = TLEAF_AVE  !Assume instantaneous TLEAF estimate for TABOVECANOPY
+        do i=1, SIZE(ZK)
+            GAMMALEAFAGE(i) = GET_GAMMA_LEAFAGE(LEAFAGEOPT, PLAI, CLAI, NDAYS, TABOVECANOPY(i), ANEW, AGRO, AMAT, AOLD)
+        end do
+
 ! Calculate emissions profile in the canopy
         EMI_OUT = 0.0_rk  ! set initial emissions profile to zero
         FLAI = 0.0_rk  ! set initial fractional FLAI (LAD) profile to zero
@@ -151,7 +179,7 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
             do i=1, SIZE(ZK)
                 if (ZK(i) .gt. 0.0 .and. ZK(i) .le. FCH) then           ! above ground level and at/below canopy top
                     FLAI(i) = ((FCLAI(i+1) - FCLAI(i)) * LAI)/MODRES    !fractional LAI in each layer converted to LAD (m2 m-3)
-                    EMI_OUT(i) = FLAI(i) * EF * GammaTLEAF_AVE(i) * GammaPPFD_AVE(i) * GAMMACO2 * CCE  ! (ug m-3 hr-1)
+                    EMI_OUT(i) = FLAI(i) * EF * GammaTLEAF_AVE(i) * GammaPPFD_AVE(i) * GAMMACO2 * CCE * GAMMALEAFAGE(i)  ! (ug m-3 hr-1)
                     EMI_OUT(i) = EMI_OUT(i) * 2.7777777777778E-13_rk    !convert emissions output to (kg m-3 s-1)
                 end if
             end do
@@ -169,7 +197,7 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
                 end if
             end do
             EMI_OUT(SIZE(ZK)) = LAI * EF * SUM(GammaTLEAF_AVE(1:LAYERS) * GammaPPFD_AVE(1:LAYERS) * &
-                VPGWT(1:LAYERS)) * GAMMACO2 * CCE   !put into top model layer (ug m-2 hr-1)
+                VPGWT(1:LAYERS)* GAMMALEAFAGE(1:LAYERS)) * GAMMACO2 * CCE   !put into top model layer (ug m-2 hr-1)
             EMI_OUT = EMI_OUT * 2.7777777777778E-13_rk    !convert emissions output to    (kg m-2 s-1)
         else if (VERT .eq. 2) then       !"MEGANv3-like": Add weighted sum of activity coefficients using normal distribution
             !across canopy layers using 5 layer numbers directly from MEGANv3
@@ -201,7 +229,7 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
                 VPGWT(i) = GAUSS(i)/sum(GAUSS(1:LAYERS))
             end do
             EMI_OUT(SIZE(ZK)) = LAI * EF * SUM(GammaTLEAF_AVE(1:LAYERS) * GammaPPFD_AVE(1:LAYERS) * &
-                VPGWT(1:LAYERS)) * GAMMACO2 * CCE   !put into top model layer (ug m-2 hr-1)
+                VPGWT(1:LAYERS)* GAMMALEAFAGE(1:LAYERS)) * GAMMACO2 * CCE    !put into top model layer (ug m-2 hr-1)
             EMI_OUT = EMI_OUT * 2.7777777777778E-13_rk    !convert emissions output to    (kg m-2 s-1)
         else if (VERT .eq. 3) then       !"MEGANv3-like": Add weighted sum of activity coefficients equally
             !across canopy layers
@@ -212,7 +240,7 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
                 VPGWT(i) = 1.0_rk/LAYERS
             end do
             EMI_OUT(SIZE(ZK)) = LAI * EF * SUM(GammaTLEAF_AVE(1:LAYERS) * GammaPPFD_AVE(1:LAYERS) * &
-                VPGWT(1:LAYERS)) * GAMMACO2 * CCE   !put into top model layer (ug m-2 hr-1)
+                VPGWT(1:LAYERS)* GAMMALEAFAGE(1:LAYERS)) * GAMMACO2 * CCE    !put into top model layer (ug m-2 hr-1)
             EMI_OUT = EMI_OUT * 2.7777777777778E-13_rk    !convert emissions output to    (kg m-2 s-1)
         else
             write(*,*)  'Wrong BIOVERT_OPT choice of ', VERT, ' in namelist...exiting'