（更新）最最最准确的中国地图（附带南海，精确到市界）

sunyue19901203

尽头的尽头 发表于 2015-3-26 22:05
可以把脚本贴出来么

1. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
   
5. load "/gfs3/wutw/sunyue/ncl/cnmap/cnmap.ncl"
   
6. ;===========================================================================================
   
7. 
   
8. begin
   
9. 
   
10. modelname = (/"Obs","MIROC4h","BCC-CSM1-1-m","EC-EARTH","MRI-CGCM3","HadGEM2-ES","IPSL-CM5A-MR",\
    
11.                "CNRM-CM5","MIROC5","INMCM4","FGOALS-s2","NorESM1-M","IPSL-CM5A-LR","GFDL-ESM2G",\
    
12.                "GFDL-ESM2M","GFDL-CM3","GISS-E2-R","BCC-CSM1-1","MIROC-ESM","MIROC-ESM-CHEM"/)
    
13. nplot = dimsizes(modelname)
    
14. 
    
15. wks   = gsn_open_wks("eps","obs_model_annualpre_perday")
    
16. plot  = new(nplot,graphic)
    
17. do i = 0, nplot-1
    
18. 
    
19.    modelfilepath    = "/gfs3/wutw/sunyue/ncl/model/cmip5/" + modelname(i) +"_00-05_annual_pre_perday.nc"
    
20.    modelfile        = addfile(modelfilepath,"r")
    
21.    prep              = modelfile->pre
    
22.    
    
23.    gsn_define_colormap(wks,"WhiteBlueGreenYellowRed")
    
24.    
    
25.    res                             = True
    
26.    res@gsnDraw                     = False
    
27.    res@gsnFrame                    = False
    
28.    
    
29.   res@cnFillOn                    = True                 ; turn on color
    
30.   res@cnLinesOn                   = False                ; no contour lines
    
31.   res@cnLineLabelsOn              = False
    
32.   res@gsnSpreadColors             = True                 ; use full color map
    
33.   res@gsnLeftString               = modelname(i)
    
34.   res@gsnLeftStringFontHeightF    = 0.03
    
35.   res@gsnAddCyclic                = False  
    
36. 
    
37.   res@cnInfoLabelOn               = False                ; Turn off info label.
    
38.   res@lbLabelBarOn                = False
    
39.   
    
40.   res@cnLevelSelectionMode        = "ManualLevels"       ; manual levels
    
41.   res@cnMinLevelValF              = 1
    
42.   res@cnMaxLevelValF              = 12
    
43. ; res@cnLevelSpacingF             = 1
    
44.   
    
45.   res@tmXBLabelsOn         = False              ; no bottom labels
    
46.   res@tmXBOn               = False              ; no bottom tickmarks
    
47.   res@tmYRLabelsOn         = False              ; no right labels
    
48.   res@tmYROn               = False              ; no right tickmarks
    
49.   res@tmYLLabelsOn         = False              ; do not draw left labels
    
50.   res@tmYLOn               = False              ; no left tickmarks
    
51.   res@tmXTLabelsOn         = False              ; do not draw top labels
    
52.   res@tmXTOn               = False              ; no top tickmarks
    
53. 
    
54.   res@tmXBTickSpacingF            = 10      
    
55.   res@tmYLTickSpacingF            = 10
    
56.   if(i%4.eq.0)
    
57.   res@tmYLOn                      = True         
    
58.   res@tmYLLabelsOn                = True
    
59.   res@tmYLLabelFontHeightF        = 0.023
    
60.   else
    
61.   res@tmYLOn                      = False         
    
62.   res@tmYLLabelsOn                = False            
    
63.   end if
    
64.   if(i.ge.16)
    
65.   res@tmXBLabelsOn                = True              
    
66.   res@tmXBOn                      = True
    
67.   res@tmXBLabelFontHeightF        = 0.023            
    
68.   else  
    
69.   res@tmXBLabelsOn                = False              
    
70.   res@tmXBOn                      = False            
    
71.   end if
    
72.   res@mpOutlineSpecifiers         = (/"China","Taiwan"/)
    
73.   res@mpNationalLineColor         = "Black"
    
74.   res@mpNationalLineThicknessF    = 1
    
75.   
    
76.   res@mpMinLatF            = 15      ; range to zoom in on
    
77.   res@mpMaxLatF            = 50
    
78.   res@mpMinLonF            = 70
    
79.   res@mpMaxLonF            = 130
    
80.   
    
81.   cnres                         =True
    
82.   cnres@china                   =False
    
83.   cnres@river                   =True
    
84.   cnres@province                =False
    
85.   cnres@nanhai                  =False
    
86.   cnres@diqu                    =False
    
87. 
    
88. 
    
89.   
    
90.   plot(i)  = gsn_csm_contour_map_ce(wks,prep,res)
    
91.   chinamap = add_china_map(wks,plot(i),cnres)
    
92.   delete(prep)
    
93. end do
    
94. 
    
95.   resp                            = True
    
96.   resp@gsnMaximize                = True
    
97.   resp@gsnPanelLabelBar           = True
    
98.   resp@lbLabelAutoStride          = True                 ; optimal labels
    
99.   resp@lbLabelFontHeightF         = 0.010
    
100.   resp@gsnPanelRowSpec            = True
     
101.   resp@gsnAttachBorderOn          = False
     
102.   resp@gsnPanelXWhiteSpacePercent = 0.25
     
103.   resp@gsnPanelYWhiteSpacePercent = 3
     
104.   resp@gsnPanelRight              = 0.98
     
105. 
     
106.   gsn_panel(wks,plot,(/4,4,4,4,4/),resp)
     
107. 
     
108. end
     

复制代码

sunyue19901203

尽头的尽头 发表于 2015-3-26 22:05
可以把脚本贴出来么

不好意思，大神，这两天没上家园，才看到，我已经把脚本贴出来了

sunyue19901203

这是最后画出来的图
C:\Users\Sun\Desktop\QOE@TE`Q~4%5DEY_~]3_EAGYB.png

MMMMMMMMMMMars

适合初学者么

MMMMMMMMMMMars

感觉 有点变形了

dtlyr

尽头的尽头 发表于 2015-3-26 22:05
可以把脚本贴出来么

支持，需要金钱，没钱了

Mr.Yes

好看的一腿啊~

tr179667208

找了半天，终于找到了

sunflower$

试一下 看效果如何 再回来感谢楼主

夏天很美