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

MeteoInfo

良辰 发表于 2017-2-13 18:42
王老师，的确是您说的问题，我在服务器上可以正常运行！但是最新版本里这个tostation这个功能改变了吗 ...

DimDataFiles类目前还没有tostation函数，可以考虑用interpn函数，参考此贴：http://bbs.06climate.com/forum.php?mod=viewthread&tid=46753

良辰

MeteoInfo 发表于 2017-2-13 22:26
DimDataFiles类目前还没有tostation函数，可以考虑用interpn函数，参考此贴：http://bbs.06climate.com/f ...

之前的贴子就用到过tostation函数，后来修改删除了吗？tostation函数和interpn函数有区别吗?影响计算结果吗？http://bbs.06climate.com/forum.p ... &extra=page%3D4

MeteoInfo

良辰 发表于 2017-2-14 04:28
之前的贴子就用到过tostation函数，后来修改删除了吗？tostation函数和interpn函数有区别吗?影响计算结果 ...

注意DimDataFile和DimDataFiles是两个不同的类，DimDataFiles是最近加的专门处理一次加载多个数据文件的类，DimDataFile类有tostation函数，但DimDataFiles类里目前没有。interpn函数是处理多维数组插值的，和tostation一样都是双线性插值算法。

良辰

MeteoInfo 发表于 2017-2-14 08:18
注意DimDataFile和DimDataFiles是两个不同的类，DimDataFiles是最近加的专门处理一次加载多个数据文件的 ...

王老师，遇到这个问题，浮点数不能插值。如何解决？
Traceback (most recent call last):
  File "<iostream>", line 143, in <module>
AttributeError: 'float' object has no attribute 'interpn'

MeteoInfo

良辰 发表于 2017-2-14 11:45
王老师，遇到这个问题，浮点数不能插值。如何解决？
Traceback (most recent call last):
  File "", l ...

多维数组才能进行插值

良辰

MeteoInfo 发表于 2017-2-14 12:38
多维数组才能进行插值

        v = vort.interpn([t,pres,lat,lon])
        ww = mpv.interpn([t,press,lat,lon])
        temp = t0.interpn([t,pres,lat,lon])
        rh = rh.interpn([t,pres,lat,lon])
是多维数组插值，我感觉问题应该出在湿度rh的高度维数比温度等多一层，会有影响吗？如有该如何修改那？谢谢王老师！维数设置代码如下：
v = fs['RELH']
rh = v[:,:,latlim,lonlim]
nx = rh.dimlen(3)
ny = rh.dimlen(2)
nz = rh.dimlen(1)
nt = rh.dimlen(0)
lat = rh.dimvalue(2)
lev = rh.dimvalue(1)
t0 = fs['TEMP'][:,:nz-1,latlim,lonlim]
uwnd = fs['UWND'][:,:nz-1,latlim,lonlim]
vwnd = fs['VWND'][:,:nz-1,latlim,lonlim]
vort = np.hcurl(uwnd, vwnd)

MeteoInfo

良辰 发表于 2017-2-14 13:19
v = vort.interpn([t,pres,lat,lon])
        ww = mpv.interpn([t,press,lat,lon])
        t ...

是哪一行代码出错？

良辰

      rh = rh.interpn([t,pres,lat,lon])

良辰

MeteoInfo 发表于 2017-2-14 13:39
是哪一行代码出错？

就是相对湿度这行！而且注释掉就可以运行了
   rh = rh.interpn([t,pres,lat,lon])
另外，王老师，如果我用到的数据涉及到跨月份该怎么办呐？咱这个最多可以加载多少个数据文件那？

MeteoInfo

良辰 发表于 2017-2-14 15:09
就是相对湿度这行！而且注释掉就可以运行了
   rh = rh.interpn([t,pres,lat,lon])
另外，王老师，如果 ...

数据文件可以加载很多，不过读取的数组太大的话可能会有内存不够的问题。