Version2 #35
Open
Version2 #35
Conversation
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
The emissions component is expanded, three new functions are included to compute etree, sigma, and gsigma, therefore sigma0, gsigma0, dersig, eland0, deland, expcost2, partfract, pbacktime0 are included as input variables. Also, a connection between the neteconomy component and this one, is created through the parameters CPRICE and PBACKTIME.
The grosseconomy component now includes the computation of AL, GA and L, which required to add a0, ga0, dela, pop0, popadj, popsym, and tsteps as input variables. The respective function are added in the first if/else statement.
The neteconomy component now considers the computation of cost1, thus the variable SIGMA is now included as input, connecting the emissions component. The function to calculate is between lines 33 and 38.
The radiativeforcing component now includes the computation of forcoth, adding fex0, fex1 and tsteps as input variables. The forcoth is an asymptotic function, which was simulated with an if/else statement.
The welfare component now includes the computation of RR, adding rr0 as an input variable.
An exception is included in the TATM function, since we follow the method in Anderson et al (2012) a3 (the damage exponent) is uniformly distributed between 1.8 and 2.2,, in the damage component TATM is elevated to a3 and a negative number to the power of a fractional number is a complex number which rises an error in a Monte Carlo simulation. So we include a3 as a parameter, only for checking.
An exception is included in the TATM function, since we follow the method in Anderson et al (2012) a3 (the damage exponent) is uniformly distributed between 1.8 and 2.2,, in the damage component TATM is elevated to a3 and a negative number to the power of a fractional number is a complex number which rises an error in a Monte Carlo simulation. So we include a3 as a parameter, only for checking.
An exception is included in the TATM function, since we follow the method in Anderson et al (2012) a3 (the damage exponent) is uniformly distributed between 1.8 and 2.2,, in the damage component TATM is elevated to a3 and a negative number to the power of a fractional number is a complex number which rises an error in a Monte Carlo simulation. So we include a3 as a parameter, only for checking.
getdice2013excelparameters() is modified to include all the input variables included in the components.
A simulation that uses the method proposed in Anderson et al (2012)
The construct function is modified to include all the parameters added in the components, and the getsimv2 function included in mcsv2 is exported as well.
This function allows to run a Monte Carlo Simulation which calculates the SCC of a given year on every experiment and creates an excel file for the SCC and can create an excel file for the inputs.
An exception is included in the TATM function, since we follow the method in Anderson et al (2012) a3 (the damage exponent) is uniformly distributed between 1.8 and 2.2,, in the damage component TATM is elevated to a3 and a negative number to the power of a fractional number is a complex number which rises an error in a Monte Carlo simulation. So we include a3 as a parameter, only for checking.
This simple code, run a Monte Carlo Simulation with 1.000 runs and save the numbers generated according to the probability functions proposed in Anderson et al (2012) and the results of the Social Cost of Carbon for the model year 2030.
Including Distributions and Dates as dependencies
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Add this suggestion to a batch that can be applied as a single commit.
This suggestion is invalid because no changes were made to the code.
Suggestions cannot be applied while the pull request is closed.
Suggestions cannot be applied while viewing a subset of changes.
Only one suggestion per line can be applied in a batch.
Add this suggestion to a batch that can be applied as a single commit.
Applying suggestions on deleted lines is not supported.
You must change the existing code in this line in order to create a valid suggestion.
Outdated suggestions cannot be applied.
This suggestion has been applied or marked resolved.
Suggestions cannot be applied from pending reviews.
Suggestions cannot be applied on multi-line comments.
Suggestions cannot be applied while the pull request is queued to merge.
Suggestion cannot be applied right now. Please check back later.
In order to increase the amount of input variables, so I can follow the method for estimating the Delta indicator created by Borgonovo(2007) and implemented in Anderson et al (2012) for the DICE model I make some changes in the components: climatedynamics, emissions, grosseconomy, neteconomy, radiativeforcing and welfare. Then I modified the getdice2013excelparameters() function in parameters.jl to find the added parameters in the same excel file and modified MimiDICE2013.jl and included these parameters in the construction of the model. Also I created 2 files, mcsv2.jl which creates a simulation for the parameters following the method proposed in Anderson et al (2012) and montecarloscc.jl which includes a function that runs a Monte Carlo Simulation and perform the SCC calculation as a post trial function, allowing prtp and eta to be uncertain parameters. I find these changes beneficial, as the model now calculates some variables that were copied from the excel thus allowing a higher degree of freedom for performing GSA.