WRF Output - Axes Trouble

I ran a 3D idealized simulation and I'm trying to plot a vertical cross
section of the winds for a specific z range.

I've created a height variable ('z') using the wrf_user_getvar and when I
check the min and max values and they match my inputs, max = 29.7 km and
min = 0.1 km.

When I try to plot the plane as u versus z though my vertical axis will
only extend to 100.

I stripped the script of anything that I thought might be influencing the
y-axis but I haven't had any luck and I'm still getting the same results.

Attached is the resulting plot and below is my script.

Any suggestions would be welcome and thanks!
Lauren Wheeler

*NCL Script:*

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;
;Plots Vertical and Horizontal winds from 3D Hill Simulation
;Feb 2014
;Modified from ;
http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/wrf_Hill2d.ncl
;
;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
;load "./WRFUserARW.ncl"

begin

  a = addfile("wrfout_d01_half-01-01_00:00:00"+".nc","r")

   type = "pdf" ;x11, ps, ncgm, pdf

  wks = gsn_open_wks(type,"Wind_XSection_test")

; Set some Basic Plot options
    res = True
    res@MainTitle = "WRF Hill 3D"
    res@InitTime = False
    res@Footer = False

    pltres = True

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

; Select time outputs for plots

  times = wrf_user_list_times(a) ; get times in the file
  ntimes = dimsizes(times) ; number of times in the file

; The plane data is plotted on

  plane = (/ 0., 50./) ; (x,y) point for vertical plane
  angle = 90.0
  pii = 3.14159
  aspect_ratio = .3 ;0.7 fits on vertical page

; This is the big loop over all of the time periods to process

  do it = 2,ntimes-1,2 ;every other time step
; do it = 1,ntimes-1 ;every time step

    time = it
    res@TimeLabel = times(it)
 ; res@AspectRatio = aspect_ratio

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Collecting data/variables from WRF output

  u_in = wrf_user_getvar(a,"ua",time) ; ua is u averaged to mass points
  v_in = wrf_user_getvar(a,"va",time) ; va is v averaged to mass points
  w_in = wrf_user_getvar(a,"wa",time) ; vertical velocity
  z_in = wrf_user_getvar(a, "z",time) ; grid point height

  u = u_in(0:119,:,:)
  v = v_in(0:119,:,:)
  w = w_in(0:119,:,:)
  z = z_in(0:119,:,:)

  u_plane = wrf_user_intrp3d( u,z,"v",plane,angle,False)
  v_plane = wrf_user_intrp3d( v,z,"v",plane,angle,False)
  w_plane = wrf_user_intrp3d( w,z,"v",plane,angle,False)
  Du_plane = (u_plane-10) ;subtracting input u0 (10 m/s) from u

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

; Horizontal Wind Perturbation Contour Fill
      opts_th = res
      opts_th@FieldTitle = u_in@description
      opts_th@cnFillOn = True
      opts_th@gsnSpreadColorEnd = -10
      opts_th@PlotOrientation = u_plane@Orientation
      contour_th = wrf_contour(a,wks,Du_plane,opts_th) ; plots
horizontal perturbation

 ; Vertical Velocity
      opts_w = res
      opts_w@FieldTitle = w_in@description
      opts_w@cnLineLabelInterval = 0
      opts_w@cnLevelSpacingF = 0.01
      contour_w = wrf_contour(a,wks, w_plane,opts_w)

      plot = wrf_overlays(a,wks,(/contour_th, contour_w/),pltres)

; ************************************************************

  end do ; end of the time loop

end

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk