# this file gives the overall sequence of R calculations for the figures used in the
# Schmidt, Mann and Rutherford comment on the McShane and Wyner (2010) paper.
# Code is adapted from McShane and Wyner (2010).
# Please report any problems or errors to Gavin Schmidt <gschmidt@giss.nasa.gov>

# Jan 2010: This is a replotting of figure 2 in the SMR paper taking account of
#   a) slightly modified (but more consistent, eiv_* reconstuctions)
#   b) variable baselines for plot (minor one year inconsistency fixed) 
#   c) Skill scores calculated using a baseline defined over the whole calibration interval
# Acknowledgements to McShane and Wyner for noticing issues (b) and (c).

### plots for fig 2: Two versions of pseudo proxy tests with two different models

# set baseline for plot
#strplotbase=1900; endplotbase=1980  
strplotbase=1856; endplotbase=1980  # align plot baseline with calibration interval

pdffile="SMR_fig2update.pdf"
pdf(pdffile, width=16, height=10)
par(mfrow=c(2,2))

panel="a) "
num=59        # 59 proxy network
src="CSM" 
source("R_pseudo_v2.txt")
panel="b) "
src="GKSS"    # choice of simulation 
source("R_pseudo_v2.txt")

num=104       # 104 proxy network
src="CSM" 
panel="c) "
source("R_pseudo_v2.txt")
src="GKSS"    # choice of simulation 
panel="d) "
source("R_pseudo_v2.txt")

dev.off()

print("Fig 2 done...")