* Solution for geoengineering

PBACK =  .005 ;
partfract(T)=1;
miu.up(T)=1;
parameters
cost1(t);
cost1(T) = (PBACK*SIGMA(T)/EXPCOST2)* ( (BACKRAT-1+ EXP (-gback* (ORD(T)-1) ) )/BACKRAT);   

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

* Definition of geoeng results

Parameters
Year(t)         Date
geoeng_y(t)
geoeng_cpc(t)
geoeng_s(t)
geoeng_indem(t)
geoeng_sigma(t)
geoeng_tatm(t)
geoeng_mat(t)
geoeng_tax(t)
geoeng_ri(t)
geoeng_rr(t)
geoeng_al(t)
geoeng_forcoth(t)
geoeng_l(t)
geoeng_etree(t)
geoeng_yy(t)
geoeng_cc(t)
geoeng_miu(t)
geoeng_wem(t)
geoeng_ri(t)
geoeng_dam(t)
geoeng_abate(t)
geoeng_mcemis(t)
geoeng_utility   ;

Year(t)         = 2005 +10*(ord(t)-1);
geoeng_y(t)=y.l(t);
geoeng_cpc(t)=cpc.l(t);
geoeng_s(t)=s.l(t)     ;
geoeng_indem(t)= e.l(t)-etree(t);;
geoeng_sigma(t)=sigma(t) ;
geoeng_tatm(t)=tatm.l(t)  ;
geoeng_mat(t)=mat.l(t)     ;
geoeng_tax(t)=-1*ee.m(t)*1000/(kk.m(t)+.00000000001)       ;
geoeng_ri(t)=ri.l(t);
geoeng_rr(t)=rr(t)   ;
geoeng_al(t)=al(t)    ;
geoeng_forcoth(t)=forcoth(t);
geoeng_l(t)=l(t);
geoeng_etree(t)=etree(t);
geoeng_yy(t)=yy.m(t)     ;
geoeng_cc(t)=cc.m(t)      ;
geoeng_miu(t)=miu.l(t)     ;
geoeng_wem(t)= e.l(t);
geoeng_ri(t)=ri.l(t)         ;
geoeng_dam(t)= damages.l(t);
geoeng_abate(t) = abatecost.l(t);
geoeng_mcemis(t)= expcost2*cost1(t)*miu.l(t)**(expcost2-1)/sigma(t)*1000;
geoeng_utility=utility.l        ;