I am trying to compute the SCC for a modified version of MimiDICE2013 but I’m not getting the right numbers.
I have added components, replaced components and optimized variables to the original MimiDICE2013 and I have used the following code to get the 5-year SCC values:
scc = DataFrame(time = Int64, scc = Float64) #empty dataframe for ii in 1:60 scc_i = MimiDICE2013.compute_scc(MyModel; year= 2005 + ii * 5, last_year = 2305, prtp = 0.03) #compute scc year_i = 2005 + ii * 5 push!(scc, [year_i, scc_i]) #fill SCC dataframe end
This is giving me some numbers but they don’t match with the results of the emissions mitigation. To give an example I’m attaching a figure that shows the emissions, temperature and SCC of the MimiDICE2013 and MyModel, lower emissions should be associated with larger SCC, right?
I know that the SCC is wrong (and not other part of the code) because I had this exact same model on matlab and I’m getting now the same results except from the SCC computation. Do you see any error or can you think about something that is causing this problem?
Thanks so much for your help in advance!