dim_avg_n_Wrap问题求助

我是小猪

我用ncep的 hgt.mon.mean.nc数据提取79-99年六月份的异常，先用dim_rmvmean_n_Wrap求异常，然后在时间维求平均，脚本如下：

1. yrStrt1 = 1979
   
2. yrLast1 = 1999
   
3. ;yrStrt2 = 2000
   
4. ;yrLast2 = 2015
   
5. 
   
6. ff=addfile("/disk3/ncep/monthly/pressure/hgt.mon.mean.nc","r")
   
7. TIME   = ff->time
   
8. YYYY   = cd_calendar(TIME,-1)/100   
   
9. 
   
10. ;1929-1999            
    
11. iYYYY1  = ind(YYYY.ge.yrStrt1 .and. YYYY.le.yrLast1)
    
12. HGT1=short2flt(ff->hgt(iYYYY1,{500},:,:))
    
13. Hgt1=HGT1(5::12,:,:)
    
14. ;ave1=dim_avg_n_Wrap(Hgt1,0)
    
15. ;anomy1=Hgt1
    
16. ;do n=0,dimsizes(Hgt1(:,0,0))-1
    
17. ;    anomy1(n,:,:)=Hgt1(n,:,:)-ave1(:,:)
    
18. ;end do
    
19. Hight = dim_rmvmean_n_Wrap(Hgt1,0)
    
20. printMinMax(Hight, False)
    
21. printVarSummary(Hight)
    
22. aanave1=dim_avg_n_Wrap(Hight,0)
    
23. ;printVarSummary(aanave1)
    
24. printMinMax(aanave1,False)

复制代码

奇怪的是，对异常检查printMinMax范围是 -217到551(m)，挺正常的，但是对异常作dim_avg_n_Wrap(或者dim_sum函数)时候，printMinMax其值竟然是
min=-2.17983e-06   max=2.17983e-06，我觉得是不是加入缺省值进去计算了啊，但是dim_avg函数不是自动忽略缺省值的嘛，这么小的数简直没法画图，不知道该怎么解决？愁啊！！ 希望大大们给点意见啊 !!在此感谢！！

井中月

你先remove了Hgt1第0维的mean，赋给了变量Hight，然后再对Hight第0维求平均，得到的结果必然接近于0啊

我是小猪

井中月 发表于 2017-7-20 08:44
你先remove了Hgt1第0维的mean，赋给了变量Hight，然后再对Hight第0维求平均，得到的结果必然接近于0啊

好像有点懂了，谢谢啦！

我是小猪

井中月 发表于 2017-7-20 08:44
你先remove了Hgt1第0维的mean，赋给了变量Hight，然后再对Hight第0维求平均，得到的结果必然接近于0啊

好像有点懂了，谢谢啦！

禾栗儿

什么意思呢，遇到了同样的问题 大佬求解答