[ncl-talk] Overlay a box and a marker over WRF field

Karin Meier-Fleischer meier-fleischer at dkrz.de
Wed May 27 04:47:23 MDT 2015


Hi,

take a look at the web page  https://www.ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml

—>
If you want to use gsn_add_polyxxxx routines (like gsn_add_polymarker <https://www.ncl.ucar.edu/Document/Graphics/Interfaces/gsn_add_polymarker.shtml>) to add primtives to a WRF plot, then you need to set:

plt_res at PanelPlot = True
plt_res at FramePlot = False
before you call wrf_map_overlays. Then, call the appropriate gsn_add_polyxxxx routine to attach the desired primitives, and finally call draw <https://www.ncl.ucar.edu/Document/Functions/Built-in/draw.shtml>(plot) and frame <https://www.ncl.ucar.edu/Document/Functions/Built-in/frame.shtml>(wks) to draw the plot with the attached primitives and advance the frame.


Bye,
Karin



> Am 27.05.2015 um 05:51 schrieb Alex Schaefer <alexander.schaefer at mines.sdsmt.edu>:
> 
> Rabah,
> 
> So if you don’t put them both both res and pltres you still end up with 2 plots?
> 
> It might also be worth trying to use the gsn_add_polymarker as opposed to using the gsn-polymarker. gsn_polymarker might have the frame internally built. The gsn_add_polymarker does return an identifier so you would do it similar to your polyline call:
> 
> pm_id = gsn_add_polymarker(wks,plot,npres at ref_lon,npres at ref_lat,mpres)
> 
> If this doesn’t work can you attach your geo_em file and script and I will see if I can see what is going on.
> 
> -Alex
> 
> 
> 
> On May 26, 2015, at 7:06 PM, Rabah Hachelaf <hachelaf at sca.uqam.ca <mailto:hachelaf at sca.uqam.ca>> wrote:
> 
>> Hi Alex ,
>> 
>> When i add 
>>  opts at gsnDraw = False
>>  opts at gsnFrame = False   
>>  or 
>>  res at gsnDraw = False
>>  res at gsnFrame = False   
>> 
>> and also 
>> 
>>  pltres at gsnDraw = False
>>  pltres at gsnFrame = False 
>> 
>> I get only grid with the box and the marker like attached. 
>> 
>> 
>> 
>> 2015-05-26 20:47 GMT-04:00 Alex Schaefer <alexander.schaefer at mines.sdsmt.edu <mailto:alexander.schaefer at mines.sdsmt.edu>>:
>> Rabah,
>> 
>> Where in your script did you add the resources suggested by Jared?
>> 
>> Since you have a bunch of different resource objects you will need to set it in each going into plots controlling plot resources. For you this should include res or opts and likely pltres.
>> 
>> -Alex
>> 
>> 
>> On May 26, 2015, at 6:36 PM, Rabah Hachelaf <hachelaf at sca.uqam.ca <mailto:hachelaf at sca.uqam.ca>> wrote:
>> 
>>> Hi Jared,
>>> 
>>> Thanks for your reponse,
>>> 
>>> But add  
>>> 
>>> 
>>> res at gsnDraw = False
>>> res at gsnFrame = False
>>> 
>>> doesn't change any things. 
>>> 
>>> Rabah
>>> 
>>> 
>>> 2015-05-26 18:34 GMT-04:00 Jared Lee <jaredlee at ucar.edu <mailto:jaredlee at ucar.edu>>:
>>> Hi Rabah,
>>> 
>>> I believe what you're missing is these two lines where you set other 'res' resources:
>>> 
>>> res at gsnDraw = False
>>> res at gsnFrame = False
>>> 
>>> By default they're both set to True, which would give you two separate images. If you set both to false you should get a single plot with your polymarker overlaid.
>>> 
>>> Hope that helps.
>>> 
>>> Jared
>>> 
>>> On Tue, May 26, 2015 at 4:21 PM, Rabah Hachelaf <hachelaf at sca.uqam.ca <mailto:hachelaf at sca.uqam.ca>> wrote:
>>> 
>>> 
>>> Hello NCL users,
>>> 
>>> I am using the script below to plot a box and a marker dot over topography field from WRF outputs,
>>> 
>>> My problem is i can not get one figure, i get two separate figures like attached. 
>>> 
>>> Any suggestion to solve this.
>>> 
>>> 
>>> 
>>> 
>>> 
>>> ;   Example script - plot terrain (using some basic WRF scripts)
>>> ;   Overwrite basic map background settings
>>> ;   November 2008
>>> 
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>>> 
>>> begin
>>> 
>>>   a = addfile("../geo_em.d02.nc <http://geo_em.d02.nc/>","r")           ; Open a file
>>>   filename = "../namelist.wps"
>>> 
>>>   npres= True
>>>   npres at ref_lat   = stringtofloat(systemfunc("grep ref_lat    " +filename+ " | cut -f2 -d'=' | cut -f1 -d','" ) )
>>>   npres at ref_lon   = stringtofloat(systemfunc("grep ref_lon    " +filename+ " | cut -f2 -d'=' | cut -f1 -d','" ) )
>>> 
>>>  type = "png"
>>>   wks = gsn_open_wks(type,"plt_geo_6")         ; Create a plot workstation
>>> 
>>> 
>>>   opts = True                                  ; Set some Basic Plot options
>>>   opts at MainTitle = "GEOGRID FIELDS"
>>>   opts at InitTime = False                        ; Do not plot time or footers
>>>   opts at Footer = False    
>>> 
>>>   ter = wrf_user_getvar(a,"HGT_M",0)           ; Read the variable to memory
>>>   lon = a->XLONG_U(0,:,:)
>>>   lat = a->XLAT_U(0,:,:)
>>>  printVarSummary(lat) 
>>>  
>>>  print (lat(0,0)+", "+lon(0,0)+",0,")
>>>  print (lat(0,60)+", "+lon(0,60)+",0,")
>>>  print (lat(59,60)+", "+lon(59,60)+",0,")
>>>  print (lat(59,0)+", "+lon(59,0)+",0")
>>>  
>>>  
>>>   res = opts                                   ; Use basic options for this field
>>>   res at cnFillOn = True                          ; Create a color fill plot
>>>   res at ContourParameters = (/ 0., 1000., 50. /) ; Set the levels
>>> 
>>>   contour = wrf_contour(a,wks,ter,res)
>>> 
>>>   pltres = True                                ; Set plot options
>>>   mpres = True                                 ; Set map options
>>>   mpres at mpGeophysicalLineColor      = "Black"  ; Overwrite basic map settings
>>>   mpres at mpGridLineColor             = "Black"
>>>   mpres at mpLimbLineColor             = "Black"
>>>   mpres at mpNationalLineColor         = "Black"
>>>   mpres at mpPerimLineColor            = "Black"
>>>   mpres at mpUSStateLineColor          = "Black"
>>>      mpres at mpProjection = "LambertConformal"
>>>        mpres at mpLambertParallel1F = 30
>>>        mpres at mpLambertParallel2F = 60 
>>>        mpres at mpLambertMeridianF = -98
>>>        mpres at mpGridSpacingF     = 0.25
>>> ;************************************************
>>> ; add the box
>>> ;************************************************
>>> ;************************************************
>>>  ypts = (/ lat(10,10), lat(10,50), lat(49,50),  lat(49,10),  lat(10,10)/)
>>>  xpts = (/ lon(10,10), lon(10,50), lon(49,50),  lon(49,10),  lon(10,10)/)
>>> ;************************************************
>>> 
>>>   resp                  = True                      ; polyline mods desired
>>>   resp at gsLineColor      = "black"                     ; color of lines
>>>   resp at gsLineThicknessF = 2.0                       ; thickness of lines
>>> ;  resp at gsLineLabelString= "test"                    ; adds a line label string
>>>    
>>> 
>>>        
>>>   plot = wrf_map_overlays(a,wks,(/contour/),pltres,mpres) ; Plot field over map background
>>>   
>>>   
>>>   dum = new(4,graphic) 
>>>   do i = 0 , 3
>>>     dum(i)=gsn_add_polyline(wks,plot,xpts(i:i+1),ypts(i:i+1),resp)      
>>>   end do 
>>> ;********************************************************* 
>>>   mpres at gsMarkerColor = "Black"
>>>   mpres at gsMarkerIndex = 16
>>>   mpres at gsMarkerSizeF = 0.005
>>>   gsn_polymarker(wks,plot,npres at ref_lon,npres at ref_lat,mpres)
>>>    
>>>   draw(plot)
>>>   frame(wks)
>>>  
>>> end
>>> 
>>> -- 
>>> ------------------------------
>>> Cordialement,
>>> Best regards,
>>> Rabah Hachelaf 
>>> 
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>>> 
>>> 
>>> 
>>> 
>>> -- 
>>> ===============================
>>> Jared A. Lee, Ph.D.
>>> Project Scientist I
>>> Research Applications Laboratory
>>> National Center for Atmospheric Research
>>> Boulder, Colorado, USA
>>> 
>>> Email: jaredlee at ucar.edu <mailto:jaredlee at ucar.edu> (w)
>>> Phone: 303.497.8485 (w)
>>> Web: https://staff.ucar.edu/users/jaredlee <https://staff.ucar.edu/users/jaredlee>
>>> ===============================
>>> 
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>>> 
>>> 
>>> 
>>> 
>>> -- 
>>> ------------------------------
>>> Cordialement,
>>> Best regards,
>>> Rabah Hachelaf 
>>> _______________________________________________
>>> ncl-talk mailing list
>>> List instructions, subscriber options, unsubscribe:
>>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk <http://mailman.ucar.edu/mailman/listinfo/ncl-talk>
>> 





-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150527/93a9860d/attachment.html 


More information about the ncl-talk mailing list