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

a292295

这是我的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'


但是出来的图区别都非常小，而且就算弄成全球数据也是很小，而文献上我看区别都很大啊，没做过合成分析，论文上要的，希望大神指点一下...我实在想不出来。