Maybe I'm missing something here, but I'm finding the documentation of wrf_user_intrp3d confusing. The docs say:
function wrf_user_intrp3d (
var3d : numeric,
vert : numeric,
plot_type : string,
loc [*] : numeric, ; up to four values
angle : numeric,
res [1] : logical
)
var3d: Data on model levels that will be interpolated. The rightmost dimensions of this array is nz x ny x nx.
vert: Array of vertical coordinate to interpolate **to**. This must either be pressure/height. Dimensions must be the same as those of var3d. If vert is pressure, the units may be either hPa or Pa.
plot_type: Plot orientation. Use "h" for horizontal plots, and "v" for vertical cross sections.
Loc: Interpolation information.
Angle: Only valid for cross sections where a plane will be plotted through a given point on the model domain. 0.0 represents a S-N cross section and 90.0 a W-E cross section.
Res: set to True if plotting a cross section from point A to point B; otherwise set to False.
But none of the call signatures in the examples match this. It seems vert should be an array of coordinates to interpolate **from**: i.e.
http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/LEVELS_INTERP/wrf_Height2.ncl
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; First get the variables we will need
tc = wrf_user_getvar(a,"tc",-1) ; T in C
u = wrf_user_getvar(a,"ua",-1) ; u averaged to mass points
v = wrf_user_getvar(a,"va",-1) ; v averaged to mass points
p = wrf_user_getvar(a, "pressure",-1) ; pressure is our vertical coordinate
z = wrf_user_getvar(a, "z",-1) ; grid point height
rh = wrf_user_getvar(a,"rh",-1) ; relative humidity
; The specific height levels that we want the data interpolated to.
; And interpolate to these levels
height_levels = (/ 250., 2000./) ; height levels to plot - in meter
nlevels = dimsizes(height_levels) ; number of height levels
p_plane = wrf_user_intrp3d( p,z,"h",height_levels,0.,False)
tc_plane = wrf_user_intrp3d(tc,z,"h",height_levels,0.,False)
rh_plane = wrf_user_intrp3d(rh,z,"h",height_levels,0.,False)
u_plane = wrf_user_intrp3d( u,z,"h",height_levels,0.,False)
v_plane = wrf_user_intrp3d( v,z,"h",height_levels,0.,False)
u_plane = u_plane*1.94386 ; kts
v_plane = v_plane*1.94386 ; kts
u_plane@units = "kts"
v_plane@units = "kts"
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Oct 15 09:28:51 2013
This archive was generated by hypermail 2.1.8 : Tue Oct 22 2013 - 10:35:27 MDT