求大神帮忙看看，有关合成分析打点的问题

fh_550

从SST里面选出12个异常年和常年平均海温做比较，信度检验用到了t检验，仿照论坛的帖子写了程序，但是图出不来，请大神帮忙看看，是不是程序有问题。以下是代码：
'sdfopen f:\m.nc''set grads off''set ccolor rainbow''define t1 = ave(sst,time=jan1965,time=dec1965)''define t2 = ave(sst,time=jan1977,time=dec1977)''define t3 = ave(sst,time=jan1979,time=dec1979)''define t4 = ave(sst,time=jan1981,time=dec1981)''define t5 = ave(sst,time=jan1986,time=dec1986)''define t6 = ave(sst,time=jan1992,time=dec1992)''define t7 = ave(sst,time=jan1997,time=dec1997)''define t8 = ave(sst,time=jan1999,time=dec1999)''define t9 = ave(sst,time=jan2012,time=dec2012)''define t10 = ave(sst,time=jan2008,time=dec2008)''define t11 = ave(sst,time=jan2009,time=dec2009)''define t12 = ave(sst,time=jan2010,time=dec2010)''define tave = ave(sst,t=1,t=648)''define ss = pow(t1-tave,2)+pow(t2-tave,2)+pow(t3-tave,2)+pow(t4-tave,2)''define ss = pow(t5-tave,2)+pow(t6-tave,2)+pow(t7-tave,2)+pow(t8-tave,2)+ss''define ss = pow(t9-tave,2)+pow(t10-tave,2)+pow(t11-tave,2)+pow(t12-tave,2)+ss''ss = ss/12''define test = tave/sqrt(ss)*sqrt(12)''set zlog on''set gxout shaded''myp=maskout(test,test*test-1.79*1.79)'*** or 'myp=maskout(test,abs(test)-1.79)''set gxout shp''set shp -pt shppt''d myp'

梦纯

这个只是生成shp文件，'set shp -pt shppt'这句可以指定你的shp文件的输出位置。

找到之后再写个gs文件画图啊！

fh_550

梦纯 发表于 2015-4-28 16:35
这个只是生成shp文件，'set shp -pt shppt'这句可以指定你的shp文件的输出位置。

找到之后再写个gs文件 ...

啊？最后那个d myp不是画图吗？这个画图我还不会，好崩溃啊

梦纯

fh_550 发表于 2015-4-28 16:38
啊？最后那个d myp不是画图吗？这个画图我还不会，好崩溃啊

你试试不就知道了么

fh_550

梦纯 发表于 2015-4-28 17:21
你试试不就知道了么

嗯，大概知道原理了，但是上面的程序，求T的数值部分是对的吗？为什么我直接画test出来只有两条0的线啊