爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3832|回复: 0

grads分年份出力数据,出图差异很小

[复制链接]

新浪微博达人勋

发表于 2017-5-2 17:11:08 | 显示全部楼层 |阅读模式
GrADS
系统平台: grads
问题截图:
问题概况: 我的grads的GS文件如下文,我吧51~06年的数据分开计算,分成鼎盛年,衰弱年还有正常年(其实也不算正常吧我觉得,还有拉尼娜),但是出图的时候风速还有水汽变化都很小,我看文献的变化都很大啊
我看过提问的智慧: 看过
自己思考时长(天): 3

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
这是我的gs文件
'reinit'
'sdfopen f:\data2\rhum.mon.mean.nc'
'sdfopen  f:\data2\uwnd.mon.mean.nc'
'sdfopen  f:\data2\vwnd.mon.mean.nc'
'enable print  f:\data2\dingsheng700.gmf'
***********************************************************
'define nrhum12 = rhum.1*0 '
'define mm=0'
tt=0
yr=1948
while(tt<=831)
*鼎盛年
if(yr=1951|yr=1953|yr=1958|yr=1963|yr=1965|yr=1969|yr=1972|yr=1977|yr=1979|yr=1982|yr=1987|yr=1991|yr=1997|yr=2002|yr=2004|yr=2006|yr=2009|yr=2015)  
*衰弱年
*if(yr=1952|yr=1954|yr=1959|yr=1964|yr=1966|yr=1970|yr=1973|yr=1978|yr=1980|yr=1983|yr=1988|yr=1992|yr=1998|yr=2003|yr=2005|yr=2007|yr=2010|yr=2016)   
*正常年
*if(yr=1955|yr=1956|yr=1957|yr=1960|yr=1961|yr=1962|yr=1967|yr=1968|yr=1971|yr=1974|yr=1975|yr=1976|yr=1981|yr=1984|yr=1985|yr=1986|yr=1989|yr=1990|yr=1993|yr=1994|yr=1995|yr=1996|yr=1999|yr=2000|yr=2001|yr=2008|yr=2011|yr=2012|yr=2013|yr=2014)   
'nrhum12 = nrhum12 + rhum.1(t='tt+6')+rhum.1(t='tt+7')+rhum.1(t='tt+8')'
'mm=mm+3'
endif
tt=tt+12
yr=yr+1
endwhile
'nrhum12 = nrhum12/mm'     ****算出夏季鼎盛年平均水汽
***********************************************************
'define nuwnd12 =uwnd.2*0 '
'define mm=0'
tt=0
yr=1948
while(tt<=831)
*鼎盛年
if(yr=1951|yr=1953|yr=1958|yr=1963|yr=1965|yr=1969|yr=1972|yr=1977|yr=1979|yr=1982|yr=1987|yr=1991|yr=1997|yr=2002|yr=2004|yr=2006|yr=2009|yr=2015)  
*衰弱年
*if(yr=1952|yr=1954|yr=1959|yr=1964|yr=1966|yr=1970|yr=1973|yr=1978|yr=1980|yr=1983|yr=1988|yr=1992|yr=1998|yr=2003|yr=2005|yr=2007|yr=2010|yr=2016)   
*正常年
*if(yr=1955|yr=1956|yr=1957|yr=1960|yr=1961|yr=1962|yr=1967|yr=1968|yr=1971|yr=1974|yr=1975|yr=1976|yr=1981|yr=1984|yr=1985|yr=1986|yr=1989|yr=1990|yr=1993|yr=1994|yr=1995|yr=1996|yr=1999|yr=2000|yr=2001|yr=2008|yr=2011|yr=2012|yr=2013|yr=2014)   
'nuwnd12 = nuwnd12 + uwnd.2(t='tt+6')+uwnd.2(t='tt+7')+uwnd.2(t='tt+8')'
'mm=mm+3'
endif
tt=tt+12
yr=yr+1
endwhile
'nuwnd12 = nuwnd12/mm'             ****算出夏季衰弱年平均u风
***********************************************************
'define nvwnd12 = vwnd.3*0 '
'define mm=0'
tt=0
yr=1948
while(tt<=831)
*鼎盛年
if(yr=1951|yr=1953|yr=1958|yr=1963|yr=1965|yr=1969|yr=1972|yr=1977|yr=1979|yr=1982|yr=1987|yr=1991|yr=1997|yr=2002|yr=2004|yr=2006|yr=2009|yr=2015)  
*衰弱年
*if(yr=1952|yr=1954|yr=1959|yr=1964|yr=1966|yr=1970|yr=1973|yr=1978|yr=1980|yr=1983|yr=1988|yr=1992|yr=1998|yr=2003|yr=2005|yr=2007|yr=2010|yr=2016)   
*正常年
*if(yr=1955|yr=1956|yr=1957|yr=1960|yr=1961|yr=1962|yr=1967|yr=1968|yr=1971|yr=1974|yr=1975|yr=1976|yr=1981|yr=1984|yr=1985|yr=1986|yr=1989|yr=1990|yr=1993|yr=1994|yr=1995|yr=1996|yr=1999|yr=2000|yr=2001|yr=2008|yr=2011|yr=2012|yr=2013|yr=2014)   
'nvwnd12 = nvwnd12 + vwnd.3(t='tt+6')+vwnd.3(t='tt+7')+vwnd.3(t='tt+8')'
'mm=mm+3'
endif
tt=tt+12
yr=yr+1
endwhile
'nvwnd12 = nvwnd12/mm'             ****算出夏季衰弱年平均v风
***************************************************************
'set grads off'
'set grid off'
'set display color white'
'set mpdset cnworld'
'set lon 70 135'
'set lat 15 55'
'set lev 700'
define_colors
'set gxout shaded'
'set clevs -15 -10 -5 0 5 10 15 20 25 30 35 40'
'set ccols   47  45  43  41 21 22 23 24 25 26 27 28 29'
'd  nrhum12-ave(rhum.1,t=37,t=828,1)'
'set gxout contour'
'set clevs -15 -10 -5 0 5 10 15 20 25 30 35 40'
'd  nrhum12-ave(rhum.1,t=37,t=828,1)'     ***距平代表水汽异常
'set gxout vector'
'set arrscl 0.8 10'
'set ccolor  3'
'set arrowhead 0.08'
'set arrlab off'
'd skip(nuwnd12,1.5,1.5);skip(nvwnd12,1.5,1.5)'        
'print'
'c'
'disable print'
'reinit'


但是出来的图区别都非常小,而且就算弄成全球数据也是很小,而文献上我看区别都很大啊,没做过合成分析,论文上要的,希望大神指点一下...我实在想不出来。
QQ图片20170502170637.jpg
QQ图片20170502170647.jpg
QQ图片20170502170700.jpg
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表