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
This archive was generated by hypermail 2.1.8 : Fri Jun 29 2012 - 15:46:14 MDT