* Definition of Gore 
 
PARTFRACT(T)=1;
 
partfract("1")=	0.25372	;
PARTFRACT("2")=.6;
PARTFRACT("3")=.7;
PARTFRACT("4")=.8;
PARTFRACT("5")=.9;
 
miu.up(T)=1;
miu.lo(T)= .95; 	

miu.fx("1")=.005	;
miu.lo("2")=.25	;
miu.lo("3")=.45	;
miu.lo("4")=.65	;
miu.lo("5")=.85	;
miu.lo("6")=.95	;
miu.lo("7")=.95	;
miu.lo("8")=.95	;


mat.up(t)= 3000;


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

Parameters
Year(t)         Date
gore_y(t)
gore_cpc(t)
gore_s(t)
gore_indem(t)
gore_sigma(t)
gore_tatm(t)
gore_mat(t)
gore_tax(t)
gore_tax_part(t)
gore_ri(t)
gore_rr(t)
gore_al(t)
gore_forcoth(t)
gore_l(t)
gore_etree(t)
gore_yy(t)
gore_cc(t)
gore_miu(t)
gore_wem(t)
gore_ri(t)
gore_dam(t)
gore_abate(t)
gore_mcemis(t)
gore_mcemis_part(t)
gore_partfract(t)
gore_utility   ;

Year(t)         = 2005 +10*(ord(t)-1);
gore_y(t)=y.l(t);
gore_cpc(t)=cpc.l(t);
gore_s(t)=s.l(t)     ;
gore_indem(t)= e.l(t)-etree(t);;
gore_sigma(t)=sigma(t) ;
gore_tatm(t)=tatm.l(t)  ;
gore_mat(t)=mat.l(t)     ;
gore_tax(t)=-1*ee.m(t)*1000/(kk.m(t)+.0000001)       ;
gore_ri(t)=ri.l(t);
gore_rr(t)=rr(t)   ;
gore_al(t)=al(t)    ;
gore_forcoth(t)=forcoth(t);
gore_l(t)=l(t);
gore_etree(t)=etree(t);
gore_yy(t)=yy.m(t)     ;
gore_cc(t)=cc.m(t)      ;
gore_miu(t)=miu.l(t)     ;
gore_wem(t)= e.l(t);
gore_ri(t)=ri.l(t)         ;
gore_dam(t)= damages.l(t);
gore_abate(t) = abatecost.l(t);
gore_mcemis(t)= expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000;
gore_mcemis_part(t)= expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000*(partfract(t)**(1-expcost2)) ;
gore_partfract(t)=partfract(t);
gore_utility=utility.l        ;