So @glarange I think to summarize there are two good options for now.
Before you employ them wait an hour or two for me to tag things and then call
pkg> up
in your environment, after which
pkg> st
should show that you have Mimi v1.2.1 and MimiIWG v 1.0.2
At this point your issues from your previous post in the forum here will be resolved, and you’ll have the ability to do one of the two following things:
- To use the deterministic version of the models, you can use the
get_marginaldamages
function to get the undercounted marginal damages vector and do with them what you want afterwards, or wrap this in a function etc.
using MimiIWG
damages = get_marginaldamages(MODEL_NAME, SCENARIO_CHOICE, gas=:CO2, year=2020, discount=0)
- To use the Monte Carlo version of the models, you can set the
save_md
flag toTRUE
in therun_scc_mcs
function with full signature
run_scc_mcs(model::model_choice;
gas::Union{Symbol, Nothing} = nothing,
trials::Int = 10000,
perturbation_years::Vector{Int} = _default_perturbation_years,
discount_rates::Vector{Float64} = _default_discount_rates,
domestic::Bool = false,
output_dir::String = nothing,
save_trials::Bool = false,
save_md::Bool = false,
tables::Bool = true)
for which the function signature states "If save_md
equals true
, then global undiscounted marginal damages from each run of the simulation will be saved in a subdirectory “output/marginal_damages”. You can then use those to do whatever you want for discounting.
- Alternatively, if you are interested in a Ramsey-type framework, in the
ramsey
branch (link) one can set the near-termdiscount_rates
, then choose parametersprtp
andeta
to discount damages according to the growth rateg
from the EMF scenarios. For now this is an experimental branch so probably less dependable, but sounds like some people are using it!
@parthum look right?