【分享】Grads批处理AWX文件CTL、GS文件分享

Crazy_Merlin

最近自己尝试用Grads批处理读取FY2E 晴空大气可降水量数据（AWX格式）。期间遇到了许多问题，但是在坛子里各位高手的帮助下，问题一一得到了解决。最后贴出我最后的CTL、gs文件和大家分享，希望大家共同学习，共同提高！

dset I:\TDDOWNLOAD\TPW\20120826\FY2E_TPW_MLT_OTG_%y4%m2%d2_%h200.AWX
title FY-2E_TPW data
options template
options yrev
undef 99
*99为无效值
xdef 1001 linear  55 0.1
ydef 1001 linear -50 0.1
zdef 1 levels 850
tdef 11 linear 00Z25aug2012 3hr
vars 1
tpw  0 -1,40,1 FY-2E_TPW (mm)
endvars

以上是CTL文件，大家在使用时只要根据自己使用的数据对(x,y,t)及路径进行修改！

'reinit'
'open I:\TDDOWNLOAD\TPW\TPW.ctl'
'set mpdraw on'
'set poli on'
'set mpdset  mres hun anh jiangx hub sic guiz'
*'draw map'
tt=1
while(tt<=11)
'set t 'tt
** http://bbs.06climate.com/forum.php?mod=viewthread&tid=6268 这是显示时间时参考的gs文件出处
     'q time'
     res = subwrd(result,3)
     year= substr(res,9,12)
      'q time'
     res = subwrd(result,3)
     month=substr(res,6,3)
     'q time'
     res = subwrd(result,3)
     day=substr(res,4,2)
     'q time'
     res = subwrd(result,3)
     hour = substr(res,1,2)
***************************************************************************************
'set gxout shaded'
'set cint 10'
'set cmin 0.0'
*'set cmax 75'
'set xlopts 1 2 0.15'
'set ylopts 1 2 0.15'
'set lon 70 135'
'set lat 15 50'
'set grads off'
'set grid on'
'set xlint 5'
'set ylint 5'
'I:\TDDOWNLOAD\TPW\rb_matlab.gs'
'd smth9(tpw)'
'I:\TDDOWNLOAD\TPW\cbar_matlab.gs 5 1 0.8'
'draw title FY-2E_TPW_MLT_OTG_'year''month''day'_'hour'00'
*'draw rec 1 1 2 2'
*'draw string 2.3 0.3 time='tt''
if(rc='0')break;endif;
*'enable print I:\TDDOWNLOAD\TPW\FY2E_TPW'year''month''day''hour'00.gmf'
'printim I:\TDDOWNLOAD\TPW\FY2E_TPW'year''month''day''hour'00.png white'
*'print'
say tt '   ' year month day hour '      completed!'
'clear'
tt=tt+1
endwhile
pull dummy
*'disable print'
return
;

以上是gs文件，这里需要改动的地方比较多，仍然要具体问题具体分析。里面要用到两个gs文件（见附件，均来自坛子），其中rb_matlab.gs已经针对晴空大气可降水量数据进行修改（范围0-75mm），大家可以根据自己的数据范围进行修改（只要修改set clevs 后面的数值）

最后附上我是用的数据包（见附件）和示例图一张。有什么问题希望各位及时提出，大家共同提高！谢谢！

言深深

辛苦，支持原创，欢迎多多交流学习工作心得

做个霸气的木头

这样的总结可以给别人提供帮助，感谢楼主啊~

mofangbao

很不错一个总结，充分利用了论坛的现有工具，这个总结一定能帮助不少人的

Crazy_Merlin

言深深 发表于 2012-8-30 11:22

                               
登录/注册后可看大图

辛苦，支持原创，欢迎多多交流学习工作心得

呵呵，谢谢支持！今后大家多多关照！

Crazy_Merlin

做个霸气的木头 发表于 2012-8-30 11:27

                               
登录/注册后可看大图

这样的总结可以给别人提供帮助，感谢楼主啊~

呵呵，我也是从论坛里学来的~谢谢版主！

Crazy_Merlin

mofangbao 发表于 2012-8-30 12:22

                               
登录/注册后可看大图

很不错一个总结，充分利用了论坛的现有工具，这个总结一定能帮助不少人的

希望大家都能分享自己的学习心得~一起提高！

sharkmo

很不错啊，正在处理这个数据要用到！！！

Crazy_Merlin

sharkmo 发表于 2012-8-30 13:18

                               
登录/注册后可看大图

很不错啊，正在处理这个数据要用到！！！

呵呵，希望能够有用~有什么问题大家相互交流！

瓶子

挺好的，先下载下来，慢慢研究，会有用的，grdas啊还是不熟悉