Re: further conversation on vertical plots

Hi Dennis,

  Thank you for your help. This was a resource was not aware of. I am
working on creating a vertical (curtain) plot using wrf_user_intrp3d and I
am running into an odd problem I could use help troubleshooting. When I
execute my script (see below) I get vertical (z - in km) values of 1-99. I
looked that the data from z and I should be going much higher than that
(model top = 4km). Any advice on where I am going wrong would be greatly
appreciated.

Thanks,
Kara

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

begin
;;;;;;;;;;;;;;
; Read in data
;;;;;;;;;;;;;;
  J = 99
  Ti = 20
  Out = "x11"
  Title = "base_unstable_64A"

  f =
addfile("/mnt/pvfs2/karay/Rath_SA/StabilitySound_xy4/"+Title+"/wrfout_d01_2006-02-23_12:00:
00.nc" , "r")
  qv = f->QVAPOR(Ti,:,:,:) ;(t,k,j,i) Water vapor
  z = wrf_user_getvar(f, "z",Ti) ; grid point height
  zk = z/1000
;************************
;Select Area of interest
;************************
  dimsqv = dimsizes (qv)

  plane = (/ dimsqv(2)/2,dimsqv(1)/2 /)
  opts = False

  qv_plane = wrf_user_intrp3d(qv,zk,"v",plane,90.0,opts)

;************************
;Plot Resources
;************************

  wks = gsn_open_wks(Out,Title)

  gsn_define_colormap(wks,"hotres") ; choose a colormap used hotres
before
  res = True ; plot mods desired
  res@lbLabelBarOn = False ; turn off individual label bars
  res@cnLineLabelsOn = False
  res@cnInfoLabelOn = False
  res@gsnDraw = False ; don't draw
  res@gsnFrame = False ; don't adnvance frame
  res@cnFillOn = True ; turn on color
  res@gsnSpreadColors = True ; use full range of colormap
  res@cnLinesOn = False ; turn off contour lines
  res@lbLabelAngleF = 45 ; angle labels
  res@lbLabelFontHeightF =.022 ; make labels larger
  res@lbTitleOn = True ; turn on title
  res@lbTitleFontHeightF = .015 ; make title smaller
  res@pmLabelBarOrthogonalPosF = .10 ; move whole thing down

;************************
;Make Plot
;************************

  plot = wrf_contour(f,wks,qv_plane,res)
  draw(plot)
  frame(wks)

end

On Fri, Jun 22, 2012 at 6:36 AM, Dennis Shea <shea@ucar.edu> wrote:

> I am not sure if this was answered.
>
> http://www.ncl.ucar.edu/**Document/Functions/WRF_arw/**
> wrf_user_intrp3d.shtml<http://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_user_intrp3d.shtml>
>
> http://www.mmm.ucar.edu/wrf/**OnLineTutorial/Graphics/NCL/**
> Examples/CROSS_SECTION/wrf_**CrossSection2.htm<http://www.mmm.ucar.edu/wrf/OnLineTutorial/Graphics/NCL/Examples/CROSS_SECTION/wrf_CrossSection2.htm>
>
>
> --
"Don't knock the weather; nine-tenths of the people couldn't start a
conversation if it didn't change once in a while." ~Kin Hubbard

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jun 26 18:20:19 2012