#'@title temp
#'@description temp
#'@return temp
#'@author Mike Vlah
#'@export
plotO2_CO2 <- function(modelResults){
par(mar=c(0,4,0.5,4))
d = modelResults$timeSeries
hour = seq(d$localHour[1], length.out=nrow(d), by=1/6)
int_hr_inds = which(hour %% 1 == 0)
# ymin = min(c(d$O2, d$CO2), na.rm=TRUE)
# ymax = max(c(d$O2, d$CO2), na.rm=TRUE)
plot(hour, d$O2, main='', xlab='',# ylim=c(ymin, ymax),
ylab='O2 conc (umol/L)', xaxt='n', type='l', col='blue')
par(new=TRUE)
plot(hour, d$pCO2, main='', xlab='',
ylab='', xaxt='n', type='l', col='salmon4', yaxt='n')
axis(4)
mtext('pCO2', 4, line=2.5)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.