爱气象，爱气象家园！ 
·设为首页 开启辅助访问 新帖子新回复热门贴精华帖

气象家园

 找回密码
 立即注册

QQ登录

只需一步，快速开始

新浪微博登陆

只需一步, 快速开始

搜索
气象家园»论坛 › 精品推荐 › MeteoInfo › MeteoInfo脚本示例：GRIB to ARL
12345678910下一页
返回列表 发新帖
楼主: MeteoInfo

MeteoInfo脚本示例：GRIB to ARL

[复制链接]
MeteoInfo

新浪微博达人勋
 楼主| 发表于 2019-3-29 15:23:19 | 显示全部楼层
搁浅 发表于 2019-3-29 14:43
王老师，你好，看了下这个程序感觉和python程序比较像，不知道可不可以用python进行调试输出呢，感觉python ...

参考这个：http://bbs.06climate.com/forum.p ... &extra=page%3D1
和这个：http://www.meteothink.org/exampl ... ctory/grib2arl.html
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

sandytiny

新浪微博达人勋
发表于 2024-3-15 19:13:20 | 显示全部楼层
王老师，您好！我使用您的grib转arl的程序，将模式输出的grib数据转为arl格式数据，但转出来的数据在meteoinfo里面打开报错，hysplit里面也用不了。通过多次测试，发现经纬度网格只要小于1000个就不会出错，请问您的代码里面有哪个参数是控制网格数的吗？还是说arl格式就是不允许太大的网格数的呢？
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

MeteoInfo

新浪微博达人勋
 楼主| 发表于 2024-3-15 22:48:35 | 显示全部楼层
sandytiny 发表于 2024-3-15 19:13
王老师，您好！我使用您的grib转arl的程序，将模式输出的grib数据转为arl格式数据，但转出来的数据在meteoi ...

你可以考虑传上来一个示例数据文件和脚本程序，我测试一下。
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

sandytiny

新浪微博达人勋
发表于 2024-3-16 17:24:10 | 显示全部楼层
MeteoInfo 发表于 2024-3-15 22:48
你可以考虑传上来一个示例数据文件和脚本程序，我测试一下。

王老师，我放了2个不同分辨率的fnl的grib2数据到网盘，麻烦王老师帮忙看看，感谢。网格数小于1000的那个grib2在我这边测试是正常转换的，大网格数的那个转出来以后meteoinfo就报错了，HYSPLIT也无法用。我这边使用的代码见附件。谢谢王老师。
数据链接：https://pan.baidu.com/s/1xuweNCixF9U96iNmWPx6ZQ?pwd=tfu8

grib2arl.py

3.01 KB, 下载次数: 0, 下载积分: 金钱 -5

密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

sandytiny

新浪微博达人勋
发表于 2024-3-16 17:26:27 | 显示全部楼层
王老师，我贴上来了

  1. #--------------------------------------------------------        
  2. # Author: Yaqiang Wang                                          
  3. # Date: 2012-12-31                                               
  4. # Purpose: Convert GRIB data to ARL data  
  5. # Note: Sample                                                   
  6. #-----------------------------------------------------------
  7. import clr
  8. from System import *
  9. from System.Collections.Generic import *
  10. clr.AddReference("MeteoInfoC.dll")
  11. from MeteoInfoC import *
  12. from MeteoInfoC.Data import *
  13. from MeteoInfoC.Data.MeteoData import *

  15. #---- Set directories
  16. dataDir = "G:\\dataprocess\\test\\grib2\"

  18. #---- Set output data file
  19. outFile = dataDir + 'fnl_20191101_06_00.arl'

  21. #---- Read a GRIB data file
  22. mydata = MeteoDataInfo()
  23. infile = dataDir + 'fnl_20191101_06_00.grib2'
  24. mydata.OpenGRIBData(infile)

  26. #---- Set output ARL data info
  27. arlDI = ARLDataInfo()

  29. #---- Set variable and level list
  30. gvars = ['U-component_of_wind@pressure','V-component_of_wind@pressure',\
  31.     'Temperature@pressure','Geopotential_height@pressure',\
  32.     'Relative_humidity@pressure','Vertical_velocity@pressure',\
  33.     'Temperature@height_above_ground','Relative_humidity@height_above_ground','U-component_of_wind@height_above_ground',\
  34.     'V-component_of_wind@height_above_ground','Pressure_reduced_to_MSL@msl',\
  35.     'Pressure@surface','Precipitation_rate@surface']
  36. avars = ['UWND','VWND','TEMP','HGTS','RELH','WWND','T02M','RH2M','U10M','V10M','MSLP','PRSS','TPP6']
  37. levels = [0,200,500,700,850,925,1000]
  38. for l in levels:
  39.     arlDI.levels.Add(l)
  40.     if l == 0:
  41.         arlDI.LevelVarList.Add(List[str](['T02M','RH2M','U10M','V10M','MSLP','PRSS','TPP6']))
  42.     else:
  43.         arlDI.LevelVarList.Add(List[str](['UWND','VWND','TEMP','HGTS','RELH','WWND']))

  45. #---- Write ARL data file
  46. arlDI.CreateDataFile(outFile)
  47. arlDI.X = mydata.GetX()
  48. arlDI.Y = mydata.GetY()
  49. variables = mydata.GetVariables()
  50. tNum = mydata.GetTimeNumber()
  51. for t in range(0, tNum):
  52.     mydata.TimeIndex = t
  53.     atime = mydata.GetTime(t)
  54.     print atime.ToString("yyyy-MM-dd HH")
  55.     aDH = arlDI.GetDataHead(mydata.ProjInfo, 'FNL1', 2)
  56.     arlDI.WriteIndexRecord(atime, aDH)
  57.     lidx = 0
  58.     for l in arlDI.levels:
  59.         print l
  60.         for v in arlDI.LevelVarList[lidx]:
  61.             vName = gvars[avars.index(v)]
  62.             print vName
  63.             if lidx == 0:
  64.                 mydata.LevelIndex = lidx
  65.             else:
  66.                 variable = mydata.GetVariable(vName)
  67.                 nlidx = variable.Levels.IndexOf(l)
  68.                 mydata.LevelIndex = nlidx
  69.             gData = mydata.GetGridData(vName)
  70.             if v == 'PRSS' or v == 'WWND' or v == 'MSLP':
  71.                 gData = gData / 100
  72.             elif v == 'TPP6':
  73.                 gData = gData / 1000
  74.             aDL = DataLabel(atime)
  75.             aDL.Level = lidx
  76.             aDL.Variable = v
  77.             aDL.Grid = 99
  78.             aDL.Forecast = 0
  79.             arlDI.WriteGridData(aDL, gData)
  80.         lidx += 1

  82. arlDI.CloseDataFile()

  84. print 'Finished!'
复制代码
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

MeteoInfo

新浪微博达人勋
 楼主| 发表于 2024-3-17 12:56:09 | 显示全部楼层
sandytiny 发表于 2024-3-16 17:26
王老师，我贴上来了

你参考这里：http://www.meteothink.org/exampl ... ctory/grib2arl.html
在MeteoInfoLab里运行代码。
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

sandytiny

新浪微博达人勋
发表于 2024-3-18 18:04:37 | 显示全部楼层
MeteoInfo 发表于 2024-3-17 12:56
你参考这里：http://www.meteothink.org/examples/meteoinfolab/trajectory/grib2arl.html
在MeteoInfoL ...

谢谢王老师，我用您的代码跑了一遍，还是生成了arl文件，但用meteoinfo打开就报错了，王老师能帮忙看看吗？代码如下：

  1. # Convert GRIB data to ARL data

  3. #---- Set data folder
  4. datadir = 'D:\grib\Test'
  5. #---- Set output data file
  6. outfn = os.path.join(datadir, 'gdas1.fnl0p25.2023010100.f00.arl')
  7. #if os.path.exists(outfn):
  8. #    os.remove(outfn)
  9. #---- Read a GRIB data file
  10. infn = os.path.join(datadir, 'gdas1.fnl0p25.2023010100.f00.grib2')
  11. print infn
  12. inf = addfile(infn)
  13. print 'GRIB data file has been opened...'
  14. #---- Set output ARL data file
  15. arlf = addfile(outfn, 'c', dtype='arl')
  16. #---- Set variable and level list
  17. gvar2d = ['Temperature_height_above_ground','Relative_humidity_height_above_ground','u-component_of_wind_height_above_ground',\
  18.     'v-component_of_wind_height_above_ground','Pressure_reduced_to_MSL_msl','Pressure_surface','Precipitation_rate_surface']
  19. gvar3d = ['u-component_of_wind_isobaric','v-component_of_wind_isobaric','Temperature_isobaric',\
  20.     'Geopotential_height_isobaric','Relative_humidity_isobaric','Vertical_velocity_geometric_isobaric']
  21. avar2d = ['T02M','RH2M','U10M','V10M','MSLP','PRSS','TPP6']
  22. avar3d = ['UWND','VWND','TEMP','HGTS','RELH','WWND']
  23. gv = inf['Geopotential_height_isobaric']
  24. nx = gv.dimlen(gv.ndim - 1)
  25. ny = gv.dimlen(gv.ndim - 2)
  26. levels = gv.dimvalue(gv.ndim - 3)[::-1]
  27. nz = len(levels)
  28. arlf.setlevels(levels)
  29. arlf.set2dvar(avar2d)
  30. for l in levels:
  31.     arlf.set3dvar(avar3d)

  33. #---- Write ARL data file
  34. arlf.setx(gv.dimvalue(gv.ndim - 1))
  35. arlf.sety(gv.dimvalue(gv.ndim - 2))
  36. tNum = inf.timenum()
  37. fhour = 0
  38. for t in range(0, tNum):
  39.     print 'Time index: ' + str(t)
  40.     atime = inf.gettime(t)
  41.     print atime.strftime('%Y-%m-%d %H:00')
  42.     dhead = arlf.getdatahead(inf.proj, 'RSMC', 2, fhour)
  43.     #Pre-write index record without checksum - will be over-write latter
  44.     arlf.writeindexrec(atime, dhead)
  45.     #Checksum list
  46.     ksumlist = []
  47.     # Write 2d variables
  48.     ksums = []
  49.     for avname,gvname in zip(avar2d, gvar2d):
  50.         print avname + ' ' + gvname
  51.         if avname == 'U10M' or avname == 'V10M':
  52.             gdata = inf[gvname][t,0,:,:]
  53.         else:
  54.             gdata = inf[gvname][t,:,:]
  55.         if avname == 'PRSS' or avname == 'MSLP':
  56.             gdata = gdata * 0.01
  57.         ksum = arlf.writedatarec(atime, 0, avname, fhour, 99, gdata)
  58.         ksums.append(ksum)
  59.     ksumlist.append(ksums)
  60.     # Write 3d variables
  61.     for lidx in range(0, nz):
  62.         ksums = []
  63.         llidx = nz - lidx - 1
  64.         print lidx
  65.         print llidx
  66.         for avname,gvname in zip(avar3d, gvar3d):
  67.             print avname + ' ' + gvname
  68.             gdata = inf[gvname][t,llidx,:,:]
  69.             if avname == 'WWND':
  70.                 gdata = gdata * 0.01
  71.             elif avname == 'SPHU':
  72.                 gdata = gdata * 1000.
  73.             ksum = arlf.writedatarec(atime, lidx + 1, avname, fhour, 99, gdata)
  74.             ksums.append(ksum)
  75.         ksumlist.append(ksums)
  76.     #Re-write index record with checksum
  77.     arlf.writeindexrec(atime, dhead, ksumlist)
  78.     fhour += 1
  79. arlf.close()
  80. print 'Finished!'
复制代码
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

MeteoInfo

新浪微博达人勋
 楼主| 发表于 2024-3-19 10:46:03 | 显示全部楼层
sandytiny 发表于 2024-3-18 18:04
谢谢王老师，我用您的代码跑了一遍，还是生成了arl文件，但用meteoinfo打开就报错了，王老师能帮忙看看吗 ...

参考这个脚本程序，不过数据文件中只有一个时次。

  1. # Convert GRIB data to ARL data

  3. #---- Set data folder
  4. datadir = 'D:/Temp/grib'
  5. #---- Set output data file
  6. outfn = os.path.join(datadir, 'gdas1.fnl0p25.2023010100.f00.arl')
  7. #if os.path.exists(outfn):
  8. #    os.remove(outfn)
  9. #---- Read a GRIB data file
  10. infn = os.path.join(datadir, 'gdas1.fnl0p25.2023010100.f00.grib2')
  11. print(infn)
  12. inf = addfile(infn)
  13. print('GRIB data file has been opened...')
  14. #---- Set output ARL data file
  15. arlf = addfile(outfn, 'c', dtype='arl')
  16. #---- Set variable and level list
  17. gvar2d = ['Temperature_height_above_ground','Relative_humidity_height_above_ground','u-component_of_wind_height_above_ground',
  18.     'v-component_of_wind_height_above_ground','Pressure_reduced_to_MSL_msl','Pressure_surface','Precipitation_rate_surface']
  19. gvar3d = ['u-component_of_wind_isobaric','v-component_of_wind_isobaric','Temperature_isobaric',
  20.     'Geopotential_height_isobaric','Relative_humidity_isobaric','Vertical_velocity_geometric_isobaric']
  21. avar2d = ['T02M','RH2M','U10M','V10M','MSLP','PRSS','TPP6']
  22. avar3d = ['UWND','VWND','TEMP','HGTS','RELH','WWND']
  23. gv = inf['Geopotential_height_isobaric']
  24. nx = gv.dimlen(gv.ndim - 1)
  25. ny = gv.dimlen(gv.ndim - 2)
  26. levels = gv.dimvalue(gv.ndim - 3)[::-1]
  27. nz = len(levels)
  28. arlf.setlevels(levels)
  29. arlf.set2dvar(avar2d)
  30. for l in levels:
  31.     arlf.set3dvar(avar3d)

  33. #---- Write ARL data file
  34. arlf.setx(gv.dimvalue(gv.ndim - 1))
  35. arlf.sety(gv.dimvalue(gv.ndim - 2))
  36. tNum = inf.timenum()
  37. fhour = 0
  38. for t in range(0, tNum):
  39.     print('Time index: ' + str(t))
  40.     atime = inf.gettime(t)
  41.     print(atime.strftime('%Y-%m-%d %H:00'))
  42.     dhead = arlf.getdatahead(inf.proj, 'RSMC', 2, fhour)
  43.     #Pre-write index record without checksum - will be over-write latter
  44.     arlf.writeindexrec(atime, dhead)
  45.     #Checksum list
  46.     ksumlist = []
  47.     # Write 2d variables
  48.     ksums = []
  49.     for avname,gvname in zip(avar2d, gvar2d):
  50.         print(avname + ' ' + gvname)
  51.         if avname in ['T02M','RH2M','U10M','V10M']:
  52.             gdata = inf[gvname][t,0,:,:]
  53.         else:
  54.             gdata = inf[gvname][t,:,:]
  55.             
  56.         if avname == 'PRSS' or avname == 'MSLP':
  57.             gdata = gdata * 0.01
  58.             
  59.         ksum = arlf.writedatarec(atime, 0, avname, fhour, 99, gdata)
  60.         ksums.append(ksum)
  61.     ksumlist.append(ksums)
  62.     # Write 3d variables
  63.     for lidx in range(0, nz):
  64.         ksums = []
  65.         llidx = nz - lidx - 1
  66.         print(lidx, llidx)
  67.         for avname,gvname in zip(avar3d, gvar3d):
  68.             print avname + ' ' + gvname
  69.             gdata = inf[gvname][t,llidx,:,:]
  70.             if avname == 'WWND':
  71.                 gdata = gdata * 0.01
  72.             elif avname == 'SPHU':
  73.                 gdata = gdata * 1000.
  74.             ksum = arlf.writedatarec(atime, lidx + 1, avname, fhour, 99, gdata)
  75.             ksums.append(ksum)
  76.         ksumlist.append(ksums)
  77.     #Re-write index record with checksum
  78.     arlf.writeindexrec(atime, dhead, ksumlist)
  79.     fhour += 1
  80. arlf.close()
  81. print('Finished!')
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

sandytiny

新浪微博达人勋
发表于 2024-3-19 11:06:58 | 显示全部楼层
MeteoInfo 发表于 2024-3-19 10:46
参考这个脚本程序，不过数据文件中只有一个时次。

谢谢王老师，我这跑了一遍，还是生成arl但就是meteoinfo打开报错，不知道是哪里出了问题
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

MeteoInfo

新浪微博达人勋
 楼主| 发表于 2024-3-19 11:23:00 | 显示全部楼层
sandytiny 发表于 2024-3-19 11:06
谢谢王老师，我这跑了一遍，还是生成arl但就是meteoinfo打开报错，不知道是哪里出了问题

我这里没什么问题，你看看你用的MeteoInfo是不是最新版本。也可以加入MeteoInfo QQ群方便交流（808554605）
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

下一页 »
12345678910下一页
返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

小黑屋|帮助中心|清除痕迹

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表，不代表气象家园立场，禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表