Still struggling with regrid2

[Daniel Victoria]

Hi Michael,

Xavier Francis helped me out and I figured out what was going on.
regrid2 was adding some more lines and columns to the output file. So,
the way I solved the issue was to export the regrided file to a binary
grid (set gxout fwrite). When I do this, GrADS reports the number of
columns/rows and the initial lat/lon. So I can then write a correct
.ctl file and export to NetCDF using lats4d. A lengthier process but
it works.

As for the details, I did sent a more detailed email before but with
the input and output dimensions (template ctl and netcdf dims). I had
no idea it could be a compilation problem with regrid2, specially
because I know nothing about fortran.
Just to answer your questions, the open line in regrid2beta.f is:

open (iunit_diag,file='udf.regrid.out',
     $       form='formatted',status='unknown')

Can't recall where I got the code from but it was about an year ago.
Is there a more recent version?

Cheers
Daniel


On 5/29/07, Michael Fiorino <Michael.Fiorino at noaa.gov> wrote:
> daniel, this looks like a regrid2 build problem.  the most annoying
> thing about fortran is the lack of a binary  i/o standard...  look at
> the 'open' line in the .f -- using direct access?  sequential?  and the
> usual plea -- if you want help, you HAVE to send details -- operating
> system, where you got the regrid code, the .ctl file you're working...
> /r mike
>
> Daniel Victoria wrote:
> > Hi all,
> >
> > I still haven't been able to regrid the gausian NCEP file to a uniform
> > 0.25 x 0.25 lat lon grid. The result I get is all messed up (attached
> > file). Does anyone knows what I'm doing wrong? Or is grads not able to
> > regrid the gaussian NCEP to uniform grid?
> >
> > Thanks
> > Daniel
> >
> > ------------------------------------------------------------------------
> >
>