* Definition of limit 2x co2 results with no damages


* Control statements

* First do the run with zero damages

aa1 = 0;
aa2 = 0 ;
aa3 = 0 ;
TATM.up(t) = 12 ;
*mat.up(t) = 1270;
mat.up(t) = 958;
solve CO2 maximizing UTILITY using nlp ;
solve CO2 maximizing UTILITY using nlp ;
solve CO2 maximizing UTILITY using nlp ;

* Get the control rates for this constraint

Parameters
miu450(t)   Result for run without damages;
miu450(t) = miu.l(t);
display miu450, mat.l;
* Then do the run with damages

aa1 = a1;
aa2 = a2 ;
aa3 = a3 ;

miu.up(t)= miu450(t);
miu.lo(t)= miu450(t);
solve CO2 maximizing UTILITY using nlp ;
solve CO2 maximizing UTILITY using nlp ;
solve CO2 maximizing UTILITY using nlp ;
* Definition of miu450 results

Parameters
Year(t)         Date
lim_450_y(t)
lim_450_cpc(t)
lim_450_s(t)
lim_450_indem(t)
lim_450_sigma(t)
lim_450_tatm(t)
lim_450_mat(t)
lim_450_tax(t)
lim_450_ri(t)
lim_450_rr(t)
lim_450_al(t)
lim_450_forcoth(t)
lim_450_l(t)
lim_450_etree(t)
lim_450_yy(t)
lim_450_cc(t)
lim_450_miu(t)
lim_450_wem(t)
lim_450_dam(t)
lim_450_abate(t)
lim_450_mcemis(t)
lim_450_utility   ;

Year(t)         = 2005 +10*(ord(t)-1);
lim_450_y(t)=y.l(t);
lim_450_cpc(t)=cpc.l(t);
lim_450_s(t)=s.l(t)     ;
lim_450_indem(t)= e.l(t)-etree(t);;
lim_450_sigma(t)=sigma(t) ;
lim_450_tatm(t)=tatm.l(t)  ;
lim_450_mat(t)=mat.l(t)     ;
lim_450_tax(t)=-1*ee.m(t)*1000/(kk.m(t)+.0000001)       ;
lim_450_ri(t)=ri.l(t);
lim_450_rr(t)=rr(t)   ;
lim_450_al(t)=al(t)    ;
lim_450_forcoth(t)=forcoth(t);
lim_450_l(t)=l(t);
lim_450_etree(t)=etree(t);
lim_450_yy(t)=yy.m(t)     ;
lim_450_cc(t)=cc.m(t) ;
lim_450_miu(t)=miu.l(t)     ;
lim_450_wem(t)= e.l(t);
lim_450_dam(t)= damages.l(t);
lim_450_abate(t) = abatecost.l(t);
lim_450_utility=utility.l        ;
lim_450_mcemis(t)= expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000;