Group 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 |
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;