# JonB Says:
November 25th, 2009 at 3:19 am
I’ve been trying to puzzle out for myself what this
code is actually trying to do. My IDL-fu is effectively
non-existant however…
If you expand ESRs original quote a little bit you get…
plot,timey,comptemp(*,3),/nodata,$
/xstyle,xrange=[1881,1994],xtitle='Year',$
/ystyle,yrange=[-3,3],ytitle='Normalised anomalies',$
; title='Northern Hemisphere temperatures, MXD and corrected
MXD'
title='Northern Hemisphere temperatures and MXD reconstruction'
;
yyy=reform(comptemp(*,2))
;mknormal,yyy,timey,refperiod=[1881,1940]
filter_cru,5.,/nan,tsin=yyy,tslow=tslow
oplot,timey,tslow,thick=5,color=22
yyy=reform(compmxd(*,2,1))
;mknormal,yyy,timey,refperiod=[1881,1940]
;
; Apply a VERY ARTIFICAL correction for decline!!
;
yrloc=[1400,findgen(19)*5.+1904]
valadj=[0.,0.,0.,0.,0.,-0.1,-0.25,-0.3,0.,-0.1,0.3,0.8,1.2,1.7,2.5,2.6,2.6,$
2.6,2.6,2.6]*0.75 ; fudge factor
if n_elements(yrloc) ne n_elements(valadj) then message,'Oooops!'
;
yearlyadj=interpol(valadj,yrloc,timey)
;
;filter_cru,5.,/nan,tsin=yyy+yearlyadj,tslow=tslow
;oplot,timey,tslow,thick=5,color=20
;
filter_cru,5.,/nan,tsin=yyy,tslow=tslow
oplot,timey,tslow,thick=5,color=21
I read this as being responsible for plotting the ‘Northern
Hemisphere temperatures and MXD reconstruction’.
Note however the commented out code. The way i’m
reading this, any graph titled ‘Northern Hemisphere
temperatures, MXD and corrected MXD’ with a thick
red line released before the brown matter hit the whirly
thing probably has cooked data. Likewise if you see a thick
blue line you might be ok.
Where do i get thick red line from?
According to the IDL Reference Guide for IDL v5.4, thick=5
means 5 times normal thickness. and color=21 links to (i
believe)
def_1color,20,color='red'
def_1color,21,color='blue'
def_1color,22,color='black'
from just above the code segment itself.
However there’s one thing that i’m not sure
about without either being able to play with IDL or seeing
the end graph. It does actually plot the uncooked data in
black (oplot,timey,tslow,thick=5,color=22) so whats the
point of showing a cooked data line along with the uncooked
dataline? Not to mention that the uncommented version apparantly
prints the same data series twice. |