【已解决】NARR资料求助~顺撒点小钱

飘渺居 · 发表于 2017-10-31 13:42:26

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由飘渺居于 2017-10-31 21:46 编辑

北美再分析资料（NARR）有两种资料（a和b），其中a为分析资料，b为预报资料。我下载了NARR-a，但是用grads处理时发现每个时刻的文件里面似乎包含了两个时刻（似乎也包含预报场？），因此批处理时时刻出错。
一个grb生成的ctl：
dset G:/data/NARR/narr-a_221_19790606_0000_000.grb
index G:/data/NARR/narr-a_221_19790606_0000_000.grb.idx
undef 9.999E+20
title G:/data/NARR/narr-a_221_19790606_0000_000.grb
*  produced by grib2ctl v0.9.12.5p39c
dtype grib 221
pdef 349 277 lcc 1 -145.5 1 1 50 50 -107 32463 32463
xdef 1111 linear -250 0.333333
ydef 247 linear 8 0.333333
tdef 1 linear 15Z01jun1979 1mo
*  z has 29 levels, for prs
zdef 29 levels
1000 975 950 925 900 875 850 825 800 775 750 725 700 650 600 550 500 450 400 350 300 275 250 225 200 175 150 125 100
vars 185
no4LFTX180_0mb  0 132,116,46080 ** 180-0 mb above gnd Best (4-layer) lifted index [K]
ACPCPsfc  0 63,1,0  ** surface Convective precipitation [kg/m^2]
ALBDOsfc  0 84,1,0  ** surface Albedo [%]
APCPsfc  0 61,1,0  ** surface Total precipitation [kg/m^2]

……
endvars
生成idx时的映射过程如下：
……
!!!!! MATCH: 372 48662074  221  1  12092 34  116 46230 48507603 48495506  13  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 373 48804920  221  1  12092 39  116 46230 48674255 48662158  11  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 374 48912130  221  1  12092  132  116 46080 48817101 48805004 8  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 375 49031220  221  1  12092  157  116 46080 48924311 48912214 9  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 376 49162188  221  1  12092  156  116 46080 49043401 49031304  10  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 377 49281278  221  1  12092 1 5 0    49174369 49162272 9  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 378 49412246  221  1  12092 84 1 0    49293459 49281362  10  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 379 49531336  221  1  12092 59 1 0    49424427 49412330 9  btim: 1979060600:00 tau:    3 dtim: 1979060603:00
!!!!! MATCH: 380 49686060  221  1  12092 61 1 0    49543517 49531420  12  btim: 1979060600:00 tau:    3 dtim: 1979060603:00
!!!!! MATCH: 381 49840784  221  1  12092  202 1 0    49698241 49686144  12  btim: 1979060600:00 tau:    3 dtim: 1979060603:00
!!!!! MATCH: 382 49995508  221  1  12092 63 1 0    49852965 49840868  12  btim: 1979060600:00 tau:    3 dtim: 1979060603:00
!!!!! MATCH: 383 50138354  221  1  12092 99 1 0    50007689 49995592  11  btim: 1979060600:00 tau:    3 dtim: 1979060603:00
……

后半部分似乎时刻增加了3个小时，但单个文件仍然都MATCH。。

而我批处理时将ctl改为如下：
dset G:/data/NARR/narr-a_221_19790606_%h200_000.grb
options template
index G:/data/NARR/narr-a_221_19790606.idx
undef 9.999E+20
title G:/data/NARR/19790606.grb
*  produced by grib2ctl v0.9.12.5p39c
dtype grib 221
pdef 349 277 lcc 1 -145.5 1 1 50 50 -107 32463 32463
xdef 1111 linear -250 0.333333
ydef 247 linear 8 0.333333
tdef 8 linear 00Z06jun1979 3hr
*  z has 29 levels, for prs
zdef 29 levels
1000 975 950 925 900 875 850 825 800 775 750 725 700 650 600 550 500 450 400 350 300 275 250 225 200 175 150 125 100
vars 185
no4LFTX180_0mb  0 132,116,46080 ** 180-0 mb above gnd Best (4-layer) lifted index [K]
ACPCPsfc  0 63,1,0  ** surface Convective precipitation [kg/m^2]

……
endvars
此时再生成idx时后半部分映射过程就会出错，如下：
……
!!!!! MATCH: 372 48662074  221  1  12092 34  116 46230 48507603 48495506  13  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 373 48804920  221  1  12092 39  116 46230 48674255 48662158  11  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 374 48912130  221  1  12092  132  116 46080 48817101 48805004 8  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 375 49031220  221  1  12092  157  116 46080 48924311 48912214 9  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 376 49162188  221  1  12092  156  116 46080 49043401 49031304  10  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 377 49281278  221  1  12092 1 5 0    49174369 49162272 9  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
!!!!! MATCH: 378 49412246  221  1  12092 84 1 0    49293459 49281362  10  btim: 1979060600:00 tau:    0 dtim: 1979060600:00
----- Time out of file limits: 379 49531336  221  1  12092 59 1 0    49424427 49412330 9  btim: 1979060600:00 tau:    3 dtim: 1979060603:00
----- Time out of file limits: 380 49686060  221  1  12092 61 1 0    49543517 49531420  12  btim: 1979060600:00 tau:    3 dtim: 1979060603:00
----- Time out of file limits: 381 49840784  221  1  12092  202 1 0    49698241 49686144  12  btim: 1979060600:00 tau:    3 dtim: 1979060603:00
----- Time out of file limits: 382 49995508  221  1  12092 63 1 0    49852965 49840868  12  btim: 1979060600:00 tau:    3 dtim: 1979060603:00
----- Time out of file limits: 383 50138354  221  1  12092 99 1 0    50007689 49995592  11  btim: 1979060600:00 tau:    3 dtim: 1979060603:00
----- Time out of file limits: 384 50189072  221  1  12092  235 1 0    50150535 50138438  11  btim: 1979060600:00 tau:    3 dtim: 1979060603:00
----- Time out of file limits: 385 50239790  221  1  12092  234 1 0    50201253 50189156  11  btim: 1979060600:00 tau:    3 dtim: 1979060603:00

……
ctl中的时间间隔改为3hr后，所有grb文件的后半部分映射都会出错。
请问有人处理过这种数据吗？一个文件如果包含了两个时刻的话，该如何批处理呢？求指点~~

飘渺居 · 发表于 2017-10-31 13:49:26

自己顶一下~

飘渺居 · 发表于 2017-10-31 13:49:55

{:5_231:}{:5_231:}

19921219 · 发表于 2017-10-31 16:09:51

虽然不太清楚，不过顶一下，希望能解决

rabin_xu · 发表于 2017-10-31 17:06:46

不懂是什么问题，帮顶

飘渺居 · 发表于 2017-10-31 21:45:56

已解决~我生成地图文件时在语句后加 -0 就可以了，即：
！gribmap -i ***.ctl -0
虽然不知道什么原因，但是问题解决啦~！有知道原因的小伙伴麻烦解释下呀

		自动登录	找回密码
密码			立即注册

[其他] 【已解决】NARR资料求助~顺撒点小钱

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

回帖奖励 +8 金钱

回帖奖励 +8 金钱


144 金钱	回复本帖可获得 8 金钱奖励! 每人限 1 次(中奖概率 80%)