set ENERGY := solar_H2 oil coal solar bio nuc elec res ind transp ;
set ENERGY_IN := solar_H2 oil coal solar bio nuc elec ;
set ENERGY_OUT := elec res ind transp ;
set FUELS := oil coal bio nuc ;
set YEARS := 2000 2010 2020 2030 2040 2050 2060 2070 2080 2090 2100 2110 2120;
set YEARS_plus := 1990 2000 2010 2020 2030 2040 2050 2060 2070 2080 2090 2100 2110 2120;
set YEARS_first := 2000;
set YEARS_minus := 1990;

param price := 
    coal 1    oil 3    bio 3    nuc 0.5 ;

param supply_pot :=
    coal 100000    oil 10000    nuc 30000    bio Infinity ;

param emis_fact :=
    oil 20    coal 25     bio 0     nuc 0 ;

#Use one of the three following climate targets scenarios
param emis_up := 
#400ppm
# 2000 8000   2010 8200   2020 7840   2030 6536   2040 4336   2050 2906   2060 1362   2070 1020   2080 1081   2090 1044   2100  940   2110 1149   2120 1379  ;
#450ppm
 2000 8000   2010 8200   2020 8826   2030 7862   2040 6517   2050 5619   2060 4966   2070 4469   2080 4055   2090 3658   2100 3276   2110 3057   2120 2877 ;
#500ppm
# 2000 8000   2010 8200   2020 8826   2030 9840   2040 10196   2050 9855   2060 8998   2070 8173   2080 7381   2090 6094   2100 5415   2110 4896   2120 4462 ;

param life_plant := 25 ;
param Msec_per_year := 31.536 ;
param t_step := 10 ;
param bio_pot:= 200 ;
param r := 0.05 ;

param veh_cost :=
    oil 0    coal 0    bio 0    solar_H2 20    elec 40    solar 0    nuc 0 ;

param limits:   ind  res  elec transp :=
    coal        0.8  0.5  1    1
    solar       1    0.5  1    1
    bio         0.5  1    1    1
    oil         1    1    1    1.3
    solar_H2    1    1    1    1 
    nuc         1    1    0.2  1
    elec        1    1    1    0.5 ;

param dem (tr): elec  res  ind  transp :=
    2000   50    50   100    80
    2010   60    60   120   100
    2020   70    75   140   120
    2030   90    90   160   140
    2040   110  105   180   160
    2050   140  120   190   180
    2060   170  120   190   210
    2070   200  110   190   230
    2080   230  100   190   260
    2090   260   90   190   290
    2100   280   90   190   300
    2110   280   90   190   300
    2120   280   90   190   300 ;

param cost_inv:   elec   res  ind   transp :=
        coal      1200    350  350   1000
        oil        600    200  200   1000
        bio       1500    400  400   1100
        solar     2500    500    0      0
        nuc       2000      0    0      0
        solar_H2  4000   3000 3000   3000
        elec         0      0    0      0 ;

param effic: elec   res  ind   transp :=
    coal      0.45   0.8   0.8   0.5
    oil       0.5    0.8   0.8   0.9
    bio       0.4    0.8   0.8   0.5
    solar     1      1     0     0
    nuc       0.33   0     0     0
    solar_H2  1      0.8   0.8   2
    elec      0      0     0     3   ;

param lf:     elec   res  ind   transp :=
    coal      0.7    0.7   0.7   0.8
    oil       0.7    0.7   0.7   0.8
    bio       0.7    0.7   0.7   0.8
    solar     0.25   0.25  0     0
    nuc       0.7    0     0     0
    solar_H2  0.25   0.25  0.25  0.25
    elec      0      0     0     0.7  ;

param init_energy:   elec   res   ind   transp :=
            coal      50     20    70      0
            oil       50     50    70    100
            bio        0      0     0      0
            solar      0      0     0      0
            nuc       28      0     0      0
            solar_H2   0      0     0      0
            elec       0      0     0      0 ;