python - Python 写入 GRIB 文件

Question

我想总结两个 grib 文件。我能够创建新变量（newfield），但我没有成功将它写入 grb 文件。当我执行“cdo sinfo newfile.grb”时，我收到消息：“不支持的文件类型”这是我的代码：

grib_dif = pygrib.open('file1.grb')

grib_dir = pygrib.open('file2.grb')

outfile = open('newfile.grb','wb')

#read grib1:
grib_dif.rewind()
dif_vals=[]
for mes1 in grib_dif:
        print mes1
        if (mes1.shortName == 'SWDIFDS_RAD'):
                fdif=mes1
        dif_vals.append(fdif.values)
        

#read grib2:
grib_dir.rewind()
dir_vals=[]
for mes2 in grib_dir:
        if (mes2.shortName == 'SWDIRS_RAD'):
                fdir=mes2
        dir_vals.append(fdir.values)

#sum up:
newfield=zeros((6,401,561),float)
for t in range(6):
         for i in range(401):
                for j in range(561):
                        newfield[t][i][j]=(dif_vals[t][i][j]+dir_vals[t][i][j])


#write output file:
outfile.write(newfield)
outfile.close()
grib_dif.close()
grib_dir.close()