请教如何把计算结果保存或转换成HYSPLIT模型需要的格式

MeteoInfo

aankai 发表于 2017-1-16 07:55
王老师，现在我可以把arl数据直接读出来算平均，怎样再存成ARL格式呢？之前是3-hr的数据，我现在算出了每 ...

是不是你用的MeteoInfo版本太低的缘故？

aankai

MeteoInfo 发表于 2017-1-16 09:24
是不是你用的MeteoInfo版本太低的缘故？

是最新的版本，那个FNL1是什么？会不会是这里的问题

aankai

MeteoInfo 发表于 2017-1-16 09:24
是不是你用的MeteoInfo版本太低的缘故？

王老师，我已经完成了平均值计算，也把结果保存成了ARL格式，但是最后data grid的信息（Pole pnt lat/lon，Tan lat/cone，Syn pnt x/y，Syn pnt lat/lon）变了，请问这个如何设置呢？（meteof[5]是2010年的原数据）

outfn = os.path.join(meteoDir, 'test1.arl')
arlf = addfile(outfn, 'c', dtype='arl')
inf = meteof[5]
avar2d = ['PRSS']
avar3d = ['UWND','VWND','HGTS','TEMP','WWND','RELH','TKEN']
avg3d = ['dyUWND','dyVWND','dyHGTS','dyTEMP','dyWWND','dyRELH','dyTKEN']
wv = inf['UWND']
nx = 185
ny = 129
levels = [1000,975,950,925,900,875,850,825,800,775,750,725,700,\
              650,600,550,500,450,400,350,300,250,200,150,100,50]
nz = len(levels)
arlf.setlevels(levels)
arlf.set2dvar(avar2d)
arlf.setx(wv.dimvalue(wv.ndim - 1))
arlf.sety(wv.dimvalue(wv.ndim - 2))
for l in levels:
    arlf.set3dvar(avar3d)


tNum = 15
fhour = 0
for t in range(0, 15):
    print 'Time index: ' + str(t)
    atime = inf.gettime(t*8)
    print atime.strftime('%Y-%m-%d')
    dhead = arlf.getdatahead(inf.proj, 'EDAS', 2, fhour)
    arlf.writeindexrec(atime, dhead)
    ksumlist = []
    ksums = []
    for avname in avar2d:
        print avname
        gdata = PRSS[t*8,:,:]
        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 = []
        print lidx
        for avname,i in zip(avar3d, range(0,7)):
            print avname,i
            gdata = meteovar[t,lidx,:,:]
            ksum = arlf.writedatarec(atime, lidx + 1, avname, fhour, 99, gdata)
            ksums.append(ksum)
        ksumlist.append(ksums)
    arlf.writeindexrec(atime, dhead, ksumlist)
    fhour += 24
arlf.close()
print 'Finished!

MeteoInfo

aankai 发表于 2017-1-24 08:03
王老师，我已经完成了平均值计算，也把结果保存成了ARL格式，但是最后data grid的信息（Pole pnt lat/lon ...

还要休假到正月15，休假期间不一定能有时间针对你的问题进行调试。

aankai

MeteoInfo 发表于 2017-2-4 18:12
还要休假到正月15，休假期间不一定能有时间针对你的问题进行调试。

好的，那就等休假过后再麻烦您解决了，谢谢

MeteoInfo

aankai 发表于 2017-2-7 00:10
好的，那就等休假过后再麻烦您解决了，谢谢

你用的是GDAS数据吗？应该是经纬度投影吧。

aankai

是EDAS40（Eta Data Assimilation System）的数据，北美地区，兰伯特投影。我只是算了平均值，所以所有的投影信息都应该跟原文件一样，就直接用inf.proj，也尝试过用projinfo重新写，但lat_0和lat_1总是变成一样的值

MeteoInfo

aankai 发表于 2017-2-11 06:08
是EDAS40（Eta Data Assimilation System）的数据，北美地区，兰伯特投影。我只是算了平均值，所以所有的投 ...

你提供几个示例文件我测试一下。

aankai

MeteoInfo 发表于 2017-2-11 10:00
你提供几个示例文件我测试一下。

数据我都是从ftp://arlftp.arlhq.noaa.gov/pub/archives/edas40下载的，以guest身份就可以链接，麻烦您了

MeteoInfo

aankai 发表于 2017-2-11 15:25
数据我都是从ftp://arlftp.arlhq.noaa.gov/pub/archives/edas40下载的，以guest身份就可以链接，麻烦您了

下载速度相当慢