2D interpolation from RUC to WRF

Hi,

I use following ncl script based on modified Dennis's codes

ftp://ftp.erc.msstate.edu/outgoing/yli/ABL/wrf-minus-ruc.ncl

to interpolate RUC grid to WRF grid

ftp://ftp.erc.msstate.edu/outgoing/yli/ABL/wrf-minus-ruc.pdf

The interpolated fields have sawtooth but the original
does not.

Would you give me any suggestions?

Yongzuo

_____________________________________
From: Dennis Shea [shea@ucar.edu]
Sent: Wednesday, June 30, 2010 12:37 PM
To: Li, Yongzuo
Cc: ncl-talk@ucar.edu
Subject: Re: 2D interpolation

Sorry.
There is no general NCL function to interpolate
from a WRF to NARR grid.

The following would be very slow

qlat = fnar->QLAT
qlon = fnar->QLON

dimq = dimsizes(qlat)
nqlat= dimq(0)
nqlon= dimq(1)

xlat = fwrf->XLAT(0,:,:)
xlon = fwrf->XLONG(0,:,:)

x1d = ndtooned( xlat )
y1d = ndtooned( xlon )

z = fwrf->Z ; whatever
dimz = dimsizes( z )
ntim = dimz(0)
klev = dimz(1)

opt = True
opt@mopt = 0 ; or 1

zNARR = new ( (/ntim,klev,nqlat,nqlon/), typeof(z), getFillValue(z) )

do nt=0,ntim-1
   do kl=0,klev-1
      z1d = ndtooned( z(nt,kl,:,:) )
      zNARR(nt,kl,:,:) = triple2grid(x1d,y1d,z1d,qlon,qlat,opt)
   end do
end do

On 6/29/10 10:13 PM, Li, Yongzuo wrote:
> Hi,
>
> In http://www.ncl.ucar.edu/Document/Functions/Built-in/rcm2rgrid.shtml,
> WRF and NARR can be interpolated to a rectilinear lat/lon grid.
>
> Can WRF grid be interpolated to NARR grid?
>
> Thanks.
>
> Yongzuo
> _______________________________________________
> 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 Wed Jun 30 23:22:48 2010