MeteoInfoLab脚本示例：同时加载多个文件

guo066600

还想请问一下王老师：MI软件的一个程序可以一次读入多少个FNL的文件？

MeteoInfo

guo066600 发表于 2020-6-13 11:14
还想请问一下王老师：MI软件的一个程序可以一次读入多少个FNL的文件？

不知道，你可以帮着测试一下

Jungle

王老师，你的示例脚本中第15行：ndata = mean(data, axis=0)  中mean函数是对前面打开的所有时次求平均吗？我打开30个文件（1981-2010），每个文件里面是每年5月的一天4次的数据，我按照你的示例代码运行，最后发现ndata是一个二维（x,y）的数据，我想分别对5月的124个时次求30年的平均，最后的数据是三维的，应该怎么做呢？

upuil

1. >>> run script...
   
2. Traceback (most recent call last):
   
3.   File "D:\read_aws_nc.py", line 10, in <module>
   
4.     fs = addfiles(fns)
   
5.   File "D:\MeteoInfo_2.4.5\MeteoInfo\pylib\mipylib\dataset\midata.py", line 78, in addfiles
   
6.     return DimDataFiles(dfs)
   
7.   File "D:\MeteoInfo_2.4.5\MeteoInfo\pylib\mipylib\dataset\dimdatafile.py", line 719, in __init__
   
8.     ftimes.append(ds.gettime(0))
   
9.   File "D:\MeteoInfo_2.4.5\MeteoInfo\pylib\mipylib\dataset\dimdatafile.py", line 265, in gettime
   
10.     t = self.dataset.getDataInfo().getTimes().get(idx)     
    
11.         at org.meteoinfo.data.meteodata.DataInfo.getTimes(DataInfo.java:156)
    
12. 
    
13.         at jdk.internal.reflect.GeneratedMethodAccessor78.invoke(Unknown Source)
    
14. 
    
15.         at java.base/jdk.internal.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43)
    
16. 
    
17.         at java.base/java.lang.reflect.Method.invoke(Method.java:564)
    
18. 
    
19. java.lang.NullPointerException: java.lang.NullPointerException: Cannot invoke "org.meteoinfo.ndarray.Dimension.getDimValue()" because "this.tDim" is null
    
20. >>>

复制代码

王老师 ，我合并多个文件的时候出现这样的错误，是由于没有时间维度吗？

代码如下

1. 
   
2. file_dir = 'E:/aws-nc-19-20/test/'
   
3. fns = []
   
4. for root, dirs, files in os.walk(file_dir):
   
5.     for file in files:  
   
6.         filedir = os.path.join(root,file)
   
7.         fns.append(filedir)
   
8. 
   
9. 
   
10. fs = addfiles(fns)

复制代码

MeteoInfo

upuil 发表于 2021-2-27 16:26
王老师 ，我合并多个文件的时候出现这样的错误，是由于没有时间维度吗？

代码如下

对，需要有时间维

guo066600

王老师，循环读入GDAS的某月w1和w2文件后，提取pblh数据，pblh = fs['PBLH'][:,:,:]
tn = pblh.dimlen(0)
后执行循环写入csv，
for i in range(tn):
    t = f.gettime(i)
    #Interpolate to station
    pblh_st = pblh[i,:,:].interpn([lat, lon])
    ws_st = ws[i,:,:].interpn([lat, lon])
    line = '%s,%.2f,%.2f' % (t.strftime('%Y-%m-%d_%H'), pblh_st, ws_st)
    print line
    of.write(line + '\n')
of.close()

报错显示只写入了w1的数据后就无法写入w2的数据了，错误提示与第一页的某位坛友的报错信息类似，应该如何解决？
抱歉，报错信息一复制进来就无法提交评论，因此只能描述了。