Re: WRF Cross-section Tick Marks Issue. Please Help.

Hi A.J.,
All questions about the wrf_* functions should be posted to the WRF-help
email list: wrfhelp@ucar.edu <mailto:wrfhelp@ucar.edu> . Those functions
were contributed by the WRF group and are supported by the folks on
WRF-help.

I tried to dig around in the WRF*.ncl files but didn't see anything that
would help, and I am not a WRF-data user. Hopefully wrf-help can assist..
Adam

On 12/17/2013 01:17 PM, A.J. Eiserloh wrote:
> Hi all,
>
> I am trying to create appropriate tick marks for my vertical WRF 2-d
> cross-section. The unit for my y-axis is kilometers. There are 100
> points in my vertical dimension in my cross section (dimension 0).
> When I try to tell NCL to put my tick marks from 0 km to 12 km, it
> seems like it is placing the tick marks at the dimension locations
> instead of the actual height data itself. How can I tell NCL to place
> the tick marks at the actual data instead of the points? I want my
> y-limit to be 0-12km.
>
>
>
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRF_contributed.ncl"
>
> begin
>
> case ="run3.1"
> cmd ="ls -1 /data2/eiserloh/data/AR/" + case + "/wrfout_d02_2012*"
> ;frost
>
> files = systemfunc(cmd) + ".nc"
> a = addfiles(files,"r")
> times = wrf_user_list_times(a) ; get times in the files
> ntimes = dimsizes(times)
>
> pt_a= wrf_user_ll_to_ij(a[0],-126.0,42.0,True)
> print(pt_a)
> pt_b= wrf_user_ll_to_ij(a[0],-121.8,35.0,True)
> print(pt_b)
> plane = new(4,float)
> plane=(/ pt_a(0),pt_b(0),pt_a(1),pt_b(1) /)
>
>
> type = "ps"
>
> FirstTime=True
>
> do it=0,1
> ;do it=0,ntimes-1
>
>
> print("Working on time: " + times(it))
> if (it.lt.10) then
> wks = gsn_open_wks(type,"cxsect_line1"+"_00"+tostring(it))
> end if
> if ((it.ge.10).and.(it.lt.100)) then
> wks = gsn_open_wks(type,"cxsect_line1"+"_0"+tostring(it))
> end if
> if (it.ge.100) then
> wks = gsn_open_wks(type,"cxsect_line1"+"_"+tostring(it))
> end if
>
> ; Extract cross section from point A to point B, as defined in "plane"
> ; And vertically interpolate to either height or p coordinates
>
> z = wrf_user_getvar(a[it],"z",0)
> z=z/1000. ;convert to km
>
> rh = wrf_user_getvar(a[it],"rh",0)
> u = wrf_user_getvar(a[it],"ua",0)
> v = wrf_user_getvar(a[it],"va",0)
> ws = ((u^2)+(v^2))^(0.5)
> ws@description = "windspeed"
> ws@units = "m/s"
>
> rh_plane = wrf_user_intrp3d(rh,z,"v",plane,0.,True) ;True if cross
> sect from ptA to ptB
> rh_plane = smth9_Wrap(rh_plane,0.5,0.2,True)
> printVarSummary(rh_plane)
>
>
> ws_plane = wrf_user_intrp3d(ws,z,"v",plane,0.,True) ;windspeed plane
> ws_plane = smth9_Wrap(ws_plane,0.5,0.2,True)
> ws_plane = smth9_Wrap(ws_plane,0.5,0.2,True)
>
> res = True
> res@gsnDraw = False
> res@gsnFrame= False
>
> ; Options for XY Plots
> opts_xy = res
> opts_xy@tiYAxisString = "Height (km)"
> opts_xy@cnMissingValPerimOn = True
> opts_xy@cnMissingValFillColor = 0
> opts_xy@cnMissingValFillPattern = 11
> opts_xy@tmYLMode = "Explicit"
> opts_xy@tmYLValues = fspan(0.0,12.0,13)
> ; Create tick marks
> opts_xy@tmYLMinorValues = fspan(0.5,11.5,12)
> opts_xy@tmYLLabels = sprintf("%.1f",fspan(0.0,12.0,13))
> ; Create labels
> opts_xy@tiXAxisFontHeightF = 0.015
> opts_xy@tiYAxisFontHeightF = 0.015
> opts_xy@tmXBMajorOutwardLengthF = 0.0
> opts_xy@tmYLMajorLengthF = 0.012
> opts_xy@tmYLLabelFontHeightF = 0.011
> opts_xy@tmXBLabelFontHeightF = 0.011
> opts_xy@PlotOrientation = ws_plane@Orientation
> opts_xy@Footer = False
> opts_xy@MainTitle = tostring(case)
> opts_xy@MainTitlePos = "Center"
> opts_xy@InitTime = False
> opts_xy@ValidTime = True
> opts_xy@TimeLabel = times(it)
>
> ; Plotting options for RH
> opts_rh = opts_xy
> opts_rh@ContourParameters = (/ 10., 90., 10. /)
> opts_rh@pmLabelBarOrthogonalPosF = -0.185
>
> opts_rh@lbLabelOffsetF = 0.02 ;determines amount of
> space between actual label and labelbar
> opts_rh@lbBoxMinorExtentF = 0.12 ; perpendicular width
> of labelbar
> opts_rh@lbLabelFontHeightF = 0.01
> opts_rh@lbTitleFontHeightF = 0.01
>
> opts_rh@cnFillOn = True
> opts_rh@cnFillColors = (/"White","White","White", \
> "White","Chartreuse","Green", \
> "Green3","Green4", \
> "ForestGreen","PaleGreen4"/)
>
> ; Plotting options for windspeed
>
> opts_ws = opts_xy
> opts_ws@ContourParameters =(/5.0,100.0,5.0/)
> opts_ws@cnFillOn = False
> opts_ws@cnLineLabelBackgroundColor = -1
> opts_ws@cnInfoLabelOn = False
> opts_ws@cnLineLabelPerimOn = False
> opts_ws@cnLineLabelDensityF = -3.0
> opts_ws@cnLineLabelFontHeightF = 0.009
> ; opts_ws@cnSmoothingOn = True
> ; opts_ws@cnSmoothingTensionF = 0.02
>
> ;Get the contour info for the rh and temp
> contour_ws = wrf_contour(a[it],wks,ws_plane,opts_ws)
> contour_rh = wrf_contour(a[it],wks,rh_plane,opts_rh)
>
>
> ; MAKE PLOTS
> pltres = True
> pltres@Footer = False
> plot = wrf_overlays(a[it],wks,(/contour_rh,contour_ws/),pltres)
>
> delete(wks)
> end do
>
> end
>
> --
> Arthur J. Eiserloh, Jr.
> San Jose State University
> Graduate Student
> Dept. of Meteorology and Climate Science
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk

-- 
______________________________________________________________
Adam Phillips                                asphilli@ucar.edu
NCAR/Climate and Global Dynamics Division       (303) 497-1726
P.O. Box 3000				
Boulder, CO 80307-3000    http://www.cgd.ucar.edu/staff/asphilli

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Dec 17 14:22:00 2013