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

MeteoInfo

良辰 发表于 2017-2-13 10:58
王老师，程序还是不通呀，能帮忙测试一下吗？

我测试过了，没有问题呀

良辰

MeteoInfo 发表于 2017-2-13 12:49
我测试过了，没有问题呀

我用的1.4r3最新版本呀。脚本我发您私信了！为啥我的运行不了哪？老师您有时间看一下私信呗！

MeteoInfo

良辰 发表于 2017-2-13 13:42
我用的1.4r3最新版本呀。脚本我发您私信了！为啥我的运行不了哪？老师您有时间看一下私信呗！

我测试用的脚本供参考：

1. datadir = r'U:\data\ARL\2016'
   
2. fns = []
   
3. for i in range(1, 5):
   
4.     fn = os.path.join(datadir, 'gdas1.jun16.w' + str(i))
   
5.     print fn
   
6.     fns.append(fn)
   
7. 
   
8. fs = addfiles(fns)
   
9. 
   
10. latlim = [10,60]
    
11. lonlim = [60,140]
    
12. v = fs['RELH']
    
13. rh = v[:,:,latlim,lonlim]
    
14. nx = rh.dimlen(3)
    
15. ny = rh.dimlen(2)
    
16. nz = rh.dimlen(1)
    
17. nt = rh.dimlen(0)
    
18. uwnd = fs['UWND'][:,:nz-1,latlim,lonlim]
    
19. vwnd = fs['VWND'][:,:nz-1,latlim,lonlim]
    
20. vort = np.hcurl(uwnd, vwnd)
    
21. prs = zeros([nt,nz,ny,nx])
    
22. prs = dim_array(prs, rh.dims)

datadir = r'U:\data\ARL\2016'<br /> fns = []<br /> for i in range(1, 5):<br />     fn = os.path.join(datadir, 'gdas1.jun16.w' + str(i))<br />     print fn<br />     fns.append(fn)<br /> <br /> fs = addfiles(fns)<br /> <br /> latlim = [10,60]<br /> lonlim = [60,140]<br /> v = fs['RELH']<br /> rh = v[:,:,latlim,lonlim]<br /> nx = rh.dimlen(3)<br /> ny = rh.dimlen(2)<br /> nz = rh.dimlen(1)<br /> nt = rh.dimlen(0)<br /> uwnd = fs['UWND'][:,:nz-1,latlim,lonlim]<br /> vwnd = fs['VWND'][:,:nz-1,latlim,lonlim]<br /> vort = np.hcurl(uwnd, vwnd)<br /> prs = zeros([nt,nz,ny,nx])<br /> prs = dim_array(prs, rh.dims)

良辰

MeteoInfo 发表于 2017-2-13 13:51
我测试用的脚本供参考：

好的，我再试一试，谢谢王老师！

良辰

MeteoInfo 发表于 2017-2-13 13:51
我测试用的脚本供参考：

王老师，用你的脚本我再次测试了一下：

1. # Set working directory
   
2. trajDir = 'D:/Temp/HYSPLIT'
   
3. datadir = r'D:\Temp\arl'
   
4. # Open meteorological data file
   
5. print 'Open meteorological data file...'
   
6. fns = []
   
7. for i in range(1, 5):
   
8.     fn = os.path.join(datadir, 'gdas1.jul09.w' + str(i))
   
9.     print fn
   
10.     fns.append(fn)
    
11. fs = addfiles(fns)
    
12. 
    
13. #print 'Meteorological file: ' + meteof
    
14. # Read data
    
15. print 'Read data...'
    
16. latlim = [10,60]
    
17. lonlim = [60,140]
    
18. v = fs['RELH']
    
19. rh = v[:,:,latlim,lonlim]
    
20. nx = rh.dimlen(3)
    
21. ny = rh.dimlen(2)
    
22. nz = rh.dimlen(1)
    
23. nt = rh.dimlen(0)
    
24. lat = rh.dimvalue(2)
    
25. lev = rh.dimvalue(1)
    
26. t0 = fs['TEMP'][:,:nz-1,latlim,lonlim]
    
27. uwnd = fs['UWND'][:,:nz-1,latlim,lonlim]
    
28. vwnd = fs['VWND'][:,:nz-1,latlim,lonlim]
    
29. vort = np.hcurl(uwnd, vwnd)
    
30. #divg = np.hdivg(uwnd, vwnd)
    
31. prs = zeros([nt,nz,ny,nx])
    
32. prs = dim_array(prs, rh.dims)

复制代码

错误提示如下：
>>> run script...
Open meteorological data file...
D:\Temp\arl\gdas1.jul09.w1
D:\Temp\arl\gdas1.jul09.w2
D:\Temp\arl\gdas1.jul09.w3
D:\Temp\arl\gdas1.jul09.w4
Read data...
Traceback (most recent call last):
  File "<iostream>", line 29, in <module>
  File "D:\MeteoInfo\pylib\mipylib\numeric\minum.py", line 1564, in hcurl
    r = ArrayMath.hcurl(u.asarray(), v.asarray(), xdim.getDimValue(), ydim.getDimValue())
        at ucar.ma2.ArrayDouble.<init>(ArrayDouble.java:103)
        at ucar.ma2.ArrayDouble$D4.<init>(ArrayDouble.java:341)
        at ucar.ma2.ArrayDouble$D4.<init>(ArrayDouble.java:333)
        at ucar.ma2.ArrayDouble.factory(ArrayDouble.java:64)
        at ucar.ma2.ArrayDouble.factory(ArrayDouble.java:49)
        at ucar.ma2.Array.factory(Array.java:110)
        at ucar.ma2.Array.factory(Array.java:92)
        at org.meteoinfo.data.ArrayMath.cdiff(ArrayMath.java:2941)
        at org.meteoinfo.data.ArrayMath.hcurl(ArrayMath.java:3055)
        at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
        at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
        at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
        at java.lang.reflect.Method.invoke(Unknown Source)

java.lang.OutOfMemoryError: java.lang.OutOfMemoryError: Java heap space
是不是因为版本问题，而且只读入了四个数据文件，我这的确有5个文件。

MeteoInfo

良辰 发表于 2017-2-13 16:08
王老师，用你的脚本我再次测试了一下：

要包括w5需要将第7行range(1, 5)改为range(1, 6)

错误信息提示内存溢出，应该是你的机器内存不够大。一次加载多个文件需要较大内存。

良辰

MeteoInfo 发表于 2017-2-13 16:19
要包括w5需要将第7行range(1, 5)改为range(1, 6)

错误信息提示内存溢出，应该是你的机器内存不够大。 ...

哦，原来是这个原因，王老师，您用的多大内存跑的呀？还是用服务器跑？

MeteoInfo

良辰 发表于 2017-2-13 16:23
哦，原来是这个原因，王老师，您用的多大内存跑的呀？还是用服务器跑？

我的机器内存8G，MeteoInfoLab.exe最多用到2G（这个例子应该够了），Java程序可以在启动的时候设置最大使用的内存。

良辰

MeteoInfo 发表于 2017-2-13 16:41
我的机器内存8G，MeteoInfoLab.exe最多用到2G（这个例子应该够了），Java程序可以在启动的时候设置最大使 ...

我的电脑是4g内存，不应该呀！我一会儿换个16g内存的服务器试一下，要是还有问题还要骚扰您!呜呜

良辰

MeteoInfo 发表于 2017-2-13 16:41
我的机器内存8G，MeteoInfoLab.exe最多用到2G（这个例子应该够了），Java程序可以在启动的时候设置最大使 ...

王老师，的确是您说的问题，我在服务器上可以正常运行！但是最新版本里这个tostation这个功能改变了吗？代码为：
142       temp = fs.tostation(tvar, lon, lat, pres, t)
143       rh = fs.tostation(rhvar, lon, lat, pres, t)
提示如下错误：
Traceback (most recent call last):
  File "<iostream>", line 142, in <module>
AttributeError: 'DimDataFiles' object has no attribute 'tostation'