GRIB转ARL数据的Groovy脚本

MeteoInfo

MeteoInfo Java版具有更好的跨平台能力，对于做模式的人来说通常会用Linux或者Unix，很适合使用Java版。之前说过由于Jython对中文的支持不好，选择了Groovy作为MeteoInfo Java版的脚本语言。这里给出一个将GRIB数据转为HYSPLIT模式用的ARL数据的Groovy脚本例子。

MeteoInfo Java版里包含了一个脚本代码编辑器，可以在MeteoInfo桌面程序的“工具 -> 脚本”菜单打开编辑器，在编辑器中可以新建、打开、编辑、保存脚本程序并运行。


也可以在命令行中运行MeteoInfo.jar时加上“-e”参数，从而在不启动MeteoInfo桌面程序的情况下直接打开脚本剪辑器：


还可以直接在命令行中运行脚本程序，不启动任何界面，适合脚本自动运行：


示例程序转成的ARL数据：


脚本代码：

1. //--------------------------------------------------------        
   
2. // Author: Yaqiang Wang                                          
   
3. // Date: 2013-8-27                                               
   
4. // Purpose: Convert GRIB data to ARL data  
   
5. // Note: Sample                                                   
   
6. //-----------------------------------------------------------
   
7. import org.meteoinfo.data.*
   
8. import org.meteoinfo.data.meteodata.*
   
9. import org.meteoinfo.data.meteodata.arl.*
   
10. 
    
11. // Set directory
    
12. dataDir = 'D:/Temp/'
    
13. 
    
14. // Set output data file
    
15. outfn = dataDir + 'arl/test.arl'
    
16. 
    
17. // Read GRIB data file
    
18. mydata = new MeteoDataInfo()
    
19. infn = dataDir + 'grib/fnl_20110416_00_00'
    
20. println infn + '\n'
    
21. mydata.openNetCDFData(infn)
    
22. println 'GRIB file has been opened...\n'
    
23. 
    
24. // Set output ARL data info
    
25. arlDI = new ARLDataInfo()
    
26. 
    
27. // Set variable and level list
    
28. gvars = ['Pressure_surface','Temperature_height_above_ground',
    
29.         'u-component_of_wind_height_above_ground','v-component_of_wind_height_above_ground',
    
30.         'Total_precipitation_surface_0_Hour_Accumulation','Geopotential_height_isobaric','Temperature_isobaric',
    
31.         'u-component_of_wind_isobaric','v-component_of_wind_isobaric','Vertical_velocity_pressure_isobaric',
    
32.         'Relative_humidity_isobaric']
    
33. avars = ['PRSS','T02M','U10M','V10M','TPP6','HGTS','TEMP','UWND','VWND','WWND','RELH']
    
34. levels = [0,1000.0,2000.0,3000.0,5000.0,7000.0,10000.0,15000.0,20000.0,25000.0,30000.0,35000.0,40000.0,45000.0,
    
35.         50000.0,55000.0,60000.0,65000.0,70000.0,75000.0,80000.0,
    
36.         85000.0,90000.0,92500.0,95000.0,97500.0,100000.0]
    
37. 
    
38. levels.each { l->
    
39.         arlDI.levels.add(l)
    
40.         if (l == 0){
    
41.                 arlDI.LevelVarList.add(['PRSS','T02M','U10M','V10M','TPP6'])               
    
42.         } else if (l < 10000) {
    
43.                 arlDI.LevelVarList.add(['HGTS','TEMP','UWND','VWND'])
    
44.         } else {
    
45.                 arlDI.LevelVarList.add(['HGTS','TEMP','UWND','VWND','WWND','RELH'])
    
46.         }
    
47. }
    
48. 
    
49. // Write ARL data file
    
50. dataInfo = mydata.getDataInfo()
    
51. arlDI.createDataFile(outfn)
    
52. arlDI.X = dataInfo.getXDimension().getValues()
    
53. arlDI.Y = dataInfo.getYDimension().getValues()
    
54. variables = dataInfo.getVariables()
    
55. tNum = dataInfo.getTimeNum()
    
56. for (t in 0..tNum-1){
    
57.         mydata.setTimeIndex(t)
    
58.         atime = dataInfo.getTimes().get(t)
    
59.         aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'FNL1', 2)
    
60.         arlDI.writeIndexRecord(atime, aDH)
    
61.         lidx = 0
    
62.         arlDI.levels.each { l->
    
63.                 println l
    
64.                 println '\n'
    
65.                 arlDI.LevelVarList[lidx].each { v->
    
66.                         vName = gvars[avars.indexOf(v)]
    
67.                         println vName + '\n'
    
68.                         if (lidx == 0){
    
69.                                 mydata.setLevelIndex(lidx)
    
70.                         } else {
    
71.                                 variable = dataInfo.getVariable(vName)
    
72.                                 nlidx = variable.getZDimension().getDimValue().indexOf (new Double(l))
    
73.                                 println 'Level index: ' + nlidx.toString() + '\n'
    
74.                                 mydata.setLevelIndex(nlidx)
    
75.                         }
    
76.                         gData = mydata.getGridData(vName)
    
77.                         if (gData == null){
    
78.                                 println 'Grid data is null!\n'
    
79.                         }
    
80.                         if (v == 'PRSS' || v == 'WWND'){
    
81.                                 gData = gData.divide(100)
    
82.                         } else if (v == 'TPP6'){
    
83.                                 gData = gData.divide(1000)
    
84.                         }
    
85.                         aDL = new DataLabel(atime)
    
86.                         aDL.setLevel(lidx)
    
87.                         aDL.setVarName(v)
    
88.                         aDL.setGrid(99)
    
89.                         aDL.setForecast(0)
    
90.                         arlDI.writeGridData(aDL, gData)
    
91.                 }
    
92.                 lidx += 1
    
93.         }
    
94. }
    
95. 
    
96. arlDI.closeDataFile()
    
97. 
    
98. println 'Finished!'
    

//--------------------------------------------------------        <br /> // Author: Yaqiang Wang                                           <br /> // Date: 2013-8-27                                               <br /> // Purpose: Convert GRIB data to ARL data  <br /> // Note: Sample                                                   <br /> //-----------------------------------------------------------<br /> import org.meteoinfo.data.*<br /> import org.meteoinfo.data.meteodata.*<br /> import org.meteoinfo.data.meteodata.arl.*<br /> <br /> // Set directory<br /> dataDir = 'D:/Temp/'<br /> <br /> // Set output data file<br /> outfn = dataDir + 'arl/test.arl'<br /> <br /> // Read GRIB data file<br /> mydata = new MeteoDataInfo()<br /> infn = dataDir + 'grib/fnl_20110416_00_00'<br /> println infn + '\n'<br /> mydata.openNetCDFData(infn)<br /> println 'GRIB file has been opened...\n'<br /> <br /> // Set output ARL data info<br /> arlDI = new ARLDataInfo()<br /> <br /> // Set variable and level list<br /> gvars = ['Pressure_surface','Temperature_height_above_ground',<br />         'u-component_of_wind_height_above_ground','v-component_of_wind_height_above_ground',<br />         'Total_precipitation_surface_0_Hour_Accumulation','Geopotential_height_isobaric','Temperature_isobaric',<br />         'u-component_of_wind_isobaric','v-component_of_wind_isobaric','Vertical_velocity_pressure_isobaric',<br />         'Relative_humidity_isobaric']<br /> avars = ['PRSS','T02M','U10M','V10M','TPP6','HGTS','TEMP','UWND','VWND','WWND','RELH']<br /> levels = [0,1000.0,2000.0,3000.0,5000.0,7000.0,10000.0,15000.0,20000.0,25000.0,30000.0,35000.0,40000.0,45000.0,<br />         50000.0,55000.0,60000.0,65000.0,70000.0,75000.0,80000.0,<br />         85000.0,90000.0,92500.0,95000.0,97500.0,100000.0]<br /> <br /> levels.each { l-><br />         arlDI.levels.add(l)<br />         if (l == 0){<br />                 arlDI.LevelVarList.add(['PRSS','T02M','U10M','V10M','TPP6'])                <br />         } else if (l < 10000) {<br />                 arlDI.LevelVarList.add(['HGTS','TEMP','UWND','VWND'])<br />         } else {<br />                 arlDI.LevelVarList.add(['HGTS','TEMP','UWND','VWND','WWND','RELH'])<br />         }<br /> }<br /> <br /> // Write ARL data file<br /> dataInfo = mydata.getDataInfo()<br /> arlDI.createDataFile(outfn)<br /> arlDI.X = dataInfo.getXDimension().getValues()<br /> arlDI.Y = dataInfo.getYDimension().getValues()<br /> variables = dataInfo.getVariables()<br /> tNum = dataInfo.getTimeNum()<br /> for (t in 0..tNum-1){<br />         mydata.setTimeIndex(t)<br />         atime = dataInfo.getTimes().get(t)<br />         aDH = arlDI.getDataHead(mydata.getProjectionInfo(), 'FNL1', 2)<br />         arlDI.writeIndexRecord(atime, aDH)<br />         lidx = 0<br />         arlDI.levels.each { l-><br />                 println l<br />                 println '\n'<br />                 arlDI.LevelVarList[lidx].each { v-><br />                         vName = gvars[avars.indexOf(v)]<br />                         println vName + '\n'<br />                         if (lidx == 0){<br />                                 mydata.setLevelIndex(lidx)<br />                         } else {<br />                                 variable = dataInfo.getVariable(vName)<br />                                 nlidx = variable.getZDimension().getDimValue().indexOf (new Double(l))<br />                                 println 'Level index: ' + nlidx.toString() + '\n'<br />                                 mydata.setLevelIndex(nlidx)<br />                         }<br />                         gData = mydata.getGridData(vName)<br />                         if (gData == null){<br />                                 println 'Grid data is null!\n'<br />                         }<br />                         if (v == 'PRSS' || v == 'WWND'){<br />                                 gData = gData.divide(100)<br />                         } else if (v == 'TPP6'){<br />                                 gData = gData.divide(1000)<br />                         }<br />                         aDL = new DataLabel(atime)<br />                         aDL.setLevel(lidx)<br />                         aDL.setVarName(v)<br />                         aDL.setGrid(99)<br />                         aDL.setForecast(0)<br />                         arlDI.writeGridData(aDL, gData)<br />                 }<br />                 lidx += 1<br />         }<br /> }<br /> <br /> arlDI.closeDataFile()<br /> <br /> println 'Finished!'<br />

随缘

继续膜拜！

La_respiration

谢谢老师

tianshui

王老师：
      我看到程序中地面数据有5个，每层数据有6个，这是ARL数据所需要的最少变量数吗？ARL的数据格式是什么？

MeteoInfo

tianshui 发表于 2013-9-3 11:26

                               
登录/注册后可看大图

王老师：
      我看到程序中地面数据有5个，每层数据有6个，这是ARL数据所需要的最少变量数吗？ARL的数据 ...

差不多吧。

ARL格式见此网页：https://hysplitbbs.arl.noaa.gov/viewtopic.php?f=25&t=48

tianshui

MeteoInfo 发表于 2013-9-3 11:56

                               
登录/注册后可看大图

差不多吧。

ARL格式见此网页：https://hysplitbbs.arl.noaa.gov/viewtopic.php?f=25&t=48

我看到PDF上第6页开始有个rotating checksum of the packed data，我在读取第一个单元的记录时，后面会有.00000 5PRSS 13 T02M  2 U10M 86 V10M 86 TPP6 84 1000.0 6HGTS189 TEMP 47 UWND201

  VWND198 WWND241 RELH 40 ，变量后面的应该就是这个吧，请问代表什么含义啊？

MeteoInfo

tianshui 发表于 2013-9-4 11:32

                               
登录/注册后可看大图

我看到PDF上第6页开始有个rotating checksum of the packed data，我在读取第一个单元的记录时，后面会有 ...

是这样的。这个应该是程序里数据格式检验生成的数字，应该对于使用没有影响。

tianshui

MeteoInfo 发表于 2013-9-4 14:26

                               
登录/注册后可看大图

是这样的。这个应该是程序里数据格式检验生成的数字，应该对于使用没有影响。

王老师可曾做过将4个字节的数据转换为ARL数据？有fortran文件吗

MeteoInfo

tianshui 发表于 2013-9-4 17:27

                               
登录/注册后可看大图

王老师可曾做过将4个字节的数据转换为ARL数据？有fortran文件吗

4个字节的数据不知道指的是什么？

tianshui

MeteoInfo 发表于 2013-9-4 20:30

                               
登录/注册后可看大图

4个字节的数据不知道指的是什么？

就是普通的数据格式