- 积分
- 2591
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-4-4
- 最后登录
- 1970-1-1
|
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'
但是出来的图区别都非常小,而且就算弄成全球数据也是很小,而文献上我看区别都很大啊,没做过合成分析,论文上要的,希望大神指点一下...我实在想不出来。
|
|