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

[Agud Pique, Enric]

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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160429/73381805/attachment.html