* put file for Optimal results

File all_d07e;
all_d07e.pc=5;
all_d07e.pw=250;
Put all_d07e;
Put / "optimal run";
Put / "year";
Loop (tearly, put year(tearly)::0);
Put / "output";
Loop (tearly, put opt_y(tearly)::3);
Put / "pccon";
Loop (tearly, put opt_cpc(tearly)::3);
Put / "savrate";
Loop (tearly, put opt_s(tearly)::4);
Put / "indem";
Loop (tearly, put opt_indem(tearly)::4);
Put / "sigma";
Loop (tearly, put opt_sigma(tearly)::4);
Put / "temp";
Loop (tearly, put opt_tatm(tearly)::3);
Put / "conc";
Loop (tearly, put opt_mat(tearly)::3);
Put / "soc cost carbon";
Loop (tearly, put opt_tax(tearly)::2);
Put / "intrate";
Loop (tearly, put opt_ri(tearly)::3);
Put / "discrate";
Loop (tearly, put opt_rr(tearly)::5);
Put / "prod";
Loop (tearly, put opt_al(tearly)::5);
Put / "exogforc";
Loop (tearly, put opt_forcoth(tearly)::3);
Put / "pop";
Loop (tearly, put opt_l(tearly)::3);
Put / "carbon tax";
Loop (tearly, put opt_mcemis(tearly)::4);
Put / "margy";
Loop (tearly, put opt_yy(tearly)::3);
Put / "margc";
Loop (tearly, put opt_cc(tearly)::5);
Put / "miu";
Loop (tearly, put opt_miu(tearly)::3);
Put / "total emissions";
Loop (tearly, put opt_wem(tearly)::3);
Put / "interest rate";
Loop (tearly, put opt_ri(tearly)::4);
Put / "damages";
Loop (tearly, put opt_dam(tearly)::3);
Put / "abatement cost";
Loop (tearly, put opt_abate(tearly)::2);
Put /"objective function";
Put opt_utility::3;

* put file for 15 period control

Put / "base run with no controls for 15 periods";
Put / "year";
Loop (tearly, put year(tearly)::0);
Put / "output";
Loop (tearly, put base25_y(tearly)::3);
Put / "pccon";
Loop (tearly, put base25_cpc(tearly)::3);
Put / "savrate";
Loop (tearly, put base25_s(tearly)::4);
Put / "indem";
Loop (tearly, put base25_indem(tearly)::4);
Put / "sigma";
Loop (tearly, put base25_sigma(tearly)::4);
Put / "temp";
Loop (tearly, put base25_tatm(tearly)::3);
Put / "conc";
Loop (tearly, put base25_mat(tearly)::3);
Put / "soc cost carbon";
Loop (tearly, put base25_tax(tearly)::2);
Put / "intrate";
Loop (tearly, put base25_ri(tearly)::3);
Put / "discrate";
Loop (tearly, put base25_rr(tearly)::5);
Put / "prod";
Loop (tearly, put base25_al(tearly)::5);
Put / "exogforc";
Loop (tearly, put base25_forcoth(tearly)::3);
Put / "pop";
Loop (tearly, put base25_l(tearly)::3);
Put / "carbon tax";
Loop (tearly, put base25_mcemis(tearly)::4);
Put / "margy";
Loop (tearly, put base25_yy(tearly)::3);
Put / "margc";
Loop (tearly, put base25_cc(tearly)::3);
Put / "miu";
Loop (tearly, put base25_miu(tearly)::3);
Put / "total emissions";
Loop (tearly, put base25_wem(tearly)::3);
Put / "interest rate";
Loop (tearly, put base25_ri(tearly)::4);
Put / "damages";
Loop (tearly, put base25_dam(tearly)::2);
Put / "abatement cost";
Loop (tearly, put base25_abate(tearly)::2);
Put /"objective function";
Put base25_utility::3;

*put file for 5 period control


Put / "base run with no controls for 5 periods";
Put / "year";
Loop (tearly, put year(tearly)::0);
Put / "output";
Loop (tearly, put base5_y(tearly)::3);
Put / "pccon";
Loop (tearly, put base5_cpc(tearly)::3);
Put / "savrate";
Loop (tearly, put base5_s(tearly)::4);
Put / "indem";
Loop (tearly, put base5_indem(tearly)::4);
Put / "sigma";
Loop (tearly, put base5_sigma(tearly)::4);
Put / "temp";
Loop (tearly, put base5_tatm(tearly)::3);
Put / "conc";
Loop (tearly, put base5_mat(tearly)::3);
Put / "soc cost carbon";
Loop (tearly, put base5_tax(tearly)::2);
Put / "intrate";
Loop (tearly, put base5_ri(tearly)::3);
Put / "discrate";
Loop (tearly, put base5_rr(tearly)::5);
Put / "prod";
Loop (tearly, put base5_al(tearly)::5);
Put / "exogforc";
Loop (tearly, put base5_forcoth(tearly)::3);
Put / "pop";
Loop (tearly, put base5_l(tearly)::3);
Put / "carbon tax";
Loop (tearly, put base5_mcemis(tearly)::4);
Put / "margy";
Loop (tearly, put base5_yy(tearly)::3);
Put / "margc";
Loop (tearly, put base5_cc(tearly)::3);
Put / "miu";
Loop (tearly, put base5_miu(tearly)::3);
Put / "total emissions";
Loop (tearly, put base5_wem(tearly)::3);
Put / "interest rate";
Loop (tearly, put base5_ri(tearly)::4);
Put / "damages";
Loop (tearly, put base5_dam(tearly)::2);
Put / "abatement cost";
Loop (tearly, put base5_abate(tearly)::2);
Put /"objective function";
Put base5_utility::3;

*put file for limit to 2 degrees

Put / "limit T increase to 2 degrees";
Put / "year";
Loop (tearly, put year(tearly)::0);
Put / "output";
Loop (tearly, put lim2deg_y(tearly)::3);
Put / "pccon";
Loop (tearly, put lim2deg_cpc(tearly)::3);
Put / "savrate";
Loop (tearly, put lim2deg_s(tearly)::4);
Put / "indem";
Loop (tearly, put lim2deg_indem(tearly)::4);
Put / "sigma";
Loop (tearly, put lim2deg_sigma(tearly)::4);
Put / "temp";
Loop (tearly, put lim2deg_tatm(tearly)::3);
Put / "conc";
Loop (tearly, put lim2deg_mat(tearly)::3);
Put / "soc cost carbon";
Loop (tearly, put lim2deg_tax(tearly)::2);
Put / "intrate";
Loop (tearly, put lim2deg_ri(tearly)::3);
Put / "discrate";
Loop (tearly, put lim2deg_rr(tearly)::5);
Put / "prod";
Loop (tearly, put lim2deg_al(tearly)::5);
Put / "exogforc";
Loop (tearly, put lim2deg_forcoth(tearly)::3);
Put / "pop";
Loop (tearly, put lim2deg_l(tearly)::3);
Put / "carbon tax";
Loop (tearly, put lim2deg_mcemis(tearly)::4);
Put / "margy";
Loop (tearly, put lim2deg_yy(tearly)::3);
Put / "margc";
Loop (tearly, put lim2deg_cc(tearly)::3);
Put / "miu";
Loop (tearly, put lim2deg_miu(tearly)::3);
Put / "total emissions";
Loop (tearly, put lim2deg_wem(tearly)::3);
Put / "interest rate";
Loop (tearly, put lim2deg_ri(tearly)::4);
Put / "damages";
Loop (tearly, put lim2deg_dam(tearly)::2);
Put / "abatement cost";
Loop (tearly, put lim2deg_abate(tearly)::2);
Put /"objective function";
Put lim2deg_utility::3;


*put file for limit to 2.5 degrees

Put / "limit T increase to 2.5 degrees";
Put / "year";
Loop (tearly, put year(tearly)::0);
Put / "output";
Loop (tearly, put lim25deg_y(tearly)::3);
Put / "pccon";
Loop (tearly, put lim25deg_cpc(tearly)::3);
Put / "savrate";
Loop (tearly, put lim25deg_s(tearly)::4);
Put / "indem";
Loop (tearly, put lim25deg_indem(tearly)::4);
Put / "sigma";
Loop (tearly, put lim25deg_sigma(tearly)::4);
Put / "temp";
Loop (tearly, put lim25deg_tatm(tearly)::3);
Put / "conc";
Loop (tearly, put lim25deg_mat(tearly)::3);
Put / "soc cost carbon";
Loop (tearly, put lim25deg_tax(tearly)::2);
Put / "intrate";
Loop (tearly, put lim25deg_ri(tearly)::3);
Put / "discrate";
Loop (tearly, put lim25deg_rr(tearly)::5);
Put / "prod";
Loop (tearly, put lim25deg_al(tearly)::5);
Put / "exogforc";
Loop (tearly, put lim25deg_forcoth(tearly)::3);
Put / "pop";
Loop (tearly, put lim25deg_l(tearly)::3);
Put / "carbon tax";
Loop (tearly, put lim25deg_mcemis(tearly)::4);
Put / "margy";
Loop (tearly, put lim25deg_yy(tearly)::3);
Put / "margc";
Loop (tearly, put lim25deg_cc(tearly)::3);
Put / "miu";
Loop (tearly, put lim25deg_miu(tearly)::3);
Put / "total emissions";
Loop (tearly, put lim25deg_wem(tearly)::3);
Put / "interest rate";
Loop (tearly, put lim25deg_ri(tearly)::4);
Put / "damages";
Loop (tearly, put lim25deg_dam(tearly)::2);
Put / "abatement cost";
Loop (tearly, put lim25deg_abate(tearly)::2);
Put /"objective function";
Put lim25deg_utility::3;

*put file for limit to 2x preindustrial co2

Put / "limit co2 to 2x preindustrial";
Put / "year";
Loop (tearly, put year(tearly)::0);
Put / "output";
Loop (tearly, put lim2xco2_y(tearly)::3);
Put / "pccon";
Loop (tearly, put lim2xco2_cpc(tearly)::3);
Put / "savrate";
Loop (tearly, put lim2xco2_s(tearly)::4);
Put / "indem";
Loop (tearly, put lim2xco2_indem(tearly)::4);
Put / "sigma";
Loop (tearly, put lim2xco2_sigma(tearly)::4);
Put / "temp";
Loop (tearly, put lim2xco2_tatm(tearly)::3);
Put / "conc";
Loop (tearly, put lim2xco2_mat(tearly)::3);
Put / "soc cost carbon";
Loop (tearly, put lim2xco2_tax(tearly)::2);
Put / "intrate";
Loop (tearly, put lim2xco2_ri(tearly)::3);
Put / "discrate";
Loop (tearly, put lim2xco2_rr(tearly)::5);
Put / "prod";
Loop (tearly, put lim2xco2_al(tearly)::5);
Put / "exogforc";
Loop (tearly, put lim2xco2_forcoth(tearly)::3);
Put / "pop";
Loop (tearly, put lim2xco2_l(tearly)::3);
Put / "carbon tax";
Loop (tearly, put lim2xco2_mcemis(tearly)::4);
Put / "margy";
Loop (tearly, put lim2xco2_yy(tearly)::3);
Put / "margc";
Loop (tearly, put lim2xco2_cc(tearly)::3);
Put / "miu";
Loop (tearly, put lim2xco2_miu(tearly)::3);
Put / "total emissions";
Loop (tearly, put lim2xco2_wem(tearly)::3);
Put / "interest rate";
Loop (tearly, put lim2xco2_ri(tearly)::4);
Put / "damages";
Loop (tearly, put lim2xco2_dam(tearly)::2);
Put / "abatement cost";
Loop (tearly, put lim2xco2_abate(tearly)::2);
Put /"objective function";
Put lim2xco2_utility::3;

* put file for hotelling results

Put / "hotelling rents run";
Put / "year";
Loop (tearly, put year(tearly)::0);
Put / "output";
Loop (tearly, put hotel_y(tearly)::3);
Put / "pccon";
Loop (tearly, put hotel_cpc(tearly)::3);
Put / "savrate";
Loop (tearly, put hotel_s(tearly)::4);
Put / "indem";
Loop (tearly, put hotel_indem(tearly)::4);
Put / "sigma";
Loop (tearly, put hotel_sigma(tearly)::4);
Put / "temp";
Loop (tearly, put hotel_tatm(tearly)::3);
Put / "conc";
Loop (tearly, put hotel_mat(tearly)::3);
Put / "soc cost carbon";
Loop (tearly, put hotel_tax(tearly)::2);
Put / "intrate";
Loop (tearly, put hotel_ri(tearly)::3);
Put / "discrate";
Loop (tearly, put hotel_rr(tearly)::5);
Put / "prod";
Loop (tearly, put hotel_al(tearly)::5);
Put / "exogforc";
Loop (tearly, put hotel_forcoth(tearly)::3);
Put / "pop";
Loop (tearly, put hotel_l(tearly)::3);
Put / "carbon tax";
Loop (tearly, put hotel_mcemis(tearly)::4);
Put / "margy";
Loop (tearly, put hotel_yy(tearly)::3);
Put / "margc";
Loop (tearly, put hotel_cc(tearly)::3);
Put / "miu";
Loop (tearly, put hotel_miu(tearly)::3);
Put / "total emissions";
Loop (tearly, put hotel_wem(tearly)::3);
Put / "interest rate";
Loop (tearly, put hotel_ri(tearly)::4);
Put / "damages";
Loop (tearly, put hotel_dam(tearly)::5);
Put / "abatement cost";
Loop (tearly, put hotel_abate(tearly)::5);
Put /"objective function";
Put hotel_utility::3;

* put file for limited participation results

Put / "limit participation run";
Put / "year";
Loop (tearly, put year(tearly)::0);
Put / "output";
Loop (tearly, put limpart_y(tearly)::3);
Put / "pccon";
Loop (tearly, put limpart_cpc(tearly)::3);
Put / "savrate";
Loop (tearly, put limpart_s(tearly)::4);
Put / "indem";
Loop (tearly, put limpart_indem(tearly)::4);
Put / "participation rate";
Loop (tearly, put limpart_partfract(tearly)::4);
Put / "temp";
Loop (tearly, put limpart_tatm(tearly)::3);
Put / "conc";
Loop (tearly, put limpart_mat(tearly)::3);
Put / "soc cost carbon";
Loop (tearly, put limpart_tax(tearly)::2);
Put / "intrate";
Loop (tearly, put limpart_ri(tearly)::3);
Put / "discrate";
Loop (tearly, put limpart_rr(tearly)::5);
Put / "prod";
Loop (tearly, put limpart_al(tearly)::5);
Put / "exogforc";
Loop (tearly, put limpart_forcoth(tearly)::3);
Put / "pop";
Loop (tearly, put limpart_l(tearly)::3);
Put / "carbon tax";
Loop (tearly, put limpart_mcemis(tearly)::4);
Put / "margy";
Loop (tearly, put limpart_yy(tearly)::3);
Put / "margc";
Loop (tearly, put limpart_cc(tearly)::3);
Put / "miu";
Loop (tearly, put limpart_miu(tearly)::3);
Put / "total emissions";
Loop (tearly, put limpart_wem(tearly)::3);
Put / "interest rate";
Loop (tearly, put limpart_ri(tearly)::4);
Put / "damages";
Loop (tearly, put limpart_dam(tearly)::5);
Put / "abatement cost";
Loop (tearly, put limpart_abate(tearly)::5);
Put /"objective function";
Put limpart_utility::3;
* put file for Optimal results


Put / "Stern controls";
Put / "year";
Loop (tearly, put year(tearly)::0);
Put / "output";
Loop (tearly, put stern_y(tearly)::3);
Put / "pccon";
Loop (tearly, put stern_cpc(tearly)::3);
Put / "savrate";
Loop (tearly, put stern_s(tearly)::4);
Put / "indem";
Loop (tearly, put stern_indem(tearly)::4);
Put / "sigma";
Loop (tearly, put stern_sigma(tearly)::4);
Put / "temp";
Loop (tearly, put stern_tatm(tearly)::3);
Put / "conc";
Loop (tearly, put stern_mat(tearly)::3);
Put / "soc cost carbon";
Loop (tearly, put stern_tax(tearly)::2);
Put / "intrate";
Loop (tearly, put stern_ri(tearly)::3);
Put / "discrate";
Loop (tearly, put stern_rr(tearly)::5);
Put / "prod";
Loop (tearly, put stern_al(tearly)::5);
Put / "exogforc";
Loop (tearly, put stern_forcoth(tearly)::3);
Put / "pop";
Loop (tearly, put stern_l(tearly)::3);
Put / "carbon tax";
Loop (tearly, put stern_mcemis(tearly)::4);
Put / "margy";
Loop (tearly, put stern_yy(tearly)::3);
Put / "margc";
Loop (tearly, put stern_cc(tearly)::3);
Put / "miu";
Loop (tearly, put stern_miu(tearly)::3);
Put / "total emissions";
Loop (tearly, put stern_wem(tearly)::3);
Put / "interest rate";
Loop (tearly, put stern_ri(tearly)::4);
Put / "damages";
Loop (tearly, put stern_dam(tearly)::5);
Put / "abatement cost";
Loop (tearly, put stern_abate(tearly)::5);
Put /"objective function";
Put stern_utility::3;

* put file for Kyoto with US


Put / "Kyoto with US";
Put / "year";
Loop (tearly, put year(tearly)::0);
Put / "output";
Loop (tearly, put kyoto_withus_y(tearly)::3);
Put / "pccon";
Loop (tearly, put kyoto_withus_cpc(tearly)::3);
Put / "savrate";
Loop (tearly, put kyoto_withus_s(tearly)::4);
Put / "indem";
Loop (tearly, put kyoto_withus_indem(tearly)::4);
Put / "participation rate";
Loop (tearly, put kyoto_withus_partfract(tearly)::4);
Put / "temp";
Loop (tearly, put kyoto_withus_tatm(tearly)::3);
Put / "conc";
Loop (tearly, put kyoto_withus_mat(tearly)::3);
Put / "soc cost carbon";
Loop (tearly, put kyoto_withus_tax(tearly)::2);
Put / "intrate";
Loop (tearly, put kyoto_withus_ri(tearly)::3);
Put / "discrate";
Loop (tearly, put kyoto_withus_rr(tearly)::5);
Put / "prod";
Loop (tearly, put kyoto_withus_al(tearly)::5);
Put / "exogforc";
Loop (tearly, put kyoto_withus_forcoth(tearly)::3);
Put / "pop";
Loop (tearly, put kyoto_withus_l(tearly)::3);
Put / "carbon tax";
Loop (tearly, put kyoto_withus_mcemis(tearly)::4);
Put / "margy";
Loop (tearly, put kyoto_withus_yy(tearly)::3);
Put / "margc";
Loop (tearly, put kyoto_withus_cc(tearly)::3);
Put / "miu";
Loop (tearly, put kyoto_withus_miu(tearly)::3);
Put / "total emissions";
Loop (tearly, put kyoto_withus_wem(tearly)::3);
Put / "interest rate";
Loop (tearly, put kyoto_withus_ri(tearly)::4);
Put / "damages";
Loop (tearly, put kyoto_withus_dam(tearly)::5);
Put / "abatement cost";
Loop (tearly, put kyoto_withus_abate(tearly)::5);
Put /"objective function";
Put kyoto_withus_utility::3;
* put file for Optimal results


Put / "Kyoto without US";
Put / "year";
Loop (tearly, put year(tearly)::0);
Put / "output";
Loop (tearly, put kyoto_nous_y(tearly)::3);
Put / "pccon";
Loop (tearly, put kyoto_nous_cpc(tearly)::3);
Put / "savrate";
Loop (tearly, put kyoto_nous_s(tearly)::4);
Put / "indem";
Loop (tearly, put kyoto_nous_indem(tearly)::4);
Put / "participation rate";
Loop (tearly, put kyoto_nous_partfract(tearly)::4);
Put / "temp";
Loop (tearly, put kyoto_nous_tatm(tearly)::3);
Put / "conc";
Loop (tearly, put kyoto_nous_mat(tearly)::3);
Put / "soc cost carbon";
Loop (tearly, put kyoto_nous_tax(tearly)::2);
Put / "intrate";
Loop (tearly, put kyoto_nous_ri(tearly)::3);
Put / "discrate";
Loop (tearly, put kyoto_nous_rr(tearly)::5);
Put / "prod";
Loop (tearly, put kyoto_nous_al(tearly)::5);
Put / "exogforc";
Loop (tearly, put kyoto_nous_forcoth(tearly)::3);
Put / "pop";
Loop (tearly, put kyoto_nous_l(tearly)::3);
Put / "carbon tax";
Loop (tearly, put kyoto_nous_mcemis(tearly)::4);
Put / "margy";
Loop (tearly, put kyoto_nous_yy(tearly)::3);
Put / "margc";
Loop (tearly, put kyoto_nous_cc(tearly)::3);
Put / "miu";
Loop (tearly, put kyoto_nous_miu(tearly)::3);
Put / "total emissions";
Loop (tearly, put kyoto_nous_wem(tearly)::3);
Put / "interest rate";
Loop (tearly, put kyoto_nous_ri(tearly)::4);
Put / "damages";
Loop (tearly, put kyoto_nous_dam(tearly)::5);
Put / "abatement cost";
Loop (tearly, put kyoto_nous_abate(tearly)::8);
Put /"objective function";
Put kyoto_nous_utility::3;
* put file for Optimal results



Put / "optimal run";
Put / "year";
Loop (tlate, put year(tlate)::0);
Put / "output";
Loop (tlate, put opt_y(tlate)::3);
Put / "pccon";
Loop (tlate, put opt_cpc(tlate)::3);
Put / "savrate";
Loop (tlate, put opt_s(tlate)::4);
Put / "indem";
Loop (tlate, put opt_indem(tlate)::4);
Put / "sigma";
Loop (tlate, put opt_sigma(tlate)::4);
Put / "temp";
Loop (tlate, put opt_tatm(tlate)::3);
Put / "conc";
Loop (tlate, put opt_mat(tlate)::3);
Put / "soc cost carbon";
Loop (tlate, put opt_tax(tlate)::2);
Put / "intrate";
Loop (tlate, put opt_ri(tlate)::3);
Put / "discrate";
Loop (tlate, put opt_rr(tlate)::5);
Put / "prod";
Loop (tlate, put opt_al(tlate)::5);
Put / "exogforc";
Loop (tlate, put opt_forcoth(tlate)::3);
Put / "pop";
Loop (tlate, put opt_l(tlate)::3);
Put / "carbon tax";
Loop (tlate, put opt_mcemis(tlate)::4);
Put / "margy";
Loop (tlate, put opt_yy(tlate)::3);
Put / "margc";
Loop (tlate, put opt_cc(tlate)::7);
Put / "miu";
Loop (tlate, put opt_miu(tlate)::3);
Put / "total emissions";
Loop (tlate, put opt_wem(tlate)::3);
Put / "interest rate";
Loop (tlate, put opt_ri(tlate)::4);
Put / "damages";
Loop (tlate, put opt_dam(tlate)::2);
Put / "abatement cost";
Loop (tlate, put opt_abate(tlate)::2);
Put /"objective function";
Put opt_utility::3;

* put file for 15 period control

Put / "base run with no controls for 15 periods";
Put / "year";
Loop (tlate, put year(tlate)::0);
Put / "output";
Loop (tlate, put base25_y(tlate)::3);
Put / "pccon";
Loop (tlate, put base25_cpc(tlate)::3);
Put / "savrate";
Loop (tlate, put base25_s(tlate)::4);
Put / "indem";
Loop (tlate, put base25_indem(tlate)::4);
Put / "sigma";
Loop (tlate, put base25_sigma(tlate)::4);
Put / "temp";
Loop (tlate, put base25_tatm(tlate)::3);
Put / "conc";
Loop (tlate, put base25_mat(tlate)::3);
Put / "soc cost carbon";
Loop (tlate, put base25_tax(tlate)::2);
Put / "intrate";
Loop (tlate, put base25_ri(tlate)::3);
Put / "discrate";
Loop (tlate, put base25_rr(tlate)::5);
Put / "prod";
Loop (tlate, put base25_al(tlate)::5);
Put / "exogforc";
Loop (tlate, put base25_forcoth(tlate)::3);
Put / "pop";
Loop (tlate, put base25_l(tlate)::3);
Put / "carbon tax";
Loop (tlate, put base25_mcemis(tlate)::4);
Put / "margy";
Loop (tlate, put base25_yy(tlate)::3);
Put / "margc";
Loop (tlate, put base25_cc(tlate)::3);
Put / "miu";
Loop (tlate, put base25_miu(tlate)::3);
Put / "total emissions";
Loop (tlate, put base25_wem(tlate)::3);
Put / "interest rate";
Loop (tlate, put base25_ri(tlate)::4);
Put / "damages";
Loop (tlate, put base25_dam(tlate)::2);
Put / "abatement cost";
Loop (tlate, put base25_abate(tlate)::2);
Put /"objective function";
Put base25_utility::3;

*put file for 5 period control


Put / "base run with no controls for 5 periods";
Put / "year";
Loop (tlate, put year(tlate)::0);
Put / "output";
Loop (tlate, put base5_y(tlate)::3);
Put / "pccon";
Loop (tlate, put base5_cpc(tlate)::3);
Put / "savrate";
Loop (tlate, put base5_s(tlate)::4);
Put / "indem";
Loop (tlate, put base5_indem(tlate)::4);
Put / "sigma";
Loop (tlate, put base5_sigma(tlate)::4);
Put / "temp";
Loop (tlate, put base5_tatm(tlate)::3);
Put / "conc";
Loop (tlate, put base5_mat(tlate)::3);
Put / "soc cost carbon";
Loop (tlate, put base5_tax(tlate)::2);
Put / "intrate";
Loop (tlate, put base5_ri(tlate)::3);
Put / "discrate";
Loop (tlate, put base5_rr(tlate)::5);
Put / "prod";
Loop (tlate, put base5_al(tlate)::5);
Put / "exogforc";
Loop (tlate, put base5_forcoth(tlate)::3);
Put / "pop";
Loop (tlate, put base5_l(tlate)::3);
Put / "carbon tax";
Loop (tlate, put base5_mcemis(tlate)::4);
Put / "margy";
Loop (tlate, put base5_yy(tlate)::3);
Put / "margc";
Loop (tlate, put base5_cc(tlate)::3);
Put / "miu";
Loop (tlate, put base5_miu(tlate)::3);
Put / "total emissions";
Loop (tlate, put base5_wem(tlate)::3);
Put / "interest rate";
Loop (tlate, put base5_ri(tlate)::4);
Put / "damages";
Loop (tlate, put base5_dam(tlate)::2);
Put / "abatement cost";
Loop (tlate, put base5_abate(tlate)::2);
Put /"objective function";
Put base5_utility::3;

*put file for limit to 2 degrees

Put / "limit T increase to 2 degrees";
Put / "year";
Loop (tlate, put year(tlate)::0);
Put / "output";
Loop (tlate, put lim2deg_y(tlate)::3);
Put / "pccon";
Loop (tlate, put lim2deg_cpc(tlate)::3);
Put / "savrate";
Loop (tlate, put lim2deg_s(tlate)::4);
Put / "indem";
Loop (tlate, put lim2deg_indem(tlate)::4);
Put / "sigma";
Loop (tlate, put lim2deg_sigma(tlate)::4);
Put / "temp";
Loop (tlate, put lim2deg_tatm(tlate)::3);
Put / "conc";
Loop (tlate, put lim2deg_mat(tlate)::3);
Put / "soc cost carbon";
Loop (tlate, put lim2deg_tax(tlate)::2);
Put / "intrate";
Loop (tlate, put lim2deg_ri(tlate)::3);
Put / "discrate";
Loop (tlate, put lim2deg_rr(tlate)::5);
Put / "prod";
Loop (tlate, put lim2deg_al(tlate)::5);
Put / "exogforc";
Loop (tlate, put lim2deg_forcoth(tlate)::3);
Put / "pop";
Loop (tlate, put lim2deg_l(tlate)::3);
Put / "carbon tax";
Loop (tlate, put lim2deg_mcemis(tlate)::4);
Put / "margy";
Loop (tlate, put lim2deg_yy(tlate)::3);
Put / "margc";
Loop (tlate, put lim2deg_cc(tlate)::3);
Put / "miu";
Loop (tlate, put lim2deg_miu(tlate)::3);
Put / "total emissions";
Loop (tlate, put lim2deg_wem(tlate)::3);
Put / "interest rate";
Loop (tlate, put lim2deg_ri(tlate)::4);
Put / "damages";
Loop (tlate, put lim2deg_dam(tlate)::2);
Put / "abatement cost";
Loop (tlate, put lim2deg_abate(tlate)::2);
Put /"objective function";
Put lim2deg_utility::3;


*put file for limit to 2.5 degrees

Put / "limit T increase to 2.5 degrees";
Put / "year";
Loop (tlate, put year(tlate)::0);
Put / "output";
Loop (tlate, put lim25deg_y(tlate)::3);
Put / "pccon";
Loop (tlate, put lim25deg_cpc(tlate)::3);
Put / "savrate";
Loop (tlate, put lim25deg_s(tlate)::4);
Put / "indem";
Loop (tlate, put lim25deg_indem(tlate)::4);
Put / "sigma";
Loop (tlate, put lim25deg_sigma(tlate)::4);
Put / "temp";
Loop (tlate, put lim25deg_tatm(tlate)::3);
Put / "conc";
Loop (tlate, put lim25deg_mat(tlate)::3);
Put / "soc cost carbon";
Loop (tlate, put lim25deg_tax(tlate)::2);
Put / "intrate";
Loop (tlate, put lim25deg_ri(tlate)::3);
Put / "discrate";
Loop (tlate, put lim25deg_rr(tlate)::5);
Put / "prod";
Loop (tlate, put lim25deg_al(tlate)::5);
Put / "exogforc";
Loop (tlate, put lim25deg_forcoth(tlate)::3);
Put / "pop";
Loop (tlate, put lim25deg_l(tlate)::3);
Put / "carbon tax";
Loop (tlate, put lim25deg_mcemis(tlate)::4);
Put / "margy";
Loop (tlate, put lim25deg_yy(tlate)::3);
Put / "margc";
Loop (tlate, put lim25deg_cc(tlate)::3);
Put / "miu";
Loop (tlate, put lim25deg_miu(tlate)::3);
Put / "total emissions";
Loop (tlate, put lim25deg_wem(tlate)::3);
Put / "interest rate";
Loop (tlate, put lim25deg_ri(tlate)::4);
Put / "damages";
Loop (tlate, put lim25deg_dam(tlate)::2);
Put / "abatement cost";
Loop (tlate, put lim25deg_abate(tlate)::2);
Put /"objective function";
Put lim25deg_utility::3;



*put file for limit to 2x preindustrial co2

Put / "limit co2 to 2x preindustrial";
Put / "year";
Loop (tlate, put year(tlate)::0);
Put / "output";
Loop (tlate, put lim2xco2_y(tlate)::3);
Put / "pccon";
Loop (tlate, put lim2xco2_cpc(tlate)::3);
Put / "savrate";
Loop (tlate, put lim2xco2_s(tlate)::4);
Put / "indem";
Loop (tlate, put lim2xco2_indem(tlate)::4);
Put / "sigma";
Loop (tlate, put lim2xco2_sigma(tlate)::4);
Put / "temp";
Loop (tlate, put lim2xco2_tatm(tlate)::3);
Put / "conc";
Loop (tlate, put lim2xco2_mat(tlate)::3);
Put / "soc cost carbon";
Loop (tlate, put lim2xco2_tax(tlate)::2);
Put / "intrate";
Loop (tlate, put lim2xco2_ri(tlate)::3);
Put / "discrate";
Loop (tlate, put lim2xco2_rr(tlate)::5);
Put / "prod";
Loop (tlate, put lim2xco2_al(tlate)::5);
Put / "exogforc";
Loop (tlate, put lim2xco2_forcoth(tlate)::3);
Put / "pop";
Loop (tlate, put lim2xco2_l(tlate)::3);
Put / "carbon tax";
Loop (tlate, put lim2xco2_mcemis(tlate)::4);
Put / "margy";
Loop (tlate, put lim2xco2_yy(tlate)::3);
Put / "margc";
Loop (tlate, put lim2xco2_cc(tlate)::3);
Put / "miu";
Loop (tlate, put lim2xco2_miu(tlate)::3);
Put / "total emissions";
Loop (tlate, put lim2xco2_wem(tlate)::3);
Put / "interest rate";
Loop (tlate, put lim2xco2_ri(tlate)::4);
Put / "damages";
Loop (tlate, put lim2xco2_dam(tlate)::2);
Put / "abatement cost";
Loop (tlate, put lim2xco2_abate(tlate)::2);
Put /"objective function";
Put lim2xco2_utility::3;



* put file for hotelling results

Put / "hotelling rents run";
Put / "year";
Loop (tlate, put year(tlate)::0);
Put / "output";
Loop (tlate, put hotel_y(tlate)::3);
Put / "pccon";
Loop (tlate, put hotel_cpc(tlate)::3);
Put / "savrate";
Loop (tlate, put hotel_s(tlate)::4);
Put / "indem";
Loop (tlate, put hotel_indem(tlate)::4);
Put / "sigma";
Loop (tlate, put hotel_sigma(tlate)::4);
Put / "temp";
Loop (tlate, put hotel_tatm(tlate)::3);
Put / "conc";
Loop (tlate, put hotel_mat(tlate)::3);
Put / "soc cost carbon";
Loop (tlate, put hotel_tax(tlate)::2);
Put / "intrate";
Loop (tlate, put hotel_ri(tlate)::3);
Put / "discrate";
Loop (tlate, put hotel_rr(tlate)::5);
Put / "prod";
Loop (tlate, put hotel_al(tlate)::5);
Put / "exogforc";
Loop (tlate, put hotel_forcoth(tlate)::3);
Put / "pop";
Loop (tlate, put hotel_l(tlate)::3);
Put / "carbon tax";
Loop (tlate, put hotel_mcemis(tlate)::4);
Put / "margy";
Loop (tlate, put hotel_yy(tlate)::3);
Put / "margc";
Loop (tlate, put hotel_cc(tlate)::3);
Put / "miu";
Loop (tlate, put hotel_miu(tlate)::3);
Put / "total emissions";
Loop (tlate, put hotel_wem(tlate)::3);
Put / "interest rate";
Loop (tlate, put hotel_ri(tlate)::4);
Put / "damages";
Loop (tlate, put hotel_dam(tlate)::5);
Put / "abatement cost";
Loop (tlate, put hotel_abate(tlate)::5);
Put /"objective function";
Put hotel_utility::3;


* put file for limited participation results

Put / "limit participation run";
Put / "year";
Loop (tlate, put year(tlate)::0);
Put / "output";
Loop (tlate, put limpart_y(tlate)::3);
Put / "pccon";
Loop (tlate, put limpart_cpc(tlate)::3);
Put / "savrate";
Loop (tlate, put limpart_s(tlate)::4);
Put / "indem";
Loop (tlate, put limpart_indem(tlate)::4);
Put / "participation rate";
Loop (tlate, put limpart_partfract(tlate)::4);
Put / "temp";
Loop (tlate, put limpart_tatm(tlate)::3);
Put / "conc";
Loop (tlate, put limpart_mat(tlate)::3);
Put / "soc cost carbon";
Loop (tlate, put limpart_tax(tlate)::2);
Put / "intrate";
Loop (tlate, put limpart_ri(tlate)::3);
Put / "discrate";
Loop (tlate, put limpart_rr(tlate)::5);
Put / "prod";
Loop (tlate, put limpart_al(tlate)::5);
Put / "exogforc";
Loop (tlate, put limpart_forcoth(tlate)::3);
Put / "pop";
Loop (tlate, put limpart_l(tlate)::3);
Put / "carbon tax";
Loop (tlate, put limpart_mcemis(tlate)::4);
Put / "margy";
Loop (tlate, put limpart_yy(tlate)::3);
Put / "margc";
Loop (tlate, put limpart_cc(tlate)::3);
Put / "miu";
Loop (tlate, put limpart_miu(tlate)::3);
Put / "total emissions";
Loop (tlate, put limpart_wem(tlate)::3);
Put / "interest rate";
Loop (tlate, put limpart_ri(tlate)::4);
Put / "damages";
Loop (tlate, put limpart_dam(tlate)::5);
Put / "abatement cost";
Loop (tlate, put limpart_abate(tlate)::5);
Put /"objective function";
Put limpart_utility::3;



* put file for original damage

Put / "Stern controls";
Put / "year";
Loop (tlate, put year(tlate)::0);
Put / "output";
Loop (tlate, put stern_y(tlate)::3);
Put / "pccon";
Loop (tlate, put stern_cpc(tlate)::3);
Put / "savrate";
Loop (tlate, put stern_s(tlate)::4);
Put / "indem";
Loop (tlate, put stern_indem(tlate)::4);
Put / "sigma";
Loop (tlate, put stern_sigma(tlate)::4);
Put / "temp";
Loop (tlate, put stern_tatm(tlate)::3);
Put / "conc";
Loop (tlate, put stern_mat(tlate)::3);
Put / "soc cost carbon";
Loop (tlate, put stern_tax(tlate)::2);
Put / "intrate";
Loop (tlate, put stern_ri(tlate)::3);
Put / "discrate";
Loop (tlate, put stern_rr(tlate)::5);
Put / "prod";
Loop (tlate, put stern_al(tlate)::5);
Put / "exogforc";
Loop (tlate, put stern_forcoth(tlate)::3);
Put / "pop";
Loop (tlate, put stern_l(tlate)::3);
Put / "carbon tax";
Loop (tlate, put stern_mcemis(tlate)::4);
Put / "margy";
Loop (tlate, put stern_yy(tlate)::3);
Put / "margc";
Loop (tlate, put stern_cc(tlate)::3);
Put / "miu";
Loop (tlate, put stern_miu(tlate)::3);
Put / "total emissions";
Loop (tlate, put stern_wem(tlate)::3);
Put / "interest rate";
Loop (tlate, put stern_ri(tlate)::4);
Put / "damages";
Loop (tlate, put stern_dam(tlate)::5);
Put / "abatement cost";
Loop (tlate, put stern_abate(tlate)::5);
Put /"objective function";
Put stern_utility::3;



* put file for Kyoto with US late

Put / "Kyoto with US";
Put / "year";
Loop (tlate, put year(tlate)::0);
Put / "output";
Loop (tlate, put kyoto_withus_y(tlate)::3);
Put / "pccon";
Loop (tlate, put kyoto_withus_cpc(tlate)::3);
Put / "savrate";
Loop (tlate, put kyoto_withus_s(tlate)::4);
Put / "indem";
Loop (tlate, put kyoto_withus_indem(tlate)::4);
Put / "participation rate";
Loop (tlate, put kyoto_withus_partfract(tlate)::4);
Put / "temp";
Loop (tlate, put kyoto_withus_tatm(tlate)::3);
Put / "conc";
Loop (tlate, put kyoto_withus_mat(tlate)::3);
Put / "soc cost carbon";
Loop (tlate, put kyoto_withus_tax(tlate)::2);
Put / "intrate";
Loop (tlate, put kyoto_withus_ri(tlate)::3);
Put / "discrate";
Loop (tlate, put kyoto_withus_rr(tlate)::5);
Put / "prod";
Loop (tlate, put kyoto_withus_al(tlate)::5);
Put / "exogforc";
Loop (tlate, put kyoto_withus_forcoth(tlate)::3);
Put / "pop";
Loop (tlate, put kyoto_withus_l(tlate)::3);
Put / "carbon tax";
Loop (tlate, put kyoto_withus_mcemis(tlate)::4);
Put / "margy";
Loop (tlate, put kyoto_withus_yy(tlate)::3);
Put / "margc";
Loop (tlate, put kyoto_withus_cc(tlate)::3);
Put / "miu";
Loop (tlate, put kyoto_withus_miu(tlate)::3);
Put / "total emissions";
Loop (tlate, put kyoto_withus_wem(tlate)::3);
Put / "interest rate";
Loop (tlate, put kyoto_withus_ri(tlate)::4);
Put / "damages";
Loop (tlate, put kyoto_withus_dam(tlate)::5);
Put / "abatement cost";
Loop (tlate, put kyoto_withus_abate(tlate)::5);
Put /"objective function";
Put kyoto_withus_utility::3;


* put file for Kyoto without US run

Put / "Kyoto without US";
Put / "year";
Loop (tlate, put year(tlate)::0);
Put / "output";
Loop (tlate, put kyoto_nous_y(tlate)::3);
Put / "pccon";
Loop (tlate, put kyoto_nous_cpc(tlate)::3);
Put / "savrate";
Loop (tlate, put kyoto_nous_s(tlate)::4);
Put / "indem";
Loop (tlate, put kyoto_nous_indem(tlate)::4);
Put / "participation rate";
Loop (tlate, put kyoto_nous_partfract(tlate)::4);
Put / "temp";
Loop (tlate, put kyoto_nous_tatm(tlate)::3);
Put / "conc";
Loop (tlate, put kyoto_nous_mat(tlate)::3);
Put / "soc cost carbon";
Loop (tlate, put kyoto_nous_tax(tlate)::2);
Put / "intrate";
Loop (tlate, put kyoto_nous_ri(tlate)::3);
Put / "discrate";
Loop (tlate, put kyoto_nous_rr(tlate)::5);
Put / "prod";
Loop (tlate, put kyoto_nous_al(tlate)::5);
Put / "exogforc";
Loop (tlate, put kyoto_nous_forcoth(tlate)::3);
Put / "pop";
Loop (tlate, put kyoto_nous_l(tlate)::3);
Put / "carbon tax";
Loop (tlate, put kyoto_nous_mcemis(tlate)::4);
Put / "margy";
Loop (tlate, put kyoto_nous_yy(tlate)::3);
Put / "margc";
Loop (tlate, put kyoto_nous_cc(tlate)::3);
Put / "miu";
Loop (tlate, put kyoto_nous_miu(tlate)::3);
Put / "total emissions";
Loop (tlate, put kyoto_nous_wem(tlate)::3);
Put / "interest rate";
Loop (tlate, put kyoto_nous_ri(tlate)::4);
Put / "damages";
Loop (tlate, put kyoto_nous_dam(tlate)::5);
Put / "abatement cost";
Loop (tlate, put kyoto_nous_abate(tlate)::5);
Put /"objective function";
Put kyoto_nous_utility::3;