Re: WRF overlays

Thanks Mary-
That is what I usually do, use the gsn routines, but I didnt know I could
use them to zoom in on an area. Are the gsn routines able to zoom in to a
portion of the domain while retaining the map?
Thanks for your help.
jimmyc

On Wed, Jun 4, 2008 at 8:03 AM, Mary Haley <haley_at_ucar.edu> wrote:

>
> Jimmy,
>
> The error message you're getting indicates something is probably
> syntactically wrong with your vort.ncl script, and not with
> wrf_overlay per se. It's hard to tell from looking at the code
> snippet, however.
>
> The WRF-NCL plotting scripts are meant to be more black-box, and not
> so much be customizable.
>
> If you want to try to use gsn_csm_xxx functions to do your plotting,
> then there are some examples that are similar to what you want to do.
> See example WRF_lc_8 at:
>
> http://www.ncl.ucar.edu/Applications/wrflc.shtml
>
> This one generates vectors and contours over a map, but you can
> substitute the "gsn_csm_vector_scalar_map" call with a call to
> "gsn_csm_contour_map_overlay" which allows you to overlay two contour
> plots.
>
> To see an actual example of using gsn_csm_contour_map_overlay, see
> example 2 at:
>
> http://www.ncl.ucar.edu/Applications/conOncon.shtml
>
> You still want to refer to the WRF_lc example for using WRF_map_c
> to set up your map resources, but you can use the above example
> for tips on how to set up your two contouring resource lists.
>
> --Mary
>
>
> On Tue, 3 Jun 2008, James Correia wrote:
>
> Hi All-
>> I am trying to overlay some variables using the wrf_contour function. I
>> have
>> completely given up on the whole zooming in feature listed at the WRF-NCL
>> website. It appears these will work in the future NCL release at least
>> thats
>> what I am guessing.
>>
>> My attempt to use wrf_overlay always results in an error. The error is:
>> fatal:syntax error: line xxx in the file vort.ncl before or near
>> wrf_overlay
>> plot = wrf_overlay
>> ______________^
>>
>> Does anybody have an example script using the 5.0.0 NCL release that can
>> accomplish:
>> 1. zooming in to a region of the domain
>> 2. overlay two variables (one with shading and one with contouring)
>> 3. maximize the plot size to fit the page.
>>
>> Thanks
>> jimmyc
>>
>>
>>
>> Due to the length of the code, I have provided my plotting procedure:
>> res = True
>> colors = (/"white","black" \ ; {fore/back}ground
>> ,"white" \
>> ,"orange" \
>> ,"red","violet" \
>> ,"blueviolet","blue" /)
>> gsn_define_colormap(wks, colors)
>> res_at_gsnMaximize = True ; maximize size
>> res_at_gsnSpreadColors = True ; use full range of
>> colormap
>> res_at_gsnRightString = "Vorticity and W"
>> res_at_gsnLeftString = "Wind speed"
>> res_at_mpGeophysicalLineThicknessF = 2.0
>> res_at_mpUSStateLineThicknessF = 2.0
>> ;************************************************
>> ; Use WRF_contributed procedure to set map resources
>> ;************************************************
>> WRF_map_c(f, res, 0) ; reads info from file
>>
>> ;************************************************
>> ; set True for native (direct) mapping (faster plotting)
>> ; set False for non-native mappings
>> ;************************************************
>> res_at_tfDoNDCOverlay = True
>>
>> ;************************************************
>> ; associate the 2-dimensional coordinates to the variable for plotting
>> ; only if not a native projection [ res_at_tfDoNDCOverlay=False ]
>> ;************************************************
>> if (.not.res_at_tfDoNDCOverlay) then
>> lat2d = f->XLAT(0,:,:)
>> lon2d = f->XLONG(0,:,:)
>>
>> u10_at_lat2d = lat2d
>> u10_at_lon2d = lon2d
>> v10_at_lat2d = lat2d
>> v10_at_lon2d = lon2d
>>
>> rainTot_at_lat2d = lat2d
>> rainTot_at_lon2d = lon2d
>> end if
>>
>>
>> nt=64
>> x1 = 120
>> x2 = 210
>> ;x1 = 120
>> ;x2 = 210
>> y1 = 80 ;-20
>> y2 = 140;-20
>> zres = True
>> zres_at_ZoomIn = True
>> zres_at_Xstart = x1
>> zres_at_Xend = x2
>> zres_at_Ystart = y1
>> zres_at_Yend = y2
>> res_at_cnFillMode = "RasterFill" ; raster
>> res_at_cnLevels = (/0.5,1,2,4,6,8,12/)
>> res_at_cnFillOn = True ; color plot desired
>> res_at_cnLinesOn = False ; turn off contour lines
>> res_at_cnLineLabelsOn = False ; turn off contour labels
>> res_at_cnLevelSelectionMode = "ExplicitLevels" ; explicit [unequal] cn
>> levels
>>
>> w = f->W(64,:,:,:)
>> wa = new((/nlevs,224,299/),float)
>> do k=0,33
>> wa(k,:,:) = (w(k,:,:) + w(k+1,:,:))/2.
>> end do
>> delete(w)
>>
>> fill = wrf_contour(f,wks,wa(4,y1:y2,x1:x2),res)
>> ;draw(fill)
>> ;frame(wks)
>>
>> resa = True
>> resa_at_cnLevels = (/-6,-4,-3,-2,-1,-0.5,0.5,1,2,3,4,6/)
>> resa_at_cnFillOn = False ; color plot desired
>> resa_at_cnLinesOn = True ; turn off contour lines
>> resa_at_cnLineLabelsOn = True ; turn off contour labels
>> resa_at_cnLevelSelectionMode = "ExplicitLevels" ; explicit [unequal] cn
>> levels
>>
>> con = wrf_contour(f,wks,vor(nt,y1:y2,x1:x2),resa)
>> plot = wrf_overlay(wks,(/fill,con/),zres)
>>
>> ;also tried using True in place of zres; same problem
>>
>> ;draw(con)
>> ;frame(wks)
>> end
>>
>>
>>
>>
>> --
>> James Correia Jr.
>> Post Doc
>> Climate Physics Group, PNNL
>>
>>

-- 
James Correia Jr.
Post Doc
Climate Physics Group, PNNL

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jun 04 2008 - 11:04:00 MDT