热带太平洋海温SST年平均趋势的程序设计

雷暴大风天气

'reinit'
'sdfopen c:\OpenGrADS\1.nc'
'set grads off'
'set grid off'
'set lev 0'
'set lat -30 30'
'set lon 120 280'
'c:\OpenGrADS\output.gs'
'set gxout shaded'
'define a1=ave(sst-273.15,t=1,t=1620,12)'
'define a2=ave(sst-273.15,t=2,t=1620,12)'
'define a3=ave(sst-273.15,t=3,t=1620,12)'
'define a4=ave(sst-273.15,t=4,t=1620,12)'
'define a5=ave(sst-273.15,t=5,t=1620,12)'
'define a6=ave(sst-273.15,t=6,t=1620,12)'
'define a7=ave(sst-273.15,t=7,t=1620,12)'
'define a8=ave(sst-273.15,t=8,t=1620,12)'
'define a9=ave(sst-273.15,t=9,t=1620,12)'
'define a10=ave(sst-273.15,t=10,t=1620,12)'
'define a11=ave(sst-273.15,t=11,t=1620,12)'
'define a12=ave(sst-273.15,t=12,t=1620,12)'
i=1
while(i<=1620)
'define s0=sst(t='i')-a1'
'define s1=sst(t='i+1')-a2'
'define s2=sst(t='i+2')-a3'
'define s3=sst(t='i+3')-a4'
'define s4=sst(t='i+4')-a5'
'define s5=sst(t='i+5')-a6'
'define s6=sst(t='i+6')-a7'
'define s7=sst(t='i+7')-a8'
'define s8=sst(t='i+8')-a9'
'define s9=sst(t='i+9')-a10'
'define s10=sst(t='i+10')-a11'
'define s11=sst(t='i+11')-a12'
'define cc=s0+s1+s2+s3+s4+s5+s6+s7+s8+s9+s10+s11'
' bb=s0 * ('i')+s1 * ('i+1')+s2 * ('i+2')+s3 * ('i+3')+s4 * ('i+4')+s5 * ('i+5')+s6 * ('i+6')+s7 * ('i+7')+s8 * ('i+8')+s9 * ('i+9')+s10 * ('i+10')+s11 * ('i+11')'
'define ee=(('i')*('i')+('i'+1)*('i'+1)+('i'+2)*('i'+2)+('i'+3)*('i'+3)+('i'+4)*('i'+4)+('i'+5)*('i'+5)+('i'+6)*('i'+6)+('i'+7)*('i'+7)+('i'+8)*('i'+8)+('i'+9)*('i'+9)+('i'+10)*('i'+10)+('i'+11)*('i'+11))-(12*'i'+66)*(12*'i'+66)/12 '
'define dd=(bb-('i'+5.5)*cc)/ee'
i=i+12
endwhile
'd dd'
'set cint 6'
'cbarn.gs'
'printim c:\OpenGrADS\0.png white'
'reinit'

趋势值是一元线性回归方程B值
这个是我自己写的gs程序，可以正常出图，由于对grads了解不是很好，感觉写的程序很繁琐。希望各位前辈指导简便方法。

Houckin

没细看，不过学到了一些东西