* Sync NCO SPA's ecf definition file updates to launch the HAFS workflow

  at T+3:05 so that the HAFS track products can be delivered on time (before T+6:05).
* Update the ecf definition file to trigger the jhafs_obs_prep job properly.
* Update ecf/scripts/prep/jhafs_atm_prep.ecf to use 1 node (instead of 2
  nodes, same as jhafs_atm_prep_mvnest.ecf)
* Update exhafs_atm_lbc.sh to check and wait for the GFS pgrb2b files as well.
* Update JHAFS_LAUNCH so that it will check and ensure HFSB will only
  run NHC storms (but not run JTWC storms).
* Update storminfo.py to deal with storms at extremely high latitudes.
* HAFS datool related bug fixes related to the atm_vi job failures for
  storms at extremely high latitudes. (From @YonghuiWeng-NOAA)
* Update submodules hafs_forecast.fd/FV3/upp to add DPT at six pressure
  levels to match HWRF grib2 products. (From @LinZhu-NOAA)

Author: BinLiu-NOAA

ecf/defs/hafs.def

@@ -4,7 +4,7 @@
 #PBS -q %QUEUE%
 #PBS -A %PROJ%-%PROJENVIR%
 #PBS -l walltime=00:30:00
-#PBS -l place=vscatter:excl,select=2:ncpus=60:mpiprocs=6:ompthreads=10
+#PBS -l place=vscatter:excl,select=1:ncpus=108:mpiprocs=18:ompthreads=6
 #PBS -l debug=true
 
 model=hafs
@@ -14,9 +14,9 @@ export cyc="%CYC%"
 %include <head.h>
 %include <envir-p1.h>
 
-export TOTAL_TASKS='12'
-export NCTSK='6'
-export OMP_THREADS='10'
+export TOTAL_TASKS='18'
+export NCTSK='18'
+export OMP_THREADS='6'
 
 export storm_num="%STORMNUM%"
 

ecf/scripts/prep/jhafs_atm_prep.ecf

@@ -103,6 +103,10 @@ else
   if [[ "$mYMDH" != "$PDY$cyc" ]]; then
     abort_reason="Message file cycle $mYMDH is not current cycle $PDY$cyc."
   fi
+  # Under NCO environment, HFSB only runs NHC/CPHC storms
+  if [[ "${RUN_ENVIR^^}" = "NCO" ]] && [[ "${RUN^^}" = "HFSB" ]] && [[ "$center" != "NHC" ]]; then
+    abort_reason="The storm in $message file is not a NHC storm for ${RUN^^} to run."
+  fi
 fi
 
 if [[ -n "$abort_reason" ]]; then

jobs/JHAFS_LAUNCH

@@ -188,6 +188,28 @@ while [ $n -le 360 ]; do
   n=$(( n+1 ))
 done
 
+# Check and wait for the pgrb2b input data if needed.
+if [ $input_type = "grib2" ] && [ $bctype = gfsgrib2ab_0p25 ]; then
+
+n=1
+while [ $n -le 360 ]; do
+  if [ -s ${INIDIR}/${CDUMP}.t${cyc}z.pgrb2b.0p25.f${FHR3} ]; then
+	while [ $(( $(date +%s) - $(stat -c %Y ${INIDIR}/${CDUMP}.t${cyc}z.pgrb2b.0p25.f${FHR3}) )) -lt 10  ]; do sleep 10; done
+    echo "${INIDIR}/${CDUMP}.t${cyc}z.pgrb2b.0p25.f${FHR3} ready, do chgres_bc"
+    break
+  else
+    echo "Either ${INIDIR}/${CDUMP}.t${cyc}z.pgrb2b.0p25.f${FHR3} not ready, sleep 10"
+    sleep 10s
+  fi
+  if [ $n -ge 360 ]; then
+    echo "FATAL ERROR: Waited too many times: $n. Exiting"
+    exit 1
+  fi
+  n=$(( n+1 ))
+done
+
+fi
+
 if [ $input_type = "grib2" ]; then
   if [ $bctype = gfsgrib2ab_0p25 ]; then
     # Use both ${CDUMP}.t${cyc}z.pgrb2.0p25.f${FHR3} and ${CDUMP}.t${cyc}z.pgrb2b.0p25.f${FHR3} files

scripts/exhafs_atm_lbc.sh

@@ -2232,7 +2232,7 @@ subroutine hafsvi_postproc(in_file, in_date, out_dir, nestdoms)
                  ke=min(iz,my_proc_id*nm+nm)
                  if ( ks >= 1 .and. ks <= iz .and. ke >= 1 .and. ke <= iz ) then
                     allocate(dat44(ix+nv-1, iy+2-nv, ke-ks+1, 1))
-                    call mpi_recv(dat44, size(dat43), mpi_real, io_proc, 200*nv+ks, comm, status, ierr)
+                    call mpi_recv(dat44, size(dat44), mpi_real, io_proc, 200*nv+ks, comm, status, ierr)
                  endif
               endif  !if ( my_proc_id == io_proc ) then
 

sorc/hafs_forecast.fd

@@ -387,8 +387,7 @@ subroutine cal_src_dst_grid_weight(grid_src, grid_dst)
   call search_nearst_grid(grid_src%grid_xt, grid_src%grid_y, lat_src, lon_src, &
                           grid_dst%grid_xt, grid_dst%grid_y, lat_dst, lon_dst, x_oini, y_oini)
 
-  write(*,'(a3, 2i6, 2f11.2,2i11, 2f11.2)')' u:', i, j, lon_dst(i,j), lat_dst(i,j), x_oini(i,j), y_oini(i,j), &
-     lon_src(x_oini(i,j),y_oini(i,j)), lat_src(x_oini(i,j),y_oini(i,j))
+  write(*,'(a3, 2i6, 2f11.2,2i11)')' u:', i, j, lon_dst(i,j), lat_dst(i,j), x_oini(i,j), y_oini(i,j)
 
   allocate(gwt%gwt_u(grid_dst%grid_xt, grid_dst%grid_y))  !grid_weight_info
   call cal_grid_weight(grid_src%grid_xt, grid_src%grid_y, lat_src, lon_src, &
@@ -422,8 +421,7 @@ subroutine cal_src_dst_grid_weight(grid_src, grid_dst)
   call search_nearst_grid(grid_src%grid_x, grid_src%grid_yt, lat_src, lon_src, &
                           grid_dst%grid_x, grid_dst%grid_yt, lat_dst, lon_dst, x_oini, y_oini)
 
-  write(*,'(a3,2i6, 2f11.2,2i11, 2f11.2)')' v:', i, j, lon_dst(i,j), lat_dst(i,j), x_oini(i,j), y_oini(i,j), &
-     lon_src(x_oini(i,j),y_oini(i,j)), lat_src(x_oini(i,j),y_oini(i,j))
+  write(*,'(a3,2i6, 2f11.2,2i11)')' v:', i, j, lon_dst(i,j), lat_dst(i,j), x_oini(i,j), y_oini(i,j)
 
   allocate(gwt%gwt_v(grid_dst%grid_x, grid_dst%grid_yt))  !grid_weight_info
   call cal_grid_weight(grid_src%grid_x, grid_src%grid_yt, lat_src, lon_src, &

sorc/hafs_tools.fd/sorc/hafs_datool/sub_hafsvi_proc.f90

@@ -606,7 +606,6 @@ def tcutil_domain_center(self,logger=None):
         if ilat >= 40: cenla=35.0
         if ilat >= 44: cenla=40.0
         if ilat >= 50: cenla=45.0
-        if ilat >= 55: cenla=50.0
         if storm_lat<0: cenla=-cenla
 
         # Decide the center longitude.