From 924deeed0136867890977565fa7bb7f9ce7e5d1b Mon Sep 17 00:00:00 2001
From: "Kate.Zhang-NOAA" <64156100+zhanglikate@users.noreply.github.com>
Date: Fri, 3 May 2024 16:18:09 -0600
Subject: [PATCH] Update to GBBPEx v005 (#55)

* Update to GBBPEx v005

* Update jobs/rocoto/prepchem.sh

* Update jobs/rocoto/prepchem.sh

* Update prepchem.sh

* Fix jobs/rocoto/prepchem.sh
---
 jobs/rocoto/prepchem.sh | 52 ++++++++++++++++++++++++++++++-----------
 1 file changed, 38 insertions(+), 14 deletions(-)

diff --git a/jobs/rocoto/prepchem.sh b/jobs/rocoto/prepchem.sh
index 762e4a4284..fcaf21ad33 100755
--- a/jobs/rocoto/prepchem.sh
+++ b/jobs/rocoto/prepchem.sh
@@ -109,12 +109,18 @@ for n in $(seq 1 6); do
       else
 
         if [ ${SYEAR} -eq 2016 -o ${emiss_date1} -ge 20230115 ];  then           ## JKH - change date
-          if [ ${emiss_date1} -ge 20230509 ];  then                              ## JKH - change date
+	  if [[ ${emiss_date1} -ge 20230509  &&  ${emiss_date1} -le 20240501 ]]; then 
             BC=GBBEPxemis-BC-${CASE}GT${n}_v4r0_${emiss_date1}.bin
             OC=GBBEPxemis-OC-${CASE}GT${n}_v4r0_${emiss_date1}.bin
             SO2=GBBEPxemis-SO2-${CASE}GT${n}_v4r0_${emiss_date1}.bin
             FRP=GBBEPxFRP-MeanFRP-${CASE}GT${n}_v4r0_${emiss_date1}.bin
             PM25=GBBEPxemis-PM25-${CASE}GT${n}_v4r0_${emiss_date1}.bin
+	  elif [[ ${emiss_date1} -ge 20240502 ]]; then
+            BC=GBBEPxemis-BC-${CASE}GT${n}_v5r0_${emiss_date1}.bin
+            OC=GBBEPxemis-OC-${CASE}GT${n}_v5r0_${emiss_date1}.bin
+            SO2=GBBEPxemis-SO2-${CASE}GT${n}_v5r0_${emiss_date1}.bin
+            FRP=GBBEPxFRP-MeanFRP-${CASE}GT${n}_v5r0_${emiss_date1}.bin
+            PM25=GBBEPxemis-PM25-${CASE}GT${n}_v5r0_${emiss_date1}.bin
           else
             BC=GBBEPx.emis_BC.003.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
             OC=GBBEPx.emis_OC.003.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
@@ -186,19 +192,37 @@ EOF
         echo "NetCDF GBBEPx File $NCGB/${SYEAR}${nmonth}${nday}/FIRE_GBBEPx_data.tile${n}.nc  exists, just link."
       else
 
-        #if [ ${SYEAR} -eq 2016 ];  then
-          BC=GBBEPxemis-BC-${CASE}GT${n}_v4r0_${SYEAR}${nmonth}${nday}.bin
-          OC=GBBEPxemis-OC-${CASE}GT${n}_v4r0_${SYEAR}${nmonth}${nday}.bin
-          PM25=GBBEPxemis-PM25-${CASE}GT${n}_v4r0_${SYEAR}${nmonth}${nday}.bin
-          SO2=GBBEPxemis-SO2-${CASE}GT${n}_v4r0_${SYEAR}${nmonth}${nday}.bin
-          FRP=GBBEPxFRP-MeanFRP-${CASE}GT${n}_v4r0_${SYEAR}${nmonth}${nday}.bin
-        #else
-        #  BC=GBBEPx.bc.${SYEAR}${nmonth}${nday}.FV3.${CASE}Grid.tile${n}.bin
-        #  OC=GBBEPx.oc.${SYEAR}${nmonth}${nday}.FV3.${CASE}Grid.tile${n}.bin
-        #  PM25=GBBEPx.pm25.${SYEAR}${nmonth}${nday}.FV3.${CASE}Grid.tile${n}.bin
-        #  SO2=GBBEPx.so2.${SYEAR}${nmonth}${nday}.FV3.${CASE}Grid.tile${n}.bin
-        #  FRP=meanFRP.${SYEAR}${nmonth}${nday}.FV3.${CASE}Grid.tile${n}.bin
-        #fi
+	 if [ ${SYEAR} -eq 2016 -o ${emiss_date1} -ge 20230115 ];  then           ## JKH - change date
+          if [[ ${emiss_date1} -ge 20230509  &&  ${emiss_date1} -le 20240501 ]];  then                              ## JKH - change date
+            BC=GBBEPxemis-BC-${CASE}GT${n}_v4r0_${emiss_date1}.bin
+            OC=GBBEPxemis-OC-${CASE}GT${n}_v4r0_${emiss_date1}.bin
+            SO2=GBBEPxemis-SO2-${CASE}GT${n}_v4r0_${emiss_date1}.bin
+            FRP=GBBEPxFRP-MeanFRP-${CASE}GT${n}_v4r0_${emiss_date1}.bin
+            PM25=GBBEPxemis-PM25-${CASE}GT${n}_v4r0_${emiss_date1}.bin
+	  elif [[ ${emiss_date1} -ge 20240502 ]]; then
+            BC=GBBEPxemis-BC-${CASE}GT${n}_v5r0_${emiss_date1}.bin
+            OC=GBBEPxemis-OC-${CASE}GT${n}_v5r0_${emiss_date1}.bin
+            SO2=GBBEPxemis-SO2-${CASE}GT${n}_v5r0_${emiss_date1}.bin
+            FRP=GBBEPxFRP-MeanFRP-${CASE}GT${n}_v5r0_${emiss_date1}.bin
+            PM25=GBBEPxemis-PM25-${CASE}GT${n}_v5r0_${emiss_date1}.bin
+          else
+            BC=GBBEPx.emis_BC.003.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
+            OC=GBBEPx.emis_OC.003.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
+            SO2=GBBEPx.emis_SO2.003.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
+            FRP=GBBEPx.FRP.003.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
+            if [ ${SYEAR} -eq 2016 ];  then
+              PM25=GBBEPx.emis_PM25.003.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
+            else
+              PM25=GBBEPx.emis_PM2.5.003.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
+            fi
+          fi
+        else
+          BC=GBBEPx.bc.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
+          OC=GBBEPx.oc.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
+          PM25=GBBEPx.pm25.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
+          SO2=GBBEPx.so2.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
+          FRP=meanFRP.${emiss_date1}.FV3.${CASE}Grid.tile${n}.bin
+        fi
 
         mkdir -p $NCGB/${SYEAR}${nmonth}${nday}
         set -ue