[GrADS] 分月输出日资料，并求月极值位置

兰溪之水

1. 
   
2. 
   
3. *****************************************************
   
4. *                                                   *
   
5. *  该 GrADS 脚本用于寻找时间序列的某个月中变量最大  *
   
6. *  值在某年该月中的某日。                           *
   
7. *                                                   *
   
8. *    Created by 兰溪之水, Feb 18, 2012.             *
   
9. *                                                   *
   
10. *****************************************************
    
11. datadir = 'H:/ocean/output/'            ;* 数据所在路径,最后的"/"不能省
    
12. datafname = '1979_2011.air.sig995.ctl'  ;* 要打开的数据文件
    
13. var = 'air'                             ;* 变量名
    
14. * 打开要处理的数据文件
    
15. 'reinit'
    
16. datafile = datadir''datafname
    
17. 'open 'datafile   ;* $$$ Note: 若为其他格式数据这里需修改
    
18. 
    
19. *####################################################
    
20. *------------- 以下无特殊情况,毋需修改 --------------
    
21. *####################################################
    
22. * 获取时间维最大时次
    
23. 'q ctlinfo'
    
24. nline = 1
    
25. while(1)
    
26.   infolin = sublin(result,nline)
    
27.   word1 = subwrd(infolin,1)      
    
28.   if(word1='tdef')           ;* 判断是否到了时间行,
    
29.     tmax = subwrd(infolin,2) ;* 是,则读取最大时次
    
30.     break                    ;* 并退出循环
    
31.   endif
    
32.   nline = nline + 1
    
33. endwhile
    
34. * 定义月份字符串
    
35. mon='JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC'
    
36. * 以月份为外循环
    
37. nmon = 1
    
38. while(nmon<=12)
    
39.   tmon = subwrd(mon,nmon) ;* 得到当前要统计的月份
    
40. *============== 对所有时次统计一遍 ==============
    
41.   maxval = -999.9
    
42.   minval = 999.9
    
43. * -----------------------------------------------------------  
    
44.   t0 = 1
    
45.   while(t0<=tmax)  ;* 全部时次为内循环
    
46.   
    
47.     'set t 't0
    
48.     'q time'
    
49.     time = subwrd(result,3)
    
50.     ttmon = substr(time,6,3)        ;* 得到当前时次所到的月份
    
51.     if(ttmon = tmon)                ;* 判断,若不是当前统计月份
    
52.       'd 'var                       ;* 则进入下一时次
    
53.       temp = subwrd(result,4)
    
54.       if(temp>maxval); maxval = temp; maxvalloc = time; endif
    
55.       if(temp<minval); minval = temp; minvalloc = time; endif
    
56.     endif
    
57.    
    
58.     t0 = t0 + 1
    
59.   endwhile
    
60. * -----------------------------------------------------------   
    
61. * 输出：  
    
62.   say ' 'tmon': '
    
63.   say ' max: 'maxval', at 'maxvalloc
    
64.   say ' min: 'minval', at 'minvalloc
    
65.   say ' -------------'
    
66.    
    
67.   nmon = nmon + 1
    
68. endwhile
    
69. 'reinit'
    
70. ;

复制代码

topmad

兰溪之水 发表于 2012-2-18 02:57

                               
登录/注册后可看大图

兰溪 你悲催了....

兰溪之水

topmad 发表于 2012-2-18 10:03

                               
登录/注册后可看大图

兰溪 你悲催了....

mofangbao

topmad 发表于 2012-2-18 10:03

                               
登录/注册后可看大图

兰溪 你悲催了....

竟然你是沙发。。

topmad

mofangbao 发表于 2012-2-18 10:05

                               
登录/注册后可看大图

竟然你是沙发。。

阴差阳错....呃

topmad

兰溪之水 发表于 2012-2-18 10:05

                               
登录/注册后可看大图

给你盖个合格的戳真难....

kc121

你所用的资料是什么？'1979_2011的日资料还是一天四个时次的资料呀。这工作是很好，可用于哪种数据没说清，没法使用。谢谢！

平流层的萝卜

顶，学习！

river

兰溪之水 发表于 2012-2-18 02:57

                               
登录/注册后可看大图

原来兰溪哥哥是grads的程序猿啊

兰溪之水

river 发表于 2013-4-5 12:24

                               
登录/注册后可看大图

原来兰溪哥哥是grads的程序猿啊

。。。。。。。。工作需要