Naming conventions imprved and README description of leafage_opt added

README.md

@@ -218,6 +218,7 @@ You can also [generate global inputs using Python (see python/global_data_proces
 | `crop_set`      | user-set real value of constant crop vegetation type heights (m) (only used if `crop_opt=1`) |
 | `co2_opt`       | user-set options for applying a CO2 inhibition factor for biogenic isoprene-only emissions using either the [Possell & Hewitt (2011)](https://doi.org/10.1111/j.1365-2486.2010.02306.x) (= `0`, default) or [Wilkinson et al. (2009)](https://doi.org/10.1111/j.1365-2486.2008.01803.x) method (= `1`). Use of option = `1` (Possell & Hewitt 2011) is especially recommended for sub-ambient CO2 concentrations.  To turn off co2 inhibition factor set `co2_opt=2`  |
 | `co2_set`       | user-set real value of atmospheric co2 concentration (ppmv) (only used if `co2_opt=0` or `co2_opt=1`) |
+| `leafage_opt`   | user-set options for applying leaf-age response to  biogenic  VOC emissions based on Guenther et al. 2006 https://doi.org/10.5194/acp-6-3181-2006, 2006. (= `0`, default), other wise assigns GAMMA_LEAGE=1 |
 | `lai_thresh`    | user-set real value of LAI threshold for contiguous canopy (m2/m2)                 |
 | `frt_thresh`    | user-set real value of forest fraction threshold for contiguous canopy             |
 | `fch_thresh`    | user-set real value of canopy height  threshold for contiguous canopy (m)          |

src/canopy_bioemi_mod.F90

@@ -8,7 +8,7 @@ module canopy_bioemi_mod
     SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
         PPFD_SHADE, TLEAF_SUN, TLEAF_SHADE, TLEAF_AVE, TEMP2, LU_OPT, &
         VTYPE, MODRES, CCE, VERT, CO2OPT, CO2SET, &
-        LEAFAGEOPT, PLAI, CLAI, NDAYS, TABOVECANOPY, &
+        LEAFAGEOPT, PASTLAI, CURRENTLAI, TSTEPLAI, TABOVECANOPY, &
         ANEW, AGRO, AMAT, AOLD, &
         EMI_IND, EMI_OUT)
 
@@ -67,9 +67,9 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
         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
+        REAL(RK),    INTENT( IN )       :: PASTLAI           ! Past LAI [cm2/cm2]
+        REAL(RK),    INTENT( IN )       :: CURRENTLAI           ! Current LAI [cm2/cm2]
+        REAL(RK),    INTENT( IN )       :: TSTEPLAI           !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)
@@ -171,7 +171,7 @@ SUBROUTINE CANOPY_BIO( ZK, FCLAI, FCH, LAI, FSUN, PPFD_SUN, &
 ! Get LEAF AGE factor
         TABOVECANOPY  = TEMP2   !TEMP2 (above air temp) for TABOVECANOPY
         !do i=1, SIZE(ZK)
-        GAMMALEAFAGE = GET_GAMMA_LEAFAGE(LEAFAGEOPT, PLAI, CLAI, NDAYS, TABOVECANOPY, ANEW, AGRO, AMAT, AOLD)
+        GAMMALEAFAGE = GET_GAMMA_LEAFAGE(LEAFAGEOPT, PASTLAI, CURRENTLAI, TSTEPLAI, TABOVECANOPY, ANEW, AGRO, AMAT, AOLD)
         !end do
 
 ! Calculate emissions profile in the canopy