Re: further conversation on vertical plots

Hi Kara,

I think the issue is with 'wrf_contour'

I don't work with WRF much so I think this will have to be
answered by the wrfhelp@ucar.edu or ncl-talk@ucar.ed communities.

D

On 6/26/12 6:20 PM, Kara Yedinak wrote:
> 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
> <http://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
> <mailto: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 21:39:43 2012