爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3429|回复: 3

[分享资料] 求指点 解码grib2数据所遇到的问题

[复制链接]

新浪微博达人勋

发表于 2016-10-14 10:12:34 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
我想处理cfsv2的grib2数据,遇到了一些问题,求指点grib2数据是从cfsv2的网上下的,预测数据
有flxf2016101100.01.2016101100.grb2;flxf.01.2016101100.201610.avrg.grib.grb2;cdas1.t00z.ipvgrbf02.grib2;cdas1.t00z.ipvgrbf01.grib2。也有对应的idx文件,经过在meteoinfo上面的显示,证明数据应该没有太大问题。
查阅论坛上兰溪和其他的人帖子,正确使用了!g2ctl命令,生成的ctl里面也有东西,用!gribmap命令,生成的txt就有问题了,


111.ctl内容如下
dset d:/gc/cfsv2/flxf2016101100.01.2016101100.grb2
index d:/gc/cfsv2/flxf2016101100.01.2016101100.grb2.idx
undef 9.999E+20
title d:/gc/cfsv2/flxf2016101100.01.2016101100.grb2
* produced by g2ctl v0.1.2
* command line options: d:/gc/cfsv2/flxf2016101100.01.2016101100.grb2
* griddef=1:0:(384 x 190):grid_template=40:winds(N/S): Gaussian grid: (384 x 190) units 1e-06 input WE:NS output WE:SN number of latitudes between pole-equator=95 #points=72960 lat 89.277000 to -89.277000 lon 0.000000 to 359.062000 :winds(N/S)


dtype grib2
* ny=190 nlat=190 lat0=-89.277000 lat1=89.277000
xdef 384 linear 0.000000 0.937498694516971
ydef 190 levels
-89.277 -88.340 -87.397 -86.454 -85.509 -84.565 -83.620 -82.676 -81.731 -80.786 -79.841 -78.897 -77.952 -77.007
-76.062 -75.117 -74.173 -73.228 -72.283 -71.338 -70.393 -69.448 -68.503 -67.559 -66.614 -65.669 -64.724 -63.779
(略略略,太多不全写了)

tdef 1 linear 00Z11oct2016 1mo
zdef 1 linear 1 1
vars 101
ACONDsfc   0,1,0   2,0,228 ** surface Aerodynamic conductance [m/s]
ALBDOsfc   0,1,0   0,19,1 ** surface Albedo [%]
CDUVBsfc   0,1,0   0,4,195 ** surface Clear sky UV-B downward solar flux [W/m^2]

(略略略,太多不全写了)
ENDVARS


This exe file was created with the evaluation version of Perl2Exe.
For more information visit http://www.indigostar.com
(The full version does not display this message with a 2 second delay.)
...


生成的1112.txt如下
Open Error:  Unknown keyword in description file
  --> The invalid description file record is:
  --> This exe file was created with the evaluation version of Perl2Exe.
  The data file was not opened.

不知道该怎么继续下去。
单独用grads打开111.ctl,也显示错误

求指导



111.bmp

flxf2016101100.01.2016101100.grb2

4.07 MB, 下载次数: 0, 下载积分: 金钱 -5

数据文件

111.ctl

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

1112.txt

202 Bytes, 下载次数: 0, 下载积分: 金钱 -5

密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-10-14 10:24:38 | 显示全部楼层
按照
yesyesnotion的帖子
手动修改了下ctl文件,用grads可以直接打开而不出错误了
dset d:/gc/cfsv2/flxf2016101100.01.2016101100.grb2
*index d:/gc/cfsv2/flxf2016101100.01.2016101100.grb2.idx
undef 9.999E+20
title d:/gc/cfsv2/flxf2016101100.01.2016101100.grb2
*dtype grib2
xdef 384 linear 0.000000 0.937498694516971
ydef 14 levels
-89.277 -88.340 -87.397 -86.454 -85.509 -84.565 -83.620 -82.676 -81.731 -80.786 -79.841 -78.897 -77.952 -77.007
tdef 1 linear 00Z11oct2016 1mo
zdef 1 linear 1 1
vars 8
ACONDsfc   0,1,0   2,0,228 ** surface Aerodynamic conductance [m/s]
ALBDOsfc   0,1,0   0,19,1 ** surface Albedo [%]
CDUVBsfc   0,1,0   0,4,195 ** surface Clear sky UV-B downward solar flux [W/m^2]
CNWATsfc   0,1,0   2,0,196 ** surface Plant Canopy Surface Water [kg/m^2]
CPRATsfc   0,1,0   0,1,196 ** surface Convective Precipitation Rate [kg/m^2/s]
CRAINsfc   0,1,0   0,1,192 ** surface Categorical Rain (yes=1; no=0) [non-dim]
CSDLFsfc   0,1,0   0,5,196 ** surface Clear Sky Downward Long Wave Flux [W/m^2]
CSDSFsfc   0,1,0   0,4,196 ** surface Clear Sky Downward Solar Flux [W/m^2]
ENDVARS

但我还不知道原理,还是不知道!g2ctl 和!bribmap处理的为什么不对
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2016-10-14 10:34:54 | 显示全部楼层
可能g2ctl和gribmap和grads版本不配套
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2016-10-14 10:43:20 | 显示全部楼层
谢谢指导
g2ctl 和gribmap运行都没有提示错误
g2ctl一开始是不行,运行就出错误,后来参照了个达人的帖子。重新写了g2ctl.pl,下载ActivePerl-5.10.0.1003-MSWin32-x86-285500,p2x-9.110-Win32重新编译的g2ctl,这样就没有报错,能用了。
gribmap这个是从哪里更新,做到和grads配套。
grads版本是2.0.a9.oga.1
以前用的是grads2.1家园版,但是环境参数怎么都弄不对,就又换回来了
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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