Re: RCM to RCM

Dear Ed,

You may want to consider using the ESMF suite of regridding routines.

See example ESMF_regrid_16.ncl at:

http://www.ncl.ucar.edu/Applications/ESMF.shtml#ex16

which regrids from a high-res curvilinear grid to a low-res curvilinear grid.

You can also see some ESMF regridding templates at:

http://www.ncl.ucar.edu/Applications/Templates/

The one you want is "ESMF_curv_to_curv.ncl.

--Mary

On Nov 6, 2013, at 9:07 PM, Yang Yang <yang123yang@yahoo.com> wrote:

> Hi There:
>
> I would like to do interpolation from a low-resolution RCM to a high-resolution RCM. I checked NCL and couldn't find a function like rcm2rcm. I had to use rcm2points with the following codes:
> do i = 1, 1000
> do j= 1, 1000
> pp(i,j) = rcm2points (tlatb,tlonb,mslpb,tlat(i,j),tlon(i,j),0)
> uu(i,j) = rcm2points (tlatb,tlonb,u10mb,tlat(i,j),tlon(i,j),0)
> vv(i,j) = rcm2points (tlatb,tlonb,v10mb,tlat(i,j),tlon(i,j),0)
> end do
> end do
> However, the codes takes a lot of CPU time. Very slow. Can someone give me some suggestions how to optimize the codes? Thanks!
>
> Ed
> _______________________________________________
> 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
Received on Thu Nov 7 08:23:31 2013