- 积分
- 370
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2011-9-23
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
数据来源于国家卫星气象中心网站,下载的是FY-2E 3小时降水估计产品(AWX格式)。因为之前下载过云顶亮温的awx文件并参照http://bbs.06climate.com/forum.php?mod=viewthread&tid=4982 用grads成功处理过,所以这次也打算这样处理,可是在用grads处理时遇到问题,请大家帮忙查找错误原因,不甚感激。
ctl文件如下:
dset H:\AWX\FY2E_PRE_003_OTG_2010%m2%d2_%h200.AWX
* dset H:\AWX\FY2E_PRE_003_OTG_20110808_0500.AWX
title 2011 twpice fy2e_pre
options template
options yrev
* yrev表示y轴反向
undef 0.0
xdef 1201 linear 45 0.1
ydef 1201 linear -60 0.1
zdef 1 linear 1000 1
tdef 1 linear 05z08Aug2011 1hr
vars 1
PRE 1 -1,40,1 Geopotential Heights (m)
endvars
gs文件如下:
'reinit'
'open h:\AWX\2011pre.ctl'
**************************
tt=1
while(tt<9)
*******控制坐标轴*******
'set xlint 5'
'set ylint 5'
'set xlopts 1 7 0.14'
'set ylopts 1 7 0.14'
*************************
'set t ' tt
'set grads off'
'set grid off'
'set lon 100 140'
'set lat 20 40'
'set gxout shaded'
'set clevs 0 10 20 30 40 50 60 '
'set ccols 9 4 13 12 7 8 2 6 '
*******添加地图********
'set mpdset cnworld '
'set map 1 1 8'
'd PRE'
'set mpdset cnriver'
'set map 0 1 8'
'draw map'
**********************
'cbarn 1 0 5 8'
'draw string 2.3 0.2 time='tt-1''
if(rc='0')break;endif;
'enable print f:\FY2E\pre_20110808.gmf'
'print'
'clear'
tt=tt+1
endwhile
pull dummy
'disable print'
'clear'
return
在grads界面中运行gs文件时报错如下:
提示出错后,画得图里面只有叠加的地图资料,并没有降水的信息。
实在不知道错误出在什么地方,请各位赐教,谢谢
|
|