日降水资料的处理出错，是资料的问题还是处理的问题

咩咩的羊

本人ncl小白，目前刚接触，老师让做一个滤波我就从NCEP Operational Analysis上面下载了日降水的资料 逐年下载的，但是我处理之后怎么都出不来图，grads阅读之后的代码显示

1. ga-> sdfopen pr_wtr.1980.nc
   
2. Scanning self-describing file:  pr_wtr.1980.nc
   
3. SDF file pr_wtr.1980.nc is open as file 1
   
4. LON set to 0 360
   
5. LAT set to -90 90
   
6. LEV set to 0 0
   
7. Time values set: 1980:1:1:0 1980:1:1:0
   
8. E set to 1 1
   
9. ga-> q ctlinfo
   
10. dset pr_wtr.1980.nc
    
11. title Once daily NMC precipitable water data (1980)
    
12. undef -9.99e+33
    
13. dtype netcdf
    
14. xdef 144 linear 0 2.5
    
15. ydef 73 linear -90 2.5
    
16. zdef 1 linear 0 1
    
17. tdef 366 linear 00Z01JAN1980 1440mn
    
18. vars 2
    
19. pr_wtr=>pr_wtr  0  t,y,x  Precitable Water
    
20. head=>head  0  t  Missing
    
21. endvars
    

复制代码

我自己的滤波作图程序：

1. fils=systemfunc("ls pr_wtr.*.nc")
   
2. f = addfiles(fils, "r")
   
3. time=f[:]->head
   
4. x = f[:]->pr_wtr(:,:,:) ;
   
5. ihp = 2 ; band pass
   
6. sigma = 1.0 ; Lanczos sigma
   
7. 
   
8. nWgt = 201 ; loose 100 each end ()
   
9. fca = 1./30. ; start freq
   
10. fcb = 1./10. ; last freq
    
11. wgt = filwgts_lanczos (nWgt, ihp, fca, fcb, sigma )
    
12. printVarSummary(wgt)
    
13. xBPF = wgt_runave_leftdim( x, wgt, 0 ) ; 10~30 day
    
14. copy_VarMeta(x, xBPF)
    
15. printVarSummary(xBPF)
    
16. 
    
17. date = cd_calendar(time&time,-2)
    
18. yrfrac  = yyyymmdd_to_yyyyfrac (date, 0)
    
19. delete(yrfrac@long_name)
    
20. delete(time@long_name)
    
21. pStrt = 19800101 ; 20 years:
    
22. pLast = 19991231
    
23. iStrt = ind(date.eq.pStrt) ; user specified dates
    
24. iLast = ind(date.eq.pLast)
    
25. delete(date)
    
26. 
    
27. wks = gsn_open_wks ("png","filters_precipitation")
    
28. gsn_define_colormap(wks, "GMT_gray")
    
29. res                      = True
    
30. res@gsnAddCyclic         = True   
    
31. res@cnFillOn             = True
    
32. res@cnLineLabelsOn       = False   
    
33. res@lbLabelBarOn         = True
    
34. res@lbLabelBarOn         = False             ; turn off individual cb's
    
35. 
    
36. res@mpMinLatF            = -90.            ; range to zoom in on
    
37. res@mpMaxLatF            = 90.
    
38. res@mpMinLonF            = 0.
    
39. res@mpMaxLonF            = 360.
    
40. res@mpCenterLonF         = 180.
    
41. 
    
42. plot = gsn_csm_contour_map(wks,xBPF(120,:,:),res)
    
43. 
    
44. 
    
45. end

复制代码

最上面的load没有加，我觉得没有必要，另外，ncl运行不报错，就是出不了图，我现在比较困惑的是数据的问题还是我自己代码的问题，求大神解答

咩咩的羊

1. Variable: x
   
2. Type: short
   
3. Total Size: 1690084800 bytes
   
4.             845042400 values
   
5. Number of Dimensions: 3
   
6. Dimensions and sizes:        [time | 7305] x [latitude | 241] x [longitude | 480]
   
7. Coordinates:
   
8.             time: [701256..876552]
   
9.             latitude: [90..-90]
   
10.             longitude: [ 0..359.25]
    
11. Number Of Attributes: 6
    
12.   scale_factor :        0.002021210590883124
    
13.   add_offset :        256.4960291896264
    
14.   _FillValue :        -32767
    
15.   missing_value :        -32767
    
16.   units :        K
    
17.   long_name :        2 metre temperature
    
18. 
    
19. Variable: wgt
    
20. Type: float
    
21. Total Size: 804 bytes
    
22.             201 values
    
23. Number of Dimensions: 1
    
24. Dimensions and sizes:        [201]
    
25. Coordinates:
    
26. Number Of Attributes: 2
    
27.   resp :        <ARRAY of 405 elements>
    
28.   freq :        <ARRAY of 405 elements>
    
29. 
    
30. Variable: xBPF
    
31. Type: float
    
32. Total Size: 3380169600 bytes
    
33.             845042400 values
    
34. Number of Dimensions: 3
    
35. Dimensions and sizes:        [time | 7305] x [latitude | 241] x [longitude | 480]
    
36. Coordinates:
    
37.             time: [701256..876552]
    
38.             latitude: [90..-90]
    
39.             longitude: [ 0..359.25]
    
40. Number Of Attributes: 9
    
41.   missing_value :        -32767
    
42.   missing_value_original :        -32767
    
43.   _FillValue :        -32767
    
44.   _FillValue_original :        -32767
    
45.   long_name :        2 metre temperature
    
46.   units :        K
    
47.   add_offset :        256.4960291896264
    
48.   scale_factor :        0.002021210590883124
    
49.   wgt_runave_op_ncl :        wgt_runave_n
    

复制代码

加一下我运行后输出的数据