Revisions for submission to *Environmental Research Letters* journal (#…

…73)

* Increase cost exponent to 3

Increased to cost exponent from 2 to 3, which brings MARGO's cost
functions more in line with other models.

Also renamed the units file.

* Improved adaptation model and costs

**Adaptation model**
Following Reviewer #1's suggestion, we improve the adaptation model somewhat following DeBruin et al's AD-DICE implementation. In particular, we assume that adaptation reduces damages proportionally, such that the damages used to compute benefits are the *residual damages*. The marginal costs of these adaptive measures increase super-linearly proportionally with the GWP. Based on suggestions from a different anonymous reviewer (from a previous submission at a different journal), we also make the cost proportional to the controlled gross damages. We tune the proportionality constant to match DeBruin's estimate for the damages that correspond to a warming of 3ºC (this reference temperature comes in as an additional input parameter in the economic submodule, Tb=3.0).

These changes are implemented in both the Diagnostic and Optimization modules.

We remove the constraint that forces adaptation to be constant in time, also following DeBruin's implementation.

**Misc**
Also fixed some bugs with the expenditure and budget objective functions (although these are not used in the paper at all).

Also improved the interface for the standard plots of the model state.

* Increment version number to reflect major (breaking) changes

* Updates tests with new results

* Updated default parameter values based on referrees

- Updated total remove costs in line with new cost function.

- Increased mitigation costs to match other estimates in 2100, rather than in 2020

- Increased climate damages according to Frankhauser review

- Increased geoengineering costs following costs of climate damages

* Fixed bugs in economic submodel and optimization

- Forgot to add Tb as paramter in some sub-module constructor methods

- Forgot to index background damages Db[i] in some places of optimization

* Increased default adaptation limit

* Bugfixes

* Back to de Bruin adaptation module

* Slightly thinner lines

* Fixed one-two dim example

* Removed legacy Tb changes

* Fixed default configuration

* Minor update

* Updated hard-coded correct results for tests

Project.toml

@@ -1,7 +1,8 @@
 name = "ClimateMARGO"
 uuid = "d3f62095-a717-45bf-aadc-ac9dfc258fa6"
 repo = "https://github.com/ClimateMARGO/ClimateMARGO.jl"
-version = "0.2.2"
+version = "0.3.0"
+
 
 [deps]
 CompilerSupportLibraries_jll = "e66e0078-7015-5450-92f7-15fbd957f2ae"
@@ -13,6 +14,7 @@ Ipopt_jll = "9cc047cb-c261-5740-88fc-0cf96f7bdcc7"
 JSON2 = "2535ab7d-5cd8-5a07-80ac-9b1792aadce3"
 JuMP = "4076af6c-e467-56ae-b986-b466b2749572"
 PyPlot = "d330b81b-6aea-500a-939a-2ce795aea3ee"
+Revise = "295af30f-e4ad-537b-8983-00126c2a3abe"
 
 [compat]
 DiffResults = "= 1.0.2"

configurations/default_params.json

@@ -1 +1,112 @@
-{"name":"default","domain":{"dt":5.0,"present_year":2020.0,"initial_year":2020.0,"final_year":2200.0},"economics":{"E0":100.0,"γ":0.02,"β":0.0022222222222222222,"ρ":0.02,"Finf":8.5,"mitigate_cost":0.034,"remove_cost":13.0,"geoeng_cost":0.126,"adapt_cost":0.18,"mitigate_init":0.0,"remove_init":0.0,"geoeng_init":0.0,"adapt_init":null,"baseline_emissions":[7.5,8.4375,9.375,10.3125,11.25,12.1875,13.125,14.0625,15.0,15.9375,16.875,17.8125,18.75,19.6875,20.625,21.5625,22.5,20.25,18.0,15.75,13.5,11.25,9.0,6.75,4.5,2.25,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0],"epsilon_cost":0.0,"extra_CO₂":[0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0,0.0]},"physics":{"c0":460.0,"T0":1.1,"a":4.977297891066924,"B":1.13,"Cd":106.0,"κ":0.73,"r":0.5}}
+{
+    "name": "default",
+    "domain": {
+        "dt": 5.0,
+        "present_year": 2020.0,
+        "initial_year": 2020.0,
+        "final_year": 2200.0
+    },
+    "economics": {
+        "E0": 100.0,
+        "γ": 0.02,
+        "β": 0.01,
+        "ρ": 0.02,
+        "Finf": 8.5,
+        "mitigate_cost": 0.055489419901509705,
+        "remove_cost": 36.97631285953187,
+        "geoeng_cost": 0.5658234787375677,
+        "adapt_cost": 0.115,
+        "mitigate_init": 0.0,
+        "remove_init": 0.0,
+        "geoeng_init": 0.0,
+        "adapt_init": 0.0,
+        "baseline_emissions": [
+            7.5,
+            8.4375,
+            9.375,
+            10.3125,
+            11.25,
+            12.1875,
+            13.125,
+            14.0625,
+            15.0,
+            15.9375,
+            16.875,
+            17.8125,
+            18.75,
+            19.6875,
+            20.625,
+            21.5625,
+            22.5,
+            20.25,
+            18.0,
+            15.75,
+            13.5,
+            11.25,
+            9.0,
+            6.75,
+            4.5,
+            2.25,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0
+        ],
+        "epsilon_cost": 0.0,
+        "extra_CO₂": [
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0,
+            0.0
+        ]
+    },
+    "physics": {
+        "c0": 460.0,
+        "T0": 1.1,
+        "a": 4.977297891066924,
+        "B": 1.13,
+        "Cd": 106.0,
+        "κ": 0.73,
+        "r": 0.5
+    }
+}