Update jobs/rocoto/prepchem.sh

zhanglikate · May 3, 2024 · 4530c46 · 4530c46
1 parent 087cf50
commit 4530c46
Showing 1 changed file with 32 additions and 14 deletions.
diff --git a/jobs/rocoto/prepchem.sh b/jobs/rocoto/prepchem.sh
@@ -109,7 +109,7 @@ 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 ] && [ ${emiss_date1} -le 20240501 ] ];  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
@@ -192,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 20230509 ]; 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