[ncl-talk] control size of domain in a map

Rick Brownrigg brownrig at ucar.edu
Fri Apr 29 14:37:08 MDT 2016


As I realized you are working with wrf mapping functions, my previous
advise on mpLimitMode might not wotk.  If that's the case, you might take a
look at the "Zoomin" resource, described in the docs for wrf_map_overlays:

http://ncl.ucar.edu/Document/Functions/WRF_arw/wrf_map_overlays.shtml

On Fri, Apr 29, 2016 at 2:32 PM, Rick Brownrigg <brownrig at ucar.edu> wrote:

> I'm not certain, but you may need to set
>
> mpres at mpLimitMode = "LatLon"
>
> as the default for that resource is "MaximalArea".  Take a look at the
> docs for what this all means:
>
> http://ncl.ucar.edu/Document/Graphics/Resources/mp.shtml#mpLimitMode
>
> Hope that helps...
> Rick
>
>
> On Fri, Apr 29, 2016 at 1:57 PM, Agud Pique, Enric <eagud.q at ccma.cat>
> wrote:
>
>>
>>
>>
>> Hello people,
>>
>> I am quite new user of NCL,
>>
>> I use WRF with map_proj = 'lambert'...when I run the following ncl
>> script, I cannot get the domain I really want. I describe the left and
>> right longitude, the low and high latitude  I wish, but I get always the
>> big domain...I supose I have something wrong in the code...
>>
>>
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
>> load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
>> begin
>> a = addfile("./wrfout_d01_2016-04-28_18:00:00","r")
>> wks = gsn_open_wks("X11","plt_Surface")
>>
>>
>> T2 = wrf_user_getvar(a,"T2",0)
>> Tc = T2-273
>> pltres = True
>> mpres = True
>> opts = True
>>
>>
>> mpres at mpMaxLatF  = 44
>> mpres at mpMinLatF  = 40
>> mpres at mpMaxLonF  = 4
>> mpres at mpMinLonF  = 0
>>
>>
>> opts at cnFillOn = True
>> opts at cnLevelSelectionMode = "ExplicitLevels"
>> opts at cnLevels   =(/ -10, -5,0,5,10,15,20,25,30,35,40/)
>> opts at cnFillColors =(/"Blue","DarkOliveGreen1", \
>>                       "DarkOliveGreen3","Chartreuse", \
>>                       "Chartreuse3","Green","ForestGreen", \
>>                       "Yellow","Orange","Red","Violet"/)
>> contour_t2 = wrf_contour(a,wks,Tc,opts)
>>
>>
>> plot= wrf_map_overlays(a,wks,(/contour_t2/),pltres ,mpres)
>>
>>
>> Thank you very much, best regards
>>
>>
>> _______________________________________________
>> ncl-talk mailing list
>> ncl-talk at ucar.edu
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160429/b6e1bafc/attachment.html 


More information about the ncl-talk mailing list