tonygardella / BUGE512

2 stars 4 forks source link

Emissions units #3

Open ashiklom opened 9 years ago

ashiklom commented 9 years ago

@tonygardella What are the units of carbon emissions in your Solow Swan model? I'm getting some weird behavior when I try to couple it with the DICE carbon cycle, and I think it may be because of a unit mismatch (DICE has emissions in GtC per year).

Also, just generally, it would help to have units in all our models to avoid confusion when we merge them.

tonygardella commented 9 years ago

They should be in Megatons, but the behavior of the series seems off to me and I haven't been able to figure out why.

tonygardella commented 9 years ago

It comes out to be gigatons of CO2, but it doesn't reflect world emissions. We would be underestimating temp perturbation in the DICE carbon model. Do you think we should use the RCP database emissions set instead?

ashiklom commented 9 years ago

I think I know the problem. As it is right now, emissions are a product of income y_net and energy intensity epsi. Both of these are large numbers (~10^6 or 10^7), so their product is something like 10^15 Megatons, which is clearly wrong. There must be a unit mismatch somewhere, or else a formula that's entered wrong.

tonygardella commented 9 years ago

I realized my old commit didn't. Check my code for emissions now.