R/plotO2_CO2.R

#'@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)
}
gholtgrieve/gassyPants documentation built on May 9, 2019, 5:02 a.m.