请问ec下发的细网格资料（后缀为*.bin）的怎么生成idx文件用grads画图呢

clare

企图处理cmcast接收的ec细网格资料（0.25x0.25分辨率的），是bin结尾的。用grib2ctl得到描述文件后，用gribmap直接对bin文件生成idx出错。是不是要先转换数据格式呢？但是描述文件中又称其为grib格式。
ctl：
dset ^W_NAFP_C_ECMF_20150102054416_P_C1D01020000010200001.bin
index ^W_NAFP_C_ECMF_20150102054416_P_C1D01020000010200001.bin.idx
undef 9.999E+20
title W_NAFP_C_ECMF_20150102054416_P_C1D01020000010200001.bin
*  produced by grib2ctl v0.9.12.5p46
dtype grib 255
options yrev
ydef 281 linear -10.000000 0.25
xdef 361 linear 60.000000 0.250000
tdef 1 linear 00Z02jan2015 1mo
*  z has 7 levels, for prs
zdef 7 levels
1000 925 850 700 500 200 100
vars 22
no10Usfc  0 165,1,0  ** surface 10 metre U wind component [m s**-1]
no10Vsfc  0 166,1,0  ** surface 10 metre V wind component [m s**-1]
no2Dsfc  0 168,1,0  ** surface 2 metre dewpoint temperature [K]
no2Tsfc  0 167,1,0  ** surface 2 metre temperature [K]
ALsfc  0 174,1,0  ** surface Albedo [(0 - 1)]
Dprs 7 155,100,0 ** (profile) Divergence [s**-1]
GHprs 7 156,100,0 ** (profile) Height [gpm]
LCCsfc  0 186,1,0  ** surface Low cloud cover [(0 - 1)]
MSLsfc  0 151,1,0  ** surface Mean sea level pressure [Pa]
PVprs 7 60,100,0 ** (profile) Potential vorticity [K m**2 kg**-1 s**-1]
Qprs 7 133,100,0 ** (profile) Specific humidity [kg kg**-1]
Rprs 7 157,100,0 ** (profile) Relative humidity [%]
RSNsfc  0 33,1,0  ** surface Snow density [kg m**-3]
SDsfc  0 141,1,0  ** surface Snow depth [m of water equivalent]
SKTsfc  0 235,1,0  ** surface Skin temperature [K]
SSTKsfc  0 34,1,0  ** surface Sea surface temperature [K]
Tprs 7 130,100,0 ** (profile) Temperature [K]
TCCsfc  0 164,1,0  ** surface Total cloud cover [(0 - 1)]
TCWVsfc  0 137,1,0  ** surface Total column water vapour [kg m**-2]
Uprs 7 131,100,0 ** (profile) U velocity [m s**-1]
Vprs 7 132,100,0 ** (profile) V velocity [m s**-1]
Wprs 7 135,100,0 ** (profile) Vertical velocity [Pa s**-1]
ENDVARS


生成的idx的映射过程提示：
grib1map error: file is not GRIB version 0 or 1, version number is 2
grib1map error: GRIB file format error (rc = 99)


同时无idx文件生出。

clare

麻烦大神指教，是不是bin不能这么直接使用？~@淡水河

balfulosa

同样求问，据说是grib1、2的混合码，直接解码不行。不太懂具体要怎们弄、求高人提示一下步骤。。

river

试试使用g2ctl 命令生成ctl

clare

river 发表于 2015-1-4 07:49
试试使用g2ctl 命令生成ctl

谢谢回复。但貌似该数据的确由grib1和grib2混合的，采用g2ctl后只得到部分变量：
ctl如下：
dset ^W_NAFP_C_ECMF_20150102054416_P_C1D01020000010200001.bin
index ^W_NAFP_C_ECMF_20150102054416_P_C1D01020000010200001.bin.idx
undef 9.999E+20
title W_NAFP_C_ECMF_20150102054416_P_C1D01020000010200001.bin
*  produced by g2ctl v0.0.4b
* griddef=1:304560:(361 x 281):grid_template=0:winds(N/S): lat-lon grid:(361 x 281) units 1e-06 input WE:NS output WE:SN res 48 lat 60.000000 to -10.000000 by 0.250000 lon 60.000000 to 150.000000 by 0.250000 #points=101441:winds(N/S)

dtype grib2
ydef 281 linear -10.000000 0.25
xdef 361 linear 60.000000 0.250000
tdef 1 linear 00Z02jan2015 1mo
* PROFILE hybrid
zdef 91 levels 1 2 4 6 7 9 10 12 14 15 17 19 20 22 24 25 27 29 30 32 34 35 37 38 40 42 43 45 46 48 49 51 52 54 55 57 58 59 61 62 63 65 66 67 69 70 71 73 74 76 77 78 80 81 82 84 85 87 88 90 91 93 94 96 97 98 100 101 103 104 105 107 108 110 111 112 114 116 117 119 120 122 124 126 128 129 131 133 134 136 137
vars 4
NLPREShy1  0,105,1   0,3,25 ** 1 hybrid level none Natural Logarithm Of Pressure in Pa [Numeric]
var0183hlev   91,105  0,1,83 ** (137 136 134 133 131.. 7 6 4 2 1) none desc [unit]
var0184hlev   91,105  0,1,84 ** (137 136 134 133 131.. 7 6 4 2 1) none desc [unit]
var0632hlev   91,105  0,6,32 ** (137 136 134 133 131.. 7 6 4 2 1) none desc [unit]
ENDVARS

river

clare 发表于 2015-1-4 11:31
谢谢回复。但貌似该数据的确由grib1和grib2混合的，采用g2ctl后只得到部分变量：
ctl如下：
dset ^W_NA ...

那把得到的两个一起利用，不就是全了么

cherven

细网格数据采用的是复杂压缩格式。搞不懂怎么解压   论坛里面有人会  不知道愿不愿意分享出来