All,<br><br>Apologies in advance if this is a duplicate post....<br>
<br>
I was wondering if anyone could lend me some insight into using
regrid2. I ended up getting a pre-compiled binary which seems to be
working ok (i.e. it doesn't seg-fault when I run it). I've basically
got 2 datasets I want to subtract from one another, one is a model grib
file, the other is a binary file.I've come up with a control file which
will allow me to view the field in the binary file, and I've already
got a functioning control file for my model grib file. The binary file
contains 1km gridded data, the model grib file contains 3km gridded
data, and I believe both are on different map projections. Here's a
portion of each of the control files....<br>
<br>
Binary control file:<br>
dset /usr/local/leins/test/us_ssmv01025SlL01T0024TTNATS2007040705DP001.dat<br>
title test<br>
#undef 65535<br>
undef 55537<br>
UNPACK 0.1 0.0<br>
options yrev<br>
xdef 6935 linear -124.729583333333 0.00833333333333333<br>
ydef 3351 linear 24.9495833333324 0.00833333333333333<br>
zdef 1 levels 1000<br>
tdef 1 linear 22jul2007 1dy<br>
<br>
<br>
Model grib file:<br>
set /usr/local/leins/test//07040700_arw_d01.GrbF02400<br>
index /usr/local/leins/test//07040700_arw_d01.GrbF02400.idx<br>
undef 9.999E+20<br>
title /usr/local/leins/test//07040700_arw_d01.GrbF02400<br>
* produced by grib2ctl v0.9.12.5p34c<br>
dtype grib 255<br>
pdef 263 250 lcc 39.447000 -87.971000 1 1 42.900000 42.900000 -83.400000 3000 30<br>
00<br>
xdef 263 linear -87.971000 0.0353176389800391<br>
ydef 250 linear 39.447000 0.0272727272727273<br>
tdef 1 linear 00Z08apr2007 3hr<br>
<br>
In short, if I attempt to subtract one from the other, GrADS complains
with "Entire Grid Undefined". Thinking that because they are on
different map projections and of different resolutions, I thought I'd
have to regrid the binary data to "dumb it down" to the model 3km
resolution. Does that sound right? So far I haven't seen a whole lot of
helpful documentation on regrid2, but I did find some help for lats4d.