GrADS中闰年判断及标准差的计算

artmunich

    前言：好些日子没发帖了。。。近来，论坛上提问数大大增多呀，真正愿意分享的却不多了，就连以前找到的好书、好资料分享也少了。。。。
    其实，一些帖子问的都是非常基础的，只要稍微搜索一下或者看一下manul什么的都能搞定。
    这种风气可能成了一种恶性循环，有些人一见大家都是来问问题的，本来愿意分享的也不敢了，怕被人拍砖挑刺。
    作为陪着论坛一路走来的老人吧，我觉得大家还是多多交流多多讨论，而不是简单的抛出个问题，让别人帮着做。或者一些隐藏的高手们呀，拿出来自己的藏家之宝吧。
    为了鼓励交流、分享、讨论等学术上的优良品质，我就抛砖引玉吧，希望大家挑刺猛拍。

    4月22日更新：源代码经过兰溪@兰溪之水，清风@mofangbao 等指点，改写成如下形式，简洁很多，特此致谢。并且取消隐藏设置，直接贴出代码了。



     事情起源于@黎大页子提供了一些数据，要用GrADS来画出数据的标准差图。数据的ctl如下，是26年9496天的逐日通量数据，希望画出逐年的标准差图：

1. dset D:\sd\26_eof.dat
   
2. title daily mean surface sensible heat flux
   
3. undef 32766
   
4. xdef 61 linear 119.5 1.0
   
5. ydef 67 linear -15.5 1.0
   
6. tdef 9496 linear 00Z01JAN1985 1dy
   
7. zdef 1 linear 0 1
   
8. vars 1
   
9. shtfl=>shtfl 0 99 eof mod of daily mean surface sensible heat
   
10. endvars
    

dset D:\sd\26_eof.dat<br /> title daily mean surface sensible heat flux<br /> undef 32766<br /> xdef 61 linear 119.5 1.0<br /> ydef 67 linear -15.5 1.0<br /> tdef 9496 linear 00Z01JAN1985 1dy<br /> zdef 1 linear 0 1<br /> vars 1<br /> shtfl=>shtfl 0 99 eof mod of daily mean surface sensible heat<br /> endvars<br />

   最后的图选一张吧，给大家看看，gs见下:

1. *********************************************************************
   
2. *  |variables uesd in this script:
   
3. *  |yr:     Every year
   
4. *  |mean:   Mean value of the variable shtfl.
   
5. *  |sd:     Standard deviation of one year.
   
6. *  |4213：  This is meant for Li.
   
7. *----------------------------------------------------------
   
8. *Author: Wayne Huai @ SYSU
   
9. *Based on Liye's data.
   
10. *Last updated: Apr 22, 2012
    
11. *Special thanks to h t t p ://b bs.06climate.com/forum.php
    
12. ******************************************
    
13. 
    
14. 'reinit'
    
15. 'open D:\sd\26_eof.ctl'
    
16. 'set mpdset hires'
    
17. yr=1985
    
18. while(yr<=2010)
    
19. 
    
20. *******************************************************
    
21. *Standard deviation will be caculated
    
22. 'set lat 0 50.5'
    
23. 'set lon 119.5 179.5'
    
24. 
    
25. *shtfl means surface latent heat flux
    
26. 'define mean=ave(shtfl,time=00Z01JAN'yr',time=00Z31DEC'yr')'
    
27. 'define sd=sqrt(ave(pow(shtfl-mean,2),time=00Z01JAN'yr',time=00Z31DEC'yr'))'
    
28. 
    
29. *********************************************************
    
30. *This is a skill to plot a shaded map with contour on it
    
31. 'set grads off'
    
32. 'set gxout shaded'
    
33. 'd sd'
    
34. 'cbarn'
    
35. 'set gxout contour'
    
36. 'set cthick 5'
    
37. 'd sd'
    
38. 'draw title Standard Deviation of 'yr' '
    
39. 'printim D:\sd\4213\'yr'.png x1000 y800 white'
    
40. 'c'
    
41. yr=yr+1
    
42. endwhile
    
43. ;

*********************************************************************<br /> *  |variables uesd in this script:<br /> *  |yr:     Every year<br /> *  |mean:   Mean value of the variable shtfl.<br /> *  |sd:     Standard deviation of one year.<br /> *  |4213：  This is meant for Li.<br /> *----------------------------------------------------------<br /> *Author: Wayne Huai @ SYSU<br /> *Based on Liye's data.<br /> *Last updated: Apr 22, 2012<br /> *Special thanks to h t t p ://b bs.06climate.com/forum.php<br /> ******************************************<br /> <br /> 'reinit'<br /> 'open D:\sd\26_eof.ctl'<br /> 'set mpdset hires'<br /> yr=1985<br /> while(yr<=2010)<br /> <br /> *******************************************************<br /> *Standard deviation will be caculated<br /> 'set lat 0 50.5'<br /> 'set lon 119.5 179.5'<br /> <br /> *shtfl means surface latent heat flux<br /> 'define mean=ave(shtfl,time=00Z01JAN'yr',time=00Z31DEC'yr')'<br /> 'define sd=sqrt(ave(pow(shtfl-mean,2),time=00Z01JAN'yr',time=00Z31DEC'yr'))'<br /> <br /> *********************************************************<br /> *This is a skill to plot a shaded map with contour on it <br /> 'set grads off'<br /> 'set gxout shaded'<br /> 'd sd'<br /> 'cbarn'<br /> 'set gxout contour'<br /> 'set cthick 5'<br /> 'd sd'<br /> 'draw title Standard Deviation of 'yr' '<br /> 'printim D:\sd\4213\'yr'.png x1000 y800 white'<br /> 'c'<br /> yr=yr+1<br /> endwhile<br /> ;

随缘

兰溪之水

A帮妹纸就有作品出了啊~

artmunich

兰溪之水 发表于 2012-4-22 16:20

                               
登录/注册后可看大图

A帮妹纸就有作品出了啊~

就算熬通宵也得废寝忘食的做你懂的

兰溪之水

刚看了A的大作，花了1个金钱，看来A花了不少力气，不想打击A，不过还是提点意见，其实你可以直接设置时间，这样可以避开闰年的判断~

artmunich

兰溪之水 发表于 2012-4-22 16:25

                               
登录/注册后可看大图

刚看了A的大作，花了1个金钱，看来A花了不少力气，不想打击A，不过还是提点意见，其实你可以直接设置时间， ...

我只听过options里面有这个设置，不过我在手册上没找到，你在documentation上看的吗？

黎大页子

兰溪之水 发表于 2012-4-22 16:25

                               
登录/注册后可看大图

刚看了A的大作，花了1个金钱，看来A花了不少力气，不想打击A，不过还是提点意见，其实你可以直接设置时间， ...

你是说用%y4的那种吗？

兰溪之水

artmunich 发表于 2012-4-22 16:32

                               
登录/注册后可看大图

我只听过options里面有这个设置，不过我在手册上没找到，你在documentation上看的吗？

类比摸索+实践证明的啦，不是有set lon .., set lat ..., set lev ...,那推理应该就有set time ... 你有看我之前发的帖子就知道我经常会这样用，看，不够给力吧，都不看兰溪哥哥的帖子的~@管理员 @审核员 @版主 @超级版主

mofangbao

当然有 set time了啊。。。A的精神可嘉啊

兰溪之水

黎大页子 发表于 2012-4-22 16:34

                               
登录/注册后可看大图

你是说用%y4的那种吗？

比如'define mean=ave(shtfl,t='idyb',t='idye')'这句
你可直接用'define mean=ave(shtfl,time=00Z01JAN'yr',time=00Z31DEC'yr')'