转换GRIB格式数据至MM5中间文件

MeteoInfo

运行MM5模式需要更大范围边界气象场的驱动，T639是目前中国气象局一个比较好的业务数据源，T639是GRIB2格式，需要转换为MM5的中间文件格式。MeteoInfo Java版支持MM5输出和中间文件格式，可以通过写脚本来方便的进行批量转换，并可以通过定时提交作业自动运行，满足业务预报的需求。

示例脚本如下：

1. #--------------------------------------------------------        
   
2. # Author: Yaqiang Wang                                          
   
3. # Date: 2014-5-8                                               
   
4. # Purpose: Convert GRIB data to MM5 intermediate data  
   
5. # Note: Sample                                                   
   
6. #-----------------------------------------------------------
   
7. from org.meteoinfo.data.meteodata import MeteoDataInfo
   
8. from org.meteoinfo.data.meteodata.mm5 import MM5IMDataInfo
   
9. from org.meteoinfo.data.meteodata import Variable
   
10. from org.meteoinfo.data.meteodata.mm5 import DataHead
    
11. import os
    
12. import datetime
    
13. 
    
14. #---- Define convert function
    
15. def convert(inDir, outDir, sdate, fhour):
    
16.   hdate = sdate + datetime.timedelta(hours=fhour)
    
17.   hdatestr = hdate.strftime('%Y-%m-%d_%H:00:00')
    
18.   print hdatestr
    
19. 
    
20.   #---- Set intput and output data file
    
21.   infn = os.path.join(inDir, 'gmf.639.' + sdate.strftime('%Y%m%d%H') + \
    
22.     '%03d.grb2'%fhour)
    
23.   print 'T639 file: ' + infn
    
24.   outfn = os.path.join(outDir, 'T639_' + sdate.strftime('%Y-%m-%d') + '_%02d'%fhour)
    
25.   print 'MM5 intermediate file: ' + outfn
    
26. 
    
27.   #---- Read a GRIB data file
    
28.   mydata = MeteoDataInfo()
    
29.   mydata.openNetCDFData(infn)
    
30.   print 'GRIB file has been opened...'
    
31. 
    
32.   #---- Set output MM5 intermediate data info
    
33.   mm5DI = MM5IMDataInfo()
    
34. 
    
35.   #---- Set variable and level list
    
36.   gvars = ['Geopotential_height_isobaric','Temperature_isobaric',\
    
37.     'u-component_of_wind_isobaric','v-component_of_wind_isobaric',\
    
38.     'Relative_humidity_isobaric','Pressure_reduced_to_MSL_msl',\
    
39.     'Temperature_surface']
    
40.   mvars = ['HGT','T','U','V','RH','PMSL','SST']
    
41.   units = ['m','K','m s{-1}','m s{-1}','%','Pa','K']
    
42.   descs = ['Height','Temperature','U','V','Relative Humidity',\
    
43.     'Sea-level Pressure','Sea-surface Temperature']
    
44. 
    
45.   #---- Write MM5 intermediate data file
    
46.   dataInfo = mydata.getDataInfo()
    
47.   mm5DI.createDataFile(outfn)
    
48.   vNum = len(gvars)
    
49.   for i in range(0, vNum):
    
50.     #Read input grid data
    
51.     varName = gvars
    
52.     print 'Variable: ' + mvars
    
53.     var = dataInfo.getVariable(varName)   
    
54.     levNum = var.getLevelNum()
    
55.     print 'Level number: ' + str(levNum)
    
56.     for l in range(0, levNum):
    
57.       mydata.setLevelIndex(l)
    
58.       gData = mydata.getGridData(varName)
    
59.       dh = DataHead()
    
60.       dh.iversion = 3
    
61.       dh.hdate = hdatestr
    
62.       dh.xfcst = 0.0
    
63.       dh.field = mvars
    
64.       dh.units = units
    
65.       dh.desc = descs
    
66.       if levNum <= 1:
    
67.         dh.level = 200100.0
    
68.       else:
    
69.         dh.level = var.getLevels().get(l)
    
70.       print 'Level: ' + str(dh.level)
    
71.       dh.idim = gData.getXNum()
    
72.       dh.jdim = gData.getYNum()
    
73.       dh.llflag = 0
    
74.       dh.startlat = gData.getYMax()
    
75.       dh.startlon = gData.getXMin()
    
76.       dh.deltalat = -gData.getYDelt()
    
77.       dh.deltalon = gData.getXDelt()
    
78.       #Write data
    
79.       mm5DI.writeDataHead(dh)
    
80.       mm5DI.writeGridData(gData)
    
81. 
    
82.   mm5DI.closeDataFile()
    
83. 
    
84. #---- Run program
    
85. #---- Set datetime
    
86. stime = datetime.datetime.now()
    
87. sdate = datetime.datetime.now()
    
88. sdate = sdate.replace(hour=0)
    
89. sdate = sdate - datetime.timedelta(days=2)
    
90. print sdate
    
91. 
    
92. #---- Set directories
    
93. inDir = "D:/Working/MIScript/baoding/rdT639"
    
94. outDir = "D:/Working/MIScript/baoding/rdT639"
    
95. 
    
96. #---- Convert data
    
97. #convert(inDir, outDir, sdate, 18)
    
98. for fhour in range(0, 109, 6):
    
99.   convert(inDir, outDir, sdate, fhour)
    
100. 
     
101. print 'Run time: %s'%(datetime.datetime.now() - stime)
     
102. print 'Finished!'

#--------------------------------------------------------        <br /> # Author: Yaqiang Wang                                           <br /> # Date: 2014-5-8                                               <br /> # Purpose: Convert GRIB data to MM5 intermediate data  <br /> # Note: Sample                                                   <br /> #-----------------------------------------------------------<br /> from org.meteoinfo.data.meteodata import MeteoDataInfo<br /> from org.meteoinfo.data.meteodata.mm5 import MM5IMDataInfo<br /> from org.meteoinfo.data.meteodata import Variable<br /> from org.meteoinfo.data.meteodata.mm5 import DataHead<br /> import os<br /> import datetime<br /> <br /> #---- Define convert function<br /> def convert(inDir, outDir, sdate, fhour):<br />   hdate = sdate + datetime.timedelta(hours=fhour)<br />   hdatestr = hdate.strftime('%Y-%m-%d_%H:00:00')<br />   print hdatestr<br /> <br />   #---- Set intput and output data file<br />   infn = os.path.join(inDir, 'gmf.639.' + sdate.strftime('%Y%m%d%H') + \<br />     '%03d.grb2'%fhour)<br />   print 'T639 file: ' + infn<br />   outfn = os.path.join(outDir, 'T639_' + sdate.strftime('%Y-%m-%d') + '_%02d'%fhour)<br />   print 'MM5 intermediate file: ' + outfn<br /> <br />   #---- Read a GRIB data file<br />   mydata = MeteoDataInfo()<br />   mydata.openNetCDFData(infn)<br />   print 'GRIB file has been opened...'<br /> <br />   #---- Set output MM5 intermediate data info<br />   mm5DI = MM5IMDataInfo()<br /> <br />   #---- Set variable and level list<br />   gvars = ['Geopotential_height_isobaric','Temperature_isobaric',\<br />     'u-component_of_wind_isobaric','v-component_of_wind_isobaric',\<br />     'Relative_humidity_isobaric','Pressure_reduced_to_MSL_msl',\<br />     'Temperature_surface']<br />   mvars = ['HGT','T','U','V','RH','PMSL','SST']<br />   units = ['m','K','m s{-1}','m s{-1}','%','Pa','K']<br />   descs = ['Height','Temperature','U','V','Relative Humidity',\<br />     'Sea-level Pressure','Sea-surface Temperature']<br /> <br />   #---- Write MM5 intermediate data file<br />   dataInfo = mydata.getDataInfo()<br />   mm5DI.createDataFile(outfn)<br />   vNum = len(gvars)<br />   for i in range(0, vNum):<br />     #Read input grid data<br />     varName = gvars<i><br />     print 'Variable: ' + mvars<i><br />     var = dataInfo.getVariable(varName)    <br />     levNum = var.getLevelNum()<br />     print 'Level number: ' + str(levNum)<br />     for l in range(0, levNum):<br />       mydata.setLevelIndex(l)<br />       gData = mydata.getGridData(varName)<br />       dh = DataHead()<br />       dh.iversion = 3<br />       dh.hdate = hdatestr<br />       dh.xfcst = 0.0<br />       dh.field = mvars<i><br />       dh.units = units<i><br />       dh.desc = descs<i><br />       if levNum <= 1:<br />         dh.level = 200100.0<br />       else:<br />         dh.level = var.getLevels().get(l)<br />       print 'Level: ' + str(dh.level)<br />       dh.idim = gData.getXNum()<br />       dh.jdim = gData.getYNum()<br />       dh.llflag = 0<br />       dh.startlat = gData.getYMax()<br />       dh.startlon = gData.getXMin()<br />       dh.deltalat = -gData.getYDelt()<br />       dh.deltalon = gData.getXDelt()<br />       #Write data<br />       mm5DI.writeDataHead(dh)<br />       mm5DI.writeGridData(gData)<br /> <br />   mm5DI.closeDataFile()<br /> <br /> #---- Run program<br /> #---- Set datetime<br /> stime = datetime.datetime.now()<br /> sdate = datetime.datetime.now()<br /> sdate = sdate.replace(hour=0)<br /> sdate = sdate - datetime.timedelta(days=2)<br /> print sdate<br /> <br /> #---- Set directories<br /> inDir = "D:/Working/MIScript/baoding/rdT639"<br /> outDir = "D:/Working/MIScript/baoding/rdT639"<br /> <br /> #---- Convert data<br /> #convert(inDir, outDir, sdate, 18)<br /> for fhour in range(0, 109, 6):<br />   convert(inDir, outDir, sdate, fhour)<br /> <br /> print 'Run time: %s'%(datetime.datetime.now() - stime)<br /> print 'Finished!'

转换前的GRIB2数据图形：



转换后的MM5中间文件数据图形：