Adding Gamma AQ, HT, LT, and HW stress factors to biogenics

README.md

@@ -284,11 +284,11 @@ Otherwise, please contact [email protected] for other GFSv16 data peri
 | `hist_opt`      | user-set option to use historically averaged short-term (24-hr) and long-term (240-hr) rolling averages for leaf temperature and PAR for biogenic emissions  (default is off i.e., `hist_opt=0`)  Note: If simulation is </= 24 hours, instantaneous values for leaf temperature and PAR will be used even if historical averaging is turned on (i.e., `hist_opt=1`).  Recommend turning on `hist_opt=1` and running at least for 25 hours to create a model spin-up, and use subsequent simulation hours for analysis. Overall a 10-day (240 hr) spinup is optimal for best analysis of biogenic emissions. |
 | `soim_opt`   | user-set options for applying soil moisture response to biogenic VOC emissions based on [Guenther et al. 2006](https://doi.org/10.5194/acp-6-3181-2006), which depends on input soil moisture at depth and the wilting point.  This includes additional PFT dependent approach for cumulative root fraction within each soil layer from [Zeng (2001)](https://doi.org/10.1175/1525-7541(2001)002<0525:GVRDFL>2.0.CO;2). (default is off i.e., `soim_opt=1`, the corresponding $\gamma$ is set to 1). If turned on (`soim_opt=0`), which is recommended, soim $\gamma$ is calculated and the prescribed 4-layer soil depths (`soild[1-4]` below) are used.  Four layers are assumed, and are based on GFS Noah and Noah-MP LSM. |
 | `soild[1-4]` | user-set real values of four level soil depths at centerpoint (cm).  Four layers are based on the GFS Noah and Noah-MP LSM, default values are `soild1=5.0`, `soild2=25.0`, `soild3=70.0`, and  `soild4=150.0`. |
-| `aq_opt`       | user-set options for applying an air quality stress factor for biogenic emissions using calculated, spatially-dependent and global GFS-based ozone W126 values (= `0`, default) or a constant user-set W126 value (= `1`). To turn off aq stress factor set `aq_opt=2`  Note:  The aq_opt should only be turned on during simulations across respective ozone season for specific region (e.g., April-October in the U.S.)|
+| `aq_opt`       | user-set options for applying an air quality stress factor for biogenic emissions using calculated, spatially-dependent and global GFS-based ozone W126 values (= `0`) or a constant user-set W126 value (= `1`). To turn off aq stress factor set `aq_opt=2` (set as default, Off).  Note:  The aq_opt should only be turned on during simulations across respective ozone season for specific region (e.g., April-October in the U.S.)|
 | `w126_set`       | user-set real value of constant ozone W126 values (ppm-hours) (only used if `aq_opt=1`) |
-| `ht_opt`       | user-set options for applying a daily high temperature stress factor for biogenic emissions using daily maximum 2-meter input temperature (= `0`, default).  This is based on MEGAN3 and it is recommended that this option is only used when turning the historical option on (i.e., `hist_opt=1`) and running longer than 1-day simulations to obtain the daily max.   To turn off ht stress factor set `ht_opt=1`  |
-| `lt_opt`       | user-set options for applying a daily low temperature stress factor for biogenic emissions using daily minimum 2-meter input temperature (= `0`, default).  This is based on MEGAN3 and it is recommended that this option is only used when turning the historical option on (i.e., `hist_opt=1`) and running longer than 1-day simulations to obtain the daily min.   To turn off lt stress factor set `lt_opt=1`  |
-| `hw_opt`       | user-set options for applying a daily high wind speed stress factor for biogenic emissions using daily maximum 10-meter input wind speed (= `0`, default).  This is based on MEGAN3 and it is recommended that this option is only used when turning the historical option on (i.e., `hist_opt=1`) and running longer than 1-day simulations to obtain the daily max.   To turn off hw stress factor set `hw_opt=1`  |
+| `ht_opt`       | user-set options for applying a daily high temperature stress factor for biogenic emissions using daily maximum 2-meter input temperature (= `0`).  This is based on MEGAN3 and it is recommended that this option is only used when turning the historical option on (i.e., `hist_opt=1`) and running longer than 1-day simulations to obtain the daily max.   To turn off ht stress factor set `ht_opt=1` (set as default, Off)  |
+| `lt_opt`       | user-set options for applying a daily low temperature stress factor for biogenic emissions using daily minimum 2-meter input temperature (= `0`).  This is based on MEGAN3 and it is recommended that this option is only used when turning the historical option on (i.e., `hist_opt=1`) and running longer than 1-day simulations to obtain the daily min.   To turn off lt stress factor set `lt_opt=1` (set as default, Off) |
+| `hw_opt`       | user-set options for applying a daily high wind speed stress factor for biogenic emissions using daily maximum 10-meter input wind speed (= `0`).  This is based on MEGAN3 and it is recommended that this option is only used when turning the historical option on (i.e., `hist_opt=1`) and running longer than 1-day simulations to obtain the daily max.   To turn off hw stress factor set `hw_opt=1` (set as default, Off) |
 
 
 **\*\*** If `modres` >> `flameh` then some error in WAF calculation will be incurred.  Suggestion is to use relative fine `modres` (at least <= 0.5 m) compared to average flame heights (e.g., ~ 1.0 m) if WAF is required.

input/namelist.canopy

@@ -83,10 +83,10 @@
   soild2      =  25.0
   soild3      =  70.0
   soild4      =  150.0
-  aq_opt      =  0
+  aq_opt      =  2
   w126_set    =  30.0
   hist_opt    =  1
-  ht_opt      =  0
-  lt_opt      =  0
-  hw_opt      =  0
+  ht_opt      =  1
+  lt_opt      =  1
+  hw_opt      =  1
 /

src/canopy_readnml.F90

@@ -344,8 +344,8 @@ SUBROUTINE canopy_readnml
 !-------------------------------------------------------------------------------
 
 !------------------------------------------------------------------------------
-! Set default integer for using air quality stress gamma for biogenic emissions (default=0; On)
-    aq_opt = 0
+! Set default integer for using air quality stress gamma for biogenic emissions (default=1; Off)
+    aq_opt = 2
 !-------------------------------------------------------------------------------
 
 !-------------------------------------------------------------------------------
@@ -354,18 +354,18 @@ SUBROUTINE canopy_readnml
 !-------------------------------------------------------------------------------
 
 !------------------------------------------------------------------------------
-! Set default integer for using high temperature stress gamma for biogenic emissions (default=0; On)
-    ht_opt = 0
+! Set default integer for using high temperature stress gamma for biogenic emissions (default=1; Off)
+    ht_opt = 1
 !-------------------------------------------------------------------------------
 
 !------------------------------------------------------------------------------
-! Set default integer for using low temperature stress gamma for biogenic emissions (default=0; On)
-    lt_opt = 0
+! Set default integer for using low temperature stress gamma for biogenic emissions (default=1; Off)
+    lt_opt = 1
 !-------------------------------------------------------------------------------
 
 !------------------------------------------------------------------------------
-! Set default integer for using high wind stress gamma for biogenic emissions (default=0; On)
-    hw_opt = 0
+! Set default integer for using high wind stress gamma for biogenic emissions (default=1; Off)
+    hw_opt = 1
 !-------------------------------------------------------------------------------
 
 !-------------------------------------------------------------------------------