yrange = range(ref_var, global_var)
par(mar=c(4.7,5,1,1), cex.axis=1.5, cex.lab=1.5)
plot(ref_var, type="l", col=1, lwd=3, ylim=yrange,
yrange = range(ref_var, global_var)
par(mar=c(4.7,5,1,1), cex.axis=1.5, cex.lab=1.5)
plot(ref_var, type="l", col=1, lwd=3, ylim=yrange,