Re: interpolate rotated grid to rectangulat

Alex,

The ESMF_regrid routines do not write the regridded output to a file.

The Rectilinear.nc and WRF_SCRIP.nc files that you see are intermediate files that ESMF needs to create in order to generate the weights files, which are HIRHAM_to_Rect.nc and WRF_to_Rect.nc.

If you want to write the regridded variables to a file, you will need to do this yourself. See the attached script, which writes a file called "example.nc".

For more information about writing data to a NetCDF file, see:

http://www.ncl.ucar.edu/Applications/o-netcdf.shtml

--Mary

On Apr 22, 2013, at 11:57 AM, igor akkerman wrote:

> Thanks Mary,
>
> The reason I said incorrect was because I can't find the netcdf file with slp, u, v data. I tried ncview'ing all the netcdf files created by the program and none of them had any data on slp, and wind. Where did the program save new interpolated data?
>
> These are the files that I checked:
> Rectilinear.nc, HIRHAM_to_Rect.nc, WRF_SCRIP.nc, WRF_to_Rect.nc
>
> Thank you again,
> Alex
>
> 2013/4/22 Mary Haley <haley@ucar.edu>
> Hi Igor,
>
> I have your files and am able to run your script. I need more information on why you think the interpolated values are not correct. Please see the attached plots, which I created by modifying your interpolate.ncl script.
>
> --Mary
>
>
>
> On Apr 19, 2013, at 4:04 PM, igor akkerman wrote:
>
> > Hi
> >
> > I have modified the code to interpolate the rotated grid onto rectangular lat-lon grid using ESMF function, the code is working but the interpolated it is not correct. I was hoping you could take a look at the files.
> >
> > I have dumped 3 files onto yout ftp server. Their names are:
> > 1) text file: hirham_rot_geo_koor.dat - text file containing info on the rotated grid (1st column - lat, 2nd column - lon) and rectangular equivalent of the grid (3d column - lat, 4th - lon)
> > 2) model output file: atm.1961-02.daily.nc
> > 3) ncl code: interpolate.ncl
> >
> > Below is some additional info on the model output grid:
> >
> > Here is a sketch of the rotated grid (rot.lat, rot.lon)
> > (24.5,-27.5) ------ (24.5,27.0)
> > | |
> > | |
> > (-25.0,-27.5) ------ (-25.0,27.0)
> > This is a regular grid with uniform resol. of 0.5.
> >
> > and in geograf. grid (geograf. lat, lon):
> > (53.8,-134.6) ------ (54.2,135.1)
> > | |
> > | |
> > (53.5,-44.7) ------ (53.85,44.2)
> >
> >
> > We have 110 (nlon) x 100 (nlat) grid points;
> > resolution is ~50 km (i.e. 0.5 deg in rotated grid; DLON=0.5,DLAT=0.5)
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
>

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk