; ; Material for volcanic cooling exercise for Reepham High School tree ring activity ; ; Get regional data (incl. ALL) plus instrumental ; restore,filename='regtemp_calibrated.allversion.idlsave' ; Gets: nyr, nreg, calregts,regname,timey,tempregts,tempnyr,temptimey nyr1=nyr ; ; Get combined HugHF AgeLF for the ALL region (overwrites nyr) ; restore,filename='combtempALL_calibrated.idlsave' ; Gets: nyr, yrmxd, prednh ; ; Adjust ALL trees to have mean of 100 and sd of 15 ; mknormal,prednh prednh=(prednh*15)+100. ; ireg=9 print,'region ',regname[ireg] ; per=[1881,1960] kmxd=where((yrmxd ge per[0]) and (yrmxd le per[1]),nkmxd) kobs=where((temptimey ge per[0]) and (temptimey le per[1]),nkobs) if nkmxd ne nkobs then message,'Ooops!' ; ; Prepare for plotting ; loadct,39 multi_plot,nrow=1,layout='large',/landscape if !d.name eq 'X' then begin win !p.font=-1 endif else begin !p.font=0 device,/helvetica,/bold,font_size=18 endelse makecolors,['magenta','red','mred','brown','orange','mdyellow','green',$ 'vlblue','mlblue','deepblue','lpurple','palepurple','black','mgrey'],indgen(14)+20 ; acoeff=linfit(prednh[kmxd],tempregts[kobs,ireg]) print,correlate(prednh[kmxd],tempregts[kobs,ireg]) ; xr=[20,140] yr=acoeff[0]+acoeff[1]*xr print,xr print,yr plot,prednh[kmxd],tempregts[kobs,ireg],psym=1,$ xtitle='Tree-ring measurement',xstyle=1,xrange=xr,$ ytitle='Summer temperature deviation (!Uo!NC)',ystyle=1,yrange=[-3,1] ; oplot,xr,yr,noclip=0,thick=5 ; end