* Definition of base5 results

* Control statements 

mat.up(t)    = 3000;
TATM.up(t)   = 12 ;

MIU.lo("1")=miuhotel("1");
MIU.lo("2")=miuhotel("2");
MIU.lo("3")=miuhotel("3");
MIU.lo("4")=miuhotel("4");
MIU.lo("5")=miuhotel("5");

MIU.up("1")=miuhotel("1");
MIU.up("2")=miuhotel("2");
MIU.up("3")=miuhotel("3");
MIU.up("4")=miuhotel("4");
MIU.up("5")=miuhotel("5");




model co2lim5 /all/
solve CO2lim5 maximizing UTILITY using nlp ;
solve CO2lim5 maximizing UTILITY using nlp ;
solve CO2lim5 maximizing UTILITY using nlp ;

* Output

Parameters
Year(t)         Date
base5_y(t)
base5_cpc(t)
base5_s(t)
base5_indem(t)
base5_sigma(t)
base5_tatm(t)
base5_mat(t)
base5_tax(t)
base5_ri(t)
base5_rr(t)
base5_al(t)
base5_forcoth(t)
base5_l(t)
base5_etree(t)
base5_yy(t)
base5_cc(t)
base5_miu(t)
base5_wem(t)
base5_ri(t)
base5_dam(t)
base5_abate(t)
base5_mcemis(t)
base5_utility   ;

Year(t)         = 2005 +10*(ord(t)-1);
base5_y(t)=y.l(t);
base5_cpc(t)=cpc.l(t);
base5_s(t)=s.l(t)     ;
base5_indem(t)= e.l(t)-etree(t);;
base5_sigma(t)=sigma(t) ;
base5_tatm(t)=tatm.l(t)  ;
base5_mat(t)=mat.l(t)     ;
base5_tax(t)=-1*ee.m(t)*1000/(kk.m(t)+.0000001)       ;
base5_ri(t)=ri.l(t);
base5_rr(t)=rr(t)   ;
base5_al(t)=al(t)    ;
base5_forcoth(t)=forcoth(t);
base5_l(t)=l(t);
base5_etree(t)=etree(t);
base5_yy(t)=yy.m(t)     ;
base5_cc(t)=cc.m(t)      ;
base5_miu(t)=miu.l(t)     ;
base5_wem(t)= e.l(t);
base5_ri(t)=ri.l(t)         ;
base5_dam(t)= damages.l(t);
base5_abate(t) = abatecost.l(t);
base5_utility=utility.l        ;
base5_mcemis(t)= expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000;


miu.up(t)= 1;
miu.lo(t)= 0;
miu.up("1")     = 0.2;