MeteoInfoLab脚本示例：GRIB转MICAPS4

MeteoInfo

GRIB数据转为MICAP4格式数据的例子：

1. fn = 'D:/Temp/grib/fnl_20110416_00_00'
   
2. f = addfile(fn)
   
3. v = f['Temperature_isobaric']
   
4. zn = v.dimlen(1)
   
5. data = v[0,zn-1,:,:]    #First level from ground
   
6. data = data - 273.16
   
7. #Save as MICAPS4 data file
   
8. outfn = 'D:/Temp/grib/grib2m4.000'
   
9. t = f.gettime(0)
   
10. data.savegrid(outfn, format='micaps4', description='Temperature', \
    
11.     date=t, hours=0, level=0, smooth=1, boldvalue=0)

fn = 'D:/Temp/grib/fnl_20110416_00_00'<br /> f = addfile(fn)<br /> v = f['Temperature_isobaric']<br /> zn = v.dimlen(1)<br /> data = v[0,zn-1,:,:]    #First level from ground<br /> data = data - 273.16<br /> #Save as MICAPS4 data file<br /> outfn = 'D:/Temp/grib/grib2m4.000'<br /> t = f.gettime(0)<br /> data.savegrid(outfn, format='micaps4', description='Temperature', \<br />     date=t, hours=0, level=0, smooth=1, boldvalue=0)

批量转换：

1. fn = 'D:/Temp/grib/Z_NAFP_C_BABJ_2016041800_P_gmf-639-2016041800003.grib1'
   
2. outdir = 'E:/Temp/test'
   
3. f = addfile(fn)
   
4. t = f.gettime(0)
   
5. vnames = ['Pressure_surface', 'Temperature_isobaric']
   
6. for vname in vnames:
   
7.     v = f[vname]
   
8.     print vname
   
9.     if v.ndim == 3:
   
10.         data = v[0,:,:]
    
11.         #Save as MICAPS4 data file
    
12.         outfn = os.path.join(outdir, 'm4_' + vname + '_surface' + '.000')        
    
13.         data.savegrid(outfn, format='micaps4', description=vname + '_1000', \
    
14.             date=t, hours=0, level=0, smooth=1, boldvalue=0)
    
15.     elif v.ndim == 4:
    
16.         #Get levels
    
17.         levels = v.dimvalue(1)
    
18.         for idx in range(0, len(levels)):
    
19.             level = levels[idx]
    
20.             print '\t' + str(level)
    
21.             data = v[0,idx,:,:]
    
22.             #Save as MICAPS4 data file
    
23.             outfn = os.path.join(outdir, 'm4_' + vname + '_' + str(int(level)) + '.000')        
    
24.             data.savegrid(outfn, format='micaps4', description=vname + '_' + \
    
25.                 str(int(level)), date=t, hours=0, level=0, smooth=1, boldvalue=0)

fn = 'D:/Temp/grib/Z_NAFP_C_BABJ_2016041800_P_gmf-639-2016041800003.grib1'<br /> outdir = 'E:/Temp/test'<br /> f = addfile(fn)<br /> t = f.gettime(0)<br /> vnames = ['Pressure_surface', 'Temperature_isobaric']<br /> for vname in vnames:<br />     v = f[vname]<br />     print vname<br />     if v.ndim == 3:<br />         data = v[0,:,:]<br />         #Save as MICAPS4 data file<br />         outfn = os.path.join(outdir, 'm4_' + vname + '_surface' + '.000')        <br />         data.savegrid(outfn, format='micaps4', description=vname + '_1000', \<br />             date=t, hours=0, level=0, smooth=1, boldvalue=0)<br />     elif v.ndim == 4:<br />         #Get levels<br />         levels = v.dimvalue(1)<br />         for idx in range(0, len(levels)):<br />             level = levels[idx]<br />             print '\t' + str(level)<br />             data = v[0,idx,:,:]<br />             #Save as MICAPS4 data file<br />             outfn = os.path.join(outdir, 'm4_' + vname + '_' + str(int(level)) + '.000')        <br />             data.savegrid(outfn, format='micaps4', description=vname + '_' + \<br />                 str(int(level)), date=t, hours=0, level=0, smooth=1, boldvalue=0)

my251394667

王老师，提示DimArray instance has no attribute 'savegrib'

MeteoInfo

my251394667 发表于 2016-5-23 19:25
王老师，提示DimArray instance has no attribute 'savegrib'

下载最新版本

my251394667

我下的是MeteoInfo_Java_1.3.2R13_Files，不是最新的吗

MeteoInfo

my251394667 发表于 2016-5-24 08:41
我下的是MeteoInfo_Java_1.3.2R13_Files，不是最新的吗

savegrid 不是 savegrib

my251394667

MeteoInfo 发表于 2016-5-24 10:29
savegrid 不是 savegrib

'NoneType' object has no attribute 'savegrid'。怎么遇到这个问题，对不起王老师，老是麻烦你

my251394667

my251394667 发表于 2016-5-24 15:53
'NoneType' object has no attribute 'savegrid'。怎么遇到这个问题，对不起王老师，老是麻烦你

王老师这个能不能一次性将grib2文件中的所有数据转成micaps4格式的？

MeteoInfo

my251394667 发表于 2016-5-24 16:02
王老师这个能不能一次性将grib2文件中的所有数据转成micaps4格式的？

可以用循环输出多个变量、时次、高度等。要不你给一个你的GRIB文件，我抽空帮你写个转换的脚本。

my251394667

MeteoInfo 发表于 2016-5-24 16:57
可以用循环输出多个变量、时次、高度等。要不你给一个你的GRIB文件，我抽空帮你写个转换的脚本。

就是这个。麻烦你啦王老师

MeteoInfo

my251394667 发表于 2016-5-24 18:33
就是这个。麻烦你啦王老师

已经在1楼给出。