为什么地图和等值线填图没法叠加

红儿老子

不知道为什么我画出来的地图和等值线填图没法叠加，ncl新手求大神看看脚本哪里不对，或者给个链接也行，多谢～

1. 
   
2. 
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
   
5. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
6. load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
   
7. 
   
8. begin
   
9. ;
   
10. ; The WRF ARW input file.  
    
11. ; This needs to have a ".nc" appended, so just do it.
    
12.   a = addfile("wrfout_d01_2008-07-22_19:00:00.nc","r")
    
13.   b = addfile("wrfout_d01_2008-07-22_20:00:00.nc","r")
    
14. 
    
15. 
    
16. ; We generate plots, but what kind do we prefer?
    
17.   ;type = "x11"
    
18. ; type = "pdf"
    
19. type = "png"
    
20. ;type="ps"
    
21. ; type = "ncgm"
    
22.   wks = gsn_open_wks(type,"plt_Precip")
    
23. 
    
24. ; Set some basic resources
    
25.   res = True
    
26.   res@MainTitle = "REAL-TIME WRF"
    
27. 
    
28. 
    
29.   pltres = True
    
30.   mpres = True
    
31. 
    
32. 
    
33. 
    
34. 
    
35. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    
36. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    
37. 
    
38. 
    
39. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    
40. ;First get the variables we will need        
    
41. 
    
42. 
    
43. 
    
44.   ; Get non-convective, convective and total precipitation
    
45.   ; Calculate tendency values   
    
46.   ;-------------------------------------------------------------------------------------------------------------------                           
    
47.     rain_exp1 = wrf_user_getvar(a,"RAINNC",0)
    
48.     rain_con1 = wrf_user_getvar(a,"RAINC",0)
    
49.     rain_exp2 = wrf_user_getvar(b,"RAINNC",0)
    
50.     rain_con2 = wrf_user_getvar(b,"RAINC",0)
    
51.     rain_exp=rain_exp2-rain_exp1
    
52.     rain_con=rain_con2-rain_con1
    
53.     rain_tot = rain_exp + rain_con
    
54.     rain_tot@description = "Total Precipitation"
    
55. 
    
56.   ;------------------------------------------------------------------------------------------------------------ -------------------
    
57. 
    
58.       ; Plotting options for Precipitation
    
59.         opts_r = res                        
    
60.         opts_r@UnitLabel            = "mm"
    
61.         opts_r@cnLevelSelectionMode = "ExplicitLevels"
    
62.         opts_r@cnLevels             = (/ .1, .2, .4, .8, 1.6, 3.2, 6.4, \
    
63.                                         12.8, 25.6, 51.2, 102.4/)
    
64.         opts_r@cnFillColors         = (/"White","White","DarkOliveGreen1", \
    
65.                                         "DarkOliveGreen3","Chartreuse", \
    
66.                                         "Chartreuse3","Green","ForestGreen", \
    
67.                                         "Yellow","Orange","Red","Violet"/)
    
68.         opts_r@cnInfoLabelOn        = False
    
69.         opts_r@cnConstFLabelOn      = False
    
70.         opts_r@cnFillOn             = True
    
71.    
    
72.         lat  = wrf_user_getvar(a,"XLAT",0)
    
73.         lon = wrf_user_getvar(a, "XLONG", 0)
    
74. ;----------------------------------------------------------------------
    
75. ; Plot full domain first.
    
76. ;----------------------------------------------------------------------
    
77.         opts_r@cnFillOn=True
    
78.         contour = wrf_contour(a,wks,rain_tot,opts_r)
    
79.      
    
80. ;----------------------------------------------------------------------
    
81. ; Plot partial domain.
    
82. ;----------------------------------------------------------------------
    
83.   opts_r@sfXArray = lon
    
84.   opts_r@sfYArray = lat
    
85.   contour = wrf_contour(a,wks,rain_tot,opts_r)
    
86. 
    
87. ;---Set special resource to indicate we are using XLAT/XLONG coordinates.
    
88.   pltres@LatLonOverlay = True
    
89. 
    
90. ;---Zoom in on map, which we can do because we're using lat/lon coordinates.
    
91.   mpres@mpLeftCornerLatF  =  28.
    
92.   mpres@mpRightCornerLatF =  34.
    
93.   mpres@mpLeftCornerLonF  = 110.
    
94.   mpres@mpRightCornerLonF = 122.
    
95. 
    
96. 
    
97. contour_tot = wrf_contour(a,wks, rain_tot, opts_r)
    
98.       
    
99. delete(opts_r)
    
100. 
     
101. 
     
102. ; MAKE PLOTS                                       
     
103. 
     
104. plot = wrf_map_overlays(a,wks,contour_tot,pltres,mpres)
     
105.       
     
106. ;--------------------------------------------------------------------------------------------------------------------------------------------
     
107. 
     
108. ;---------------------------------------------------
     
109. ;         add China map
     
110. ;----------------------------------------------------
     
111. 
     
112.   map=wrf_map(wks, a, mpres)
     
113.   shp_name1="/usr/local/NCL/lib/ncarg/nclscripts/cnmap/cnmap.shp"
     
114.   lnres1=True
     
115.   lnres1@gsLineColor      = "black"
     
116.   lnres1@gsLineThicknessF = 2.0   
     
117.   poly1=gsn_add_shapefile_polylines(wks, map, shp_name1,lnres1)
     
118. 
     
119. 
     
120.   shp_name2="/usr/local/NCL/lib/ncarg/nclscripts/cnmap/river.shp"
     
121.   lnres2=True
     
122.   lnres2@gsLineThicknessF = 2.0      
     
123.   lnres2@gsLineColor = "blue"
     
124.    poly2=gsn_add_shapefile_polylines(wks, map, shp_name2, lnres2)
     
125. 
     
126.   shp_name3="/usr/local/NCL/lib/ncarg/nclscripts/cnmap/cnhimap.shp"
     
127.   lnres3=True
     
128.   lnres3@gsLineThicknessF = 2.0      
     
129.   lnres3@gsLineColor = "black"
     
130.    poly3=gsn_add_shapefile_polylines(wks, map, shp_name3, lnres3)
     
131.    draw(map)
     
132.    frame(wks)
     
133. 
     
134. 
     
135. 
     
136. end
     

复制代码

红儿老子

是后面加地图的时候出了问题。已解决～

aa22556688

ncl两幅图叠加需要overlay吧

jiangjunxia

楼主解决方案没贴出来啊

Eudora99

pltres = True 后加
pltres@FramePlot = False   ; do not frame plot - will do this manually later

李子

{:5_231:}{:5_231:}{:5_231:}{:5_231:}

李子

{:lxm_25:}{:lxm_25:}{:lxm_25:}{:lxm_25:}

李子

{:lxm_26:}{:lxm_26:}{:lxm_26:}{:lxm_26:}{:lxm_26:}

I彭彭Ed

按楼主的方法改了，是叠加在一起了，可是图的大小不一样
请问怎么改变大小啊

wawatou83

占位…………