* Definition of Kyoto runs with US

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");

partfract(T)=1; 

partfract("1")=	0.25372	; 
partfract("2")=	0.45647	;
partfract("3")=	0.41182	;
partfract("4")=	0.37981	;
partfract("5")=	0.35337	;
partfract("6")=	0.33010	;
partfract("7")=	0.30974	;
partfract("8")=	0.29203	;
partfract("9")=	0.27667	;
partfract("10")=	0.26338	;
partfract("11")=	0.25192	;
partfract("12")=	0.24206	;
 

 		
 
miu.up("2")=	0.090171958	;
miu.up("3")=	0.093733295	;
miu.up("4")=	0.091635228	;
miu.up("5")=	0.089109129	;
miu.up("6")=	0.087445789	;
miu.up("7")=	0.086385877	;
miu.up("8")=	0.085760818	;
miu.up("9")=	0.085460649	;
miu.up("10")=	0.085412828	;
miu.up("11")=	0.085568495	;
miu.up("12")=	0.085893697	;
 


 
miu.lo("2")=	0.090171958	;
miu.lo("3")=	0.093733295	;
miu.lo("4")=	0.091635228	;
miu.lo("5")=	0.089109129	;
miu.lo("6")=	0.087445789	;
miu.lo("7")=	0.086385877	;
miu.lo("8")=	0.085760818	;
miu.lo("9")=	0.085460649	;
miu.lo("10")=	0.085412828	;
miu.lo("11")=	0.085568495	;
miu.lo("12")=	0.085893697	;
 



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

Parameters
Year(t)         Date
kyoto_withus_y(t)
kyoto_withus_cpc(t)
kyoto_withus_s(t)
kyoto_withus_indem(t)
kyoto_withus_sigma(t)
kyoto_withus_tatm(t)
kyoto_withus_mat(t)
kyoto_withus_tax(t)
kyoto_withus_tax_part(t)
kyoto_withus_ri(t)
kyoto_withus_rr(t)
kyoto_withus_al(t)
kyoto_withus_forcoth(t)
kyoto_withus_l(t)
kyoto_withus_etree(t)
kyoto_withus_yy(t)
kyoto_withus_cc(t)
kyoto_withus_miu(t)
kyoto_withus_wem(t)
kyoto_withus_ri(t)
kyoto_withus_dam(t)
kyoto_withus_abate(t)
kyoto_withus_mcemis(t)
kyoto_withus_mcemis_part(t)
kyoto_withus_partfract(t)
kyoto_withus_utility   ;

Year(t)         = 2005 +10*(ord(t)-1);
kyoto_withus_y(t)=y.l(t);
kyoto_withus_cpc(t)=cpc.l(t);
kyoto_withus_s(t)=s.l(t)     ;
kyoto_withus_indem(t)= e.l(t)-etree(t);;
kyoto_withus_sigma(t)=sigma(t) ;
kyoto_withus_tatm(t)=tatm.l(t)  ;
kyoto_withus_mat(t)=mat.l(t)     ;
kyoto_withus_tax(t)=-1*ee.m(t)*1000/(kk.m(t)+.0000001)       ;
kyoto_withus_ri(t)=ri.l(t);
kyoto_withus_rr(t)=rr(t)   ;
kyoto_withus_al(t)=al(t)    ;
kyoto_withus_forcoth(t)=forcoth(t);
kyoto_withus_l(t)=l(t);
kyoto_withus_etree(t)=etree(t);
kyoto_withus_yy(t)=yy.m(t)     ;
kyoto_withus_cc(t)=cc.m(t)      ;
kyoto_withus_miu(t)=miu.l(t)     ;
kyoto_withus_wem(t)= e.l(t);
kyoto_withus_ri(t)=ri.l(t)         ;
kyoto_withus_dam(t)= damages.l(t);
kyoto_withus_abate(t) = abatecost.l(t);
kyoto_withus_mcemis(t)= expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000;
kyoto_withus_mcemis_part(t)= expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000*(partfract(t)**(1-expcost2)) ;
kyoto_withus_partfract(t)=partfract(t);
kyoto_withus_utility=utility.l        ;