Halocarbon Concentrations & RF

[kdorheim]

While comparing Fair and Hector halocarbon RFs we've run into a few questionable elements of Hector's halocarbon implementation. The following things require further attention/investigation.

1. Units of emissMol are incorrect, it looks like we should be in units of gigga mol but that is not the case. The molar mass is being read in from the ini file in grams/mol
2. What are the units on the AtmosphereDryAirConstant? And why is it being multiplied by .1?
3. On the RHS of the update atmospheric concentration equation should the Ha be the concentration at the previous time step? Or should it be the pre-industrial concentration? This is a mismatch with the equation in the documentation manuscripts.
4. What should the units for rho be?
5. Confirm that the halocarbons with non-zero preindustrial concentrations are using the correct values? We seem to be off

---

I've attached the halocarbon concentrations from AR6 and RCMIP if helpful down the line...

RCMIP_haloconc_hector.csv

AR6_Fig2_6_data_hector.csv

[kdorheim]

3. On the RHS of the update atmospheric concentration equation should the Ha be the concentration at the previous time step? Or should it be the pre-industrial concentration? This is a mismatch with the equation in the documentation manuscripts.

Yes, the Ha should be the previous year's concentration, whereas in some of our documentation, it states that it is the concentration at 1750.

[kdorheim]

Add the halocarbon RF to the outputstream csv file