[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