Hi-
I'm trying to use a custom 12km grid in order to match previous work, and I'm finding that the INPUT/sfc_data.nc file covers a slightly different domain than the output grid, so we are getting fill values where the two diverge in extent.
See the attached image. On the left is the sfc_data.nc file generated by SRW. On the right is the vegetation type field from the output. You can see where fill values are present in the output along the north and west edges of the grid. You can also see part of Cuba in sfc_data.nc, but the output area does not extend that far south. It almost seems as if SRW is using different params to generate sfc_data.nc than it is to specify the output grid. But perhaps I have my grid params specified incorrectly? I used the RRFS_CONUS_13km as a template and modified it to match the domain I need, but it seems that even for an improperly-specified grid the input and output domains should match. Any help would be greatly appreciated.
My grid params are copied below.
Thank you!
-Paddy McCarthy.
#
#-----------------------------------------------------------------------
#
# The NOAA Fire domain with ~12km cells.
#
#-----------------------------------------------------------------------
#
"NOAA_Fire_12km")
GRID_GEN_METHOD="ESGgrid"
ESGgrid_LON_CTR="-97.0"
ESGgrid_LAT_CTR="40.0"
ESGgrid_DELX="12000.0"
ESGgrid_DELY="12000.0"
ESGgrid_NX="481"
ESGgrid_NY="300"
ESGgrid_PAZI="0.0"
ESGgrid_WIDE_HALO_WIDTH="6"
DT_ATMOS="${DT_ATMOS:-40}"
LAYOUT_X="${LAYOUT_X:-5}"
LAYOUT_Y="${LAYOUT_Y:-2}"
BLOCKSIZE="${BLOCKSIZE:-40}"
if [ "$QUILTING" = "TRUE" ]; then
WRTCMP_write_groups="1"
WRTCMP_write_tasks_per_group="2"
WRTCMP_output_grid="lambert_conformal"
WRTCMP_cen_lon="${ESGgrid_LON_CTR}"
WRTCMP_cen_lat="${ESGgrid_LAT_CTR}"
WRTCMP_stdlat1="${ESGgrid_LAT_CTR}"
WRTCMP_stdlat2="${ESGgrid_LAT_CTR}"
WRTCMP_nx="476"
WRTCMP_ny="295"
WRTCMP_lon_lwr_left="-125.0"
WRTCMP_lat_lwr_left="22.0"
WRTCMP_dx="${ESGgrid_DELX}"
WRTCMP_dy="${ESGgrid_DELY}"
fi
;;
- 52 views