From 716a9caebf6020b244f6d61c58b59c983f980aae Mon Sep 17 00:00:00 2001
From: "Judy.K.Henderson" <Judy.K.Henderson@noaa.gov>
Date: Fri, 29 Mar 2024 15:01:29 +0000
Subject: [PATCH]  Add changes for running pygraf under Rocky

  - call remapgrib.sh due to wgrib2 change
  - current updated checkout is under /scratch1/BMC/gsd-fv3/exp/pygraf
---
 FV3GFSwfm/rt_v17p8_ugwpv1_c3_mynn/runcmds |   2 +-
 jobs/rocoto/remapgrib.sh                  | 110 ++++++++++++++++++++++
 2 files changed, 111 insertions(+), 1 deletion(-)
 create mode 100755 jobs/rocoto/remapgrib.sh

diff --git a/FV3GFSwfm/rt_v17p8_ugwpv1_c3_mynn/runcmds b/FV3GFSwfm/rt_v17p8_ugwpv1_c3_mynn/runcmds
index aa43700fbc..52f4d12be1 100644
--- a/FV3GFSwfm/rt_v17p8_ugwpv1_c3_mynn/runcmds
+++ b/FV3GFSwfm/rt_v17p8_ugwpv1_c3_mynn/runcmds
@@ -6,7 +6,7 @@ rocotostat   -w ${HOME}/UFS-CAMsuite_dev1/FV3GFSwfm/rt_v17p8_ugwpv1_c3_mynn/v17p
 
 
 
-PyGraf workflow
+PyGraf workflow  (working as of 03/29/24)
 ===============
 rocotorun -w ${HOME}/UFS-CAMsuite_dev1/FV3GFSwfm/rt_v17p8_ugwpv1_c3_mynn/rt_pygraf_global_ugwpv1_c3_mynn.xml -d ${HOME}/rt_dbfiles/rt_pygraf_global_ugwpv1_c3_mynn.db
 rocotostat -w ${HOME}/UFS-CAMsuite_dev1/FV3GFSwfm/rt_v17p8_ugwpv1_c3_mynn/rt_pygraf_global_ugwpv1_c3_mynn.xml -d ${HOME}/rt_dbfiles/rt_pygraf_global_ugwpv1_c3_mynn.db -c  `date --date='4 days ago' +%Y%m%d0000`: | m
diff --git a/jobs/rocoto/remapgrib.sh b/jobs/rocoto/remapgrib.sh
new file mode 100755
index 0000000000..7c01ab1c1b
--- /dev/null
+++ b/jobs/rocoto/remapgrib.sh
@@ -0,0 +1,110 @@
+#!/bin/bash --login
+
+# this file remaps 1/4 degree lat/lon grids to other grids
+#    236 201 244 130 224 242
+
+# initialize
+module load gnu/13.2.0 intel/2023.2.0 netcdf/4.7.0 wgrib2/3.1.2_ncep
+module list
+ECHO=echo
+MKDIR=mkdir
+LN=ln
+echo entering remapgrib.sh....
+echo "****************************"
+echo "ROTDIR      = ${ROTDIR}"
+echo "CDUMP       = ${CDUMP}"
+echo "COMPONENT   = ${COMPONENT}"
+echo "yyyymmdd    = ${yyyymmdd}"
+echo "hh          = ${hh}"
+echo "fcst        = ${fcst}"
+echo "GRID_NAMES  = ${GRID_NAMES}"
+echo
+
+echo `which wgrib2`
+
+#########################
+# Grid Definitions
+#########################
+# Full domain
+#---------
+# 32 km
+export grid_specs_221="lambert:253:50.000000 214.500000:349:32463.000000 1.000000:277:32463.000000"
+
+#---------
+# CONUS
+#---------
+# 13 km
+export grid_specs_130="lambert:265:25.000000 233.862000:451:13545.000000 16.281000:337:13545.000000"
+# 20 km
+export grid_specs_252="lambert:265:25.000000 233.862000:301:20318.000000 16.281000:225:20318.000000"
+# 40 km
+export grid_specs_236="lambert:265:25.000000 233.862000:151:40635.000000 16.281000:113:40635.000000"
+
+#---------
+# Alaska
+#---------
+export grid_specs_242="nps:225:60.000000 187.000000:553:11250.000000 30.000000:425:11250.000000"
+
+#---------
+# Hawaii
+#---------
+export grid_specs_243="latlon 190.0:126:0.400 10.000:101:0.400"
+
+#---------
+# Puerto Rico
+#---------
+export grid_specs_200="lambert:253:50.000000 285.720000:108:16232.000000 16.201000:94:16232.000000"
+
+#---------
+#HRRRE
+#---------
+export grid_specs_999="lambert:253:50.000000 227.500000:675:13545.000000 7.500000:500:13545.000000"
+
+#---------
+# North Polar Stereographic
+#---------
+export grid_specs_201="nps:-105:60.000000 -150.000000:259:94512.000000 -20.826000:259:94512.000000"
+
+#---------
+# South Polar Stereographic
+#---------
+export grid_specs_224="sps:75:-60.000000 120.000000:257:95250.000000 20.826000:257:95250.000000"
+
+#---------
+# North Atlantic
+#---------
+export grid_specs_244="latlon 261.750:275:0.25 0.250:203:0.25"
+#########################
+
+# make post directory if doesn't exist
+postDir=${ROTDIR}/${CDUMP}.${yyyymmdd}/${hh}/products/${COMPONENT}/grib2/0p25/post
+echo "postDir:  $postDir"
+if [ ! -d ${postDir} ]
+then
+  echo "creating ${postDir} "
+  mkdir -p ${postDir}
+fi
+
+# parse out domain
+grids=$(echo $GRID_NAMES|sed 's/D/ /g')
+src_gribfile=gfs.t${hh}z.pgrb2.0p25.f${fcst}
+
+# loop through each domain and remap 1/4 degree grib2 file to new grid
+grids=$(echo $GRID_NAMES|sed 's/D/ /g')
+for grid in $grids
+do
+    tgt_gribfile_dir=${postDir}/${grid}
+    ${MKDIR} -p ${tgt_gribfile_dir}
+    tgt_gribfile=${tgt_gribfile_dir}/${src_gribfile}
+    ${ECHO} "Processing grids for grid ${grid}"
+    eval grid_specs=\${grid_specs_${grid}}
+    wgrib2 ${postDir}/../${src_gribfile} -set_grib_type c3 -new_grid_winds grid \
+      -new_grid_interpolation bilinear \
+      -new_grid ${grid_specs} ${tgt_gribfile}
+done
+
+# make links for full domain
+tgt_gribfile_dir=${postDir}/full
+${MKDIR} -p ${tgt_gribfile_dir}
+cd ${tgt_gribfile_dir}
+${LN} -fs ../../${src_gribfile}