Turn DustDead emissions on by default and set OFFLINE_DUST to false

Randall Martin advised not using offline dust emissions until further notice, as they are buggy. They will be regenerated following Meng et al. (2021, submitted). In the meantime, online dust emissions will be used by default. For GEOS-Chem Classic, mass tuning factors will be modified by createRunDir.sh based on met field and resolution. For GCHP, the mass tuning factors will be modified by runConfig.sh. Signed-off-by: Melissa Sulprizio <mpayer@seas.harvard.edu>
geoschem · Mar 23, 2021 · eb97d0b · eb97d0b
1 parent f679ef9
commit eb97d0b
Show file tree

Hide file tree

Showing 4 changed files with 32 additions and 27 deletions.
diff --git a/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.aerosol b/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.aerosol
@@ -92,15 +92,16 @@ Warnings:                    1
     --> QFED2                  :       false    # 2000-2020
     --> GFAS                   :       false    # 2003-2021
 # ----- OFFLINE EMISSIONS -------------------------------------
-# To use online emissions instead set the
-# offline emissions to 'false' and the
-# corresponding HEMCO extension to 'on':
+# To use online emissions instead set the offline emissions to
+# 'false' and the corresponding HEMCO extension to 'on':
 #   OFFLINE_DUST        - DustDead or DustGinoux
 #   OFFLINE_BIOGENICVOC - MEGAN
 #   OFFLINE_SEASALT     - SeaSalt
 #   OFFLINE_SOILNOX     - SoilNOx
+# NOTE: The offline dust emissions are considered to be buggy.
+#  Until further notice online dust emissions are recommended.
 #--------------------------------------------------------------
-    --> OFFLINE_DUST           :       true     # 1980-2021
+    --> OFFLINE_DUST           :       false    # 1980-2021
     --> OFFLINE_BIOGENICVOC    :       true     # 1980-2021
     --> OFFLINE_SEASALT        :       true     # 1980-2021
     --> OFFLINE_SOILNOX        :       true     # 1980-2021
@@ -114,7 +115,7 @@ Warnings:                    1
 # -----------------------------------------------------------------------------
 100     Custom                 : off   -
 101     SeaFlux                : on    DMS
-105     DustDead               : off   DST1/DST2/DST3/DST4
+105     DustDead               : on    DST1/DST2/DST3/DST4
     --> Mass tuning factor     :       -999.0e0
 106     DustGinoux             : off   DST1/DST2/DST3/DST4
 107     SeaSalt                : off   SALA/SALC/SALACL/SALCCL/SALAAL/SALCAL/BrSALA/BrSALC/MOPO/MOPI

diff --git a/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.fullchem b/run/GCClassic/HEMCO_Config.rc.templates/HEMCO_Config.rc.fullchem
@@ -96,17 +96,18 @@ Warnings:                    1
     --> QFED2                  :       false    # 2000-2020
     --> GFAS                   :       false    # 2003-2021
 # ----- OFFLINE EMISSIONS -------------------------------------
-# To use online emissions instead set the
-# offline emissions to 'false' and the
-# corresponding HEMCO extension to 'on':
+# To use online emissions instead set the offline emissions to
+# 'false' and the corresponding HEMCO extension to 'on':
 #   OFFLINE_DUST        - DustDead or DustGinoux
 #   OFFLINE_BIOGENICVOC - MEGAN
 #    - Note: MEGAN must always be on for certain
 #            species not in OFFLINE_BIOGENICVOC
 #   OFFLINE_SEASALT     - SeaSalt
 #   OFFLINE_SOILNOX     - SoilNOx
+# NOTE: The offline dust emissions are considered to be buggy.
+#  Until further notice online dust emissions are recommended.
 #--------------------------------------------------------------
-    --> OFFLINE_DUST           :       true     # 1980-2021
+    --> OFFLINE_DUST           :       false    # 1980-2021
     --> OFFLINE_BIOGENICVOC    :       true     # 1980-2021
     --> OFFLINE_SEASALT        :       true     # 1980-2021
     --> OFFLINE_SOILNOX        :       true     # 1980-2021
@@ -134,7 +135,7 @@ Warnings:                    1
     --> CDF table              :       $ROOT/LIGHTNOX/v2014-07/light_dist.ott2010.dat
 104     SoilNOx                : off   NO
     --> Use fertilizer NOx     :       true
-105     DustDead               : off   DST1/DST2/DST3/DST4
+105     DustDead               : on    DST1/DST2/DST3/DST4
     --> Mass tuning factor     :       -999.0e0
 106     DustGinoux             : off   DST1/DST2/DST3/DST4
 107     SeaSalt                : off   SALA/SALC/SALACL/SALCCL/SALAAL/SALCAL/BrSALA/BrSALC/MOPO/MOPI

diff --git a/run/GCClassic/createRunDir.sh b/run/GCClassic/createRunDir.sh
@@ -278,7 +278,7 @@ while [ "${valid_met}" -eq 0 ]; do
 	met_cn_year='2015'
 	pressure_unit='Pa '
 	pressure_scale='0.01'
-	dust_sf='3.86e-4'
+	offline_dust_sf='3.86e-4'
     elif [[ ${met_num} = "2" ]]; then
 	met_name='GEOSFP'
 	met_name_lc="geosfp"
@@ -291,7 +291,7 @@ while [ "${valid_met}" -eq 0 ]; do
 	met_cn_year='2011'
 	pressure_unit='hPa'
 	pressure_scale='1.0 '
-	dust_sf='6.42e-5'
+	offline_dust_sf='6.42e-5'
     else
 	valid_met=0
 	printf "Invalid meteorology option. Try again.\n"
@@ -591,7 +591,7 @@ sed_ie "s|{MET_DIR}|${met_dir}|"          HEMCO_Config.rc
 sed_ie "s|{NATIVE_RES}|${met_native}|"    HEMCO_Config.rc
 sed_ie "s|{LATRES}|${met_latres}|"        HEMCO_Config.rc
 sed_ie "s|{LONRES}|${met_lonres}|"        HEMCO_Config.rc
-sed_ie "s|{DUST_SF}|${dust_sf}|"          HEMCO_Config.rc
+sed_ie "s|{DUST_SF}|${offline_dust_sf}|"  HEMCO_Config.rc
 
 # Special handling for start/end date based on simulation so that
 # start year/month/day matches default initial restart file.
@@ -630,26 +630,29 @@ fi
 # Modify input files for benchmark that are specific to GEOS-Chem Classic
 if [[ "x${sim_extra_option}" == "xbenchmark" ]]; then
     replace_colon_sep_val "Use GC classic timers?"   T    input.geos
-    if [[ "x${met_name}" == "xGEOSFP" && "x${grid_res}" == "x4x5" ]]; then
-	replace_colon_sep_val "--> Mass tuning factor" 8.3286e-4 HEMCO_Config.rc
-    fi
-    if [[ "x${met_name}" == "xMERRA2" && "x${grid_res}" == "x4x5" ]]; then
-	replace_colon_sep_val "--> Mass tuning factor" 7.8533e-4 HEMCO_Config.rc
-    fi
 fi
 
 # Modify input files for TOMAS that are specific to GEOS-Chem Classic
-# NOTE: Also use the same dust tuning factors as for the benchmark
-# for the Tomas_DustDead extension (mps, bmy, 3/11/12)
 if [[ ${sim_extra_option} =~ "TOMAS" ]]; then
     replace_colon_sep_val "Tran/conv timestep [sec]" 1800 input.geos
     replace_colon_sep_val "Chem/emis timestep [sec]" 3600 input.geos
+fi
+
+# Set online dust emission mass tuning factor according to met field
+# and resolution. These values were obtained from hcox_dustdead_mod.F90.
+if [[ "x${sim_name}" == "xfullchem" || "x${sim_name}" == "xaerosol" ]]; then
     if [[ "x${met_name}" == "xGEOSFP" && "x${grid_res}" == "x4x5" ]]; then
 	replace_colon_sep_val "--> Mass tuning factor" 8.3286e-4 HEMCO_Config.rc
     fi
+    if [[ "x${met_name}" == "xGEOSFP" && "x${grid_res}" == "x2x25" ]]; then
+	replace_colon_sep_val "--> Mass tuning factor" 5.0416e-4 HEMCO_Config.rc
+    fi
     if [[ "x${met_name}" == "xMERRA2" && "x${grid_res}" == "x4x5" ]]; then
 	replace_colon_sep_val "--> Mass tuning factor" 7.8533e-4 HEMCO_Config.rc
     fi
+    if [[ "x${met_name}" == "xMERRA2" && "x${grid_res}" == "x2x25" ]]; then
+	replace_colon_sep_val "--> Mass tuning factor" 4.7586e-4 HEMCO_Config.rc
+    fi
 fi
 
 # Modify input files for troposphere-only chemistry grids

diff --git a/run/GCHP/HEMCO_Config.rc.templates/HEMCO_Config.rc.fullchem b/run/GCHP/HEMCO_Config.rc.templates/HEMCO_Config.rc.fullchem
@@ -96,17 +96,18 @@ Warnings:                    1
     --> QFED2                  :       false    # 2000-2020
     --> GFAS                   :       false    # 2003-2021
 # ----- OFFLINE EMISSIONS -------------------------------------
-# To use online emissions instead set the
-# offline emissions to 'false' and the
-# corresponding HEMCO extension to 'on':
+# To use online emissions instead set the offline emissions to
+# 'false' and the corresponding HEMCO extension to 'on':
 #   OFFLINE_DUST        - DustDead or DustGinoux
 #   OFFLINE_BIOGENICVOC - MEGAN
 #    - Note: MEGAN must always be on for certain
 #            species not in OFFLINE_BIOGENICVOC
 #   OFFLINE_SEASALT     - SeaSalt
 #   OFFLINE_SOILNOX     - SoilNOx
+# NOTE: The offline dust emissions are considered to be buggy.
+#  Until further notice online dust emissions are recommended.
 #--------------------------------------------------------------
-    --> OFFLINE_DUST           :       true     # 1980-2021
+    --> OFFLINE_DUST           :       false   # 1980-2021
     --> OFFLINE_BIOGENICVOC    :       true     # 1980-2021
     --> OFFLINE_SEASALT        :       true     # 1980-2021
     --> OFFLINE_SOILNOX        :       true     # 1980-2021
@@ -134,8 +135,7 @@ Warnings:                    1
     --> CDF table              :       $ROOT/LIGHTNOX/v2014-07/light_dist.ott2010.dat
 104     SoilNOx                : off   NO
     --> Use fertilizer NOx     :       true
-105     DustDead               : off   DST1/DST2/DST3/DST4
-# Mass tuning factor auto-updated with resolution-dependent value set in runConfig.sh
+105     DustDead               : on    DST1/DST2/DST3/DST4
     --> Mass tuning factor     :       -999.0e0
 106     DustGinoux             : off   DST1/DST2/DST3/DST4
 107     SeaSalt                : off   SALA/SALC/SALACL/SALCCL/SALAAL/SALCAL/BrSALA/BrSALC/MOPO/MOPI