[ncl-talk] masking only over Myanmar
김가은
gaeun.kim0921 at gmail.com
Mon Sep 5 19:11:22 MDT 2016
Hi, me again. always thanks for your big help.
I'm trying to present data only over Myanmar territory
like this example :
https://www.ncl.ucar.edu/Applications/Images/mask_4_lg.png
and this is what I've done:
;---- Source for plot (Myanmar)
fill_specs = (/"water", "Land"/)
mask_specs = (/"Myanmar"/)
res3 = res
res3 at mpDataSetName = "Earth..4"
res3 at mpFillBoundarySets = "NoBoundaries"
res3 at mpFillAreaSpecifiers = fill_specs
res3 at mpSpecifiedFillColors = (/7,2/)
res3 at mpAreaMaskingOn = 1
res3 at mpMaskAreaSpecifiers = mask_specs
res3 at mpMinLonF = 90 ;123 ;25 Horn of Africa
res3 at mpMaxLonF = 105 ;132 ;53
res3 at mpMinLatF = 0 ;32 ;-10
res3 at mpMaxLatF = 30 ;45 ;20
;---- Source for Panel
pres3 = pres
pres3 at txString = "Total Precipitation of ICON (40km) over Myanmar"
;---- Draw panel plot (Myanmar)
res3 at gsnCenterString = "SEP 2016"
plot_3(0) = gsn_csm_contour_map_ce(wks,Total_p(0,:,:),res3)
res3 at gsnCenterString = "OCT 2016"
plot_3(1) = gsn_csm_contour_map_ce(wks,Total_p(1,:,:),res3)
res3 at gsnCenterString = "NOV 2016"
plot_3(2) = gsn_csm_contour_map_ce(wks,Total_p(2,:,:),res3)
res3 at gsnCenterString = "SON 2016"
plot_3(3) = gsn_csm_contour_map_ce(wks,Total_p(3,:,:),res3)
gsn_panel(wks,plot_3,(/2,2/),pres3)
the result is ...
warning:MapV40DHUpdateDrawList: invalid boundary specification string:
"myanmar"
warning:MapV40DHUpdateDrawList: invalid boundary specification string:
"myanmar"
warning:MapV40DHUpdateDrawList: invalid boundary specification string:
"myanmar"
warning:MapV40DHUpdateDrawList: invalid boundary specification string:
"myanmar"
could anyone give me an advice?
Regards,
Gaeun
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160906/4e90502a/attachment.html
More information about the ncl-talk
mailing list