Skip to content

Commit

Permalink
Merge pull request NOAA-EMC#192 from JeffBeck-NOAA/community_develop
Browse files Browse the repository at this point in the history 
Mods to update real-time data on Jet and update pre-defined Alaska domains/GFDLgrid option for GSD 3-km CONUS domain
  • Loading branch information
@gsketefian
gsketefian authored Apr 14, 2020
2 parents c1bd40a + 6243029 commit 0e35e2b
Show file tree
Hide file tree
Showing 4 changed files with 151 additions and 39 deletions.
57 changes: 45 additions & 12 deletions ush/get_extrn_mdl_file_dir_info.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -331,7 +331,8 @@ fi
#-----------------------------------------------------------------------
#
if [ "${extrn_mdl_name}" = "RAPX" ] || \
[ "${extrn_mdl_name}" = "HRRRX" ]; then
[ "${extrn_mdl_name}" = "HRRRX" ] || \
[ "${extrn_mdl_name}" = "FV3GFS" -a "${MACHINE}" = "JET" ]; then
#
# Get the Julian day-of-year of the starting date and time of the exter-
# nal model run.
Expand Down Expand Up @@ -388,7 +389,11 @@ fi

# fns=( "atm" "sfc" "nst" )
fns=( "atm" "sfc" )
prefix="gfs.t${hh}z."
if [ "${MACHINE}" = "JET" ]; then
prefix="${yy}${ddd}${hh}00.gfs.t${hh}z."
else
prefix="gfs.t${hh}z."
fi
fns=( "${fns[@]/#/$prefix}" )
suffix="anl.nemsio"
fns=( "${fns[@]/%/$suffix}" )
Expand All @@ -405,11 +410,19 @@ fi
;;

"RAPX")
fns=( "${yy}${ddd}${hh}${mn}${fcst_hh}${fcst_mn}" )
if [ "${MACHINE}" = "JET" ]; then
fns=( "wrfnat_130_${fcst_hh}.grib2" )
else
fns=( "${yy}${ddd}${hh}${mn}${fcst_hh}${fcst_mn}" )
fi
;;

"HRRRX")
fns=( "${yy}${ddd}${hh}${mn}${fcst_hh}${fcst_mn}" )
if [ "${MACHINE}" = "JET" ]; then
fns=( "wrfnat_hrconus_${fcst_hh}.grib2" )
else
fns=( "${yy}${ddd}${hh}${mn}${fcst_hh}${fcst_mn}" )
fi
;;

*)
Expand Down Expand Up @@ -447,7 +460,11 @@ bination of external model (extrn_mdl_name) and analysis or forecast
"FV3GFS")
if [ "${fv3gfs_file_fmt}" = "nemsio" ]; then
fcst_hhh=( $( printf "%03d " "${lbc_update_fhrs[@]}" ) )
prefix="gfs.t${hh}z.atmf"
if [ "${MACHINE}" = "JET" ]; then
prefix="${yy}${ddd}${hh}00.gfs.t${hh}z.atmf"
else
prefix="gfs.t${hh}z.atmf"
fi
fns=( "${fcst_hhh[@]/#/$prefix}" )
suffix=".nemsio"
fns=( "${fns[@]/%/$suffix}" )
Expand All @@ -460,17 +477,33 @@ bination of external model (extrn_mdl_name) and analysis or forecast

"RAPX")
fcst_hh=( $( printf "%02d " "${lbc_update_fhrs[@]}" ) )
prefix="${yy}${ddd}${hh}${mn}"
if [ "${MACHINE}" = "JET" ]; then
prefix="wrfnat_130_"
else
prefix="${yy}${ddd}${hh}${mn}"
fi
fns=( "${fcst_hh[@]/#/$prefix}" )
suffix="${fcst_mn}"
if [ "${MACHINE}" = "JET" ]; then
suffix=".grib2"
else
suffix="${fcst_mn}"
fi
fns=( "${fns[@]/%/$suffix}" )
;;

"HRRRX")
fcst_hh=( $( printf "%02d " "${lbc_update_fhrs[@]}" ) )
prefix="${yy}${ddd}${hh}${mn}"
if [ "${MACHINE}" = "JET" ]; then
prefix="wrfnat_hrconus_"
else
prefix="${yy}${ddd}${hh}${mn}"
fi
fns=( "${fcst_hh[@]/#/$prefix}" )
suffix="${fcst_mn}"
if [ "${MACHINE}" = "JET" ]; then
suffix=".grib2"
else
suffix="${fcst_mn}"
fi
fns=( "${fns[@]/%/$suffix}" )
;;

Expand Down Expand Up @@ -554,7 +587,7 @@ has not been specified for this external model and machine combination:
sysdir="$sysbasedir/gfs.${yyyymmdd}/${hh}"
;;
"JET")
sysdir="$sysbasedir/${yyyymmdd}"
sysdir="$sysbasedir"
;;
"ODIN")
sysdir="$sysbasedir/${yyyymmdd}"
Expand Down Expand Up @@ -585,7 +618,7 @@ has not been specified for this external model and machine combination:
sysdir="$sysbasedir"
;;
"JET")
sysdir="$sysbasedir"
sysdir="$sysbasedir/${yyyymmdd}${hh}/postprd"
;;
"ODIN")
sysdir="$sysbasedir"
Expand Down Expand Up @@ -616,7 +649,7 @@ has not been specified for this external model and machine combination:
sysdir="$sysbasedir"
;;
"JET")
sysdir="$sysbasedir"
sysdir="$sysbasedir/${yyyymmdd}${hh}/postprd"
;;
"ODIN")
sysdir="$sysbasedir"
Expand Down
12 changes: 6 additions & 6 deletions ush/set_extrn_mdl_params.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -91,7 +91,7 @@ this machine and external model combination:
EXTRN_MDL_FILES_SYSBASEDIR_ICS="/scratch1/NCEPDEV/rstprod/com/gfs/prod"
;;
"JET")
EXTRN_MDL_FILES_SYSBASEDIR_ICS="/lfs3/projects/hpc-wof1/ywang/regional_fv3/gfs"
EXTRN_MDL_FILES_SYSBASEDIR_ICS="/public/data/grids/gfs/nemsio"
;;
"ODIN")
EXTRN_MDL_FILES_SYSBASEDIR_ICS="/scratch/ywang/test_runs/FV3_regional/gfs"
Expand All @@ -117,7 +117,7 @@ this machine and external model combination:
EXTRN_MDL_FILES_SYSBASEDIR_ICS="/scratch2/BMC/public/data/gsd/rap/full/wrfnat"
;;
"JET")
EXTRN_MDL_FILES_SYSBASEDIR_ICS="/public/data/grids/rap/full/wrfnat/grib2"
EXTRN_MDL_FILES_SYSBASEDIR_ICS="/misc/whome/rtrr/rap"
;;
*)
print_err_msg_exit "\
Expand All @@ -137,7 +137,7 @@ this machine and external model combination:
EXTRN_MDL_FILES_SYSBASEDIR_ICS="/scratch2/BMC/public/data/gsd/hrrr/conus/wrfnat"
;;
"JET")
EXTRN_MDL_FILES_SYSBASEDIR_ICS="/public/data/grids/hrrr/conus/wrfnat/grib2"
EXTRN_MDL_FILES_SYSBASEDIR_ICS="/misc/whome/rtrr/hrrr"
;;
*)
print_err_msg_exit "\
Expand Down Expand Up @@ -241,7 +241,7 @@ this machine and external model combination:
EXTRN_MDL_FILES_SYSBASEDIR_LBCS="/scratch1/NCEPDEV/rstprod/com/gfs/prod"
;;
"JET")
EXTRN_MDL_FILES_SYSBASEDIR_LBCS="/lfs3/projects/hpc-wof1/ywang/regional_fv3/gfs"
EXTRN_MDL_FILES_SYSBASEDIR_LBCS="/public/data/grids/gfs/nemsio"
;;
"ODIN")
EXTRN_MDL_FILES_SYSBASEDIR_LBCS="/scratch/ywang/test_runs/FV3_regional/gfs"
Expand All @@ -268,7 +268,7 @@ this machine and external model combination:
EXTRN_MDL_FILES_SYSBASEDIR_LBCS="/scratch2/BMC/public/data/gsd/rap/full/wrfnat"
;;
"JET")
EXTRN_MDL_FILES_SYSBASEDIR_LBCS="/public/data/grids/rap/full/wrfnat/grib2"
EXTRN_MDL_FILES_SYSBASEDIR_LBCS="/misc/whome/rtrr/rap"
;;
*)
print_err_msg_exit "\
Expand All @@ -290,7 +290,7 @@ this machine and external model combination:
EXTRN_MDL_FILES_SYSBASEDIR_LBCS="/scratch2/BMC/public/data/gsd/hrrr/conus/wrfnat"
;;
"JET")
EXTRN_MDL_FILES_SYSBASEDIR_LBCS="/public/data/grids/hrrr/conus/wrfnat/grib2"
EXTRN_MDL_FILES_SYSBASEDIR_LBCS="/misc/whome/rtrr/hrrr"
;;

*)
Expand Down
120 changes: 99 additions & 21 deletions ush/set_predef_grid_params.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -434,11 +434,53 @@ predefined domain:

if [ "${GRID_GEN_METHOD}" = "GFDLgrid" ]; then

print_err_msg_exit "\
The parameters for a \"${GRID_GEN_METHOD}\" type grid have not yet been specified for this
predefined domain:
PREDEF_GRID_NAME = \"${PREDEF_GRID_NAME}\"
GRID_GEN_METHOD = \"${GRID_GEN_METHOD}\""
GFDLgrid_LON_T6_CTR=-97.5
GFDLgrid_LAT_T6_CTR=38.5
GFDLgrid_STRETCH_FAC=1.0001
GFDLgrid_RES="96"
GFDLgrid_REFINE_RATIO=36

#num_margin_#num_margin_cells_T6_left=9
#GFDLgrid_ISTART_OF_RGNL_DOM_ON_T6G=$(( num_margin_cells_T6_left + 1 ))
GFDLgrid_ISTART_OF_RGNL_DOM_ON_T6G="26"

#num_margin_cells_T6_right=9
#GFDLgrid_IEND_OF_RGNL_DOM_ON_T6G=$(( GFDLgrid_RES - num_margin_cells_T6_right ))
GFDLgrid_IEND_OF_RGNL_DOM_ON_T6G="71"

#num_margin_cells_T6_bottom=9
#GFDLgrid_JSTART_OF_RGNL_DOM_ON_T6G=$(( num_margin_cells_T6_bottom + 1 ))
GFDLgrid_JSTART_OF_RGNL_DOM_ON_T6G="36"

#num_margin_cells_T6_top=9
#GFDLgrid_JEND_OF_RGNL_DOM_ON_T6G=$(( GFDLgrid_RES - num_margin_cells_T6_top ))
GFDLgrid_JEND_OF_RGNL_DOM_ON_T6G="61"

GFDLgrid_USE_GFDLgrid_RES_IN_FILENAMES="FALSE"

DT_ATMOS="40"

LAYOUT_X="36"
LAYOUT_Y="24"
BLOCKSIZE="26"

QUILTING="TRUE"

if [ "$QUILTING" = "TRUE" ]; then
WRTCMP_write_groups="1"
WRTCMP_write_tasks_per_group=$(( 1*LAYOUT_Y ))
WRTCMP_output_grid="lambert_conformal"
WRTCMP_cen_lon="-97.5"
WRTCMP_cen_lat="38.5"
WRTCMP_stdlat1="38.5"
WRTCMP_stdlat2="38.5"
WRTCMP_nx="1738"
WRTCMP_ny="974"
WRTCMP_lon_lwr_left="-122.21414225"
WRTCMP_lat_lwr_left="22.41403305"
WRTCMP_dx="3000.0"
WRTCMP_dy="3000.0"
fi

elif [ "${GRID_GEN_METHOD}" = "JPgrid" ]; then

Expand Down Expand Up @@ -800,26 +842,62 @@ predefined domain:
#
#-----------------------------------------------------------------------
#
"GSD_HRRR_AK_3km")
"GSD_RRFSAK_3km")

if [ "${GRID_GEN_METHOD}" = "GFDLgrid" ]; then
if [ "${GRID_GEN_METHOD}" = "GFDLgrid" ]; then

GFDLgrid_LON_T6_CTR=-160.8
GFDLgrid_LAT_T6_CTR=63.0
GFDLgrid_STRETCH_FAC=1.161
GFDLgrid_RES="768"
GFDLgrid_REFINE_RATIO=4

num_margin_cells_T6_left=204
GFDLgrid_ISTART_OF_RGNL_DOM_ON_T6G=$(( num_margin_cells_T6_left + 1 ))

num_margin_cells_T6_right=204
GFDLgrid_IEND_OF_RGNL_DOM_ON_T6G=$(( GFDLgrid_RES - num_margin_cells_T6_right ))

num_margin_cells_T6_bottom=249
GFDLgrid_JSTART_OF_RGNL_DOM_ON_T6G=$(( num_margin_cells_T6_bottom + 1 ))

num_margin_cells_T6_top=249
GFDLgrid_JEND_OF_RGNL_DOM_ON_T6G=$(( GFDLgrid_RES - num_margin_cells_T6_top ))

GFDLgrid_USE_GFDLgrid_RES_IN_FILENAMES="FALSE"

DT_ATMOS="18"

LAYOUT_X="24"
LAYOUT_Y="24"
BLOCKSIZE="15"

if [ "$QUILTING" = "TRUE" ]; then
WRTCMP_write_groups="1"
WRTCMP_write_tasks_per_group="2"
WRTCMP_output_grid="lambert_conformal"
WRTCMP_cen_lon="${GFDLgrid_LON_T6_CTR}"
WRTCMP_cen_lat="${GFDLgrid_LAT_T6_CTR}"
WRTCMP_stdlat1="${GFDLgrid_LAT_T6_CTR}"
WRTCMP_stdlat2="${GFDLgrid_LAT_T6_CTR}"
WRTCMP_nx="1320"
WRTCMP_ny="950"
WRTCMP_lon_lwr_left="173.734"
WRTCMP_lat_lwr_left="46.740347"
WRTCMP_dx="3000.0"
WRTCMP_dy="3000.0"
fi

print_err_msg_exit "\
The parameters for a \"${GRID_GEN_METHOD}\" type grid have not yet been specified for this
predefined domain:
PREDEF_GRID_NAME = \"${PREDEF_GRID_NAME}\"
GRID_GEN_METHOD = \"${GRID_GEN_METHOD}\"
"
elif [ "${GRID_GEN_METHOD}" = "JPgrid" ]; then

JPgrid_LON_CTR=-163.5 #HRRR-AK is -163.5
JPgrid_LAT_CTR=62.8 #HRRR-AK is 60.8
JPgrid_LON_CTR=-161.5
JPgrid_LAT_CTR=63.0

JPgrid_DELX="3000.0"
JPgrid_DELY="3000.0"

JPgrid_NX=1230 #HRRR-AK is 1300
JPgrid_NY=850 #HRRR-AK is 920
JPgrid_NX=1380
JPgrid_NY=1020

JPgrid_WIDE_HALO_WIDTH=6

Expand All @@ -837,10 +915,10 @@ predefined domain:
WRTCMP_cen_lat="${JPgrid_LAT_CTR}"
WRTCMP_stdlat1="${JPgrid_LAT_CTR}"
WRTCMP_stdlat2="${JPgrid_LAT_CTR}"
WRTCMP_nx="1169"
WRTCMP_ny="762"
WRTCMP_lon_lwr_left="172.0"
WRTCMP_lat_lwr_left="49.0"
WRTCMP_nx="1320"
WRTCMP_ny="950"
WRTCMP_lon_lwr_left="173.734"
WRTCMP_lat_lwr_left="46.740347"
WRTCMP_dx="${JPgrid_DELX}"
WRTCMP_dy="${JPgrid_DELY}"
fi
Expand Down
1 change: 1 addition & 0 deletions ush/valid_param_vals.sh
View file
Original file line number Diff line number Diff line change
Expand Up @@ -16,6 +16,7 @@ valid_vals_PREDEF_GRID_NAME=( \
"GSD_HRRR3km" \
"GSD_HRRR13km" \
"GSD_HRRR25km" \
"GSD_RRFSAK_3km" \
"GSD_RAP13km" )
valid_vals_EMC_GRID_NAME=("ak" "conus" "conus_c96" "conus_orig" "guam" "hi" "pr")
valid_vals_USE_CCPP=("TRUE" "true" "YES" "yes" "FALSE" "false" "NO" "no")
Expand Down

0 comments on commit 0e35e2b

Please sign in to comment.