wrf_user_interp3d documentation

From: <sam.hawkins_at_nyahnyahspammersnyahnyah>
Date: Tue Oct 15 2013 - 09:28:29 MDT

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.

; 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:
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