# Re: rcm2points/WRF

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Tue Jun 07 2011 - 11:29:17 MDT

On 6/3/11 3:32 PM, Nikola Marjanovic wrote:
> Hi,
>
> What is the best way to interpolate WRF data to a specific point in my
> domain at specified elevations. My approach has been to use rcm2points
> to interpolate data in 3 spatial dimensions (time, nz, ny, nx) to a
> specified latitude and longitude. Does rcm2points interpolate at each nz
> point independently in 2D (without using the levels above or below)

It interpolates in 2D only.

or
> does it interpolate at each nz level using data points at nz+1 and nz-1?
> If it only interpolates in 2D, is there a function that does this in 3D
> (I think that would be useful in complex terrain)?

You could use rcm2points and then
http://www.ncl.ucar.edu/Document/Functions/Built-in/int2p_n.shtml
to interpolate in the vertical. [Ignore the text about pressure levels.]
>
> What I then do is extract the Full Model Height of the domain by using
> wrf_user_getvar(a,"z",-1) for all times to get z(time,nz,ny,nx). I use
> rcm2points again to get the elevations at the specified lat lon location
> and use this output to interpolate to my specified elevations by using
> wrf_interp_1d(dom, z, z_spec). I create z_spec by using the function
> onedtond to create a repetitive array with my specified elevations at
> each time. (dom has dimensions [time, nz])
>
> Is this the best way to go about doing this or are there better options?

==
I think your approach is the best way to do it.

>
> Thanks,
> Nikola
> _______________________________________________
> 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 Tue Jun 7 11:29:25 2011

This archive was generated by hypermail 2.1.8 : Thu Jun 09 2011 - 14:56:48 MDT