关于Python提取美国GFS资料降水量的一点经验分享

haozl07

    最近借助Python和wgrib提取了美国GFS资料的总降水量，用于其他计算和作图。琢磨尝试花费了一点时间，终于成功提取，现将关键代码分享出来。
    import os
    convert = "wgrib2\wgrib2 "+filename_gfs+" -match '(:APCP:surface:0-"+str(lday)+")' -text GFSrain1.txt"
    os.system("%s" %(convert))
    关于wgrib工具，我的使用习惯是在*.py文件的所在位置新建wgrib文件夹，然后将相关的库、可执行文件放置其中，代码体现为“wgrib2\wgrib2”。
    filename_gfs为GFS文件名，lday为1、2、3…这样的数字，代码中“:APCPsurface:0-3”即指模式预报产品的3天累积降水量（这时对应的GFS文件应为72小时预报文件）。如果要计算48-72h预报区间的24小时累积降水量，则是用72小时累积降水量减去48小时累积降水量。
    后续的具体读取代码如下：

    gfsr1 = []
    fgfs1 = open('GFSrain1.txt','r')
    fgfs1.readline()
    for line in fgfs1:
        word = line.split()
        gfsr1.append(word[0])
    fgfs1.close()
    gfsrain1 = np.array(gfsr1)
    gfsrain1.resize(361, 720)
    gfsrain1 = gfsrain1.astype(float)
    这一段是代码效率并不高，如果有大神路过，还需指点一二。

puck66

请教楼主， import os  这句是何用意？

一大碗年糕

{:5_213:}{:5_213:}

haozl07

puck66 发表于 2021-05-23 11:45
请教楼主， import os  这句是何用意？

调用电脑系统命令用的

来自安卓APP客户端

飞哥哥

在网上找到的别人的程序也是这样的 不知道有没有什么更高级的办法

werewolf

飞哥哥 发表于 2021-5-25 14:33
在网上找到的别人的程序也是这样的 不知道有没有什么更高级的办法

用wgrib2或grib_copy提取降水后，直接读grib2文件。可以参考这里：
https://mp.weixin.qq.com/s/xfCvgQ7tS7E7hNFnUnqUZQ

如果喜欢可以关注公众号，不定期推送一些原创程序