CO2 and stuff notes

[bpbond]

I am re-reading the Abramoff 2018 and 2022 papers...

• [ ] In the 2022 paper, they did a type of parameter sensitivity exercise! Huh. See Figure S2. But your analysis goes much further.
• [ ] CO2: when I run the Millenial example RMarkdown (so not Leeya's work) the CO2 outputs looks like what we'd expect:

  So why is it producing the funky behavior we're seeing?!? Hmm.

[bpbond]

AH HA! @stephpenn1

From the model_tutorial.Rmd file:

#Define function which plots the modeled CO2 emissions over time
##CO2 emissions are cumulative, so we have to take the difference between timesteps to plot the CO2 emission rate
plot.co2 <- function(output){
  plot(output[1:(dim(output)[1]-1),"time"], 
       diff(output[,"CO2"]), 
       xlab="time (days)", 
       ylab="co2 (gC m-2 d-1)", 
       col=1, 
       type="l")
}

Who the heck decided that CO2 emissions should be cumulative, and why is this not documented in the papers?!? 😐

Anyway, now we know. You can use the diff function (see above) to calculate the differences between timesteps, and that will be the output we care about.

[stephpenn1]

Nice investigative work! I will use that diff function for our output!