Skip to content

Group bioparm_mod

Modules > bioparm_mod

Module for biogenic emission factors and parameters.

Modules

Type Name
module Input Variables
Get biogenic emission factors and parameters from MEGAN2.1.

Public Attributes

Type Name
real(rk), dimension((ug/m2 hr) real(rk), intent(out) ::ldf !> light-dependent fraction real(rk), intent(out) ::beta !> empirical coefficient for temperature dependence of light-independent fraction real(rk), intent(out) ::ct1 !> out activation energy(kj/mol) real(rk), intent(out) ::ceo !> out empirical coefficient real(rk), intent(out) ::anew, agro, amat, aold !> empirical factors or coefficients for:growing, mature, and old/senescing foliage, as per table 4 of guenther et al., 2012 real(rk), intent(out) ::roota, rootb !> coefficients a and b used for pft dependent cumulative root depth fraction[m-1] real(rk), intent(out) ::caq !> coefficient for poor air quality stress real(rk), intent(out) ::taq !> threshold for poor air quality stress(ppm-hours) real(rk), intent(out) ::dtaq !> delta threshold for poor air quality stress(ppm-hours) real(rk), intent(out) ::cht !> coefficient for high temperature stress real(rk), intent(out) ::tht !> threshold for high temperature stress(k) real(rk), intent(out) ::dtht !> delta threshold high temperature stress(k) real(rk), intent(out) ::clt !> coefficient for low temperature stress real(rk), intent(out) ::tlt !> threshold for low temperature stress(k) real(rk), intent(out) ::dtlt !> delta threshold low temperature stress(k) real(rk), intent(out) ::chw !> coefficient for high wind stress real(rk), intent(out) ::thw !> threshold for high wind stress(m/s) real(rk), intent(out) ::dthw !> delta threshold high wind stress(m/s) !> } !> \defgroup bioparm_local_vars local variables !! \brief local variables for parameter assignment !! { real(rk) ::ef1, ef2, ef3, ef4, ef5, ef6, ef7 !> plant emission factors(ef)(ug/m2 hr) real(rk) ::ef8, ef9, ef10, ef11, ef12, ef13 !> plant emission factors(ef)(ug/m2 hr) real(rk) ::ef14, ef15 !> plant emission factors(ef)(ug/m2 hr) !> } !> \defgroup bioparm_isop_params isoprene parameters !! \brief plant-dependent emission capacity factors for isoprene from tables 2-3 of guenther et al.(2012) !! { !> \brief needleleaf evergreen temperate tree isoprene ef(μg/m²/hr) real(rk), parameter ::ef1_isop= 600.0_rk !> \brief needleleaf evergreen boreal tree isoprene ef(μg/m²/hr) real(rk), parameter ::ef2_isop= 3000.0_rk !> \brief needleleaf deciduous boreal tree isoprene ef(μg/m²/hr) real(rk), parameter ::ef3_isop= 1.0_rk !> \brief broadleaf evergreen tropical tree isoprene ef(μg/m²/hr) real(rk), parameter ::ef4_isop= 7000.0_rk !> \brief broadleaf evergreen temperate tree isoprene ef(μg/m²/hr) real(rk), parameter ::ef5_isop= 10000.0_rk !> \brief broadleaf deciduous tropical tree isoprene ef(μg/m²/hr) real(rk), parameter ::ef6_isop= 7000.0_rk !> \brief broadleaf deciduous temperate tree isoprene ef(μg/m²/hr) real(rk), parameter ::ef7_isop= 10000.0_rk !> \brief broadleaf deciduous boreal tree isoprene ef(μg/m²/hr) real(rk), parameter ::ef8_isop= 11000.0_rk !> \brief broadleaf evergreen temperate shrub isoprene ef(μg/m²/hr) real(rk), parameter ::ef9_isop= 2000.0_rk !> \brief broadleaf deciduous temperate shrub isoprene ef(μg/m²/hr) real(rk), parameter ::ef10_isop= 4000.0_rk !> \brief broadleaf deciduous boreal shrub isoprene ef(μg/m²/hr) real(rk), parameter ::ef11_isop= 4000.0_rk !> \brief arctic c3 grass isoprene ef(μg/m²/hr) real(rk), parameter ::ef12_isop= 1600.0_rk !> \brief cool c3 grass isoprene ef(μg/m²/hr) real(rk), parameter ::ef13_isop= 800.0_rk !> \brief warm c4 grass isoprene ef(μg/m²/hr) real(rk), parameter ::ef14_isop= 200.0_rk !> \brief crop1 isoprene ef(μg/m²/hr) real(rk), parameter ::ef15_isop= 1.0_rk !> \brief isoprene leaf age factor for new foliage(table 4 of guenther et al., 2012) real(rk), parameter ::anew_isop=0.05_rk !> \brief isoprene leaf age factor for growing foliage(table 4 of guenther et al., 2012) real(rk), parameter ::agro_isop=0.6_rk !> \brief isoprene leaf age factor for mature foliage(table 4 of guenther et al., 2012) real(rk), parameter ::amat_isop=1.0_rk !> \brief isoprene leaf age factor for old/senescing foliage(table 4 of guenther et al., 2012) real(rk), parameter ::aold_isop=0.9_rk !> } ! plant-dependent emissions capacity/factors(efs) for myrcene(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_myrc= 70.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_myrc= 70.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_myrc= 60.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_myrc= 80.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_myrc= 30.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_myrc= 80.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_myrc= 30.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_myrc= 30.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_myrc= 30.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_myrc= 50.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_myrc= 30.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_myrc= 0.3_rk ! arctic c3 grass real(rk), parameter ::ef13_myrc= 0.3_rk ! cool c3 grass real(rk), parameter ::ef14_myrc= 0.3_rk ! warm c4 grass real(rk), parameter ::ef15_myrc= 0.3_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for myrcene as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_myrc=2.0_rk real(rk), parameter ::agro_myrc=1.8_rk real(rk), parameter ::amat_myrc=1.0_rk real(rk), parameter ::aold_myrc=1.05_rk ! plant-dependent emissions capacity/factors(efs) for sabinene(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_sabi= 70.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_sabi= 70.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_sabi= 40.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_sabi= 80.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_sabi= 50.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_sabi= 80.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_sabi= 50.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_sabi= 50.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_sabi= 50.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_sabi= 70.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_sabi= 50.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_sabi= 0.7_rk ! arctic c3 grass real(rk), parameter ::ef13_sabi= 0.7_rk ! cool c3 grass real(rk), parameter ::ef14_sabi= 0.7_rk ! warm c4 grass real(rk), parameter ::ef15_sabi= 0.7_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for sabinene as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_sabi=2.0_rk real(rk), parameter ::agro_sabi=1.8_rk real(rk), parameter ::amat_sabi=1.0_rk real(rk), parameter ::aold_sabi=1.05_rk ! plant-dependent emissions capacity/factors(efs) for limonene(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_limo= 100.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_limo= 100.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_limo= 130.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_limo= 80.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_limo= 80.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_limo= 80.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_limo= 80.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_limo= 80.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_limo= 60.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_limo= 100.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_limo= 60.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_limo= 0.7_rk ! arctic c3 grass real(rk), parameter ::ef13_limo= 0.7_rk ! cool c3 grass real(rk), parameter ::ef14_limo= 0.7_rk ! warm c4 grass real(rk), parameter ::ef15_limo= 0.7_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for limonene as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_limo=2.0_rk real(rk), parameter ::agro_limo=1.8_rk real(rk), parameter ::amat_limo=1.0_rk real(rk), parameter ::aold_limo=1.05_rk ! plant-dependent emissions capacity/factors(efs) for 3-carene(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_care= 160.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_care= 160.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_care= 80.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_care= 40.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_care= 30.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_care= 40.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_care= 30.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_care= 30.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_care= 30.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_care= 100.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_care= 30.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_care= 0.3_rk ! arctic c3 grass real(rk), parameter ::ef13_care= 0.3_rk ! cool c3 grass real(rk), parameter ::ef14_care= 0.3_rk ! warm c4 grass real(rk), parameter ::ef15_care= 0.3_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for 3-carene as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_care=2.0_rk real(rk), parameter ::agro_care=1.8_rk real(rk), parameter ::amat_care=1.0_rk real(rk), parameter ::aold_care=1.05_rk ! plant-dependent emissions capacity/factors(efs) for t-beta-ocimene(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_ocim= 70.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_ocim= 70.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_ocim= 60.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_ocim= 150.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_ocim= 120.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_ocim= 150.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_ocim= 120.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_ocim= 120.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_ocim= 90.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_ocim= 150.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_ocim= 90.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_ocim= 2.0_rk ! arctic c3 grass real(rk), parameter ::ef13_ocim= 2.0_rk ! cool c3 grass real(rk), parameter ::ef14_ocim= 2.0_rk ! warm c4 grass real(rk), parameter ::ef15_ocim= 2.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for t-beta-ocimene as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_ocim=2.0_rk real(rk), parameter ::agro_ocim=1.8_rk real(rk), parameter ::amat_ocim=1.0_rk real(rk), parameter ::aold_ocim=1.05_rk ! plant-dependent emissions capacity/factors(efs) for beta-pinene(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_bpin= 300.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_bpin= 300.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_bpin= 200.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_bpin= 120.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_bpin= 130.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_bpin= 120.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_bpin= 130.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_bpin= 130.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_bpin= 100.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_bpin= 150.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_bpin= 100.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_bpin= 1.5_rk ! arctic c3 grass real(rk), parameter ::ef13_bpin= 1.5_rk ! cool c3 grass real(rk), parameter ::ef14_bpin= 1.5_rk ! warm c4 grass real(rk), parameter ::ef15_bpin= 1.5_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for beta-pinene as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_bpin=2.0_rk real(rk), parameter ::agro_bpin=1.8_rk real(rk), parameter ::amat_bpin=1.0_rk real(rk), parameter ::aold_bpin=1.05_rk ! plant-dependent emissions capacity/factors(efs) for alpha-pinene(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_apin= 500.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_apin= 500.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_apin= 510.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_apin= 600.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_apin= 400.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_apin= 600.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_apin= 400.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_apin= 400.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_apin= 200.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_apin= 300.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_apin= 200.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_apin= 2.0_rk ! arctic c3 grass real(rk), parameter ::ef13_apin= 2.0_rk ! cool c3 grass real(rk), parameter ::ef14_apin= 2.0_rk ! warm c4 grass real(rk), parameter ::ef15_apin= 2.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for alpha-pinene as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_apin=2.0_rk real(rk), parameter ::agro_apin=1.8_rk real(rk), parameter ::amat_apin=1.0_rk real(rk), parameter ::aold_apin=1.05_rk ! plant-dependent emissions capacity/factors(efs) for other monoterpenes(tables 2-3 of guenther et al., 2012)(ug/m2 hr) ! ! other monoterpenes category(34 compounds): see table 1 of guenther et al.(2012) real(rk), parameter ::ef1_mono= 180.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_mono= 180.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_mono= 170.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_mono= 150.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_mono= 150.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_mono= 150.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_mono= 150.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_mono= 150.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_mono= 110.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_mono= 200.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_mono= 110.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_mono= 5.0_rk ! arctic c3 grass real(rk), parameter ::ef13_mono= 5.0_rk ! cool c3 grass real(rk), parameter ::ef14_mono= 5.0_rk ! warm c4 grass real(rk), parameter ::ef15_mono= 5.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for other monoterpenes as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_mono=2.0_rk real(rk), parameter ::agro_mono=1.8_rk real(rk), parameter ::amat_mono=1.0_rk real(rk), parameter ::aold_mono=1.05_rk ! plant-dependent emissions capacity/factors(efs) for alpha-farnesene(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_farn= 40.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_farn= 40.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_farn= 40.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_farn= 60.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_farn= 40.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_farn= 60.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_farn= 40.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_farn= 40.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_farn= 40.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_farn= 40.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_farn= 40.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_farn= 3.0_rk ! arctic c3 grass real(rk), parameter ::ef13_farn= 3.0_rk ! cool c3 grass real(rk), parameter ::ef14_farn= 3.0_rk ! warm c4 grass real(rk), parameter ::ef15_farn= 4.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for alpha-farnesene as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_farn=0.4_rk real(rk), parameter ::agro_farn=0.6_rk real(rk), parameter ::amat_farn=1.0_rk real(rk), parameter ::aold_farn=0.95_rk ! plant-dependent emissions capacity/factors(efs) for beta-caryophyllene(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_cary= 80.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_cary= 80.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_cary= 80.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_cary= 60.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_cary= 40.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_cary= 60.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_cary= 40.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_cary= 40.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_cary= 50.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_cary= 50.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_cary= 50.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_cary= 1.0_rk ! arctic c3 grass real(rk), parameter ::ef13_cary= 1.0_rk ! cool c3 grass real(rk), parameter ::ef14_cary= 1.0_rk ! warm c4 grass real(rk), parameter ::ef15_cary= 4.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for beta-caryophyllene as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_cary=0.4_rk real(rk), parameter ::agro_cary=0.6_rk real(rk), parameter ::amat_cary=1.0_rk real(rk), parameter ::aold_cary=0.95_rk ! plant-dependent emissions capacity/factors(efs) for other sesquieterpenes(tables 2-3 of guenther et al., 2012)(ug/m2 hr) ! other sesquiterpenes category(30 compounds): see table 1 of guenther et al.(2012) real(rk), parameter ::ef1_sesq= 120.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_sesq= 120.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_sesq= 120.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_sesq= 120.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_sesq= 100.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_sesq= 120.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_sesq= 100.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_sesq= 100.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_sesq= 100.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_sesq= 100.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_sesq= 100.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_sesq= 1.0_rk ! arctic c3 grass real(rk), parameter ::ef13_sesq= 1.0_rk ! cool c3 grass real(rk), parameter ::ef14_sesq= 1.0_rk ! warm c4 grass real(rk), parameter ::ef15_sesq= 1.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for other sesquieterpenes as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_sesq=0.4_rk real(rk), parameter ::agro_sesq=0.6_rk real(rk), parameter ::amat_sesq=1.0_rk real(rk), parameter ::aold_sesq=0.95_rk ! plant-dependent emissions capacity/factors(efs) for 232-mbo(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_mbol= 700.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_mbol= 60.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_mbol= 0.01_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_mbol= 0.01_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_mbol= 0.01_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_mbol= 0.01_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_mbol= 0.01_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_mbol= 2.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_mbol= 0.01_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_mbol= 0.01_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_mbol= 0.01_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_mbol= 0.01_rk ! arctic c3 grass real(rk), parameter ::ef13_mbol= 0.01_rk ! cool c3 grass real(rk), parameter ::ef14_mbol= 0.01_rk ! warm c4 grass real(rk), parameter ::ef15_mbol= 0.01_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for 232-mbo as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_mbol=0.05_rk real(rk), parameter ::agro_mbol=0.6_rk real(rk), parameter ::amat_mbol=1.0_rk real(rk), parameter ::aold_mbol=0.9_rk ! plant-dependent emissions capacity/factors(efs) for methanol(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_meth= 900.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_meth= 900.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_meth= 900.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_meth= 500.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_meth= 900.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_meth= 500.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_meth= 900.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_meth= 900.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_meth= 900.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_meth= 900.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_meth= 900.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_meth= 500.0_rk ! arctic c3 grass real(rk), parameter ::ef13_meth= 500.0_rk ! cool c3 grass real(rk), parameter ::ef14_meth= 500.0_rk ! warm c4 grass real(rk), parameter ::ef15_meth= 900.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for methanol as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_meth=3.5_rk real(rk), parameter ::agro_meth=3.0_rk real(rk), parameter ::amat_meth=1.0_rk real(rk), parameter ::aold_meth=1.2_rk ! plant-dependent emissions capacity/factors(efs) for acetone(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_acet= 240.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_acet= 240.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_acet= 240.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_acet= 240.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_acet= 240.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_acet= 240.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_acet= 240.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_acet= 240.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_acet= 240.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_acet= 240.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_acet= 240.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_acet= 80.0_rk ! arctic c3 grass real(rk), parameter ::ef13_acet= 80.0_rk ! cool c3 grass real(rk), parameter ::ef14_acet= 80.0_rk ! warm c4 grass real(rk), parameter ::ef15_acet= 80.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for acetone as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_acet=1.0_rk real(rk), parameter ::agro_acet=1.0_rk real(rk), parameter ::amat_acet=1.0_rk real(rk), parameter ::aold_acet=1.0_rk ! plant-dependent emissions capacity/factors(efs) for carbon monoxide(tables 2-3 of guenther et al., 2012)(ug/m2 hr) real(rk), parameter ::ef1_co= 600.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_co= 600.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_co= 600.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_co= 600.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_co= 600.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_co= 600.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_co= 600.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_co= 600.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_co= 600.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_co= 600.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_co= 600.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_co= 600.0_rk ! arctic c3 grass real(rk), parameter ::ef13_co= 600.0_rk ! cool c3 grass real(rk), parameter ::ef14_co= 600.0_rk ! warm c4 grass real(rk), parameter ::ef15_co= 600.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for carbon monoxide as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_co=1.0_rk real(rk), parameter ::agro_co=1.0_rk real(rk), parameter ::amat_co=1.0_rk real(rk), parameter ::aold_co=1.0_rk ! plant-dependent emissions capacity/factors(efs) for bidi voc species(tables 2-3 of guenther et al., 2012)(ug/m2 hr) ! bidirectional voc(5 compounds):see table 1 of guenther et al.(2012) real(rk), parameter ::ef1_bvoc= 500.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_bvoc= 500.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_bvoc= 500.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_bvoc= 500.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_bvoc= 500.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_bvoc= 500.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_bvoc= 500.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_bvoc= 500.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_bvoc= 500.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_bvoc= 500.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_bvoc= 500.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_bvoc= 80.0_rk ! arctic c3 grass real(rk), parameter ::ef13_bvoc= 80.0_rk ! cool c3 grass real(rk), parameter ::ef14_bvoc= 80.0_rk ! warm c4 grass real(rk), parameter ::ef15_bvoc= 80.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for bidi voc species as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_bvoc=1.0_rk real(rk), parameter ::agro_bvoc=1.0_rk real(rk), parameter ::amat_bvoc=1.0_rk real(rk), parameter ::aold_bvoc=1.0_rk ! plant-dependent emissions capacity/factors(efs) for stress vocs(tables 2-3 of guenther et al., 2012)(ug/m2 hr) ! stress voc(15 compounds):see table 1 of guenther et al.(2012) real(rk), parameter ::ef1_svoc= 300.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_svoc= 300.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_svoc= 300.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_svoc= 300.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_svoc= 300.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_svoc= 300.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_svoc= 300.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_svoc= 300.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_svoc= 300.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_svoc= 300.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_svoc= 300.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_svoc= 300.0_rk ! arctic c3 grass real(rk), parameter ::ef13_svoc= 300.0_rk ! cool c3 grass real(rk), parameter ::ef14_svoc= 300.0_rk ! warm c4 grass real(rk), parameter ::ef15_svoc= 300.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for stress vocs as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_svoc=1.0_rk real(rk), parameter ::agro_svoc=1.0_rk real(rk), parameter ::amat_svoc=1.0_rk real(rk), parameter ::aold_svoc=1.0_rk ! plant-dependent emissions capacity/factors(efs) for other vocs(tables 2-3 of guenther et al., 2012)(ug/m2 hr) ! other voc(49 compounds):see table 1 of guenther et al.(2012) real(rk), parameter ::ef1_ovoc= 140.0_rk ! needleleaf evergreen temperate tree real(rk), parameter ::ef2_ovoc= 140.0_rk ! needleleaf evergreen boreal tree real(rk), parameter ::ef3_ovoc= 140.0_rk ! needleleaf deciduous boreal tree real(rk), parameter ::ef4_ovoc= 140.0_rk ! broadleaf evergreen tropical tree real(rk), parameter ::ef5_ovoc= 140.0_rk ! broadleaf evergreen temperate tree real(rk), parameter ::ef6_ovoc= 140.0_rk ! broadleaf deciduous tropical tree real(rk), parameter ::ef7_ovoc= 140.0_rk ! broadleaf deciduous temperate tree real(rk), parameter ::ef8_ovoc= 140.0_rk ! broadleaf deciduous boreal tree real(rk), parameter ::ef9_ovoc= 140.0_rk ! broadleaf evergreen temperate shrub real(rk), parameter ::ef10_ovoc= 140.0_rk ! broadleaf deciduous temperate shrub real(rk), parameter ::ef11_ovoc= 140.0_rk ! broadleaf deciduous boreal shrub real(rk), parameter ::ef12_ovoc= 140.0_rk ! arctic c3 grass real(rk), parameter ::ef13_ovoc= 140.0_rk ! cool c3 grass real(rk), parameter ::ef14_ovoc= 140.0_rk ! warm c4 grass real(rk), parameter ::ef15_ovoc= 140.0_rk ! crop1 !empirical factors or coefficients for:growing, mature, and old/senescing foliage, for other vocs as per table 4 of guenther et al., 2012 real(rk), parameter ::anew_ovoc=1.0_rk real(rk), parameter ::agro_ovoc=1.0_rk real(rk), parameter ::amat_ovoc=1.0_rk real(rk), parameter ::aold_ovoc=1.0_rk ! species-dependent parameterized canopy model parameters(table 4 of guenther et al., 2012) real(rk), parameter ::ldf_isop= 1.0_rk !light-dependent fraction real(rk), parameter ::beta_isop= 0.13_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_isop= 95.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_isop= 2.0_rk !empirical coefficient real(rk), parameter ::caq_isop= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_isop= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_isop= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_isop= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_isop= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_isop= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_isop= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_isop= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_isop= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_isop= 1.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_isop= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_isop= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_myrc= 0.6_rk !light-dependent fraction real(rk), parameter ::beta_myrc= 0.1_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_myrc= 80.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_myrc= 1.83_rk !empirical coefficient real(rk), parameter ::caq_myrc= 5.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_myrc= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_myrc= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_myrc= 5.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_myrc= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_myrc= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_myrc= 5.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_myrc= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_myrc= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_myrc= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_myrc= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_myrc= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_sabi= 0.6_rk !light-dependent fraction real(rk), parameter ::beta_sabi= 0.1_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_sabi= 80.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_sabi= 1.83_rk !empirical coefficient real(rk), parameter ::caq_sabi= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_sabi= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_sabi= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_sabi= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_sabi= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_sabi= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_sabi= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_sabi= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_sabi= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_sabi= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_sabi= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_sabi= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_limo= 0.2_rk !light-dependent fraction real(rk), parameter ::beta_limo= 0.1_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_limo= 80.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_limo= 1.83_rk !empirical coefficient real(rk), parameter ::caq_limo= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_limo= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_limo= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_limo= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_limo= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_limo= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_limo= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_limo= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_limo= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_limo= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_limo= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_limo= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_care= 0.2_rk !light-dependent fraction real(rk), parameter ::beta_care= 0.1_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_care= 80.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_care= 1.83_rk !empirical coefficient real(rk), parameter ::caq_care= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_care= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_care= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_care= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_care= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_care= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_care= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_care= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_care= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_care= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_care= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_care= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_ocim= 0.8_rk !light-dependent fraction real(rk), parameter ::beta_ocim= 0.1_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_ocim= 80.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_ocim= 1.83_rk !empirical coefficient real(rk), parameter ::caq_ocim= 5.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_ocim= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_ocim= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_ocim= 5.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_ocim= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_ocim= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_ocim= 5.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_ocim= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_ocim= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_ocim= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_ocim= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_ocim= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_bpin= 0.2_rk !light-dependent fraction real(rk), parameter ::beta_bpin= 0.1_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_bpin= 80.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_bpin= 1.83_rk !empirical coefficient real(rk), parameter ::caq_bpin= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_bpin= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_bpin= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_bpin= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_bpin= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_bpin= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_bpin= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_bpin= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_bpin= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_bpin= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_bpin= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_bpin= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_apin= 0.6_rk !light-dependent fraction real(rk), parameter ::beta_apin= 0.1_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_apin= 80.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_apin= 1.83_rk !empirical coefficient real(rk), parameter ::caq_apin= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_apin= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_apin= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_apin= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_apin= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_apin= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_apin= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_apin= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_apin= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_apin= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_apin= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_apin= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_mono= 0.4_rk !light-dependent fraction real(rk), parameter ::beta_mono= 0.1_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_mono= 80.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_mono= 1.83_rk !empirical coefficient real(rk), parameter ::caq_mono= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_mono= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_mono= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_mono= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_mono= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_mono= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_mono= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_mono= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_mono= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_mono= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_mono= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_mono= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_farn= 0.5_rk !light-dependent fraction real(rk), parameter ::beta_farn= 0.17_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_farn= 130.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_farn= 2.37_rk !empirical coefficient real(rk), parameter ::caq_farn= 5.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_farn= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_farn= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_farn= 5.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_farn= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_farn= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_farn= 5.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_farn= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_farn= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_farn= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_farn= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_farn= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_cary= 0.5_rk !light-dependent fraction real(rk), parameter ::beta_cary= 0.17_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_cary= 130.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_cary= 2.37_rk !empirical coefficient real(rk), parameter ::caq_cary= 5.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_cary= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_cary= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_cary= 5.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_cary= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_cary= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_cary= 5.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_cary= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_cary= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_cary= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_cary= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_cary= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_sesq= 0.5_rk !light-dependent fraction real(rk), parameter ::beta_sesq= 0.17_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_sesq= 130.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_sesq= 2.37_rk !empirical coefficient real(rk), parameter ::caq_sesq= 5.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_sesq= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_sesq= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_sesq= 5.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_sesq= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_sesq= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_sesq= 5.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_sesq= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_sesq= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_sesq= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_sesq= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_sesq= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_mbol= 1.0_rk !light-dependent fraction real(rk), parameter ::beta_mbol= 0.13_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_mbol= 95.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_mbol= 2.0_rk !empirical coefficient real(rk), parameter ::caq_mbol= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_mbol= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_mbol= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_mbol= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_mbol= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_mbol= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_mbol= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_mbol= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_mbol= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_mbol= 1.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_mbol= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_mbol= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_meth= 0.8_rk !light-dependent fraction real(rk), parameter ::beta_meth= 0.08_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_meth= 60.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_meth= 1.6_rk !empirical coefficient real(rk), parameter ::caq_meth= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_meth= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_meth= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_meth= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_meth= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_meth= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_meth= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_meth= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_meth= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_meth= 1.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_meth= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_meth= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_acet= 0.2_rk !light-dependent fraction real(rk), parameter ::beta_acet= 0.1_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_acet= 80.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_acet= 1.83_rk !empirical coefficient real(rk), parameter ::caq_acet= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_acet= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_acet= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_acet= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_acet= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_acet= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_acet= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_acet= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_acet= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_acet= 1.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_acet= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_acet= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_co= 1.0_rk !light-dependent fraction real(rk), parameter ::beta_co= 0.08_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_co= 60.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_co= 1.6_rk !empirical coefficient real(rk), parameter ::caq_co= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_co= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_co= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_co= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_co= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_co= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_co= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_co= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_co= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_co= 1.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_co= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_co= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_bvoc= 0.8_rk !light-dependent fraction real(rk), parameter ::beta_bvoc= 0.13_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_bvoc= 95.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_bvoc= 2.0_rk !empirical coefficient real(rk), parameter ::caq_bvoc= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_bvoc= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_bvoc= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_bvoc= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_bvoc= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_bvoc= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_bvoc= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_bvoc= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_bvoc= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_bvoc= 1.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_bvoc= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_bvoc= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_svoc= 0.8_rk !light-dependent fraction real(rk), parameter ::beta_svoc= 0.1_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_svoc= 80.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_svoc= 1.83_rk !empirical coefficient real(rk), parameter ::caq_svoc= 5.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_svoc= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_svoc= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_svoc= 5.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_svoc= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_svoc= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_svoc= 5.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_svoc= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_svoc= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_svoc= 5.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_svoc= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_svoc= 8.0_rk !delta threshold for high wind stress(m/s) real(rk), parameter ::ldf_ovoc= 0.2_rk !light-dependent fraction real(rk), parameter ::beta_ovoc= 0.1_rk !empirical coefficient for temperature dependence of light-independent fraction real(rk), parameter ::ct1_ovoc= 80.0_rk !activation energy(kj/mol) real(rk), parameter ::ceo_ovoc= 1.83_rk !empirical coefficient real(rk), parameter ::caq_ovoc= 1.0_rk **coefficient** for poor air quality stress real(rk), parameter ::taq_ovoc= 20.0_rk !threshold for poor air quality stress(ppm-hours) real(rk), parameter ::dtaq_ovoc= 30.0_rk !delta threshold for poor air quality stress(ppm-hours) real(rk), parameter ::cht_ovoc= 1.0_rk **coefficient** for high temperature stress real(rk), parameter ::tht_ovoc= 313.15_rk !threshold for high temperature stress(k) real(rk), parameter ::dtht_ovoc= 8.0_rk !delta threshold for high temperature stress(k) real(rk), parameter ::clt_ovoc= 1.0_rk **coefficient** for low temperature stress real(rk), parameter ::tlt_ovoc= 283.15_rk !threshold for low temperature stress(k) real(rk), parameter ::dtlt_ovoc= 8.0_rk !delta threshold for low temperature stress(k) real(rk), parameter ::chw_ovoc= 1.0_rk **coefficient** for high wind stress real(rk), parameter ::thw_ovoc= 12.0_rk !threshold for high wind stress(m/s) real(rk), parameter ::dthw_ovoc= 8.0_rk !delta threshold for high wind stress(m/s) ! set tree and species dependent coefficients if(emi_ind .eq. 1) then ldf=ldf_isop beta=beta_isop ct1=ct1_isop ceo=ceo_isop ef1=ef1_isop ef2=ef2_isop ef3=ef3_isop ef4=ef4_isop ef5=ef5_isop ef6=ef6_isop ef7=ef7_isop ef8=ef8_isop ef9=ef9_isop ef10=ef10_isop ef11=ef11_isop ef12=ef12_isop ef13=ef13_isop ef14=ef14_isop ef15=ef15_isop anew=anew_isop agro=agro_isop amat=amat_isop aold=aold_isop caq=caq_isop taq=taq_isop dtaq=dtaq_isop cht=cht_isop tht=tht_isop dtht=dtht_isop clt=clt_isop tlt=tlt_isop dtlt=dtlt_isop chw=chw_isop thw=thw_isop dthw=dthw_isop else if(emi_ind .eq. 2) then ldf=ldf_myrc beta=beta_myrc ct1=ct1_myrc ceo=ceo_myrc ef1=ef1_myrc ef2=ef2_myrc ef3=ef3_myrc ef4=ef4_myrc ef5=ef5_myrc ef6=ef6_myrc ef7=ef7_myrc ef8=ef8_myrc ef9=ef9_myrc ef10=ef10_myrc ef11=ef11_myrc ef12=ef12_myrc ef13=ef13_myrc ef14=ef14_myrc ef15=ef15_myrc anew=anew_myrc agro=agro_myrc amat=amat_myrc aold=aold_myrc caq=caq_myrc taq=taq_myrc dtaq=dtaq_myrc cht=cht_myrc tht=tht_myrc dtht=dtht_myrc clt=clt_myrc tlt=tlt_myrc dtlt=dtlt_myrc chw=chw_myrc thw=thw_myrc dthw=dthw_myrc else if(emi_ind .eq. 3) then ldf=ldf_sabi beta=beta_sabi ct1=ct1_sabi ceo=ceo_sabi ef1=ef1_sabi ef2=ef2_sabi ef3=ef3_sabi ef4=ef4_sabi ef5=ef5_sabi ef6=ef6_sabi ef7=ef7_sabi ef8=ef8_sabi ef9=ef9_sabi ef10=ef10_sabi ef11=ef11_sabi ef12=ef12_sabi ef13=ef13_sabi ef14=ef14_sabi ef15=ef15_sabi anew=anew_sabi agro=agro_sabi amat=amat_sabi aold=aold_sabi caq=caq_sabi taq=taq_sabi dtaq=dtaq_sabi cht=cht_sabi tht=tht_sabi dtht=dtht_sabi clt=clt_sabi tlt=tlt_sabi dtlt=dtlt_sabi chw=chw_sabi thw=thw_sabi dthw=dthw_sabi else if(emi_ind .eq. 4) then ldf=ldf_limo beta=beta_limo ct1=ct1_limo ceo=ceo_limo ef1=ef1_limo ef2=ef2_limo ef3=ef3_limo ef4=ef4_limo ef5=ef5_limo ef6=ef6_limo ef7=ef7_limo ef8=ef8_limo ef9=ef9_limo ef10=ef10_limo ef11=ef11_limo ef12=ef12_limo ef13=ef13_limo ef14=ef14_limo ef15=ef15_limo anew=anew_limo agro=agro_limo amat=amat_limo aold=aold_limo caq=caq_limo taq=taq_limo dtaq=dtaq_limo cht=cht_limo tht=tht_limo dtht=dtht_limo clt=clt_limo tlt=tlt_limo dtlt=dtlt_limo chw=chw_limo thw=thw_limo dthw=dthw_limo else if(emi_ind .eq. 5) then ldf=ldf_care beta=beta_care ct1=ct1_care ceo=ceo_care ef1=ef1_care ef2=ef2_care ef3=ef3_care ef4=ef4_care ef5=ef5_care ef6=ef6_care ef7=ef7_care ef8=ef8_care ef9=ef9_care ef10=ef10_care ef11=ef11_care ef12=ef12_care ef13=ef13_care ef14=ef14_care ef15=ef15_care anew=anew_care agro=agro_care amat=amat_care aold=aold_care caq=caq_care taq=taq_care dtaq=dtaq_care cht=cht_care tht=tht_care dtht=dtht_care clt=clt_care tlt=tlt_care dtlt=dtlt_care chw=chw_care thw=thw_care dthw=dthw_care else if(emi_ind .eq. 6) then ldf=ldf_ocim beta=beta_ocim ct1=ct1_ocim ceo=ceo_ocim ef1=ef1_ocim ef2=ef2_ocim ef3=ef3_ocim ef4=ef4_ocim ef5=ef5_ocim ef6=ef6_ocim ef7=ef7_ocim ef8=ef8_ocim ef9=ef9_ocim ef10=ef10_ocim ef11=ef11_ocim ef12=ef12_ocim ef13=ef13_ocim ef14=ef14_ocim ef15=ef15_ocim anew=anew_ocim agro=agro_ocim amat=amat_ocim aold=aold_ocim caq=caq_ocim taq=taq_ocim dtaq=dtaq_ocim cht=cht_ocim tht=tht_ocim dtht=dtht_ocim clt=clt_ocim tlt=tlt_ocim dtlt=dtlt_ocim chw=chw_ocim thw=thw_ocim dthw=dthw_ocim else if(emi_ind .eq. 7) then ldf=ldf_bpin beta=beta_bpin ct1=ct1_bpin ceo=ceo_bpin ef1=ef1_bpin ef2=ef2_bpin ef3=ef3_bpin ef4=ef4_bpin ef5=ef5_bpin ef6=ef6_bpin ef7=ef7_bpin ef8=ef8_bpin ef9=ef9_bpin ef10=ef10_bpin ef11=ef11_bpin ef12=ef12_bpin ef13=ef13_bpin ef14=ef14_bpin ef15=ef15_bpin anew=anew_bpin agro=agro_bpin amat=amat_bpin aold=aold_bpin caq=caq_bpin taq=taq_bpin dtaq=dtaq_bpin cht=cht_bpin tht=tht_bpin dtht=dtht_bpin clt=clt_bpin tlt=tlt_bpin dtlt=dtlt_bpin chw=chw_bpin thw=thw_bpin dthw=dthw_bpin else if(emi_ind .eq. 8) then ldf=ldf_apin beta=beta_apin ct1=ct1_apin ceo=ceo_apin ef1=ef1_apin ef2=ef2_apin ef3=ef3_apin ef4=ef4_apin ef5=ef5_apin ef6=ef6_apin ef7=ef7_apin ef8=ef8_apin ef9=ef9_apin ef10=ef10_apin ef11=ef11_apin ef12=ef12_apin ef13=ef13_apin ef14=ef14_apin ef15=ef15_apin anew=anew_apin agro=agro_apin amat=amat_apin aold=aold_apin caq=caq_apin taq=taq_apin dtaq=dtaq_apin cht=cht_apin tht=tht_apin dtht=dtht_apin clt=clt_apin tlt=tlt_apin dtlt=dtlt_apin chw=chw_apin thw=thw_apin dthw=dthw_apin else if(emi_ind .eq. 9) then ldf=ldf_mono beta=beta_mono ct1=ct1_mono ceo=ceo_mono ef1=ef1_mono ef2=ef2_mono ef3=ef3_mono ef4=ef4_mono ef5=ef5_mono ef6=ef6_mono ef7=ef7_mono ef8=ef8_mono ef9=ef9_mono ef10=ef10_mono ef11=ef11_mono ef12=ef12_mono ef13=ef13_mono ef14=ef14_mono ef15=ef15_mono anew=anew_mono agro=agro_mono amat=amat_mono aold=aold_mono caq=caq_mono taq=taq_mono dtaq=dtaq_mono cht=cht_mono tht=tht_mono dtht=dtht_mono clt=clt_mono tlt=tlt_mono dtlt=dtlt_mono chw=chw_mono thw=thw_mono dthw=dthw_mono else if(emi_ind .eq. 10) then ldf=ldf_farn beta=beta_farn ct1=ct1_farn ceo=ceo_farn ef1=ef1_farn ef2=ef2_farn ef3=ef3_farn ef4=ef4_farn ef5=ef5_farn ef6=ef6_farn ef7=ef7_farn ef8=ef8_farn ef9=ef9_farn ef10=ef10_farn ef11=ef11_farn ef12=ef12_farn ef13=ef13_farn ef14=ef14_farn ef15=ef15_farn anew=anew_farn agro=agro_farn amat=amat_farn aold=aold_farn caq=caq_farn taq=taq_farn dtaq=dtaq_farn cht=cht_farn tht=tht_farn dtht=dtht_farn clt=clt_farn tlt=tlt_farn dtlt=dtlt_farn chw=chw_farn thw=thw_farn dthw=dthw_farn else if(emi_ind .eq. 11) then ldf=ldf_cary beta=beta_cary ct1=ct1_cary ceo=ceo_cary ef1=ef1_cary ef2=ef2_cary ef3=ef3_cary ef4=ef4_cary ef5=ef5_cary ef6=ef6_cary ef7=ef7_cary ef8=ef8_cary ef9=ef9_cary ef10=ef10_cary ef11=ef11_cary ef12=ef12_cary ef13=ef13_cary ef14=ef14_cary ef15=ef15_cary anew=anew_cary agro=agro_cary amat=amat_cary aold=aold_cary caq=caq_cary taq=taq_cary dtaq=dtaq_cary cht=cht_cary tht=tht_cary dtht=dtht_cary clt=clt_cary tlt=tlt_cary dtlt=dtlt_cary chw=chw_cary thw=thw_cary dthw=dthw_cary else if(emi_ind .eq. 12) then ldf=ldf_sesq beta=beta_sesq ct1=ct1_sesq ceo=ceo_sesq ef1=ef1_sesq ef2=ef2_sesq ef3=ef3_sesq ef4=ef4_sesq ef5=ef5_sesq ef6=ef6_sesq ef7=ef7_sesq ef8=ef8_sesq ef9=ef9_sesq ef10=ef10_sesq ef11=ef11_sesq ef12=ef12_sesq ef13=ef13_sesq ef14=ef14_sesq ef15=ef15_sesq anew=anew_sesq agro=agro_sesq amat=amat_sesq aold=aold_sesq caq=caq_sesq taq=taq_sesq dtaq=dtaq_sesq cht=cht_sesq tht=tht_sesq dtht=dtht_sesq clt=clt_sesq tlt=tlt_sesq dtlt=dtlt_sesq chw=chw_sesq thw=thw_sesq dthw=dthw_sesq else if(emi_ind .eq. 13) then ldf=ldf_mbol beta=beta_mbol ct1=ct1_mbol ceo=ceo_mbol ef1=ef1_mbol ef2=ef2_mbol ef3=ef3_mbol ef4=ef4_mbol ef5=ef5_mbol ef6=ef6_mbol ef7=ef7_mbol ef8=ef8_mbol ef9=ef9_mbol ef10=ef10_mbol ef11=ef11_mbol ef12=ef12_mbol ef13=ef13_mbol ef14=ef14_mbol ef15=ef15_mbol anew=anew_mbol agro=agro_mbol amat=amat_mbol aold=aold_mbol caq=caq_mbol taq=taq_mbol dtaq=dtaq_mbol cht=cht_mbol tht=tht_mbol dtht=dtht_mbol clt=clt_mbol tlt=tlt_mbol dtlt=dtlt_mbol chw=chw_mbol thw=thw_mbol dthw=dthw_mbol else if(emi_ind .eq. 14) then ldf=ldf_meth beta=beta_meth ct1=ct1_meth ceo=ceo_meth ef1=ef1_meth ef2=ef2_meth ef3=ef3_meth ef4=ef4_meth ef5=ef5_meth ef6=ef6_meth ef7=ef7_meth ef8=ef8_meth ef9=ef9_meth ef10=ef10_meth ef11=ef11_meth ef12=ef12_meth ef13=ef13_meth ef14=ef14_meth ef15=ef15_meth anew=anew_meth agro=agro_meth amat=amat_meth aold=aold_meth caq=caq_meth taq=taq_meth dtaq=dtaq_meth cht=cht_meth tht=tht_meth dtht=dtht_meth clt=clt_meth tlt=tlt_meth dtlt=dtlt_meth chw=chw_meth thw=thw_meth dthw=dthw_meth else if(emi_ind .eq. 15) then ldf=ldf_acet beta=beta_acet ct1=ct1_acet ceo=ceo_acet ef1=ef1_acet ef2=ef2_acet ef3=ef3_acet ef4=ef4_acet ef5=ef5_acet ef6=ef6_acet ef7=ef7_acet ef8=ef8_acet ef9=ef9_acet ef10=ef10_acet ef11=ef11_acet ef12=ef12_acet ef13=ef13_acet ef14=ef14_acet ef15=ef15_acet anew=anew_acet agro=agro_acet amat=amat_acet aold=aold_acet caq=caq_acet taq=taq_acet dtaq=dtaq_acet cht=cht_acet tht=tht_acet dtht=dtht_acet clt=clt_acet tlt=tlt_acet dtlt=dtlt_acet chw=chw_acet thw=thw_acet dthw=dthw_acet else if(emi_ind .eq. 16) then ldf=ldf_co beta=beta_co ct1=ct1_co ceo=ceo_co ef1=ef1_co ef2=ef2_co ef3=ef3_co ef4=ef4_co ef5=ef5_co ef6=ef6_co ef7=ef7_co ef8=ef8_co ef9=ef9_co ef10=ef10_co ef11=ef11_co ef12=ef12_co ef13=ef13_co ef14=ef14_co ef15=ef15_co anew=anew_co agro=agro_co amat=amat_co aold=aold_co caq=caq_co taq=taq_co dtaq=dtaq_co cht=cht_co tht=tht_co dtht=dtht_co clt=clt_co tlt=tlt_co dtlt=dtlt_co chw=chw_co thw=thw_co dthw=dthw_co else if(emi_ind .eq. 17) then ldf=ldf_bvoc beta=beta_bvoc ct1=ct1_bvoc ceo=ceo_bvoc ef1=ef1_bvoc ef2=ef2_bvoc ef3=ef3_bvoc ef4=ef4_bvoc ef5=ef5_bvoc ef6=ef6_bvoc ef7=ef7_bvoc ef8=ef8_bvoc ef9=ef9_bvoc ef10=ef10_bvoc ef11=ef11_bvoc ef12=ef12_bvoc ef13=ef13_bvoc ef14=ef14_bvoc ef15=ef15_bvoc anew=anew_bvoc agro=agro_bvoc amat=amat_bvoc aold=aold_bvoc caq=caq_bvoc taq=taq_bvoc dtaq=dtaq_bvoc cht=cht_bvoc tht=tht_bvoc dtht=dtht_bvoc clt=clt_bvoc tlt=tlt_bvoc dtlt=dtlt_bvoc chw=chw_bvoc thw=thw_bvoc dthw=dthw_bvoc else if(emi_ind .eq. 18) then ldf=ldf_svoc beta=beta_svoc ct1=ct1_svoc ceo=ceo_svoc ef1=ef1_svoc ef2=ef2_svoc ef3=ef3_svoc ef4=ef4_svoc ef5=ef5_svoc ef6=ef6_svoc ef7=ef7_svoc ef8=ef8_svoc ef9=ef9_svoc ef10=ef10_svoc ef11=ef11_svoc ef12=ef12_svoc ef13=ef13_svoc ef14=ef14_svoc ef15=ef15_svoc anew=anew_svoc agro=agro_svoc amat=amat_svoc aold=aold_svoc caq=caq_svoc taq=taq_svoc dtaq=dtaq_svoc cht=cht_svoc tht=tht_svoc dtht=dtht_svoc clt=clt_svoc tlt=tlt_svoc dtlt=dtlt_svoc chw=chw_svoc thw=thw_svoc dthw=dthw_svoc else ! emi_ind=19 ldf=ldf_ovoc beta=beta_ovoc ct1=ct1_ovoc ceo=ceo_ovoc ef1=ef1_ovoc ef2=ef2_ovoc ef3=ef3_ovoc ef4=ef4_ovoc ef5=ef5_ovoc ef6=ef6_ovoc ef7=ef7_ovoc ef8=ef8_ovoc ef9=ef9_ovoc ef10=ef10_ovoc ef11=ef11_ovoc ef12=ef12_ovoc ef13=ef13_ovoc ef14=ef14_ovoc ef15=ef15_ovoc anew=anew_ovoc agro=agro_ovoc amat=amat_ovoc aold=aold_ovoc caq=caq_ovoc taq=taq_ovoc dtaq=dtaq_ovoc cht=cht_ovoc tht=tht_ovoc dtht=dtht_ovoc clt=clt_ovoc tlt=tlt_ovoc dtlt=dtlt_ovoc chw=chw_ovoc thw=thw_ovoc dthw=dthw_ovoc end if if(lu_opt .eq. 0 .or. lu_opt .eq. 1) then !viirs or modis lu types ! simple megan(table 3 in guenther et al., 2012) pft to viirs/modis vtype mapping if(vtype .eq. 1) then !viirs cat 1 evergreen needleleaf !--> average needleleaf evergreen temperate tree and needleleaf evergreen boreal tree ef=(ef1+ef2)/2.0_rk **set** pft dependent a and b coefficients for cumulative root depth fraction(zeng 2001) !see table 2 for igpb classification at: https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml roota=6.706_rk rootb=2.175_rk else if(vtype .eq. 2) then !viirs/modis cat 2 evergreen broadleaf !--> average broadleaf evergreen tropical tree and broadleaf evergreen temperate tree ef=(ef4+ef5)/2.0_rk **set** pft dependent a and b coefficients for cumulative root depth fraction(zeng 2001) !see table 2 for igpb classification at: https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml roota=7.344_rk rootb=1.303_rk else if(vtype .eq. 3) then !viirs/modis cat 3 deciduous needleaf !--> average needleleaf deciduous boreal tree ef=ef3 **set** pft dependent a and b coefficients for cumulative root depth fraction(zeng 2001) !see table 2 for igpb classification at: https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml roota=7.066_rk rootb=1.953_rk else if(vtype .eq. 4) then !viirs/modis cat 4 deciduous broadleaf !--> average broadleaf deciduous tropical tree, broadleaf deciduous temperate tree, ! and broadleaf deciduous boreal tree ef=(ef6+ef7+ef8)/3.0_rk **set** pft dependent a and b coefficients for cumulative root depth fraction(zeng 2001) !see table 2 for igpb classification at: https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml roota=5.990_rk rootb=1.955_rk else if(vtype .eq. 5) then !viirs/modis cat 5 mixed forests !--> avearge of all above ef1-ef8 pfts. ef=(ef1+ef2+ef3+ef4+ef5+ef6+ef7+ef8)/8.0_rk **set** pft dependent a and b coefficients for cumulative root depth fraction(zeng 2001) !see table 2 for igpb classification at: https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml roota=4.453_rk rootb=1.631_rk else if(vtype .ge. 6 .and. vtype .le. 7) then !viirs/modis cat 6-7 closed/open shrublands !--> avearge broadleaf evergreen temperate shrub, broadleaf deciduous temperate shrub, ! and broadleaf deciduous boreal shrub ef=(ef9+ef10+ef11)/3.0_rk **set** pft dependent a and b coefficients for cumulative root depth fraction(zeng 2001) !see table 2 for igpb classification at: https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml roota=(6.326_rk+7.718_rk)/2.0_rk rootb=(1.567_rk+1.262_rk)/2.0_rk else if(vtype .ge. 8 .and. vtype .le. 11) then !viirs/modis cat 8-10 savannas and grasslands !--> avearge arctic c3 grass, cool c3 grass, warm c4 grass), intent(out) ef
real(rk), intent(out) out

Public Attributes Documentation

variable ef

real(rk), dimension ((ug/m2 hr)
        real(rk),    intent( out )      :: ldf             !> light-dependent fraction
        real(rk),    intent( out )      :: beta            !> empirical coefficient for temperature dependence of light-independent fraction
        real(rk),    intent( out )      :: ct1             !> out activation energy (kj/mol)
        real(rk),    intent( out )      :: ceo             !> out empirical coefficient
        real(rk),    intent( out )      :: anew, agro, amat, aold   !> empirical factors or coefficients for: growing, mature, and old/senescing foliage, as per table 4 of guenther et al., 2012
        real(rk),    intent( out )      :: roota, rootb    !> coefficients a and b used for pft dependent cumulative root depth fraction [m-1]
        real(rk),    intent( out )      :: caq             !> coefficient for poor air quality stress
        real(rk),    intent( out )      :: taq             !> threshold for poor air quality stress (ppm-hours)
        real(rk),    intent( out )      :: dtaq            !> delta threshold for poor air quality stress (ppm-hours)
        real(rk),    intent( out )      :: cht             !> coefficient for high temperature stress
        real(rk),    intent( out )      :: tht             !> threshold for high temperature stress (k)
        real(rk),    intent( out )      :: dtht            !> delta threshold high temperature stress (k)
        real(rk),    intent( out )      :: clt             !> coefficient for low temperature stress
        real(rk),    intent( out )      :: tlt             !> threshold for low temperature stress (k)
        real(rk),    intent( out )      :: dtlt            !> delta threshold low temperature stress (k)
        real(rk),    intent( out )      :: chw             !> coefficient for high wind stress
        real(rk),    intent( out )      :: thw             !> threshold for high wind stress (m/s)
        real(rk),    intent( out )      :: dthw            !> delta threshold high wind stress (m/s)
!> \}

!> \defgroup bioparm_local_vars local variables
!! \brief local variables for parameter assignment
!! \{
        real(rk) :: ef1,ef2,ef3,ef4,ef5,ef6,ef7    !> plant emission factors (ef) (ug/m2 hr)
        real(rk) :: ef8,ef9,ef10,ef11,ef12,ef13    !> plant emission factors (ef) (ug/m2 hr)
        real(rk) :: ef14,ef15                      !> plant emission factors (ef) (ug/m2 hr)
!> \}

!> \defgroup bioparm_isop_params isoprene parameters
!! \brief plant-dependent emission capacity factors for isoprene from tables 2-3 of guenther et al. (2012)
!! \{

        !> \brief needleleaf evergreen temperate tree isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef1_isop    =  600.0_rk
        !> \brief needleleaf evergreen boreal tree isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef2_isop    =  3000.0_rk
        !> \brief needleleaf deciduous boreal tree isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef3_isop    =  1.0_rk
        !> \brief broadleaf evergreen tropical tree isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef4_isop    =  7000.0_rk
        !> \brief broadleaf evergreen temperate tree isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef5_isop    =  10000.0_rk
        !> \brief broadleaf deciduous tropical tree isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef6_isop    =  7000.0_rk
        !> \brief broadleaf deciduous temperate tree isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef7_isop    =  10000.0_rk
        !> \brief broadleaf deciduous boreal tree isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef8_isop    =  11000.0_rk
        !> \brief broadleaf evergreen temperate shrub isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef9_isop    =  2000.0_rk
        !> \brief broadleaf deciduous temperate shrub isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef10_isop   =  4000.0_rk
        !> \brief broadleaf deciduous boreal shrub isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef11_isop   =  4000.0_rk
        !> \brief arctic c3 grass isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef12_isop   =  1600.0_rk
        !> \brief cool c3 grass isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef13_isop   =  800.0_rk
        !> \brief warm c4 grass isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef14_isop   =  200.0_rk
        !> \brief crop1 isoprene ef (μg/m²/hr)
        real(rk),          parameter     :: ef15_isop   =  1.0_rk

        !> \brief isoprene leaf age factor for new foliage (table 4 of guenther et al., 2012)
        real(rk),          parameter     :: anew_isop  = 0.05_rk
        !> \brief isoprene leaf age factor for growing foliage (table 4 of guenther et al., 2012)
        real(rk),          parameter     :: agro_isop  = 0.6_rk
        !> \brief isoprene leaf age factor for mature foliage (table 4 of guenther et al., 2012)
        real(rk),          parameter     :: amat_isop  = 1.0_rk
        !> \brief isoprene leaf age factor for old/senescing foliage (table 4 of guenther et al., 2012)
        real(rk),          parameter     :: aold_isop  = 0.9_rk

!> \}

! plant-dependent emissions capacity/factors (efs) for myrcene (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_myrc    =  70.0_rk      ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_myrc    =  70.0_rk      ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_myrc    =  60.0_rk      ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_myrc    =  80.0_rk      ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_myrc    =  30.0_rk      ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_myrc    =  80.0_rk      ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_myrc    =  30.0_rk      ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_myrc    =  30.0_rk      ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_myrc    =  30.0_rk      ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_myrc   =  50.0_rk      ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_myrc   =  30.0_rk      ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_myrc   =  0.3_rk       ! arctic c3 grass
        real(rk),          parameter     :: ef13_myrc   =  0.3_rk       ! cool c3 grass
        real(rk),          parameter     :: ef14_myrc   =  0.3_rk       ! warm c4 grass
        real(rk),          parameter     :: ef15_myrc   =  0.3_rk       ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for myrcene as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_myrc  = 2.0_rk
        real(rk),          parameter     :: agro_myrc  = 1.8_rk
        real(rk),          parameter     :: amat_myrc  = 1.0_rk
        real(rk),          parameter     :: aold_myrc  = 1.05_rk

! plant-dependent emissions capacity/factors (efs) for sabinene (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_sabi    =  70.0_rk      ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_sabi    =  70.0_rk      ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_sabi    =  40.0_rk      ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_sabi    =  80.0_rk      ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_sabi    =  50.0_rk      ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_sabi    =  80.0_rk      ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_sabi    =  50.0_rk      ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_sabi    =  50.0_rk      ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_sabi    =  50.0_rk      ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_sabi   =  70.0_rk      ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_sabi   =  50.0_rk      ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_sabi   =  0.7_rk       ! arctic c3 grass
        real(rk),          parameter     :: ef13_sabi   =  0.7_rk       ! cool c3 grass
        real(rk),          parameter     :: ef14_sabi   =  0.7_rk       ! warm c4 grass
        real(rk),          parameter     :: ef15_sabi   =  0.7_rk       ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for sabinene as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_sabi  = 2.0_rk
        real(rk),          parameter     :: agro_sabi  = 1.8_rk
        real(rk),          parameter     :: amat_sabi  = 1.0_rk
        real(rk),          parameter     :: aold_sabi  = 1.05_rk

! plant-dependent emissions capacity/factors (efs) for limonene (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_limo    =  100.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_limo    =  100.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_limo    =  130.0_rk     ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_limo    =  80.0_rk      ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_limo    =  80.0_rk      ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_limo    =  80.0_rk      ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_limo    =  80.0_rk      ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_limo    =  80.0_rk      ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_limo    =  60.0_rk      ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_limo   =  100.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_limo   =  60.0_rk      ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_limo   =  0.7_rk       ! arctic c3 grass
        real(rk),          parameter     :: ef13_limo   =  0.7_rk       ! cool c3 grass
        real(rk),          parameter     :: ef14_limo   =  0.7_rk       ! warm c4 grass
        real(rk),          parameter     :: ef15_limo   =  0.7_rk       ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for limonene as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_limo  = 2.0_rk
        real(rk),          parameter     :: agro_limo  = 1.8_rk
        real(rk),          parameter     :: amat_limo  = 1.0_rk
        real(rk),          parameter     :: aold_limo  = 1.05_rk

! plant-dependent emissions capacity/factors (efs) for 3-carene (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_care    =  160.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_care    =  160.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_care    =  80.0_rk      ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_care    =  40.0_rk      ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_care    =  30.0_rk      ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_care    =  40.0_rk      ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_care    =  30.0_rk      ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_care    =  30.0_rk      ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_care    =  30.0_rk      ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_care   =  100.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_care   =  30.0_rk      ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_care   =  0.3_rk       ! arctic c3 grass
        real(rk),          parameter     :: ef13_care   =  0.3_rk       ! cool c3 grass
        real(rk),          parameter     :: ef14_care   =  0.3_rk       ! warm c4 grass
        real(rk),          parameter     :: ef15_care   =  0.3_rk       ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for 3-carene as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_care  = 2.0_rk
        real(rk),          parameter     :: agro_care  = 1.8_rk
        real(rk),          parameter     :: amat_care  = 1.0_rk
        real(rk),          parameter     :: aold_care  = 1.05_rk

! plant-dependent emissions capacity/factors (efs) for t-beta-ocimene (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_ocim    =  70.0_rk      ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_ocim    =  70.0_rk      ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_ocim    =  60.0_rk      ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_ocim    =  150.0_rk     ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_ocim    =  120.0_rk     ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_ocim    =  150.0_rk     ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_ocim    =  120.0_rk     ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_ocim    =  120.0_rk     ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_ocim    =  90.0_rk      ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_ocim   =  150.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_ocim   =  90.0_rk      ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_ocim   =  2.0_rk       ! arctic c3 grass
        real(rk),          parameter     :: ef13_ocim   =  2.0_rk       ! cool c3 grass
        real(rk),          parameter     :: ef14_ocim   =  2.0_rk       ! warm c4 grass
        real(rk),          parameter     :: ef15_ocim   =  2.0_rk       ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for t-beta-ocimene as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_ocim  = 2.0_rk
        real(rk),          parameter     :: agro_ocim  = 1.8_rk
        real(rk),          parameter     :: amat_ocim  = 1.0_rk
        real(rk),          parameter     :: aold_ocim  = 1.05_rk

! plant-dependent emissions capacity/factors (efs) for beta-pinene (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_bpin    =  300.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_bpin    =  300.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_bpin    =  200.0_rk     ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_bpin    =  120.0_rk     ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_bpin    =  130.0_rk     ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_bpin    =  120.0_rk     ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_bpin    =  130.0_rk     ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_bpin    =  130.0_rk     ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_bpin    =  100.0_rk     ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_bpin   =  150.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_bpin   =  100.0_rk     ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_bpin   =  1.5_rk       ! arctic c3 grass
        real(rk),          parameter     :: ef13_bpin   =  1.5_rk       ! cool c3 grass
        real(rk),          parameter     :: ef14_bpin   =  1.5_rk       ! warm c4 grass
        real(rk),          parameter     :: ef15_bpin   =  1.5_rk       ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for beta-pinene as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_bpin  = 2.0_rk
        real(rk),          parameter     :: agro_bpin  = 1.8_rk
        real(rk),          parameter     :: amat_bpin  = 1.0_rk
        real(rk),          parameter     :: aold_bpin  = 1.05_rk


! plant-dependent emissions capacity/factors (efs) for alpha-pinene (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_apin    =  500.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_apin    =  500.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_apin    =  510.0_rk     ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_apin    =  600.0_rk     ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_apin    =  400.0_rk     ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_apin    =  600.0_rk     ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_apin    =  400.0_rk     ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_apin    =  400.0_rk     ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_apin    =  200.0_rk     ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_apin   =  300.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_apin   =  200.0_rk     ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_apin   =  2.0_rk       ! arctic c3 grass
        real(rk),          parameter     :: ef13_apin   =  2.0_rk       ! cool c3 grass
        real(rk),          parameter     :: ef14_apin   =  2.0_rk       ! warm c4 grass
        real(rk),          parameter     :: ef15_apin   =  2.0_rk       ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for alpha-pinene as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_apin  = 2.0_rk
        real(rk),          parameter     :: agro_apin  = 1.8_rk
        real(rk),          parameter     :: amat_apin  = 1.0_rk
        real(rk),          parameter     :: aold_apin  = 1.05_rk

! plant-dependent emissions capacity/factors (efs) for other monoterpenes (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
! ! other monoterpenes category (34 compounds):  see table 1 of guenther et al. (2012)
        real(rk),          parameter     :: ef1_mono    =  180.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_mono    =  180.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_mono    =  170.0_rk     ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_mono    =  150.0_rk     ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_mono    =  150.0_rk     ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_mono    =  150.0_rk     ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_mono    =  150.0_rk     ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_mono    =  150.0_rk     ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_mono    =  110.0_rk     ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_mono   =  200.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_mono   =  110.0_rk     ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_mono   =  5.0_rk       ! arctic c3 grass
        real(rk),          parameter     :: ef13_mono   =  5.0_rk       ! cool c3 grass
        real(rk),          parameter     :: ef14_mono   =  5.0_rk       ! warm c4 grass
        real(rk),          parameter     :: ef15_mono   =  5.0_rk       ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for other monoterpenes as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_mono  = 2.0_rk
        real(rk),          parameter     :: agro_mono  = 1.8_rk
        real(rk),          parameter     :: amat_mono  = 1.0_rk
        real(rk),          parameter     :: aold_mono  = 1.05_rk

! plant-dependent emissions capacity/factors (efs) for alpha-farnesene (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_farn    =  40.0_rk      ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_farn    =  40.0_rk      ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_farn    =  40.0_rk      ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_farn    =  60.0_rk      ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_farn    =  40.0_rk      ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_farn    =  60.0_rk      ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_farn    =  40.0_rk      ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_farn    =  40.0_rk      ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_farn    =  40.0_rk      ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_farn   =  40.0_rk      ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_farn   =  40.0_rk      ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_farn   =  3.0_rk       ! arctic c3 grass
        real(rk),          parameter     :: ef13_farn   =  3.0_rk       ! cool c3 grass
        real(rk),          parameter     :: ef14_farn   =  3.0_rk       ! warm c4 grass
        real(rk),          parameter     :: ef15_farn   =  4.0_rk       ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for alpha-farnesene as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_farn  = 0.4_rk
        real(rk),          parameter     :: agro_farn  = 0.6_rk
        real(rk),          parameter     :: amat_farn  = 1.0_rk
        real(rk),          parameter     :: aold_farn  = 0.95_rk

! plant-dependent emissions capacity/factors (efs) for beta-caryophyllene (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_cary    =  80.0_rk      ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_cary    =  80.0_rk      ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_cary    =  80.0_rk      ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_cary    =  60.0_rk      ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_cary    =  40.0_rk      ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_cary    =  60.0_rk      ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_cary    =  40.0_rk      ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_cary    =  40.0_rk      ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_cary    =  50.0_rk      ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_cary   =  50.0_rk      ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_cary   =  50.0_rk      ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_cary   =  1.0_rk       ! arctic c3 grass
        real(rk),          parameter     :: ef13_cary   =  1.0_rk       ! cool c3 grass
        real(rk),          parameter     :: ef14_cary   =  1.0_rk       ! warm c4 grass
        real(rk),          parameter     :: ef15_cary   =  4.0_rk       ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for beta-caryophyllene as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_cary  = 0.4_rk
        real(rk),          parameter     :: agro_cary  = 0.6_rk
        real(rk),          parameter     :: amat_cary  = 1.0_rk
        real(rk),          parameter     :: aold_cary  = 0.95_rk

! plant-dependent emissions capacity/factors (efs) for other sesquieterpenes (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
! other sesquiterpenes category (30 compounds):  see table 1 of guenther et al. (2012)
        real(rk),          parameter     :: ef1_sesq    =  120.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_sesq    =  120.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_sesq    =  120.0_rk     ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_sesq    =  120.0_rk     ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_sesq    =  100.0_rk     ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_sesq    =  120.0_rk     ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_sesq    =  100.0_rk     ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_sesq    =  100.0_rk     ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_sesq    =  100.0_rk     ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_sesq   =  100.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_sesq   =  100.0_rk     ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_sesq   =  1.0_rk       ! arctic c3 grass
        real(rk),          parameter     :: ef13_sesq   =  1.0_rk       ! cool c3 grass
        real(rk),          parameter     :: ef14_sesq   =  1.0_rk       ! warm c4 grass
        real(rk),          parameter     :: ef15_sesq   =  1.0_rk       ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for other sesquieterpenes as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_sesq  = 0.4_rk
        real(rk),          parameter     :: agro_sesq  = 0.6_rk
        real(rk),          parameter     :: amat_sesq  = 1.0_rk
        real(rk),          parameter     :: aold_sesq  = 0.95_rk

! plant-dependent emissions capacity/factors (efs) for 232-mbo (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_mbol    =  700.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_mbol    =  60.0_rk      ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_mbol    =  0.01_rk      ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_mbol    =  0.01_rk      ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_mbol    =  0.01_rk      ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_mbol    =  0.01_rk      ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_mbol    =  0.01_rk      ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_mbol    =  2.0_rk       ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_mbol    =  0.01_rk      ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_mbol   =  0.01_rk      ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_mbol   =  0.01_rk      ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_mbol   =  0.01_rk      ! arctic c3 grass
        real(rk),          parameter     :: ef13_mbol   =  0.01_rk      ! cool c3 grass
        real(rk),          parameter     :: ef14_mbol   =  0.01_rk      ! warm c4 grass
        real(rk),          parameter     :: ef15_mbol   =  0.01_rk      ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for 232-mbo as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_mbol  = 0.05_rk
        real(rk),          parameter     :: agro_mbol  = 0.6_rk
        real(rk),          parameter     :: amat_mbol  = 1.0_rk
        real(rk),          parameter     :: aold_mbol  = 0.9_rk

! plant-dependent emissions capacity/factors (efs) for methanol (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_meth    =  900.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_meth    =  900.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_meth    =  900.0_rk     ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_meth    =  500.0_rk     ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_meth    =  900.0_rk     ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_meth    =  500.0_rk     ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_meth    =  900.0_rk     ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_meth    =  900.0_rk     ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_meth    =  900.0_rk     ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_meth   =  900.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_meth   =  900.0_rk     ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_meth   =  500.0_rk     ! arctic c3 grass
        real(rk),          parameter     :: ef13_meth   =  500.0_rk     ! cool c3 grass
        real(rk),          parameter     :: ef14_meth   =  500.0_rk     ! warm c4 grass
        real(rk),          parameter     :: ef15_meth   =  900.0_rk     ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for methanol as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_meth  = 3.5_rk
        real(rk),          parameter     :: agro_meth  = 3.0_rk
        real(rk),          parameter     :: amat_meth  = 1.0_rk
        real(rk),          parameter     :: aold_meth  = 1.2_rk

! plant-dependent emissions capacity/factors (efs) for acetone (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_acet    =  240.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_acet    =  240.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_acet    =  240.0_rk     ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_acet    =  240.0_rk     ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_acet    =  240.0_rk     ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_acet    =  240.0_rk     ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_acet    =  240.0_rk     ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_acet    =  240.0_rk     ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_acet    =  240.0_rk     ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_acet   =  240.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_acet   =  240.0_rk     ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_acet   =  80.0_rk      ! arctic c3 grass
        real(rk),          parameter     :: ef13_acet   =  80.0_rk      ! cool c3 grass
        real(rk),          parameter     :: ef14_acet   =  80.0_rk      ! warm c4 grass
        real(rk),          parameter     :: ef15_acet   =  80.0_rk      ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for acetone as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_acet  = 1.0_rk
        real(rk),          parameter     :: agro_acet  = 1.0_rk
        real(rk),          parameter     :: amat_acet  = 1.0_rk
        real(rk),          parameter     :: aold_acet  = 1.0_rk

! plant-dependent emissions capacity/factors (efs) for carbon monoxide (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
        real(rk),          parameter     :: ef1_co      =  600.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_co      =  600.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_co      =  600.0_rk     ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_co      =  600.0_rk     ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_co      =  600.0_rk     ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_co      =  600.0_rk     ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_co      =  600.0_rk     ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_co      =  600.0_rk     ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_co      =  600.0_rk     ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_co     =  600.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_co     =  600.0_rk     ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_co     =  600.0_rk     ! arctic c3 grass
        real(rk),          parameter     :: ef13_co     =  600.0_rk     ! cool c3 grass
        real(rk),          parameter     :: ef14_co     =  600.0_rk     ! warm c4 grass
        real(rk),          parameter     :: ef15_co     =  600.0_rk     ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for carbon monoxide as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_co  = 1.0_rk
        real(rk),          parameter     :: agro_co  = 1.0_rk
        real(rk),          parameter     :: amat_co  = 1.0_rk
        real(rk),          parameter     :: aold_co  = 1.0_rk

! plant-dependent emissions capacity/factors (efs) for bidi voc species (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
! bidirectional voc (5 compounds): see table 1 of guenther et al. (2012)
        real(rk),          parameter     :: ef1_bvoc    =  500.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_bvoc    =  500.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_bvoc    =  500.0_rk     ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_bvoc    =  500.0_rk     ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_bvoc    =  500.0_rk     ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_bvoc    =  500.0_rk     ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_bvoc    =  500.0_rk     ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_bvoc    =  500.0_rk     ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_bvoc    =  500.0_rk     ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_bvoc   =  500.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_bvoc   =  500.0_rk     ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_bvoc   =  80.0_rk      ! arctic c3 grass
        real(rk),          parameter     :: ef13_bvoc   =  80.0_rk      ! cool c3 grass
        real(rk),          parameter     :: ef14_bvoc   =  80.0_rk      ! warm c4 grass
        real(rk),          parameter     :: ef15_bvoc   =  80.0_rk      ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for bidi voc species as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_bvoc  = 1.0_rk
        real(rk),          parameter     :: agro_bvoc  = 1.0_rk
        real(rk),          parameter     :: amat_bvoc  = 1.0_rk
        real(rk),          parameter     :: aold_bvoc  = 1.0_rk

! plant-dependent emissions capacity/factors (efs) for stress vocs (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
! stress voc (15 compounds): see table 1 of guenther et al. (2012)
        real(rk),          parameter     :: ef1_svoc    =  300.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_svoc    =  300.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_svoc    =  300.0_rk     ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_svoc    =  300.0_rk     ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_svoc    =  300.0_rk     ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_svoc    =  300.0_rk     ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_svoc    =  300.0_rk     ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_svoc    =  300.0_rk     ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_svoc    =  300.0_rk     ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_svoc   =  300.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_svoc   =  300.0_rk     ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_svoc   =  300.0_rk     ! arctic c3 grass
        real(rk),          parameter     :: ef13_svoc   =  300.0_rk     ! cool c3 grass
        real(rk),          parameter     :: ef14_svoc   =  300.0_rk     ! warm c4 grass
        real(rk),          parameter     :: ef15_svoc   =  300.0_rk     ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for stress vocs as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_svoc  = 1.0_rk
        real(rk),          parameter     :: agro_svoc  = 1.0_rk
        real(rk),          parameter     :: amat_svoc  = 1.0_rk
        real(rk),          parameter     :: aold_svoc  = 1.0_rk

! plant-dependent emissions capacity/factors (efs) for other vocs (tables 2-3 of guenther et al., 2012) (ug/m2 hr)
! other voc (49 compounds): see table 1 of guenther et al. (2012)
        real(rk),          parameter     :: ef1_ovoc    =  140.0_rk     ! needleleaf evergreen temperate tree
        real(rk),          parameter     :: ef2_ovoc    =  140.0_rk     ! needleleaf evergreen boreal tree
        real(rk),          parameter     :: ef3_ovoc    =  140.0_rk     ! needleleaf deciduous boreal tree
        real(rk),          parameter     :: ef4_ovoc    =  140.0_rk     ! broadleaf evergreen tropical tree
        real(rk),          parameter     :: ef5_ovoc    =  140.0_rk     ! broadleaf evergreen temperate tree
        real(rk),          parameter     :: ef6_ovoc    =  140.0_rk     ! broadleaf deciduous tropical tree
        real(rk),          parameter     :: ef7_ovoc    =  140.0_rk     ! broadleaf deciduous temperate tree
        real(rk),          parameter     :: ef8_ovoc    =  140.0_rk     ! broadleaf deciduous boreal tree
        real(rk),          parameter     :: ef9_ovoc    =  140.0_rk     ! broadleaf evergreen temperate shrub
        real(rk),          parameter     :: ef10_ovoc   =  140.0_rk     ! broadleaf deciduous temperate shrub
        real(rk),          parameter     :: ef11_ovoc   =  140.0_rk     ! broadleaf deciduous boreal shrub
        real(rk),          parameter     :: ef12_ovoc   =  140.0_rk     ! arctic c3 grass
        real(rk),          parameter     :: ef13_ovoc   =  140.0_rk     ! cool c3 grass
        real(rk),          parameter     :: ef14_ovoc   =  140.0_rk     ! warm c4 grass
        real(rk),          parameter     :: ef15_ovoc   =  140.0_rk     ! crop1

!empirical factors or coefficients for: growing, mature, and old/senescing foliage, for other vocs as per table 4 of guenther et al., 2012
        real(rk),          parameter     :: anew_ovoc  = 1.0_rk
        real(rk),          parameter     :: agro_ovoc  = 1.0_rk
        real(rk),          parameter     :: amat_ovoc  = 1.0_rk
        real(rk),          parameter     :: aold_ovoc  = 1.0_rk

! species-dependent parameterized canopy model parameters (table 4 of guenther et al., 2012)
        real(rk),          parameter     :: ldf_isop         =  1.0_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_isop        =  0.13_rk    !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_isop         =  95.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_isop         =  2.0_rk     !empirical coefficient
        real(rk),          parameter     :: caq_isop         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_isop         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_isop        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_isop         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_isop         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_isop        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_isop         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_isop         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_isop        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_isop         =  1.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_isop         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_isop        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_myrc         =  0.6_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_myrc        =  0.1_rk     !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_myrc         =  80.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_myrc         =  1.83_rk    !empirical coefficient
        real(rk),          parameter     :: caq_myrc         =  5.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_myrc         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_myrc        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_myrc         =  5.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_myrc         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_myrc        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_myrc         =  5.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_myrc         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_myrc        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_myrc         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_myrc         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_myrc        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_sabi         =  0.6_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_sabi        =  0.1_rk     !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_sabi         =  80.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_sabi         =  1.83_rk    !empirical coefficient
        real(rk),          parameter     :: caq_sabi         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_sabi         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_sabi        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_sabi         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_sabi         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_sabi        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_sabi         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_sabi         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_sabi        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_sabi         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_sabi         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_sabi        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_limo         =  0.2_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_limo        =  0.1_rk     !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_limo         =  80.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_limo         =  1.83_rk    !empirical coefficient
        real(rk),          parameter     :: caq_limo         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_limo         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_limo        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_limo         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_limo         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_limo        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_limo         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_limo         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_limo        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_limo         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_limo         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_limo        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_care         =  0.2_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_care        =  0.1_rk     !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_care         =  80.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_care         =  1.83_rk    !empirical coefficient
        real(rk),          parameter     :: caq_care         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_care         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_care        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_care         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_care         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_care        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_care         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_care         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_care        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_care         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_care         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_care        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_ocim         =  0.8_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_ocim        =  0.1_rk     !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_ocim         =  80.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_ocim         =  1.83_rk    !empirical coefficient
        real(rk),          parameter     :: caq_ocim         =  5.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_ocim         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_ocim        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_ocim         =  5.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_ocim         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_ocim        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_ocim         =  5.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_ocim         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_ocim        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_ocim         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_ocim         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_ocim        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_bpin         =  0.2_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_bpin        =  0.1_rk     !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_bpin         =  80.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_bpin         =  1.83_rk    !empirical coefficient
        real(rk),          parameter     :: caq_bpin         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_bpin         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_bpin        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_bpin         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_bpin         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_bpin        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_bpin         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_bpin         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_bpin        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_bpin         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_bpin         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_bpin        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_apin         =  0.6_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_apin        =  0.1_rk     !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_apin         =  80.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_apin         =  1.83_rk    !empirical coefficient
        real(rk),          parameter     :: caq_apin         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_apin         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_apin        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_apin         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_apin         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_apin        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_apin         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_apin         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_apin        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_apin         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_apin         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_apin        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_mono         =  0.4_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_mono        =  0.1_rk     !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_mono         =  80.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_mono         =  1.83_rk    !empirical coefficient
        real(rk),          parameter     :: caq_mono         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_mono         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_mono        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_mono         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_mono         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_mono        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_mono         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_mono         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_mono        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_mono         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_mono         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_mono        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_farn         =  0.5_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_farn        =  0.17_rk    !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_farn         =  130.0_rk   !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_farn         =  2.37_rk    !empirical coefficient
        real(rk),          parameter     :: caq_farn         =  5.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_farn         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_farn        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_farn         =  5.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_farn         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_farn        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_farn         =  5.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_farn         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_farn        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_farn         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_farn         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_farn        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_cary         =  0.5_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_cary        =  0.17_rk    !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_cary         =  130.0_rk   !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_cary         =  2.37_rk    !empirical coefficient
        real(rk),          parameter     :: caq_cary         =  5.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_cary         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_cary        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_cary         =  5.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_cary         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_cary        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_cary         =  5.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_cary         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_cary        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_cary         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_cary         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_cary        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_sesq         =  0.5_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_sesq        =  0.17_rk    !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_sesq         =  130.0_rk   !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_sesq         =  2.37_rk    !empirical coefficient
        real(rk),          parameter     :: caq_sesq         =  5.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_sesq         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_sesq        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_sesq         =  5.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_sesq         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_sesq        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_sesq         =  5.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_sesq         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_sesq        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_sesq         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_sesq         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_sesq        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_mbol         =  1.0_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_mbol        =  0.13_rk    !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_mbol         =  95.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_mbol         =  2.0_rk     !empirical coefficient
        real(rk),          parameter     :: caq_mbol         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_mbol         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_mbol        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_mbol         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_mbol         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_mbol        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_mbol         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_mbol         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_mbol        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_mbol         =  1.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_mbol         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_mbol        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_meth         =  0.8_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_meth        =  0.08_rk    !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_meth         =  60.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_meth         =  1.6_rk     !empirical coefficient
        real(rk),          parameter     :: caq_meth         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_meth         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_meth        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_meth         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_meth         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_meth        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_meth         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_meth         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_meth        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_meth         =  1.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_meth         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_meth        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_acet         =  0.2_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_acet        =  0.1_rk     !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_acet         =  80.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_acet         =  1.83_rk    !empirical coefficient
        real(rk),          parameter     :: caq_acet         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_acet         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_acet        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_acet         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_acet         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_acet        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_acet         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_acet         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_acet        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_acet         =  1.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_acet         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_acet        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_co           =  1.0_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_co          =  0.08_rk    !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_co           =  60.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_co           =  1.6_rk     !empirical coefficient
        real(rk),          parameter     :: caq_co           =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_co           =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_co          =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_co           =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_co           =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_co          =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_co           =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_co           =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_co          =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_co           =  1.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_co           =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_co          =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_bvoc         =  0.8_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_bvoc        =  0.13_rk    !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_bvoc         =  95.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_bvoc         =  2.0_rk     !empirical coefficient
        real(rk),          parameter     :: caq_bvoc         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_bvoc         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_bvoc        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_bvoc         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_bvoc         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_bvoc        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_bvoc         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_bvoc         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_bvoc        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_bvoc         =  1.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_bvoc         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_bvoc        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_svoc         =  0.8_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_svoc        =  0.1_rk     !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_svoc         =  80.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_svoc         =  1.83_rk    !empirical coefficient
        real(rk),          parameter     :: caq_svoc         =  5.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_svoc         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_svoc        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_svoc         =  5.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_svoc         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_svoc        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_svoc         =  5.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_svoc         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_svoc        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_svoc         =  5.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_svoc         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_svoc        =  8.0_rk     !delta threshold for high wind stress (m/s)

        real(rk),          parameter     :: ldf_ovoc         =  0.2_rk     !light-dependent fraction
        real(rk),          parameter     :: beta_ovoc        =  0.1_rk     !empirical coefficient for temperature dependence of light-independent fraction
        real(rk),          parameter     :: ct1_ovoc         =  80.0_rk    !activation energy (kj/mol)
        real(rk),          parameter     :: ceo_ovoc         =  1.83_rk    !empirical coefficient
        real(rk),          parameter     :: caq_ovoc         =  1.0_rk     !coefficient for poor air quality stress
        real(rk),          parameter     :: taq_ovoc         =  20.0_rk    !threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: dtaq_ovoc        =  30.0_rk    !delta threshold for poor air quality stress (ppm-hours)
        real(rk),          parameter     :: cht_ovoc         =  1.0_rk     !coefficient for high temperature stress
        real(rk),          parameter     :: tht_ovoc         =  313.15_rk  !threshold for high temperature stress (k)
        real(rk),          parameter     :: dtht_ovoc        =  8.0_rk     !delta threshold for high temperature stress (k)
        real(rk),          parameter     :: clt_ovoc         =  1.0_rk     !coefficient for low temperature stress
        real(rk),          parameter     :: tlt_ovoc         =  283.15_rk  !threshold for low temperature stress (k)
        real(rk),          parameter     :: dtlt_ovoc        =  8.0_rk     !delta threshold for low temperature stress (k)
        real(rk),          parameter     :: chw_ovoc         =  1.0_rk     !coefficient for high wind stress
        real(rk),          parameter     :: thw_ovoc         =  12.0_rk    !threshold for high wind stress (m/s)
        real(rk),          parameter     :: dthw_ovoc        =  8.0_rk     !delta threshold for high wind stress (m/s)

! set tree and species dependent coefficients
        if (emi_ind .eq. 1 ) then
            ldf  = ldf_isop
            beta = beta_isop
            ct1  = ct1_isop
            ceo  = ceo_isop
            ef1  = ef1_isop
            ef2  = ef2_isop
            ef3  = ef3_isop
            ef4  = ef4_isop
            ef5  = ef5_isop
            ef6  = ef6_isop
            ef7  = ef7_isop
            ef8  = ef8_isop
            ef9  = ef9_isop
            ef10 = ef10_isop
            ef11 = ef11_isop
            ef12 = ef12_isop
            ef13 = ef13_isop
            ef14 = ef14_isop
            ef15 = ef15_isop
            anew  = anew_isop
            agro  = agro_isop
            amat  = amat_isop
            aold  = aold_isop
            caq   = caq_isop
            taq   = taq_isop
            dtaq  = dtaq_isop
            cht   = cht_isop
            tht   = tht_isop
            dtht  = dtht_isop
            clt   = clt_isop
            tlt   = tlt_isop
            dtlt  = dtlt_isop
            chw   = chw_isop
            thw   = thw_isop
            dthw  = dthw_isop
        else if (emi_ind .eq. 2 ) then
            ldf  = ldf_myrc
            beta = beta_myrc
            ct1 = ct1_myrc
            ceo = ceo_myrc
            ef1  = ef1_myrc
            ef2  = ef2_myrc
            ef3  = ef3_myrc
            ef4  = ef4_myrc
            ef5  = ef5_myrc
            ef6  = ef6_myrc
            ef7  = ef7_myrc
            ef8  = ef8_myrc
            ef9  = ef9_myrc
            ef10 = ef10_myrc
            ef11 = ef11_myrc
            ef12 = ef12_myrc
            ef13 = ef13_myrc
            ef14 = ef14_myrc
            ef15 = ef15_myrc
            anew  = anew_myrc
            agro  = agro_myrc
            amat  = amat_myrc
            aold  = aold_myrc
            caq   = caq_myrc
            taq   = taq_myrc
            dtaq  = dtaq_myrc
            cht   = cht_myrc
            tht   = tht_myrc
            dtht  = dtht_myrc
            clt   = clt_myrc
            tlt   = tlt_myrc
            dtlt  = dtlt_myrc
            chw   = chw_myrc
            thw   = thw_myrc
            dthw  = dthw_myrc
        else if (emi_ind .eq. 3 ) then
            ldf  = ldf_sabi
            beta = beta_sabi
            ct1 = ct1_sabi
            ceo = ceo_sabi
            ef1  = ef1_sabi
            ef2  = ef2_sabi
            ef3  = ef3_sabi
            ef4  = ef4_sabi
            ef5  = ef5_sabi
            ef6  = ef6_sabi
            ef7  = ef7_sabi
            ef8  = ef8_sabi
            ef9  = ef9_sabi
            ef10 = ef10_sabi
            ef11 = ef11_sabi
            ef12 = ef12_sabi
            ef13 = ef13_sabi
            ef14 = ef14_sabi
            ef15 = ef15_sabi
            anew  = anew_sabi
            agro  = agro_sabi
            amat  = amat_sabi
            aold  = aold_sabi
            caq   = caq_sabi
            taq   = taq_sabi
            dtaq  = dtaq_sabi
            cht   = cht_sabi
            tht   = tht_sabi
            dtht  = dtht_sabi
            clt   = clt_sabi
            tlt   = tlt_sabi
            dtlt  = dtlt_sabi
            chw   = chw_sabi
            thw   = thw_sabi
            dthw  = dthw_sabi
        else if (emi_ind .eq. 4 ) then
            ldf  = ldf_limo
            beta = beta_limo
            ct1 = ct1_limo
            ceo = ceo_limo
            ef1  = ef1_limo
            ef2  = ef2_limo
            ef3  = ef3_limo
            ef4  = ef4_limo
            ef5  = ef5_limo
            ef6  = ef6_limo
            ef7  = ef7_limo
            ef8  = ef8_limo
            ef9  = ef9_limo
            ef10 = ef10_limo
            ef11 = ef11_limo
            ef12 = ef12_limo
            ef13 = ef13_limo
            ef14 = ef14_limo
            ef15 = ef15_limo
            anew  = anew_limo
            agro  = agro_limo
            amat  = amat_limo
            aold  = aold_limo
            caq   = caq_limo
            taq   = taq_limo
            dtaq  = dtaq_limo
            cht   = cht_limo
            tht   = tht_limo
            dtht  = dtht_limo
            clt   = clt_limo
            tlt   = tlt_limo
            dtlt  = dtlt_limo
            chw   = chw_limo
            thw   = thw_limo
            dthw  = dthw_limo
        else if (emi_ind .eq. 5 ) then
            ldf  = ldf_care
            beta = beta_care
            ct1 = ct1_care
            ceo = ceo_care
            ef1  = ef1_care
            ef2  = ef2_care
            ef3  = ef3_care
            ef4  = ef4_care
            ef5  = ef5_care
            ef6  = ef6_care
            ef7  = ef7_care
            ef8  = ef8_care
            ef9  = ef9_care
            ef10 = ef10_care
            ef11 = ef11_care
            ef12 = ef12_care
            ef13 = ef13_care
            ef14 = ef14_care
            ef15 = ef15_care
            anew  = anew_care
            agro  = agro_care
            amat  = amat_care
            aold  = aold_care
            caq   = caq_care
            taq   = taq_care
            dtaq  = dtaq_care
            cht   = cht_care
            tht   = tht_care
            dtht  = dtht_care
            clt   = clt_care
            tlt   = tlt_care
            dtlt  = dtlt_care
            chw   = chw_care
            thw   = thw_care
            dthw  = dthw_care
        else if (emi_ind .eq. 6 ) then
            ldf  = ldf_ocim
            beta = beta_ocim
            ct1 = ct1_ocim
            ceo = ceo_ocim
            ef1  = ef1_ocim
            ef2  = ef2_ocim
            ef3  = ef3_ocim
            ef4  = ef4_ocim
            ef5  = ef5_ocim
            ef6  = ef6_ocim
            ef7  = ef7_ocim
            ef8  = ef8_ocim
            ef9  = ef9_ocim
            ef10 = ef10_ocim
            ef11 = ef11_ocim
            ef12 = ef12_ocim
            ef13 = ef13_ocim
            ef14 = ef14_ocim
            ef15 = ef15_ocim
            anew  = anew_ocim
            agro  = agro_ocim
            amat  = amat_ocim
            aold  = aold_ocim
            caq   = caq_ocim
            taq   = taq_ocim
            dtaq  = dtaq_ocim
            cht   = cht_ocim
            tht   = tht_ocim
            dtht  = dtht_ocim
            clt   = clt_ocim
            tlt   = tlt_ocim
            dtlt  = dtlt_ocim
            chw   = chw_ocim
            thw   = thw_ocim
            dthw  = dthw_ocim
        else if (emi_ind .eq. 7 ) then
            ldf  = ldf_bpin
            beta = beta_bpin
            ct1 = ct1_bpin
            ceo = ceo_bpin
            ef1  = ef1_bpin
            ef2  = ef2_bpin
            ef3  = ef3_bpin
            ef4  = ef4_bpin
            ef5  = ef5_bpin
            ef6  = ef6_bpin
            ef7  = ef7_bpin
            ef8  = ef8_bpin
            ef9  = ef9_bpin
            ef10 = ef10_bpin
            ef11 = ef11_bpin
            ef12 = ef12_bpin
            ef13 = ef13_bpin
            ef14 = ef14_bpin
            ef15 = ef15_bpin
            anew  = anew_bpin
            agro  = agro_bpin
            amat  = amat_bpin
            aold  = aold_bpin
            caq   = caq_bpin
            taq   = taq_bpin
            dtaq  = dtaq_bpin
            cht   = cht_bpin
            tht   = tht_bpin
            dtht  = dtht_bpin
            clt   = clt_bpin
            tlt   = tlt_bpin
            dtlt  = dtlt_bpin
            chw   = chw_bpin
            thw   = thw_bpin
            dthw  = dthw_bpin
        else if (emi_ind .eq. 8 ) then
            ldf  = ldf_apin
            beta = beta_apin
            ct1 = ct1_apin
            ceo = ceo_apin
            ef1  = ef1_apin
            ef2  = ef2_apin
            ef3  = ef3_apin
            ef4  = ef4_apin
            ef5  = ef5_apin
            ef6  = ef6_apin
            ef7  = ef7_apin
            ef8  = ef8_apin
            ef9  = ef9_apin
            ef10 = ef10_apin
            ef11 = ef11_apin
            ef12 = ef12_apin
            ef13 = ef13_apin
            ef14 = ef14_apin
            ef15 = ef15_apin
            anew  = anew_apin
            agro  = agro_apin
            amat  = amat_apin
            aold  = aold_apin
            caq   = caq_apin
            taq   = taq_apin
            dtaq  = dtaq_apin
            cht   = cht_apin
            tht   = tht_apin
            dtht  = dtht_apin
            clt   = clt_apin
            tlt   = tlt_apin
            dtlt  = dtlt_apin
            chw   = chw_apin
            thw   = thw_apin
            dthw  = dthw_apin
        else if (emi_ind .eq. 9 ) then
            ldf  = ldf_mono
            beta = beta_mono
            ct1 = ct1_mono
            ceo = ceo_mono
            ef1  = ef1_mono
            ef2  = ef2_mono
            ef3  = ef3_mono
            ef4  = ef4_mono
            ef5  = ef5_mono
            ef6  = ef6_mono
            ef7  = ef7_mono
            ef8  = ef8_mono
            ef9  = ef9_mono
            ef10 = ef10_mono
            ef11 = ef11_mono
            ef12 = ef12_mono
            ef13 = ef13_mono
            ef14 = ef14_mono
            ef15 = ef15_mono
            anew  = anew_mono
            agro  = agro_mono
            amat  = amat_mono
            aold  = aold_mono
            caq   = caq_mono
            taq   = taq_mono
            dtaq  = dtaq_mono
            cht   = cht_mono
            tht   = tht_mono
            dtht  = dtht_mono
            clt   = clt_mono
            tlt   = tlt_mono
            dtlt  = dtlt_mono
            chw   = chw_mono
            thw   = thw_mono
            dthw  = dthw_mono
        else if (emi_ind .eq. 10 ) then
            ldf  = ldf_farn
            beta = beta_farn
            ct1 = ct1_farn
            ceo = ceo_farn
            ef1  = ef1_farn
            ef2  = ef2_farn
            ef3  = ef3_farn
            ef4  = ef4_farn
            ef5  = ef5_farn
            ef6  = ef6_farn
            ef7  = ef7_farn
            ef8  = ef8_farn
            ef9  = ef9_farn
            ef10 = ef10_farn
            ef11 = ef11_farn
            ef12 = ef12_farn
            ef13 = ef13_farn
            ef14 = ef14_farn
            ef15 = ef15_farn
            anew  = anew_farn
            agro  = agro_farn
            amat  = amat_farn
            aold  = aold_farn
            caq   = caq_farn
            taq   = taq_farn
            dtaq  = dtaq_farn
            cht   = cht_farn
            tht   = tht_farn
            dtht  = dtht_farn
            clt   = clt_farn
            tlt   = tlt_farn
            dtlt  = dtlt_farn
            chw   = chw_farn
            thw   = thw_farn
            dthw  = dthw_farn
        else if (emi_ind .eq. 11 ) then
            ldf  = ldf_cary
            beta = beta_cary
            ct1 = ct1_cary
            ceo = ceo_cary
            ef1  = ef1_cary
            ef2  = ef2_cary
            ef3  = ef3_cary
            ef4  = ef4_cary
            ef5  = ef5_cary
            ef6  = ef6_cary
            ef7  = ef7_cary
            ef8  = ef8_cary
            ef9  = ef9_cary
            ef10 = ef10_cary
            ef11 = ef11_cary
            ef12 = ef12_cary
            ef13 = ef13_cary
            ef14 = ef14_cary
            ef15 = ef15_cary
            anew  = anew_cary
            agro  = agro_cary
            amat  = amat_cary
            aold  = aold_cary
            caq   = caq_cary
            taq   = taq_cary
            dtaq  = dtaq_cary
            cht   = cht_cary
            tht   = tht_cary
            dtht  = dtht_cary
            clt   = clt_cary
            tlt   = tlt_cary
            dtlt  = dtlt_cary
            chw   = chw_cary
            thw   = thw_cary
            dthw  = dthw_cary
        else if (emi_ind .eq. 12 ) then
            ldf  = ldf_sesq
            beta = beta_sesq
            ct1 = ct1_sesq
            ceo = ceo_sesq
            ef1  = ef1_sesq
            ef2  = ef2_sesq
            ef3  = ef3_sesq
            ef4  = ef4_sesq
            ef5  = ef5_sesq
            ef6  = ef6_sesq
            ef7  = ef7_sesq
            ef8  = ef8_sesq
            ef9  = ef9_sesq
            ef10 = ef10_sesq
            ef11 = ef11_sesq
            ef12 = ef12_sesq
            ef13 = ef13_sesq
            ef14 = ef14_sesq
            ef15 = ef15_sesq
            anew  = anew_sesq
            agro  = agro_sesq
            amat  = amat_sesq
            aold  = aold_sesq
            caq   = caq_sesq
            taq   = taq_sesq
            dtaq  = dtaq_sesq
            cht   = cht_sesq
            tht   = tht_sesq
            dtht  = dtht_sesq
            clt   = clt_sesq
            tlt   = tlt_sesq
            dtlt  = dtlt_sesq
            chw   = chw_sesq
            thw   = thw_sesq
            dthw  = dthw_sesq
        else if (emi_ind .eq. 13 ) then
            ldf  = ldf_mbol
            beta = beta_mbol
            ct1 = ct1_mbol
            ceo = ceo_mbol
            ef1  = ef1_mbol
            ef2  = ef2_mbol
            ef3  = ef3_mbol
            ef4  = ef4_mbol
            ef5  = ef5_mbol
            ef6  = ef6_mbol
            ef7  = ef7_mbol
            ef8  = ef8_mbol
            ef9  = ef9_mbol
            ef10 = ef10_mbol
            ef11 = ef11_mbol
            ef12 = ef12_mbol
            ef13 = ef13_mbol
            ef14 = ef14_mbol
            ef15 = ef15_mbol
            anew  = anew_mbol
            agro  = agro_mbol
            amat  = amat_mbol
            aold  = aold_mbol
            caq   = caq_mbol
            taq   = taq_mbol
            dtaq  = dtaq_mbol
            cht   = cht_mbol
            tht   = tht_mbol
            dtht  = dtht_mbol
            clt   = clt_mbol
            tlt   = tlt_mbol
            dtlt  = dtlt_mbol
            chw   = chw_mbol
            thw   = thw_mbol
            dthw  = dthw_mbol
        else if (emi_ind .eq. 14 ) then
            ldf  = ldf_meth
            beta = beta_meth
            ct1 = ct1_meth
            ceo = ceo_meth
            ef1  = ef1_meth
            ef2  = ef2_meth
            ef3  = ef3_meth
            ef4  = ef4_meth
            ef5  = ef5_meth
            ef6  = ef6_meth
            ef7  = ef7_meth
            ef8  = ef8_meth
            ef9  = ef9_meth
            ef10 = ef10_meth
            ef11 = ef11_meth
            ef12 = ef12_meth
            ef13 = ef13_meth
            ef14 = ef14_meth
            ef15 = ef15_meth
            anew  = anew_meth
            agro  = agro_meth
            amat  = amat_meth
            aold  = aold_meth
            caq   = caq_meth
            taq   = taq_meth
            dtaq  = dtaq_meth
            cht   = cht_meth
            tht   = tht_meth
            dtht  = dtht_meth
            clt   = clt_meth
            tlt   = tlt_meth
            dtlt  = dtlt_meth
            chw   = chw_meth
            thw   = thw_meth
            dthw  = dthw_meth
        else if (emi_ind .eq. 15 ) then
            ldf  = ldf_acet
            beta = beta_acet
            ct1 = ct1_acet
            ceo = ceo_acet
            ef1  = ef1_acet
            ef2  = ef2_acet
            ef3  = ef3_acet
            ef4  = ef4_acet
            ef5  = ef5_acet
            ef6  = ef6_acet
            ef7  = ef7_acet
            ef8  = ef8_acet
            ef9  = ef9_acet
            ef10 = ef10_acet
            ef11 = ef11_acet
            ef12 = ef12_acet
            ef13 = ef13_acet
            ef14 = ef14_acet
            ef15 = ef15_acet
            anew  = anew_acet
            agro  = agro_acet
            amat  = amat_acet
            aold  = aold_acet
            caq   = caq_acet
            taq   = taq_acet
            dtaq  = dtaq_acet
            cht   = cht_acet
            tht   = tht_acet
            dtht  = dtht_acet
            clt   = clt_acet
            tlt   = tlt_acet
            dtlt  = dtlt_acet
            chw   = chw_acet
            thw   = thw_acet
            dthw  = dthw_acet
        else if (emi_ind .eq. 16 ) then
            ldf  = ldf_co
            beta = beta_co
            ct1 = ct1_co
            ceo = ceo_co
            ef1  = ef1_co
            ef2  = ef2_co
            ef3  = ef3_co
            ef4  = ef4_co
            ef5  = ef5_co
            ef6  = ef6_co
            ef7  = ef7_co
            ef8  = ef8_co
            ef9  = ef9_co
            ef10 = ef10_co
            ef11 = ef11_co
            ef12 = ef12_co
            ef13 = ef13_co
            ef14 = ef14_co
            ef15 = ef15_co
            anew  = anew_co
            agro  = agro_co
            amat  = amat_co
            aold  = aold_co
            caq   = caq_co
            taq   = taq_co
            dtaq  = dtaq_co
            cht   = cht_co
            tht   = tht_co
            dtht  = dtht_co
            clt   = clt_co
            tlt   = tlt_co
            dtlt  = dtlt_co
            chw   = chw_co
            thw   = thw_co
            dthw  = dthw_co
        else if (emi_ind .eq. 17 ) then
            ldf  = ldf_bvoc
            beta = beta_bvoc
            ct1 = ct1_bvoc
            ceo = ceo_bvoc
            ef1  = ef1_bvoc
            ef2  = ef2_bvoc
            ef3  = ef3_bvoc
            ef4  = ef4_bvoc
            ef5  = ef5_bvoc
            ef6  = ef6_bvoc
            ef7  = ef7_bvoc
            ef8  = ef8_bvoc
            ef9  = ef9_bvoc
            ef10 = ef10_bvoc
            ef11 = ef11_bvoc
            ef12 = ef12_bvoc
            ef13 = ef13_bvoc
            ef14 = ef14_bvoc
            ef15 = ef15_bvoc
            anew  = anew_bvoc
            agro  = agro_bvoc
            amat  = amat_bvoc
            aold  = aold_bvoc
            caq   = caq_bvoc
            taq   = taq_bvoc
            dtaq  = dtaq_bvoc
            cht   = cht_bvoc
            tht   = tht_bvoc
            dtht  = dtht_bvoc
            clt   = clt_bvoc
            tlt   = tlt_bvoc
            dtlt  = dtlt_bvoc
            chw   = chw_bvoc
            thw   = thw_bvoc
            dthw  = dthw_bvoc
        else if (emi_ind .eq. 18 ) then
            ldf  = ldf_svoc
            beta = beta_svoc
            ct1 = ct1_svoc
            ceo = ceo_svoc
            ef1  = ef1_svoc
            ef2  = ef2_svoc
            ef3  = ef3_svoc
            ef4  = ef4_svoc
            ef5  = ef5_svoc
            ef6  = ef6_svoc
            ef7  = ef7_svoc
            ef8  = ef8_svoc
            ef9  = ef9_svoc
            ef10 = ef10_svoc
            ef11 = ef11_svoc
            ef12 = ef12_svoc
            ef13 = ef13_svoc
            ef14 = ef14_svoc
            ef15 = ef15_svoc
            anew  = anew_svoc
            agro  = agro_svoc
            amat  = amat_svoc
            aold  = aold_svoc
            caq   = caq_svoc
            taq   = taq_svoc
            dtaq  = dtaq_svoc
            cht   = cht_svoc
            tht   = tht_svoc
            dtht  = dtht_svoc
            clt   = clt_svoc
            tlt   = tlt_svoc
            dtlt  = dtlt_svoc
            chw   = chw_svoc
            thw   = thw_svoc
            dthw  = dthw_svoc
        else   ! emi_ind = 19
            ldf  = ldf_ovoc
            beta = beta_ovoc
            ct1 = ct1_ovoc
            ceo = ceo_ovoc
            ef1  = ef1_ovoc
            ef2  = ef2_ovoc
            ef3  = ef3_ovoc
            ef4  = ef4_ovoc
            ef5  = ef5_ovoc
            ef6  = ef6_ovoc
            ef7  = ef7_ovoc
            ef8  = ef8_ovoc
            ef9  = ef9_ovoc
            ef10 = ef10_ovoc
            ef11 = ef11_ovoc
            ef12 = ef12_ovoc
            ef13 = ef13_ovoc
            ef14 = ef14_ovoc
            ef15 = ef15_ovoc
            anew  = anew_ovoc
            agro  = agro_ovoc
            amat  = amat_ovoc
            aold  = aold_ovoc
            caq   = caq_ovoc
            taq   = taq_ovoc
            dtaq  = dtaq_ovoc
            cht   = cht_ovoc
            tht   = tht_ovoc
            dtht  = dtht_ovoc
            clt   = clt_ovoc
            tlt   = tlt_ovoc
            dtlt  = dtlt_ovoc
            chw   = chw_ovoc
            thw   = thw_ovoc
            dthw  = dthw_ovoc
        end if

        if (lu_opt .eq. 0 .or. lu_opt .eq. 1) then !viirs or modis  lu types

! simple megan (table 3 in guenther et al., 2012) pft to viirs/modis vtype mapping
            if (vtype .eq. 1) then !viirs cat 1 evergreen needleleaf
                !--> average needleleaf evergreen temperate tree and needleleaf evergreen boreal tree

                ef = (ef1+ef2)/2.0_rk
                !set pft dependent a and b coefficients for cumulative root depth fraction (zeng 2001)
                !see table 2 for igpb classification at:  https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml
                roota = 6.706_rk
                rootb = 2.175_rk

            else if (vtype .eq. 2) then !viirs/modis cat 2 evergreen broadleaf
                !--> average broadleaf evergreen tropical tree and broadleaf evergreen temperate tree

                ef = (ef4+ef5)/2.0_rk
                !set pft dependent a and b coefficients for cumulative root depth fraction (zeng 2001)
                !see table 2 for igpb classification at:  https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml
                roota = 7.344_rk
                rootb = 1.303_rk

            else if (vtype .eq. 3) then !viirs/modis cat 3 deciduous needleaf
                !--> average needleleaf deciduous boreal tree

                ef = ef3
                !set pft dependent a and b coefficients for cumulative root depth fraction (zeng 2001)
                !see table 2 for igpb classification at:  https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml
                roota = 7.066_rk
                rootb = 1.953_rk

            else if (vtype .eq. 4) then !viirs/modis cat 4 deciduous broadleaf
                !--> average broadleaf deciduous tropical tree, broadleaf deciduous temperate tree,
                ! and broadleaf deciduous boreal tree

                ef = (ef6+ef7+ef8)/3.0_rk
                !set pft dependent a and b coefficients for cumulative root depth fraction (zeng 2001)
                !see table 2 for igpb classification at:  https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml
                roota = 5.990_rk
                rootb = 1.955_rk

            else if (vtype .eq. 5) then !viirs/modis cat 5 mixed forests
                !--> avearge of all above ef1-ef8 pfts.

                ef = (ef1+ef2+ef3+ef4+ef5+ef6+ef7+ef8)/8.0_rk
                !set pft dependent a and b coefficients for cumulative root depth fraction (zeng 2001)
                !see table 2 for igpb classification at:  https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml
                roota = 4.453_rk
                rootb = 1.631_rk

            else if (vtype .ge. 6 .and. vtype .le. 7) then !viirs/modis cat 6-7 closed/open shrublands
                !--> avearge broadleaf evergreen temperate shrub, broadleaf deciduous temperate shrub,
                ! and broadleaf deciduous boreal shrub

                ef = (ef9+ef10+ef11)/3.0_rk
                !set pft dependent a and b coefficients for cumulative root depth fraction (zeng 2001)
                !see table 2 for igpb classification at:  https://journals.ametsoc.org/view/journals/hydr/2/5/1525-7541_2001_002_0525_gvrdfl_2_0_co_2.xml
                roota = (6.326_rk + 7.718_rk)/2.0_rk
                rootb = (1.567_rk + 1.262_rk)/2.0_rk

            else if (vtype .ge. 8 .and. vtype .le. 11) then !viirs/modis cat 8-10 savannas and grasslands
                !--> avearge arctic c3 grass, cool c3 grass, warm c4 grass), intent(out) ef;

variable out

real(rk), intent(out) out;