-
Notifications
You must be signed in to change notification settings - Fork 4
/
mosaic_health_exe.gms
179 lines (120 loc) · 7.11 KB
/
mosaic_health_exe.gms
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
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
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
*health
urbanization(t,c) = ((global_param_health("urbanization_param_1")*sqrt(y_pc(t,c)))
+ (global_param_health("urbanization_param_2")*sqrt(P_dens(t,c))))/
(1 + global_param_health("urbanization_param_1")*sqrt(y_pc(t,c))
+ (global_param_health("urbanization_param_2")*sqrt(P_dens(t,c))));
*MORBIDITY AND MORTALITY
value_life(t,c) = (value_life_baseline_USA_2010 / Y_gross("2010","USA")) * Y_net(t,c);
*DIARRHEA
change_diarrhea_rate(t,c) = (country_rate(c,"diarrhea") / 100) *
(y_pc(t,c) / y_pc("2010",c))**
global_param_health("income_elas_diarrhea_mort") *
(temp_r(t,c) / (tempC_preind*Country_Tol_regional_temp(c)))**
global_param_health("diarrhea_mort_response_warming");
deaths_diarrhea(t,c) = change_diarrhea_rate(t,c) * pop(t,c) / 1000;
damages_diarrhea(t,c) = deaths_diarrhea(t,c) * value_life(t,c);
*VECTOR BORNE
change_malaria_rate(t,c) = (country_rate(c,"malaria") * 10)
* global_param_health("malaria_param") *
temp_change_national_from_2010(t,c)**
global_param_health("deg_nonlinearity_mortality_in_warming") *
(y_pc(t,c) / y_pc("2010",c))**
global_param_health("income_elas_vector_mort");
deaths_malaria(t,c) = change_malaria_rate(t,c) * pop(t,c) / 1000000;
damages_malaria(t,c) = deaths_malaria(t,c) * value_life(t,c);
change_dengue_rate(t,c) = (country_rate(c,"dengue") * 10)
* global_param_health("dengue_param") *
temp_change_national_from_2010(t,c)**
global_param_health("deg_nonlinearity_mortality_in_warming") *
(y_pc(t,c) / y_pc("2010",c))**
global_param_health("income_elas_vector_mort");
deaths_dengue(t,c) = change_dengue_rate(t,c) * pop(t,c) / 1000000;
damages_dengue(t,c) = deaths_dengue(t,c) * value_life(t,c);
change_schistosomiasis_rate(t,c) = (country_rate(c,"schistosomiasis") * 10)
* global_param_health("schistosomiasis_param") *
temp_change_national_from_2010(t,c)**
global_param_health("deg_nonlinearity_mortality_in_warming") *
(y_pc(t,c) / y_pc("2010",c))**
global_param_health("income_elas_vector_mort");
deaths_schistosomiasis(t,c) = change_schistosomiasis_rate(t,c) * pop(t,c) / 1000000;
damages_schistosomiasis(t,c) = deaths_schistosomiasis(t,c) * value_life(t,c);
change_vector_borne_rate(t,c) = change_malaria_rate(t,c) + change_dengue_rate(t,c)
+ change_schistosomiasis_rate(t,c);
deaths_vector_borne(t,c) = deaths_malaria(t,c) + deaths_dengue(t,c) + deaths_schistosomiasis(t,c);
damages_vector_borne(t,c) = damages_malaria(t,c) + damages_dengue(t,c) + damages_schistosomiasis(t,c);
*CARDIOVASCULAR
change_cardiovascular_heat_over65_rate(t,c) = Tol_heat_over65_1(c) *
temp_change_national_from_2010(t,c) + Tol_heat_over65_2(c) *
temp_change_national_from_2010(t,c)** 2 ;
deaths_cardiovascular_heat_over65(t,c) = (change_cardiovascular_heat_over65_rate(t,c) *
(pop(t,c) * (1-pop_under65(t,c))) / 100000) * urbanization(t,c) ;
damages_cardiovascular_heat_over65(t,c) = deaths_cardiovascular_heat_over65(t,c) * value_life(t,c);
change_cardiovascular_heat_under65_rate(t,c) = Tol_heat_under65_1(c) *
temp_change_national_from_2010(t,c) + Tol_heat_under65_2(c) *
temp_change_national_from_2010(t,c)** 2 ;
deaths_cardiovascular_heat_under65(t,c) = (change_cardiovascular_heat_under65_rate(t,c) *
(pop(t,c) * pop_under65(t,c)) / 100000) * urbanization(t,c) ;
damages_cardiovascular_heat_under65(t,c) = deaths_cardiovascular_heat_under65(t,c) * value_life(t,c);
change_cardiovascular_heat_rate(t,c) = damages_cardiovascular_heat_under65(t,c) + damages_cardiovascular_heat_over65(t,c);
deaths_cardiovascular_heat(t,c) = deaths_cardiovascular_heat_over65(t,c) + deaths_cardiovascular_heat_under65(t,c) ;
damages_cardiovascular_heat(t,c) = deaths_cardiovascular_heat(t,c) * value_life(t,c);
change_cardiovascular_cold_over65_rate(t,c) = Tol_cold_over65_1(c) *
temp_change_national_from_2010(t,c) + Tol_cold_over65_2(c) *
temp_change_national_from_2010(t,c)**2 ;
deaths_cardiovascular_cold_over65(t,c) = change_cardiovascular_cold_over65_rate(t,c) *
(pop(t,c) * (1 - pop_under65(t,c))) / 100000 ;
damages_cardiovascular_cold_over65(t,c) = deaths_cardiovascular_cold_over65(t,c) * value_life(t,c);
change_cardiovascular_cold_under65_rate(t,c) = Tol_cold_under65_1(c) *
temp_change_national_from_2010(t,c) + Tol_cold_under65_2(c) *
temp_change_national_from_2010(t,c)**2 ;
deaths_cardiovascular_cold_under65(t,c) = change_cardiovascular_cold_under65_rate(t,c) *
(pop(t,c) * pop_under65(t,c)) / 100000 ;
damages_cardiovascular_cold_under65(t,c) = deaths_cardiovascular_cold_under65(t,c) * value_life(t,c);
change_cardiovascular_cold_rate(t,c) = damages_cardiovascular_cold_under65(t,c) + damages_cardiovascular_cold_over65(t,c);
deaths_cardiovascular_cold(t,c) = deaths_cardiovascular_cold_over65(t,c) + deaths_cardiovascular_cold_under65(t,c);
damages_cardiovascular_cold(t,c) = deaths_cardiovascular_cold(t,c) * value_life(t,c);
change_cardiovascular_rate(t,c) = change_cardiovascular_heat_rate(t,c) + change_cardiovascular_cold_rate(t,c);
deaths_cardiovascular(t,c) = deaths_cardiovascular_heat(t,c) + deaths_cardiovascular_cold(t,c);
damages_cardiovascular(t,c) = deaths_cardiovascular(t,c) * value_life(t,c);
*RESPIRATORY
change_respiratory_rate(t,c)= Tol_respiratory_1(c) *
temp_change_national_from_2010(t,c) + Tol_respiratory_2(c) *
temp_change_national_from_2010(t,c)** 2 ;
deaths_respiratory(t,c) = change_respiratory_rate(t,c) * pop(t,c) / 100000;
damages_respiratory(t,c) = deaths_respiratory(t,c) * value_life(t,c);
*TOTAL DAMAGES LOST DUE TO CLIMATE CHANGE
total_deaths(t,c) = deaths_diarrhea(t,c) + deaths_vector_borne(t,c) +
deaths_cardiovascular(t,c) + deaths_respiratory(t,c);
total_deaths_percentage(t,c) = total_deaths(t,c) / pop(t,c) * 100;
total_health_damages(t,c) = total_deaths(t,c) * value_life(t,c);
total_health_damages_percentage(t,c) = total_health_damages(t,c) / Y_gross(t,c) * 100;
*DISPLAYS GROUPED BY SECTION
$ontext
*GENERAL
display urbanization;
display pop_under65;
display value_life;
display total_deaths;
display total_deaths_pop_fraction;
display total_health_damages;
display total_health_damages_percentage;
*DIARRHEA
display deaths_diarrhea, change_diarrhea_rate;
*VECTOR BORNE
display deaths_vector_borne;
display deaths_malaria, change_malaria_rate;
display deaths_dengue, change_dengue_rate;
display deaths_schistosomiasis, change_schistosomiasis_rate;
*CARDIOVASCULAR
display deaths_cardiovascular;
display deaths_cardiovascular_heat;
display deaths_cardiovascular_heat_over65, change_cardiovascular_hot_over65_rate;
display deaths_cardiovascular_heat_under65, change_cardiovascular_hot_under65_rate;
display deaths_cardiovascular_cold;
display deaths_cardiovascular_cold_over65, change_cardiovascular_cold_over65_rate;
display deaths_cardiovascular_cold_under65, change_cardiovascular_cold_under65_rate;
display deaths_cardiovascular_heat;
display deaths_cardiovascular_cold;
*RESPIRATORY
display deaths_respiratory, change_respiratory_rate;
$offtext