* Definition of strong Kyoto runs
 

partfract(T)=1;
			
partfract("1")=	0.25372	;
partfract("2")=	0.45766	;
partfract("3")=	0.62126	;
partfract("4")=	0.74175	;
partfract("5")=	0.81769	;
partfract("6")=	0.85136	;
partfract("7")=	0.85946	;
partfract("8")=	0.85839	;
partfract("9")=	0.85749	;
partfract("10")=	0.85673	;
partfract("11")=	0.85607	;
partfract("12")=	0.85549	;
partfract("13")=	0.85497	;
partfract("14")=	0.85450	;
partfract("15")=	0.85406	; 
 		
 
miu.up("2")=	0.105192821	;
miu.up("3")=	0.184429898	;
miu.up("4")=	0.286196434	;
miu.up("5")=	0.374194761	;
miu.up("6")=	0.425011942	;
miu.up("7")=	0.459838558	;
miu.up("8")=	0.489337979	;
miu.up("9")=	0.515386224	;
miu.up("10")=	0.538304033	;
miu.up("11")=	0.558427907	;
miu.up("12")=	0.576085216	;
miu.up("13")=	0.591580687	;
miu.up("14")=	0.605190184	;
miu.up("15")=	0.617158912	;
 
miu.lo("2")=	0.105192821	;
miu.lo("3")=	0.184429898	;
miu.lo("4")=	0.286196434	;
miu.lo("5")=	0.374194761	;
miu.lo("6")=	0.425011942	;
miu.lo("7")=	0.459838558	;
miu.lo("8")=	0.489337979	;
miu.lo("9")=	0.515386224	;
miu.lo("10")=	0.538304033	;
miu.lo("11")=	0.558427907	;
miu.lo("12")=	0.576085216	;
miu.lo("13")=	0.591580687	;
miu.lo("14")=	0.605190184	;
miu.lo("15")=	0.617158912	;

* Upper bound from Hotelling


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

Parameters
Year(t)         Date
kyoto_strong_y(t)
kyoto_strong_cpc(t)
kyoto_strong_s(t)
kyoto_strong_indem(t)
kyoto_strong_sigma(t)
kyoto_strong_tatm(t)
kyoto_strong_mat(t)
kyoto_strong_tax(t)
kyoto_strong_tax_part(t)
kyoto_strong_ri(t)
kyoto_strong_rr(t)
kyoto_strong_al(t)
kyoto_strong_forcoth(t)
kyoto_strong_l(t)
kyoto_strong_etree(t)
kyoto_strong_yy(t)
kyoto_strong_cc(t)
kyoto_strong_miu(t)
kyoto_strong_wem(t)
kyoto_strong_ri(t)
kyoto_strong_dam(t)
kyoto_strong_abate(t)
kyoto_strong_mcemis(t)
kyoto_strong_mcemis_part(t)
kyoto_strong_partfract(t)
kyoto_strong_utility   ;

Year(t)         = 2005 +10*(ord(t)-1);
kyoto_strong_y(t)=y.l(t);
kyoto_strong_cpc(t)=cpc.l(t);
kyoto_strong_s(t)=s.l(t)     ;
kyoto_strong_indem(t)= e.l(t)-etree(t);;
kyoto_strong_sigma(t)=sigma(t) ;
kyoto_strong_tatm(t)=tatm.l(t)  ;
kyoto_strong_mat(t)=mat.l(t)     ;
kyoto_strong_tax(t)=-1*ee.m(t)*1000/(kk.m(t)+.0000001)       ;
kyoto_strong_ri(t)=ri.l(t);
kyoto_strong_rr(t)=rr(t)   ;
kyoto_strong_al(t)=al(t)    ;
kyoto_strong_forcoth(t)=forcoth(t);
kyoto_strong_l(t)=l(t);
kyoto_strong_etree(t)=etree(t);
kyoto_strong_yy(t)=yy.m(t)     ;
kyoto_strong_cc(t)=cc.m(t)      ;
kyoto_strong_miu(t)=miu.l(t)     ;
kyoto_strong_wem(t)= e.l(t);
kyoto_strong_ri(t)=ri.l(t)         ;
kyoto_strong_dam(t)= damages.l(t);
kyoto_strong_abate(t) = abatecost.l(t);
kyoto_strong_mcemis(t)= expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000;
kyoto_strong_mcemis_part(t)= expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000*(partfract(t)**(1-expcost2)) ;
kyoto_strong_partfract(t)=partfract(t);
kyoto_strong_utility=utility.l        ;