SOI index, 1876-2004. By William M. Connolley.

Red: monthly averages Blue: 12 month running mean

Data from ftp://ftp.bom.gov.au/anon/home/ncc/www/sco/soi/soiplaintext.html.

IDL script:

a=readfromfile('soi.txt')
a1=a(1:12,*)
y=a(0,*)
a2=avg(a1,0)
n=n_elements(y)
a3=smooth(reform(a1,n*12),13)

y1=a1
for i=0,n-1 do y1(i*12+indgen(12))=y(i)+indgen(12)/12.

device,set_res=500*[1.5,1]

plot,/xs,y1,a1,yr=25*[-1,1],/nodata
oplot,y1,a1,col=2
oplot,y1,a3,th=3,col=3

oplot,makerange(y1),[0,0]

已授權您依據自由軟體基金會發行的無固定段落、封面文字和封底文字GNU自由文件授權條款1.2版或任意後續版本，對本檔進行複製、傳播和／或修改。該協議的副本列在GNU自由文件授權條款中。http://www.gnu.org/copyleft/fdl.htmlGFDLGNU Free Documentation Licensetruetrue

	此檔案採用共享創意 姓名標示-相同方式分享 3.0 未在地化版本授權條款。
	您可以自由： 分享 – 複製、發佈和傳播本作品 重新修改 – 創作演繹作品 惟需遵照下列條件： 姓名標示 – 您必須指名出正確的製作者，和提供授權條款的連結，以及表示是否有對內容上做出變更。您可以用任何合理的方式來行動，但不得以任何方式表明授權條款是對您許可或是由您所使用。 相同方式分享 – 如果您利用本素材進行再混合、轉換或創作，您必須基於如同原先的相同或兼容的條款，來分布您的貢獻成品。
	已新增授權條款標題至此檔案，作為GFDL授權更新的一部份。http://creativecommons.org/licenses/by-sa/3.0/CC BY-SA 3.0Creative Commons Attribution-Share Alike 3.0truetrue