* Definition of limited participation results

 PARTFRACT1 = .15;  
 PARTFRACT2 = .2;  
 PARTFRACT21 =.5; 
partfract(t) = partfract21;
PARTFRACT(T)$(ord(T)<25) = Partfract21 + (PARTFRACT2-Partfract21)*exp(-DPARTFRACT*(ORD(T)-2));
partfract("1")= PARTFRACT1;

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

Parameters
Year(t)         Date
limpart_y(t)
limpart_cpc(t)
limpart_s(t)
limpart_indem(t)
limpart_sigma(t)
limpart_tatm(t)
limpart_mat(t)
limpart_tax(t)
limpart_ri(t)
limpart_rr(t)
limpart_al(t)
limpart_forcoth(t)
limpart_l(t)
limpart_etree(t)
limpart_yy(t)
limpart_cc(t)
limpart_miu(t)
limpart_wem(t)
limpart_ri(t)
limpart_dam(t)
limpart_abate(t)
limpart_mcemis(t)
limpart_partfract(t)
limpart_utility   ;

Year(t)         = 2005 +10*(ord(t)-1);
limpart_y(t)=y.l(t);
limpart_cpc(t)=cpc.l(t);
limpart_s(t)=s.l(t)     ;
limpart_indem(t)= e.l(t)-etree(t);;
limpart_sigma(t)=sigma(t) ;
limpart_tatm(t)=tatm.l(t)  ;
limpart_mat(t)=mat.l(t)     ;
limpart_tax(t)=-1*ee.m(t)*1000/(kk.m(t)+.0000001)       ;
limpart_ri(t)=ri.l(t);
limpart_rr(t)=rr(t)   ;
limpart_al(t)=al(t)    ;
limpart_forcoth(t)=forcoth(t);
limpart_l(t)=l(t);
limpart_etree(t)=etree(t);
limpart_yy(t)=yy.m(t)     ;
limpart_cc(t)=cc.m(t)      ;
limpart_miu(t)=miu.l(t)     ;
limpart_wem(t)= e.l(t);
limpart_ri(t)=ri.l(t)         ;
limpart_dam(t)= damages.l(t);
limpart_abate(t) = abatecost.l(t);
limpart_partfract(t)=partfract(t);
limpart_utility=utility.l        ;
limpart_mcemis(t)= expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000;