Re: Help: How to add line and makers on the Zoomed domain of WRF output

  Thank you, Adam,
However, I tried the two ways you suggest, neither of them works.
Thank you again.

Kekuan

On 9/17/2010 2:00 PM, ncl-talk-request@ucar.edu wrote:
> Send ncl-talk mailing list submissions to
> ncl-talk@ucar.edu
>
> To subscribe or unsubscribe via the World Wide Web, visit
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> or, via email, send a message with subject or body 'help' to
> ncl-talk-request@ucar.edu
>
> You can reach the person managing the list at
> ncl-talk-owner@ucar.edu
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of ncl-talk digest..."
>
>
> Today's Topics:
>
> 1. Help: How to add line and makers on the Zoomed domain of WRF
> output (Kekuan Chu)
> 2. Re: Help: How to add line and makers on the Zoomed domain of
> WRF output (Adam Phillips)
> 3. Re: Setting contour levels using a min/max contour level and
> a spacing (Mary Haley)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 17 Sep 2010 08:53:41 -0400
> From: Kekuan Chu<kchu@marine.usf.edu>
> Subject: Help: How to add line and makers on the Zoomed
> domain of WRF output
> To: ncl-talk@ucar.edu
> Message-ID:<4C9364D5.9040308@marine.usf.edu>
> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>
> Hello everyone:
> I tried to plot a small area of the WRF domain, and then add some
> line and makers on it.
> However, the fields, line and makers were drawn on separated graphs.
> How to plot them on one graph?
> Thank you very much.
>
> The script is:
> ; Refer to wrf_Zoom.ncl etc.
>
> 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"
>
> begin
> ;********************************************************
> ;
> xp = new((/6,13/),float) ; allocate memory
> yp = new((/6,13/),float)
>
> ;001 cntl
> yp(0,:) =
> (/17.0,17.8,18.5,19.2,19.7,20.1,20.5,20.9,20.8,21.0,21.3,21.6,22.1/)
> xp(0,:) =
> (/125.1,125.1,125.0,124.6,124.3,124.3,124.4,124.2,124.1,124.5,124.6,124.5,124.5/)
>
> ;********************************************************
>
> ; open file
> a = addfile("./wrfinput.nc","r")
>
> ; We generate plots, but what kind do we prefer?
> wks = gsn_open_wks("x11","plt_Zoom")
>
> ; Set the color map
> gsn_define_colormap(wks,"ViBlGrWhYeOrRe"); select color map
>
> ; Set some Basic Plot options
> res = True
> res@MainTitle = "Typhoon"
> pltres = True
> mpres1 = True
>
> var = wrf_user_getvar(a,"U",0)
>
> dims = dimsizes(var)
> ; define the zoomed domain
> x_start = 22
> x_end = 60
> y_start = 0
> y_end = 40
>
> mpres1@ZoomIn = True
> mpres1@Xstart = x_start
> mpres1@Ystart = y_start
> mpres1@Xend = x_end
> mpres1@Yend = y_end
>
> var_zoom = var(:,y_start:y_end,x_start:x_end)
> var_zoom@description = var@description
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> opts = res
>
> opts@cnFillOn = True ; turn on color
> opts@cnLinesOn = False ; turn contour lines
> opts@gsnSpreadColors = True ; use entire color map
> opts@lbLabelAutoStride = True ; let NCL determine
> label spacing
>
> opts@cnLevelSelectionMode = "ManualLevels" ; manual
> contour levels
> opts@cnMinLevelValF = -6. ; set min
> contour level
> opts@cnMaxLevelValF = 6. ; set max
> contour level
> opts@cnLevelSpacingF = 1. ; contour interval
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> ;plot var in zoom area
> contour = wrf_contour(a,wks,tc_zoom(var,:,:),opts)
> plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres1)
>
> ; Add trajectory lines.
> pres = True ; polyline resources
> pres@gsLineThicknessF = 3.0 ; line thickness
> pres@gsLineColor = "black"
> line = gsn_add_polyline(wks,plot,xp(0,:),yp(0,:),pres) ; draw
> the traj
>
> ; Add markers to the trajectories.
> mres = True ; marker resources for best track
> mres@gsMarkerIndex = 16 ; marker style (filled circle)
> mres@gsMarkerSizeF = 8.0 ; marker size
> mres@gsMarkerColor = "black" ; maker color
> markers = gsn_add_polymarker(wks,plot,xp(0,:),yp(0,:),mres)
> draw(plot)
> frame(wks)
>
> ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
>
> end
>

-- 
Kekuan Chu
College of Marine Science,
University of South Florida
727-553-3358(office)
727-238-1462(cell)
School of Atmospheric Sciences,
Key Laboratory of Mesoscale Severe Weather (LMSWE)
Nanjing University, Nanjing, 210093
P. R. China
kkchu@nju.edu.cn
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk