ncl画wrf转出来的数据只有一天的图

空城不空

我所处理的数据是三天的，然后用官网上现成的脚本来画图（ http://www2.mmm.ucar.edu/wrf/OnL ... ECIP/wrf_Precip.htm ）

画出来的图只有第一天的，我用这个脚本去画WRF官网转下来的数据也是只有一天的。
(0)        Working on time: 2017-07-12_00:00:00
(0)        Working on time: 2017-07-12_06:00:00
(0)        Working on time: 2017-07-12_12:00:00
我不太清楚哪里还需要更改。。。
下面这个是ncl脚本

1. ;   Example script to produce plots for a WRF real-data run,
   
2. ;   with the ARW coordinate dynamics option.
   
3. 
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
5. load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
   
6. 
   
7. begin
   
8. ;
   
9. ; The WRF ARW input file.  
   
10. ; This needs to have a ".nc" appended, so just do it.
    
11.   a = addfile("wrfout_d01_2017-07-12_00_00_00","r")
    
12. 
    
13. 
    
14. ; We generate plots, but what kind do we prefer?
    
15. ; type = "x11"
    
16.   type = "pdf"
    
17. ; type = "ps"
    
18. ; type = "ncgm"
    
19.   wks = gsn_open_wks(type,"plt_Precip")
    
20. 
    
21. 
    
22. ; Set some basic resources
    
23.   res = True
    
24.   res@MainTitle = "Heavy Rainfall"
    
25. 
    
26.   pltres = True
    
27.   mpres = True
    
28.   mpres@mpGeophysicalLineColor = "Black"
    
29.   mpres@mpNationalLineColor    = "Black"
    
30.   mpres@mpUSStateLineColor     = "Black"
    
31.   mpres@mpGridLineColor        = "Black"
    
32.   mpres@mpLimbLineColor        = "Black"
    
33.   mpres@mpPerimLineColor       = "Black"
    
34. 
    
35. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    
36. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    
37. 
    
38. ; What times and how many time steps are in the data set?
    
39.   FirstTime = True
    
40.   times = wrf_user_getvar(a,"times",-1)  ; get all times in the file
    
41.   ntimes = dimsizes(times)         ; number of times in the file
    
42. 
    
43. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    
44. 
    
45.   do it = 0,ntimes-1,2             ; TIME LOOP
    
46. 
    
47.     print("Working on time: " + times(it) )
    
48.     if (FirstTime) then            ; Save some times for tracking tendencies
    
49.       times_sav = times(it)
    
50.     end if
    
51.     res@TimeLabel = times(it)   ; Set Valid time to use on plots
    
52. 
    
53. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    
54. ; First get the variables we will need        
    
55. 
    
56.     slp = wrf_user_getvar(a,"slp",it)  ; slp
    
57.       wrf_smooth_2d( slp, 3 )            ; smooth slp
    
58. 
    
59.   ; Get non-convective, convective and total precipitation
    
60.   ; Calculate tendency values                              
    
61.     rain_exp = wrf_user_getvar(a,"RAINNC",it)
    
62.     rain_con = wrf_user_getvar(a,"RAINC",it)
    
63.     rain_tot = rain_exp + rain_con
    
64.     rain_tot@description = "Total Precipitation"
    
65. 
    
66.     if( FirstTime ) then
    
67.       if ( it .eq. 0 ) then
    
68.         rain_exp_save = rain_exp
    
69.         rain_con_save = rain_con
    
70.         rain_tot_save = rain_tot
    
71.       else
    
72.         rain_exp_save = wrf_user_getvar(a,"RAINNC",it-1)
    
73.         rain_con_save = wrf_user_getvar(a,"RAINC",it-1)
    
74.         rain_tot_save = rain_exp_save + rain_con_save
    
75.         FirstTime = False
    
76.         times_sav = times(it-1)
    
77.       end if
    
78.     end if
    
79. 
    
80.     rain_exp_tend = rain_exp - rain_exp_save
    
81.     rain_con_tend = rain_con - rain_con_save
    
82.     rain_tot_tend = rain_tot - rain_tot_save
    
83.     rain_exp_tend@description = "Explicit Precipitation Tendency"
    
84.     rain_con_tend@description = "Param  Precipitation Tendency"
    
85.     rain_tot_tend@description = "Precipitation Tendency"
    
86. 
    
87.   ; Bookkeeping, just to allow the tendency at the next time step
    
88.     rain_exp_save = rain_exp
    
89.     rain_con_save = rain_con
    
90.     rain_tot_save = rain_tot
    
91. 
    
92. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
    
93. 
    
94.     if( .not. FirstTime ) then     ; We will skip the first time
    
95. 
    
96.       ; Plotting options for Sea Level Pressure
    
97.         opts_psl = res         
    
98.         opts_psl@ContourParameters = (/ 900., 1100., 2. /)
    
99.         opts_psl@cnLineColor       = "Blue"
    
100.         opts_psl@cnInfoLabelOn     = False
     
101.         opts_psl@cnLineLabelFontHeightF = 0.01
     
102.         opts_psl@cnLineLabelPerimOn = False
     
103.         opts_psl@gsnContourLineThicknessesScale = 1.5
     
104.         contour_psl = wrf_contour(a,wks,slp,opts_psl)
     
105.         delete(opts_psl)
     
106.    
     
107. 
     
108.       ; Plotting options for Precipitation
     
109.         opts_r = res                        
     
110.         opts_r@UnitLabel            = "mm"
     
111.         opts_r@cnLevelSelectionMode = "ExplicitLevels"
     
112.         opts_r@cnLevels             = (/ .1, .2, .4, .8, 1.6, 3.2, 6.4, \
     
113.                                         12.8, 25.6, 51.2, 102.4/)
     
114.         opts_r@cnFillColors         = (/"White","White","DarkOliveGreen1", \
     
115.                                         "DarkOliveGreen3","Chartreuse", \
     
116.                                         "Chartreuse3","Green","ForestGreen", \
     
117.                                         "Yellow","Orange","Red","Violet"/)
     
118.         opts_r@cnInfoLabelOn        = False
     
119.         opts_r@cnConstFLabelOn      = False
     
120.         opts_r@cnFillOn             = True
     
121.    
     
122. 
     
123.       ; Total Precipitation (color fill)
     
124.         contour_tot = wrf_contour(a,wks, rain_tot, opts_r)
     
125.    
     
126.       ; Precipitation Tendencies
     
127.         opts_r@SubFieldTitle = "from " + times_sav + " to " + times(it)
     
128.    
     
129.         contour_tend = wrf_contour(a,wks, rain_tot_tend,opts_r) ; total (color)
     
130.         contour_res = wrf_contour(a,wks,rain_exp_tend,opts_r)   ; exp (color)
     
131.         opts_r@cnFillOn = False
     
132.         opts_r@cnLineColor = "Red4"
     
133.         contour_prm = wrf_contour(a,wks,rain_con_tend,opts_r)   ; con (red lines)
     
134.         delete(opts_r)
     
135. 
     
136. 
     
137. 
     
138.       ; MAKE PLOTS                                       
     
139. 
     
140.         ; Total Precipitation
     
141.           plot = wrf_map_overlays(a,wks,contour_tot,pltres,mpres)
     
142. 
     
143.         ; Total Precipitation Tendency + SLP
     
144.           plot = wrf_map_overlays(a,wks,(/contour_tend,contour_psl/),pltres,mpres)
     
145. 
     
146.         ; Non-Convective and Convective Precipiation Tendencies
     
147.           plot = wrf_map_overlays(a,wks,(/contour_res,contour_prm/),pltres,mpres)
     
148. 
     
149.     end if    ; END IF FOR SKIPPING FIRST TIME
     
150. 
     
151. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
     
152. 
     
153.     times_sav = times(it)
     
154.     FirstTime = False
     
155.   end do        ; END OF TIME LOOP
     
156. 
     
157. end

复制代码

空城不空

有没有大神哇

weiguang1233

出图你先用x11看一下，是不是被覆盖了

1649518749

检查数据吧