Date: Mon Apr 22 2013 - 10:10:15 MDT
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
