* Definition of base_350yr results

* Control statements
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.lo("6")=miuhotel("6");
MIU.lo("7")=miuhotel("7");
MIU.lo("8")=miuhotel("8");
MIU.lo("9")=miuhotel("9");
MIU.lo("10")=miuhotel("10");
MIU.lo("11")=miuhotel("11");
MIU.lo("12")=miuhotel("12");
MIU.lo("13")=miuhotel("13");
MIU.lo("14")=miuhotel("14");
MIU.lo("15")=miuhotel("15");
MIU.lo("16")=miuhotel("16");
MIU.lo("17")=miuhotel("17");
MIU.lo("18")=miuhotel("18");
MIU.lo("19")=miuhotel("19");
MIU.lo("20")=miuhotel("20");
MIU.lo("21")=miuhotel("21");
MIU.lo("22")=miuhotel("22");
MIU.lo("23")=miuhotel("23");
MIU.lo("24")=miuhotel("24");
MIU.lo("25")=miuhotel("25");
MIU.lo("26")=miuhotel("26");
MIU.lo("27")=miuhotel("27");
MIU.lo("28")=miuhotel("28");
MIU.lo("29")=miuhotel("29");
MIU.lo("30")=miuhotel("30");
MIU.lo("31")=miuhotel("31");
MIU.lo("32")=miuhotel("32");
MIU.lo("33")=miuhotel("33");
MIU.lo("34")=miuhotel("34");
MIU.lo("35")=miuhotel("35");

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");
MIU.up("6")=miuhotel("6");
MIU.up("7")=miuhotel("7");
MIU.up("8")=miuhotel("8");
MIU.up("9")=miuhotel("9");
MIU.up("10")=miuhotel("10");
MIU.up("11")=miuhotel("11");
MIU.up("12")=miuhotel("12");
MIU.up("13")=miuhotel("13");
MIU.up("14")=miuhotel("14");
MIU.up("15")=miuhotel("15");
MIU.up("16")=miuhotel("16");
MIU.up("17")=miuhotel("17");
MIU.up("18")=miuhotel("18");
MIU.up("19")=miuhotel("19");
MIU.up("20")=miuhotel("20");
MIU.up("21")=miuhotel("21");
MIU.up("22")=miuhotel("22");
MIU.up("23")=miuhotel("23");
MIU.up("24")=miuhotel("24");
MIU.up("25")=miuhotel("25");
MIU.up("26")=miuhotel("26");
MIU.up("27")=miuhotel("27");
MIU.up("28")=miuhotel("28");
MIU.up("29")=miuhotel("29");
MIU.up("30")=miuhotel("30");
MIU.up("31")=miuhotel("31");
MIU.up("32")=miuhotel("32");
MIU.up("33")=miuhotel("33");
MIU.up("34")=miuhotel("34");
MIU.up("35")=miuhotel("35");

solve CO2 maximizing UTILITY using nlp ;
solve CO2 maximizing UTILITY using nlp ;
solve CO2 maximizing UTILITY using nlp ;

*Output
Parameters
Year(t)         Date
base35_y(t)
base35_cpc(t)
base35_s(t)
base35_indem(t)
base35_sigma(t)
base35_tatm(t)
base35_mat(t)
base35_tax(t)
base35_ri(t)
base35_rr(t)
base35_al(t)
base35_forcoth(t)
base35_l(t)
base35_etree(t)
base35_yy(t)
base35_cc(t)
base35_miu(t)
base35_wem(t)
base35_ri(t)
base35_dam(t)
base35_abate(t)
base35_mcemis(t)
base35_mcemis(t)
base35_utility   ;

Year(t)         = 2005 +10*(ord(t)-1);
base35_y(t)=y.l(t);
base35_cpc(t)=cpc.l(t);
base35_s(t)=s.l(t)     ;
base35_indem(t)= e.l(t)-etree(t);;
base35_sigma(t)=sigma(t) ;
base35_tatm(t)=tatm.l(t)  ;
base35_mat(t)=mat.l(t)     ;
base35_tax(t)=-1*ee.m(t)*1000/(kk.m(t)+.00000000001)       ;
base35_ri(t)=ri.l(t);
base35_rr(t)=rr(t)   ;
base35_al(t)=al(t)    ;
base35_forcoth(t)=forcoth(t);
base35_l(t)=l(t);
base35_etree(t)=etree(t);
base35_yy(t)=yy.m(t)     ;
base35_cc(t)=cc.m(t)      ;
base35_miu(t)=miu.l(t)     ;
base35_wem(t)= e.l(t);
base35_ri(t)=ri.l(t)         ;
base35_dam(t)= damages.l(t);
base35_abate(t) = abatecost.l(t);
base35_mcemis(t)= expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000;
base35_utility=utility.l        ;
base35_mcemis(t) = expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000;