Data for: An Integrated Biophysical and Economic Modeling Framework for Long-Term Sustainability Analysis: the HARMONEY Model

These files are the two simulation codes for the Human And Resources with MONEY (HARMONEY) model. These codes, programmed in the R language, were used to produce the simulation results in the main paper and supplemental information of "An Integrated Biophysical and Economic Modeling Framework for Long-Term Sustainability Analysis: the HARMONEY Model". The codes also produce several plots in the R plot window.

HARMONEY_SetPriceByConstantMarkup: This simulates the results that assume a constant price markup, as defined in the paper.
HARMONEY_SetPriceByVariableMarkup: This simulates the results that assume a "variable" price markup, as defined in the paper.

When running the code, it might not simulate to the end time for all combinations of parameter choices. This occurs for parameters that lead to system "collapse" when the net power flow from the extraction sector reaches zero. It is best to change the parameter values slightly from those listed as being used in the simulations of the paper. If you have trouble seeing results, reduce the simulation time by reducing the "to" value in this line of code:
times <- seq(from = 0, to = 160, by = 0.1);

The investment parameters that were changed to explore the "high" and "low" investment cases are the following (for "kappa 0" and "kappa 1"):
assign("inv_kappa0.1",1.0)
assign("inv_kappa1.1",1.5)

Other discussion and simulations in the paper were informed by varying the following parameters, as listed in the code:
assign("lambda_h.max1",1.0e2) # Nature carrying capacity before "tcritical.lambda_h"
assign("delta_y.1",1.2e-02) # Depletion (production) factor before "tcritical.lambda_h" (1/[Capital * yr])
assign("delta1",0.03) # Depreciation rate of capital (%/yr)
assign("v1",1.5) # "Goods capital / gross physical output" ratio
assign("a_gg",0.1) # I-O coefficient: goods consumption of goods sector per unit of goods sector output
assign("a_ge",0.2) # I-O coefficient: goods consumption of extraction sector per unit of extraction sector output
assign("eta.Kextract_max",0.16) # Resources required to operate K.extract (and the maximum value if modeling a decline in K.extract after "tcritical.Inv_Cum_extract")
assign("eta.Kgoods_max",eta.Kextract_max) # Resources required to operate K.goods (and the maximum value if modeling a decline in K.goods after "tcritical.Inv_Cum_goods")
assign("nature_per_unit_physical_good",0.1) # The quantity of nature needed to become a unit of capital
assign("markup.goods",0.07) # fraction to mark up costs to define price, goods sector
assign("markup.extract",markup.goods) # fraction to mark up costs to define price, extraction sector
assign("S",8e-2); # Subsistence salary per capita (consumption per person when not in famine in nature/(person * yr))
assign("rho.extract",0.010); # Threshold nature consumption per capita (consumption per person below which there is "famine": nature/(person * yr))