Skip to content

Documentation of model components

The MIMOSA model consists of several sub-modules, called model components. Each component is made up of multiple variables (global or regional), parameters and equations (called constraints).

  • Economics

    Cobb-Douglas production function, investments, consumption and capital stock

    Read more

  • Welfare

    Welfare and the utility function

    Read more

  • Emissions and temperature

    • Baseline emissions, regional mitigated emissions and cumulative emissions
    • Temperature module
    • Constraints on emissions like inertia

    Read more

  • Mitigation

    Mitigation costs, Marginal Abatement Cost curve and technological learning

    Read more

  • Sea-level rise

    Determines the level of global sea-level rise

    Read more

  • Damages

    Temperature related damages and sea-level rise related damages

    Read more

  • Effort-sharing

    Optional Various effort-sharing regimes to distribute the mitigation (and sometimes damage) costs

    Read more

  • Emission trading

    Optional Allows for global trade of emission reductions

    Read more

  • Financial transfers

    Optional Financial transfer schemes like a global damage cost pool

    Read more

Building the concrete model

Builds the concrete model for MIMOSA by combining all components. Some components are optional. In the parameters, different variants of some components can be chosen. The components are:

Source code in mimosa/model_builder.py 
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
def create_model(
    inputs: ModelInputs,
) -> Tuple[ConcreteModel, List[Equation]]:
    """
    ## Building the concrete model

    Builds the concrete model for MIMOSA by combining all components. Some components are optional. In the
    parameters, different variants of some components can be chosen. The components are:

    - [`damage module`](#param__model%20structure.damage%20module): The damage module to use
    - [`emissiontrade module`](#param__model%20structure.emissiontrade%20module): The emission trading module to use
    - [`financialtransfer module`](#param__model%20structure.financialtransfer%20module): The financial transfer module to use
    - [`welfare module`](#param__model%20structure.welfare%20module): The welfare module to use
    - [`objective module`](#param__model%20structure.objective%20module): The objective module to use
    - [`effortsharing module`](#param__model%20structure.effortsharing%20module): The effort-sharing module to use

    """
    m = create_base_model(inputs)

    # Each constraint will be put in this list,
    # then added to the model at the end of this file.
    constraints = []

    # Emissions, temperature and sea-level rise
    constraints.extend(emissions.get_constraints(m, inputs))
    constraints.extend(sealevelrise.get_constraints(m, inputs))

    # Damage and mitigation costs
    constraints.extend(damages.get_constraints(m, inputs))
    constraints.extend(mitigation.get_constraints(m, inputs))

    # Emission trading, financial transfers and effort sharing
    constraints.extend(emissiontrade.get_constraints(m, inputs))
    constraints.extend(financialtransfer.get_constraints(m, inputs))
    constraints.extend(effortsharing.get_constraints(m, inputs))

    # Production, consumption and welfare
    constraints.extend(cobbdouglas.get_constraints(m, inputs))
    constraints.extend(welfare.get_constraints(m, inputs))

    # Objective of optimisation
    model_objective, objective_constraints = objective.get_constraints(m, inputs)
    constraints.extend(objective_constraints)

    ######################
    # Add constraints to concrete model
    ######################

    for constraint in constraints:
        add_constraint(m, constraint.to_pyomo_constraint(m), constraint.name)

    ######################
    # Get all equations (not constraints) for simulation mode
    ######################
    equations = [eq for eq in constraints if isinstance(eq, Equation)]

    m.objective = model_objective

    return m, equations