From 0e009ddf1ba05718b50570d23ec8eebc4088cf03 Mon Sep 17 00:00:00 2001
From: Yuanjian Zhang <yuanjian.z@wustl.edu>
Date: Thu, 10 Oct 2024 23:15:20 -0500
Subject: [PATCH 1/3] Update dust scale factors for MERRA-2 GEOS-FP and GEOS-IT

A companion update to HEMCO pull request #279: use meterology Ustar for dust emission
---
 run/GCClassic/createRunDir.sh | 12 ++++++------
 1 file changed, 6 insertions(+), 6 deletions(-)

diff --git a/run/GCClassic/createRunDir.sh b/run/GCClassic/createRunDir.sh
index 39d34e84b..f385e193a 100755
--- a/run/GCClassic/createRunDir.sh
+++ b/run/GCClassic/createRunDir.sh
@@ -688,17 +688,17 @@ else
     # Use GEOS-FP values as placeholders for GEOS-IT until parameters derived
     if [[ "x${sim_name}" == "xfullchem" || "x${sim_name}" == "xaerosol" ]]; then
 	if [[ "x${met}" == "xgeosfp" && "x${grid_res}" == "x4x5" ]]; then
-	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='8.3286e-4'\n"
+	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='4.8632e-5'\n"
 	elif [[ "x${met}" == "xgeosfp" && "x${grid_res}" == "x2x25" ]]; then
-	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='5.0416e-4'\n"
+	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='3.8197e-5'\n"
 	elif [[ "x${met}" == "xmerra2" && "x${grid_res}" == "x4x5" ]]; then
-	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='7.8533e-4'\n"
+	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='5.6659e-5'\n"
 	elif [[ "x${met}" == "xmerra2" && "x${grid_res}" == "x2x25" ]]; then
-	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='4.7586e-4'\n"
+	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='4.5412e-5'\n"
 	elif [[ "x${met}" == "xgeosit" && "x${grid_res}" == "x4x5" ]]; then
-	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='8.3286e-4'\n"
+	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='3.8656e-5'\n"
 	elif [[ "x${met}" == "xgeosit" && "x${grid_res}" == "x2x25" ]]; then
-	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='5.0416e-4'\n"
+	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='3.1132e-5'\n"
 	else
 	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='-999.0e0'\n"
 	fi

From 7a47b9a9116bba6ba34df59db3cae00461975d5b Mon Sep 17 00:00:00 2001
From: Yuanjian Zhang <yuanjian.z@wustl.edu>
Date: Thu, 24 Oct 2024 15:05:25 -0500
Subject: [PATCH 2/3] Update CHANGELOG.md

---
 CHANGELOG.md | 3 +++
 1 file changed, 3 insertions(+)

diff --git a/CHANGELOG.md b/CHANGELOG.md
index 860e97617..efcec9a89 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -1,3 +1,6 @@
+### Added
+- Added dust scale factors for MERRA-2, GEOS-IT, and GEOS-FP when using USTAR for Dust DEAD extension
+
 # Changelog
 
 This file documents all notable changes to the GEOS-Chem repository starting in version 14.0.0, including all GEOS-Chem Classic and GCHP run directory updates.

From dc34515acb356ccdbec516776471d3c9e7faba26 Mon Sep 17 00:00:00 2001
From: Yuanjian Zhang <yuanjian.z@wustl.edu>
Date: Thu, 24 Oct 2024 15:13:38 -0500
Subject: [PATCH 3/3] Add comment on GCAP2 scale factors

---
 run/GCClassic/createRunDir.sh | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/run/GCClassic/createRunDir.sh b/run/GCClassic/createRunDir.sh
index f385e193a..b8ebf4384 100755
--- a/run/GCClassic/createRunDir.sh
+++ b/run/GCClassic/createRunDir.sh
@@ -658,6 +658,8 @@ fi
 #-----------------------------------------------------------------
 
 if [[ ${met} = "ModelE2.1" ]]; then
+    # Current scale factors are based on U10 and V10, while the latest implementation of DEAD is using USTAR
+    # May cause 10 times larger dust emissions than properly scaled
     if [[ "$runid" == "E213f10aF40oQ40nudge" ]]; then
         if [[ "$grid_res" ==  "4x5" ]]; then
 	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='0.00474046'\n"
@@ -685,7 +687,6 @@ if [[ ${met} = "ModelE2.1" ]]; then
     fi
 else
     RUNDIR_VARS+="RUNDIR_GISS_RES='not_used'\n"
-    # Use GEOS-FP values as placeholders for GEOS-IT until parameters derived
     if [[ "x${sim_name}" == "xfullchem" || "x${sim_name}" == "xaerosol" ]]; then
 	if [[ "x${met}" == "xgeosfp" && "x${grid_res}" == "x4x5" ]]; then
 	    RUNDIR_VARS+="RUNDIR_DUSTDEAD_TF='4.8632e-5'\n"