转arl

长野YY

我参照家园帖子把grib格式的数据转为ARL，可总是报错run script...Output ARL file: F:/2024(1)GRIB/202401.arl
F:/2024(1)GRIB/hebing.grib
Traceback (most recent call last):
  File "<iostream>", line 14, in <module>
NameError: name 'addfile' is not defined
以下是我要转换的脚本，求各位大佬指点
from org.meteoinfo.data.meteodata import MeteoDataInfo
from org.meteoinfo.data.meteodata.arl import ARLDataInfo
from org.meteoinfo.data.meteodata.arl import DataLabel
import os

#---- Set data folder
datadir = r'F:/2024(1)GRIB/'
#---- Set output ARL data file
outfn = os.path.join(datadir, '202401.arl')
print 'Output ARL file:', outfn
#---- Read a GRIB data file
infn = os.path.join(datadir, 'hebing.grib')
print infn
inf = addfile(infn)
print 'GRIB data file has been opened...'
#---- Set output ARL data file
arlf = addfile(outfn, 'c', dtype='arl')
#---- Set variable and level list
gvar2d = ['Surface_pressure_surface','2_metre_temperature_surface','10_metre_U_wind_component_surface',\
    '10_metre_V_wind_component_surface'
gvar3d = ['Geopotential_isobaric','Temperature_isobaric','Vertical_velocity_isobaric',\
    'U-component_of_wind_isobaric','V-component_of_wind_isobaric','Relative_humidity_isobaric'
avar2d = ['PRSS','T02M','U10M','V10M'
avar3d = ['HGTS','TEMP','WWND','UWND','VWND','SPHU'
gv = inf['Geopotential_isobaric'
nx = gv.dimlen(gv.ndim - 1)
ny = gv.dimlen(gv.ndim - 2)
levels = gv.dimvalue(gv.ndim - 3)[::-1
nz = len(levels)
arlf.set2dvar(avar2d)
arlf.setlevels(levels)
for l in levels:
    arlf.set3dvar(avar3d)

#---- Write ARL data file
arlf.setx(gv.dimvalue(gv.ndim - 1))
arlf.sety(gv.dimvalue(gv.ndim - 2))
tNum = inf.timenum()
fhour = 0
for t in range(0, tNum):
    print 'Time index: ' + str(t)
    atime = inf.gettime(t)
    print atime.strftime('%Y-%m-%d %H:00')
    dhead = arlf.getdatahead(inf.proj, 'RSMC', 2, fhour)
    #Pre-write index record without checksum - will be over-write latter
    arlf.writeindexrec(atime, dhead)
    #Checksum list
    ksumlist = [
    # Write 2d variables
    ksums = [
    for avname,gvname in zip(avar2d, gvar2d):
        print avname + ' ' + gvname
        gdata = inf[gvname[t,:,:
        if avname == 'PRSS':
            gdata = gdata * 0.01
        ksum = arlf.writedatarec(atime, 0, avname, fhour, 99, gdata)
        ksums.append(ksum)
    ksumlist.append(ksums)
    # Write 3d variables
    for lidx in range(0, nz):
        ksums = [
        llidx = nz - lidx - 1
        print lidx
        print llidx
        for avname,gvname in zip(avar3d, gvar3d):
            print avname + ' ' + gvname
            gdata = inf[gvname[t,llidx,:,:
            if avname == 'WWND':
                gdata = gdata * 0.01
            elif avname == 'SPHU':
                gdata = gdata * 1000.
            ksum = arlf.writedatarec(atime, lidx + 1, avname, fhour, 99, gdata)
            ksums.append(ksum)
        ksumlist.append(ksums)
    #Re-write index record with checksum
    arlf.writeindexrec(atime, dhead, ksumlist)
    fhour += 1
arlf.close()
print 'Finished!'

MeteoInfo

参考这里：http://www.meteothink.org/exampl ... ctory/grib2arl.html

长野YY

我就是参照这个帖子在meteoinfomap里面的脚本编辑器运行的，可是错误还是run script...
F:/2024Temp/hebing.grib
Traceback (most recent call last):
  File "<iostream>", line 12, in <module>
NameError: name 'addfile' is not defined

MeteoInfo

长野YY 发表于 2025-4-23 15:43
我就是参照这个帖子在meteoinfomap里面的脚本编辑器运行的，可是错误还是run script...
F:/2024Temp/hebin ...

在MeteoInfoLab里运行

长野YY

MeteoInfo 发表于 2025-4-23 22:00
在MeteoInfoLab里运行

老师我就是在meteoinlab里面的脚本编辑器运行的

长野YY

长野YY 发表于 2025-4-24 10:16
老师我就是在meteoinlab里面的脚本编辑器运行的

>>> run script...
F:/2024Temp/hebing.grib
Traceback (most recent call last):
  File "F:\?????2?.py", line 15, in <module>
    inf = addfile(infn)
  File "F:\MeteoInfo_4.0\MeteoInfo\pylib\mipylib\dataset\midata.py", line 122, in addfile
    if datainfo.getDataType() == MeteoDataType.RADAR:
AttributeError: 'NoneType' object has no attribute 'getDataType'

MeteoInfo

长野YY 发表于 2025-4-24 10:22
>>> run script...
F:/2024Temp/hebing.grib
Traceback (most recent call last):

数据格式的问题吧，程序没有识别出来，建议用MeteoInfo最新版本（如果你用的不是最新）试试。还有问题的话可以发上来一个实例数据文件，我这边测试一下。