Re: Grid to grid interpolation (Correia, James)

I have found it tolerably fast to use rcm2points in a loop over the x & y dimensions, calculating values for one gridcell at a time:

do j = 0,ny-1
  do i = 0,nx-1
    lon = wrflon2d(j,i)
    lat = wrflat2d(j,i)
    temp = rcm2points(narrlat2d,narrlon2d,indata,lat,lon,1)
    outdata(:,j,i) = (/temp/)
    end do
end do

(I can't remember whether I decided to try looping to avoid memory problems with certain big files, or because it seemed easier than massaging the output back into the shape I wanted, or both, but it works.)

I *think* rcm2points uses the same inverse distance weighting as linint2; I have no idea how it performs with regard to position errors.

Cheers,

--Seth McGinnis

>
>Hi-
>Does anyone have advice on interpolating from the NARR grid to a WRF grid
>(both have 2d lat lon arrays)? What method would be pretty fast?
>
>Triple2grid won't do because I am concerned about even minor position errors
>as I plan to do some validating of the wrf model dataset.
>Linint2 does interpolation but not grid to grid.
>
>Thanks
>
>James Correia Jr., PhD
>Climate Physics Group
>Post. Doc.
>Pacific Northwest National Lab

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun May 17 2009 - 17:59:40 MDT