[ncl-talk] masking only over Myanmar

Karin Meier-Fleischer meier-fleischer at dkrz.de
Tue Sep 6 00:20:57 MDT 2016


Hi Gaeun,

you have to set the mpDataBaseVersion resource to Ncarg4_1 or the alias MediumRes when 
accessing the MapPlot database Earth..4.  See 
http://www.ncl.ucar.edu/Document/HLUs/Classes/MapPlotData4_1_earth_4.shtml

Insert the following line to your script 

  res3 at mpDataBaseVersion = "MediumRes"


Bye,
Karin


> Am 06.09.2016 um 03:11 schrieb 김가은 <gaeun.kim0921 at gmail.com>:
> 
> 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 <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
> 
> _______________________________________________
> 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/20160906/172fe79a/attachment.html 


More information about the ncl-talk mailing list