GRIB1批量转换为ARL格式

zhangyw

HYSPLIT单机版自带的转换工具是否可以实现批量转换数据格式呢？看说明书应该是可以的啊，但是一次只能选中一个文件，求解答~~~

zhangyw

看了下运行日志，好像设置时间范围，输入文件选起始文件，会自动读取后面的文件进行转换的~~

龍_彳亍_天_下

请问楼主现在问题解决了嘛

zhangyw

龍_彳亍_天_下 发表于 2021-2-27 17:38
请问楼主现在问题解决了嘛

用Linux下的脚本进行转换，win似乎不行

龍_彳亍_天_下

zhangyw 发表于 2021-4-6 14:45
用Linux下的脚本进行转换，win似乎不行

我现在已经解决啦，用的王老师的程序，编写的python脚本，可以实现windows下的批量转格式，脚本我找一下

龍_彳亍_天_下

import os
import datetime
import time

start = '19810102'
end = '20201231'


datestart = datetime.datetime.strptime(start, '%Y%m%d')
dateend = datetime.datetime.strptime(end, '%Y%m%d')
a = os.system(r"cd D:\coldwave_data\ERA\arl")
print(a)
data_list = list()
while datestart <= dateend:
    time_format = datestart.strftime('%Y%m%d')
    print(time_format)
    grib2arl_str_1 = 'C:\\hysplit\\exec\\grib2arl.exe -iD:\\coldwave_data\\ERA\height2p5\\'
    grib2arl_str_2 = '.grib -p0 -k37 -sD:\\coldwave_data\\ERA\\surface2p5\\'
    grib2arl_str_3 = '.grib -cD:\\coldwave_data\\ERA\\GEO2P5.grib'

    grib2arl_str = str(grib2arl_str_1) + str(time_format) + str(grib2arl_str_2) + str(time_format) + str(grib2arl_str_3)
   
    b = os.popen(grib2arl_str)
    print(b)
    time.sleep(15)

    if os.path.exists('DATA.ARL') == True:
       c = os.system(r"rename DATA.ARL "+time_format+".arl")
       print(c)
       time.sleep(3)
    d = os.system(r"del DATA.ARL")
    time.sleep(3)     


    datestart += datetime.timedelta(days=1)
说明：
下载的数据包括高空层数据、地面数据和GEO数据

zhangyw

龍_彳亍_天_下 发表于 2021-4-6 15:45
import os
import datetime
import time

好哒，谢谢！

melodyz

龍_彳亍_天_下 发表于 2021-4-6 15:45
import os
import datetime
import time

您好，请问GEO数据是在哪里下载的呀？我在ERA5的数据中只看到了pressure level和single level的，但是不知道去哪找invariant的数据

天雨

melodyz 发表于 2023-07-01 15:58
您好，请问GEO数据是在哪里下载的呀？我在ERA5的数据中只看到了pressure level和single level的，但是不知道去哪找invariant的数据

你是否知道了ERA5这个数据？我还是没找到

                               
登录/注册后可看大图

来自苹果APP客户端

zhang714

melodyz 发表于 2023-7-1 15:58
您好，请问GEO数据是在哪里下载的呀？我在ERA5的数据中只看到了pressure level和single level的，但是不 ...

请问您找到GEO数据了吗