I have two SST grids of differing dimensions, one low-res grid which
covers my entire model domain (GAREA=15;-110;50;-50), and one high-res
which covers only a small subset of my model domain (GAREA=25;-90;40;-67).
Ultimately what I'd like to do is: overlay the high-res grid on the low
res-grid, perform a smoothing function, save the overlay as a new grid,
convert the overlay grid to grib so that it can be used as input into WRF.
I was thinking that I could just use gdbiint to interpolate the high-res
grid to the same grid as the low-res SSTs, hoping that GEMPAK would just
write missing values in for the area for which the high-res grid does not
cover. But when I do this my high-res grid has SST values ranging from
negative infinity to infinity. So either I've done something wrong...or
GEMPAK expects the grid to cover the same GAREA as the grid it is being
I know I could use fortran to expand the high-res grid with blank
columns/rows so that it covers my entire model domain. I am really hoping
I can avoid this though, as it would be a bit painstaking with my limited
fortran background to ensure that I'm adding the correct number of grid
points in the correct places.
I figure since I can do the overlay visually in GEMPAK, that there must be
a way to write out the results to a new grid...as it seems like such a
simple thing to do. Is there a program I'm overlooking? Can someone help
me out with this?