Dear GAMS experts,
My model has some problem.
I attached my code here, if you have any idea, then please let me know.
Thank you very much in advance.
Sincerley yours,
Redwood
SETS T TIME PERIODS /05/
TFIRST(T) FIRST PERIOD
TLAST(T) LAST PERIOD
N Regions /U, C/;
SCALARS
** constant
ep Exchange rate / 1.0 /
Gbar1 Limits of carbon emissions in country 1 / 4880 /
r1 Interest rate in country 1 / 0.03 /
L1 Labor force in country 1 / 157.5 /
Gbar2 Limits of carbon emissions in country 2 / 3410 /
r2 Interest rate in country 2 / 0.05 /
L2 Labor force in country 2 / 799.5 /
Pw World price of energy / 1 /
GHGs Target of greenhouse gas emissions / 2748390 /
delta1 Discount rate / 0.02 /
gam1 Coefficient of utility 1 / 0.8 /
eta11 Elasticity of MGS in country 1 / 0.609 /
eta12 Elasticity of EQ in country 1 / 0.1644 /
ag1 Clean capital in country 1 / 1.675 /
alpha1 Target green house gas emissions / 15.9341253 /
rhoq1 Target green house gas emissions / 0.01974 /
beta11 Elasticity of capital in county 1 / 0.25 /
beta12 Elasticity of labor force in country 1 / 0.60/
beta13 Elasticity of energy in country 1 / 0.15 /
g11 Emissions coefficient in country 1 / 2.352 /
g12 Scale factor of emissions in country 1 / -0.1869 /
xi1 Target green house gas emissions / 0.8898 /
rhom1 Target green house gas emissions / 0.004 /
delta2 Discount rate / 0.04 /
gam2 Coefficient of utility 2 / 0.6509 /
eta21 Elasticity of MGS in country 2 / 0.5 /
eta22 Elasticity of EQ in country 2 / 0.09488 /
ag2 Clean capital in country 2 / -1.239 /
alpha2 Target green house gas emissions / 3.257025325 /
rhoq2 Target green house gas emissions / 0.023270 /
beta21 Elasticity of capital in country 2 / 0.217928 /
beta22 Elasticity of labor force in country 2 / 0.667832 /
beta23 Elasticity of energy in country 2 / 0.114240 /
g21 Emissions coefficient in country 2 / 0.5398 /
g22 Scale factor of emissions in country 2 / -2.28 /
xi2 Target green house gas emissions / 0.15 /
rhom2 Target green house gas emissions / 0.0453 /
d Greenhouse gas emissions discout facot / 0.003 /
psi1 Elasticity of offsets price / 0.01 /
Ku0 Initial dirty capital stock in country 1 / 34867.22787 /
Kc0 Initial dirty capital stock in country 2 / 19511.9685 /
Kgu0 Initial clean capital stock in country 1 / 34.90213 /
Kgc0 Initial clean capital stock in country 2 / 19.5315 /
Pg0 Initial price of offset / 5 /
GHG0 Initial greenhouse gas emissions / 3026090 /
Cu0 Initial consumtion in country 1 / 12.8 /
Cc0 Initial consumtion in country 2 /2.84/
Igu0 Initial clean investment in country 1 / 0.0000214 /
Igc0 Initial clean investment in country 2 / 0.002712081 /
Eu0 Energy consumption in country 1 /2220/
Ec0 Energy consumption in country 2 /2420/
Qu0 GDP in country 1 /14400/
Qc0 GDP in country 2 /5930/
mu1 Scale factor of consumption in country 1 /0.8562/
mu2 Scale factor of consumption in country 2 /0.4911/
MGSu0 Initial imports in country 1 /2660/
MGSc0 Initial imports in country 2 /1580/
PARAMETERS ZU(T) DISCOUNT FACTOR
ZC(T) DISCOUNT FACTOR
FU(T)
FC(T);
TFIRST(T) = YES$(ORD(T) EQ 1);
TLAST(T) = YES$(ORD(T) EQ CARD(T));
ZU(T) = EXP(rhoq1(ORD(T)-1));
ZC(T) = EXP(rhoq2(ORD(T)-1));
FU(T) = EXP(-rhom1(ORD(T)-1));
FC(T) = EXP(-rhom2*(ORD(T)-1));
Parameter BigM/150000/;
VARIABLES
-
Control variables
*C(T,N) Consumption
MGS(T,N) Imports
Ig(T,N) Clean investment
E(T,N) Energy -
State variables
K(T,N) Dirty capital stocks
Kg(T,N) Clean capital stocks
Pg(T) Price of offsets
GHG(T) The concentration of carbon dioxide (CO2) in Earth’s atmosphere -
Other variables
Op(T) Offsets
Q(T,N) GDP(Output)
Utility1
Utility2
Dummy(T);
Positive Variable
- Costate varialbes
lamd(T,N) Costate variable for K
qg(T,N) Costate variable for Kg
qp(T,N) Costate variable for Pg
qe(T,N) Costate variable for GHG
;
*Binary Variable
*Y(T);
EQUATIONS
*EEu(T,N)
*EEc(T,N)
*EEEu(T,N)
*EEEc(T,N)
QQu(T,N) Output in country 1
QQc(T,N) Output in country 2
OOp1(T) Offsets 1
OOp2(T) Offsets 2
*OOp3(T) Offsets 3
*OOp4(T) Offsets 4
*OOp5(T)
dHK11(T,N)
dHK12(T,N)
dHKg11(T,N)
dHKg12(T,N)
dHPg1(T)
dHGHG1(T)
dHIg1(T,N)
dHMGS1(T,N)
dHE11(T,N)
dHE12(T,N)
####### FOC for country 2 #########
dHK21(T,N)
dHK22(T,N)
dHKg21(T,N)
dHKg22(T,N)
dHPg2(T)
dHGHG2(T)
dHIg2(T,N)
dHMGS2(T,N)
dHE21(T,N)
dHE22(T,N)
Ku(T,N) Dirty capital stocks in country 1
TKu(T,N) Termial condition for dirty capital stocks in country 1
Kc(T,N) Dirty capital stocks in country 2
TKc(T,N) Termial condition for dirty capital stocks in country 2
KKg(T,N) Clean capital stocks
TKg(T,N) Terminal condition for clean capital stocks
PPg(T) Price of offsets
TPg(T) Terminal condition for price of offsets
GGHG(T) The concentration of carbon dioxide (CO2) in Earth’s atmosphere
TGHG(T) Terminal condition for the concentration of CO2
;
EEu(T,N)… E(T,“U”) =e= beta13Q(T,“U”)/(Pw+Pg(T)Op(T)/E(T, “U”));
EEc(T,N)… E(T,“C”) =e= beta23Q(T,“C”)/Pw;
EEEu(T,N)… beta13Q(T,“U”)/(Pw+Pg(T)Op(T)/E(T, “U”)) =l= (Gbar1 + Op(T))/(g11FU(T)-g12(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1);
EEEc(T,N)… beta23Q(T,“C”)/Pw =l= (Gbar2)/(g21FC(T)-g22*(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))xi2);
QQu(T,N)… Q(T,“U”) =e= ZU(T)alpha1((K(T, “U”)+ Kg(T, “U”))beta11)*(L1beta12)(E(T, “U”)**beta13);
QQc(T,N)… Q(T,“C”) =e= ZC(T)alpha2((K(T, “C”)+ Kg(T, “C”))**beta21)(L2beta22)(E(T, “C”)**beta23);
OOp1(T)… Op(T) =e= (g11FU(T)-g12*(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))*xi1)E(T, “U”) - Gbar1$(Gbar2 - (g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2)E(T, “C”) + Gbar1 - (g11FU(T)-g12(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))*xi1)E(T, “U”) >= 0);
OOp2(T)… Op(T) =e= Gbar2 - (g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2)E(T, “C”)$(Gbar2 - (g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2)E(T, “C”) + Gbar1 - (g11FU(T)-g12(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1)*E(T, “U”) = 0);
dHK12(T,N)… - lamd(T+1, “U”) + (1+delta1)lamd(T, “U”) =e= ((mu1Q(T,“U”)*MGS(T, “U”)*eta11(GHGs/GHG(T))**eta12)**gam1)*beta11/(K(T, “U”)+ Kg(T, “U”))
- lamd(T,“U”)((1-mu1)beta11Q(T,“U”)/(K(T, “U”)+ Kg(T, “U”))- xi1g12*Pg(T)E(T,“U”)(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1/(K(T,“U”)+ Kg(T,“U”)))
- qp(T,“U”)psi1Pg(T)(g11FU(T)beta13**2Q(T,“U”)/((Pw + Pg(T)Op(T)/E(T,“U”))(K(T,“U”)+ Kg(T,“U”)))
- g12beta13Q(T,“U”)(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1)(g11 - beta11)/(Pw + Pg(T)Op(T)/E(T,“U”))(K(T,“U”)+ Kg(T,“U”))/sqrt(Gbar1*Gbar2)
- qe(T,“U”)xi1g12E(T,“U”)(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1/(K(T,“U”)+ Kg(T,“U”))
$(Gbar2 - (g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2)E(T, “C”) + Gbar1 - (g11FU(T)-g12(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1)*E(T, “U”) = 0);
dHKg12(T,N)… - qg(T+1, “U”) + (1+delta1)qg(T, “U”) =e= ((mu1Q(T,“U”)*MGS(T, “U”)*eta11(GHGs/GHG(T))**eta12)**gam1)*beta11/(K(T, “U”)+ Kg(T, “U”)) - lamd(T,“U”)((1-mu1)beta11Q(T,“U”)/(K(T, “U”)+ Kg(T, “U”)) + ag1(Ig(T,“U”)/Kg(T,“U”))**2/2)
- qp(T,“U”)psi1Pg(T)(g11FU(T)beta13**2Q(T,“U”)/((Pw + Pg(T)Op(T)/E(T,“U”))(K(T,“U”)+ Kg(T,“U”)))
- g12beta13Q(T,“U”)(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))*xi1(xi1K(T,“U”)/Kg(T,“U”) + beta11Q(T,“U”))/((Pw + Pg(T)Op(T)/E(T,“U”))(K(T,“U”)+Kg(T,“U”))))/sqrt(Gbar1Gbar2)
- qe(T,“U”)xi1g12E(T,“U”)K(T,“U”)(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi2/(K(T,“U”)+ Kg(T,“U”))
$(Gbar2 - (g21FC(T)-g22*(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2)E(T, “C”) + Gbar1 - (g11FU(T)-g12(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1)E(T, “U”) = 0);
dHE12(T,N)… beta13(mu1*Q(T,“U”)*MGS(T, “U”)*eta11(GHGs/GHG(T))**eta12)**gam1/E(T, “U”) =e= - lamd(T,“U”)*(1-mu1)beta13Q(T,“U”)/E(T, “U”)
- qp(T,“U”)psi1Pg(T)(1/sqrt(Gbar1Gbar2))(g11FU(T)-g12*(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1)beta13(beta13*Q(T,“U”)/E(T,“U”) + Q(T,“U”)*Pg(T)*Op(T)/((Pw + Pg(T)*Op(T)/E(T,“U”))*E(T,“U”)**2))/(Pw + Pg(T)*Op(T)/E(T,“U”))
- qe(T,“C”)(g11FU(T)-g12*(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1)
$(Gbar2 - (g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2)E(T, “C”) + Gbar1 - (g11FU(T)-g12(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1)*E(T, “U”) = 0);
dHK22(T,N)… - lamd(T+1, “C”) + (1-delta2)lamd(T, “C”) =e= ((mu2Q(T,“C”)*MGS(T, “C”)*eta21(GHGs/GHG(T))**eta22)**gam2)*beta21/(K(T, “C”)+ Kg(T, “C”))
- lamd(T,“C”)*((1-mu2)beta21Q(T,“C”)/(K(T, “C”)+ Kg(T, “C”)))
- qp(T,“C”)psi1Pg(T)beta23Q(T,“C”)(beta21g21FC(T) + (xi2 - beta21)g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))**xi2)/((K(T,“C”)+ Kg(T,“C”))Pwsqrt(Gbar1Gbar2))
- qe(T,“C”)xi2g22E(T,“C”)(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))**xi2/(K(T,“C”)+ Kg(T,“C”))
$(Gbar2 - (g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2)E(T, “C”) + Gbar1 - (g11FU(T)-g12(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1)*E(T, “U”) = 0);
dHKg22(T,N)… - qg(T+1, “C”) + (1-delta2)qg(T, “C”) =e= ((mu2Q(T,“C”)*MGS(T, “C”)*eta21(GHGs/GHG(T))**eta22)**gam2)*beta21/(K(T, “C”)+ Kg(T, “C”)) - lamd(T,“C”)((1-mu2)beta21Q(T,“C”)/(K(T, “C”)+ Kg(T, “C”)) + ag2(Ig(T,“C”)/Kg(T,“C”))**2/2)
- qp(T,“C”)psi1Pg(T)beta23Q(T,“C”)(beta21g21FC(T) - g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2(xi2K(T,“C”) + beta21Kg(T,“C”))/Kg(T,“C”))/((K(T,“C”)+ Kg(T,“C”))Pwsqrt(Gbar1*Gbar2))
- qe(T,“C”)xi2g22*E(T,“C”)K(T,“C”)(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2/((K(T,“C”)+ Kg(T,“C”))Kg(T,“C”))
$(Gbar2 - (g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2)E(T, “C”) + Gbar1 - (g11FU(T)-g12(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1)*E(T, “U”) = 0);
dHE22(T,N)… beta23*(mu2Q(T,“C”)MGS(T, “C”)*eta21(GHGs/GHG(T))**eta22)**gam2/E(T, “C”) =e=
-lamd(T,“C”)((1-mu2)beta23Q(T,“C”)/E(T, “C”))
-qp(T,“C”)psi1Pg(T)beta23**2Q(T,“C”)(g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))**xi2)/(PwE(T,“C”)sqrt(Gbar1Gbar2))
-qe(T,“C”)(g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))**xi2)
$(Gbar2 - (g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2)E(T, “C”) + Gbar1 - (g11FU(T)-g12(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))xi1)*E(T, “U”) < 0);
Ku(T+1,N)… K(T+1, “U”) =e= (1-mu1)ZU(T)alpha1((K(T, “U”)+ Kg(T, “U”))**beta11)(L1beta12)*(E(T, “U”)**beta13)
- MGS(T,“U”)ep - Ig(T, “U”)(1+ag1Ig(T, “U”)/(2Kg(T, “U”)))
- Pg(T)Op(T)+ K(T, “U”);
TKu(TLAST,N)… 0 =l= (1-mu1)ZU(TLAST)alpha1((K(TLAST, “U”)+ Kg(TLAST, “U”))**beta11)(L1**beta12)(E(TLAST, “U”)**beta13) - MGS(TLAST,“U”)ep - Ig(TLAST, “U”)(1+ag1Ig(TLAST, “U”)/(2Kg(TLAST, “U”)))
- Pg(TLAST)Op(TLAST)+ K(TLAST, “U”);
Kc(T+1,N)… K(T+1, “C”) =e= (1-mu2)ZC(T)alpha2((K(T, “C”)+ Kg(T, “C”))**beta21)(L2**beta22)(E(T, “C”)**beta23) - MGS(T,“C”)/ep - Ig(T, “C”)(1+ag2Ig(T, “C”)/(2*Kg(T, “C”)))
- Pg(T)Op(T)/ep + K(T, “C”);
TKc(TLAST,N)… 0 =l= (1-mu2)ZC(TLAST)alpha2((K(TLAST, “C”)+ Kg(TLAST, “C”))**beta21)(L2**beta22)(E(TLAST, “C”)**beta23)
- MGS(TLAST,“C”)ep - Ig(TLAST, “C”)(1+ag2Ig(TLAST, “C”)/(2Kg(TLAST, “C”)))
- Pg(TLAST)Op(TLAST)/ep + K(TLAST, “C”);
KKg(T+1,N)… Kg(T+1,N) =e= Ig(T,N)+ Kg(T,N);
TKg(TLAST,N)… 0 =l= Ig(TLAST,N) + Kg(TLAST,N);
PPg(T+1)… Pg(T+1) =e= psi1Pg(T)
((g11FU(T)- g12*(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1)beta13Q(T,“U”)/(Pw+Pg(T)*Op(T)/E(T, “U”)) - Gbar1
- Gbar2 + (g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))*xi2)beta23Q(T,“C”)/Pw)
/sqrt(Gbar1Gbar2) + Pg(T);
TPg(TLAST)… 0 =l= psi1Pg(TLAST)
((g11FU(TLAST)-g12(Kg(TLAST,“U”)/(K(TLAST,“U”)+ Kg(TLAST,“U”)))**xi1)beta13Q(TLAST,“U”)/(Pw+Pg(TLAST)*Op(TLAST)/E(TLAST, “U”))
- (g21FC(TLAST)-g22(Kg(TLAST,“C”)/(K(TLAST,“C”)+ Kg(TLAST,“C”)))**xi2)beta23Q(TLAST,“C”)/Pw - Gbar1 - Gbar2)
/sqrt(Gbar1Gbar2) + Pg(TLAST);
GGHG(T+1)… GHG(T+1) =e= - dGHG(T) - (g11FU(T)-g12(Kg(T,“U”)/(K(T,“U”)+ Kg(T,“U”)))**xi1)*E(T, “U”)
- (g21FC(T)-g22(Kg(T,“C”)/(K(T,“C”)+ Kg(T,“C”)))**xi2)*E(T, “C”)
- GHG(T);
TGHG(TLAST)… 0 =l= - d*GHG(TLAST) - (g11FU(TLAST)-g12(Kg(TLAST,“U”)/(K(TLAST,“U”)+ Kg(TLAST,“U”)))**xi1)*E(TLAST, “U”)
- (g21FC(TLAST)-g22(Kg(TLAST,“C”)/(K(TLAST,“C”)+ Kg(TLAST,“C”)))**xi2)*E(TLAST, “C”)
- GHG(TLAST);
MODEL Nash /all/;
K.LO(T,“U”) = 13150;
K.LO(T,“C”) = 19300;
Kg.LO(T,N) = 1;
Pg.LO(T) = 0.1; Pg.up(T) = 100;
GHG.lo(T) = 2000000; GHG.up(T) = 3500000;
*C.LO(T,N)= 0.01; C.UP(T,N) = 150;
MGS.LO(T,N)= 500; MGS.UP(T,N) = 10000;
Ig.LO(T,N)= 10; Ig.UP(T,N) = 10000;
E.LO(T,N)= 2000; E.up(T,N)= 2500;
Op.lO(T) = 0; Op.up(T) = 10000;
K.fx(TFIRST, “U”) = Ku0;
K.fx(TFIRST, “C”) = Kc0;
Kg.fx(TFIRST, “U”) = Kgu0;
Kg.fx(TFIRST, “C”) = Kgc0;
Pg.fx(TFIRST) = Pg0;
GHG.fx(TFIRST) = GHG0;
*C.fx(TFIRST, “U”) = Cu0;
*C.fx(TFIRST, “C”) = Cc0;
MGS.fx(TFIRST, “U”) = MGSu0;
MGS.fx(TFIRST, “C”) = MGSc0;
Ig.fx(TFIRST, “U”) = Igu0;
Ig.fx(TFIRST, “C”) = Igc0;
E.fx(TFIRST, “U”) = Eu0;
E.fx(TFIRST, “C”) = Ec0;
*Q.fx(TFIRST, “U”) = Qu0;
*Q.fx(TFIRST, “C”) = Qc0;
lamd.fx(TLAST,N) = 0;
qg.fx(TLAST,N) = 0;
qp.fx(TLAST,N) = 0;
qe.fx(TLAST,N) = 0;
SOLVE Nash USING MCP;
–
To unsubscribe from this group and stop receiving emails from it, send an email to gamsworld+unsubscribe@googlegroups.com.
To post to this group, send email to gamsworld@googlegroups.com.
Visit this group at http://groups.google.com/group/gamsworld.
For more options, visit https://groups.google.com/groups/opt_out.