登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
前言:好些日子没发帖了。。。近来,论坛上提问数大大增多呀,真正愿意分享的却不多了,就连以前找到的好书、好资料分享也少了。。。。
其实,一些帖子问的都是非常基础的,只要稍微搜索一下或者看一下manul什么的都能搞定。
这种风气可能成了一种恶性循环,有些人一见大家都是来问问题的,本来愿意分享的也不敢了,怕被人拍砖挑刺。
作为陪着论坛一路走来的老人吧,我觉得大家还是多多交流多多讨论,而不是简单的抛出个问题,让别人帮着做。或者一些隐藏的高手们呀,拿出来自己的藏家之宝吧。
为了鼓励交流、分享、讨论等学术上的优良品质,我就抛砖引玉吧,希望大家挑刺猛拍。
4月22日更新:源代码经过兰溪@兰溪之水,清风@mofangbao 等指点,改写成如下形式,简洁很多,特此致谢。并且取消隐藏设置,直接贴出代码了。
事情起源于@黎大页子提供了一些数据,要用GrADS来画出数据的标准差图。数据的ctl如下,是26年9496天的逐日通量数据,希望画出逐年的标准差图:
- dset D:\sd\26_eof.dat
- title daily mean surface sensible heat flux
- undef 32766
- xdef 61 linear 119.5 1.0
- ydef 67 linear -15.5 1.0
- tdef 9496 linear 00Z01JAN1985 1dy
- zdef 1 linear 0 1
- vars 1
- shtfl=>shtfl 0 99 eof mod of daily mean surface sensible heat
- endvars
最后的图选一张吧,给大家看看,gs见下:
- *********************************************************************
- * |variables uesd in this script:
- * |yr: Every year
- * |mean: Mean value of the variable shtfl.
- * |sd: Standard deviation of one year.
- * |4213: This is meant for Li.
- *----------------------------------------------------------
- *Author: Wayne Huai @ SYSU
- *Based on Liye's data.
- *Last updated: Apr 22, 2012
- *Special thanks to h t t p ://b bs.06climate.com/forum.php
- ******************************************
- 'reinit'
- 'open D:\sd\26_eof.ctl'
- 'set mpdset hires'
- yr=1985
- while(yr<=2010)
- *******************************************************
- *Standard deviation will be caculated
- 'set lat 0 50.5'
- 'set lon 119.5 179.5'
- *shtfl means surface latent heat flux
- 'define mean=ave(shtfl,time=00Z01JAN'yr',time=00Z31DEC'yr')'
- 'define sd=sqrt(ave(pow(shtfl-mean,2),time=00Z01JAN'yr',time=00Z31DEC'yr'))'
- *********************************************************
- *This is a skill to plot a shaded map with contour on it
- 'set grads off'
- 'set gxout shaded'
- 'd sd'
- 'cbarn'
- 'set gxout contour'
- 'set cthick 5'
- 'd sd'
- 'draw title Standard Deviation of 'yr' '
- 'printim D:\sd\4213\'yr'.png x1000 y800 white'
- 'c'
- yr=yr+1
- endwhile
- ;
|