search

ashiklom / hector-rcmip

Reduced Complexity Model Intercomparison using the Hector model (jgcri.github.io/hector)
Other
1 stars 2 forks source link

Fix N2O emission conversion #8

Closed kdorheim closed 4 years ago

kdorheim commented 4 years ago

We were running into issues with a cold bias in the Hector temperature results. It turns out that it was related to the way that the N2O emissions were being converted.

This has PR corrects that and add some tests.

Before this PR is merged I would like MattN to do some diagnostic runs. And check out the global mean temp, total radiative forcing, and N2O emissions.

  1. Hector rcp from ini file
  2. Hector rcp from ini file + N2O emissions from this branch
mnichol3 commented 4 years ago

This plot shows the default Hector, RCMIP, & corrected RCMIP RCP60 N2O emissions on the bottom axis and the global temp (Tgav) that results from a Hector run using each version of the N2O emissions. The corrected N2O emissions match perfectly with the default Hector emissions, as do the resulting global temperatures.

corrected_rcmip_n2o_ems.pdf

The code used to create the plots can be found here

ashiklom commented 4 years ago

Wow, that is a big difference. My bad 😬 Nice job finding and fixing it, both of you!