XXX="WCT"
Directory_Output<-"/Volumes/shared_data/paulb/Results/Results/Skills/" #Directory where you want your output to be stocked. 

AnnEns<-read.csv(paste(Directory_Output,XXX,"_","ensemble","_mean_Anomalies_Continuous_Skill_Scores_3_Month_Predictions.csv",sep=""),sep=",")
MSE_AnnEns<-AnnEns[seq(4,nrow(AnnEns),by=7),4]
AnnSing<-read.csv(paste(Directory_Output,XXX,"_","single","_mean_Anomalies_Continuous_Skill_Scores_3_Month_Predictions.csv",sep=""),sep=",")
MSE_AnnSing<-AnnSing[seq(4,nrow(AnnSing),by=7),4]

plot.new()
ymax<-max(cbind(MSE_AnnSing,MSE_AnnEns))
plot.window(c(1,12), c(0,ymax))
axis(1, at=seq(1,12,by=1))
axis(2, at=seq(1,ymax,by=floor(ymax/10)))
box()
lines(MSE_AnnEns, type="o", col=2, pch=1)
lines(MSE_AnnSing, type="o", col=1, pch=1)
title(main=paste("MSE for seasonial lead time and three month anomalies predictions by the ",XXX," model"), col.main="black", font.main=1)
title(xlab="Lead Time, from Mars to February")
title(ylab=paste("MSE")) 
legend("topright", c("Single output","Averaged ensemble output"), cex=0.8, col=c(1,2), pch=c(1,1), horiz=TRUE);