From b27a01db1516f9aa229e75957885ffa2125d31d6 Mon Sep 17 00:00:00 2001
From: "fanglin.yang" <fanglin.yang@noaa.gov>
Date: Wed, 3 Jun 2020 16:58:29 +0000
Subject: [PATCH] modified:   parm/config/config.fcst modified:  
 scripts/exglobal_fcst_nemsfv3gfs.sh

GFS forecast restrat capability from a breakpoint is no longer working
with IAU turned on.  This function has been overhauled to make it more
general and works for cases with and without IAU
---
 parm/config/config.fcst             | 16 +++++++++++++++-
 scripts/exglobal_fcst_nemsfv3gfs.sh | 23 ++++++++++++++++-------
 2 files changed, 31 insertions(+), 8 deletions(-)

diff --git a/parm/config/config.fcst b/parm/config/config.fcst
index a26729bd08..8954b6df5d 100755
--- a/parm/config/config.fcst
+++ b/parm/config/config.fcst
@@ -214,8 +214,22 @@ elif [[ "$CDUMP" == "gfs" ]] ; then # GFS cycle specific parameters
     fi
 
     # Write gfs restart files to rerun fcst from any break point         
-     export restart_interval=${restart_interval_gfs:-0}
+    export restart_interval_gfs=${restart_interval_gfs:-0}
+    if [ $restart_interval_gfs -le 0 ]; then
+        export restart_interval=0
+    else
+        rst_list=""
+        IAU_OFFSET=${IAU_OFFSET:-0}
+        [[ $DOIAU = "NO" ]] && export IAU_OFFSET=0
+        xfh=$((restart_interval_gfs+(IAU_OFFSET/2)))
+        while [ $xfh -le $FHMAX_GFS ]; do
+          rst_list="$rst_list $xfh"
+          xfh=$((xfh+restart_interval_gfs))
+        done
+        export restart_interval="$rst_list"
+    fi
 
+    
     # Choose coupling with wave
     if [ $DO_WAVE = "YES" -a "$WAVE_CDUMP" != "gdas" ]; then 
         export cplwav=".true."
diff --git a/scripts/exglobal_fcst_nemsfv3gfs.sh b/scripts/exglobal_fcst_nemsfv3gfs.sh
index e290a50998..6f60bbc51d 100755
--- a/scripts/exglobal_fcst_nemsfv3gfs.sh
+++ b/scripts/exglobal_fcst_nemsfv3gfs.sh
@@ -20,6 +20,7 @@
 # 2019-03-21  Fanglin Yang   Add restart capability for running gfs fcst from a break point.
 # 2019-12-12  Henrique Alves Added wave model blocks for coupled run
 # 2020-01-31  Henrique Alves Added IAU capability for wave component
+# 2020-06-02  Fanglin Yang   restore restart capability when IAU is turned on.                     
 #
 # $Id$
 #
@@ -160,21 +161,22 @@ fi
 RERUN="NO"
 filecount=$(find $RSTDIR_TMP -type f | wc -l) 
 if [ $CDUMP = "gfs" -a $rst_invt1 -gt 0 -a $FHMAX -gt $rst_invt1 -a $filecount -gt 10 ]; then
-    SDATE=$($NDATE +$FHMAX $CDATE)
-    EDATE=$($NDATE +$rst_invt1 $CDATE)
-    while [ $SDATE -gt $EDATE ]; do
+    reverse=$(echo "${restart_interval[@]} " | tac -s ' ')
+    for xfh in $reverse ; do
+        yfh=$((xfh-(IAU_OFFSET/2)))
+        SDATE=$($NDATE +$yfh $CDATE)
         PDYS=$(echo $SDATE | cut -c1-8)
         cycs=$(echo $SDATE | cut -c9-10)
         flag1=$RSTDIR_TMP/${PDYS}.${cycs}0000.coupler.res
         flag2=$RSTDIR_TMP/coupler.res
         if [ -s $flag1 ]; then
+            CDATE_RST=$SDATE          
+            [[ $RERUN = "YES" ]] && break
             mv $flag1 ${flag1}.old
             if [ -s $flag2 ]; then mv $flag2 ${flag2}.old ;fi
             RERUN="YES"
-            CDATE_RST=$($NDATE -$rst_invt1 $SDATE)
-            break
+            [[ $xfh = $rst_invt1 ]] && RERUN="NO"
         fi 
-        SDATE=$($NDATE -$rst_invt1 $SDATE)
     done
 fi
 
@@ -314,6 +316,13 @@ EOF
       file2=$(echo $file2 | cut -d. -f3-) 
       $NLN $file $DATA/INPUT/$file2
     done
+   
+    if [ $DOIAU = "YES" ]; then
+      IAUFHRS=-1         
+      IAU_DELTHRS=0
+      IAU_INC_FILES="''"
+    fi
+
 
   fi
 #.............................
@@ -762,7 +771,7 @@ start_day:               ${tPDY:6:2}
 start_hour:              ${tcyc}
 start_minute:            0
 start_second:            0
-fhrot:                   ${IAU_FHROT}
+fhrot:                   ${IAU_FHROT:-0}
 nhours_fcst:             $FHMAX
 RUN_CONTINUE:            ${RUN_CONTINUE:-".false."}
 ENS_SPS:                 ${ENS_SPS:-".false."}