hi , I would like to create a model about climate economics.and I wrote codes as below.But they are not working properly and system gives error as below.
Appreciate if someone can help me. thanks
var L A k c q omega Tat;
varexo e;
parameters alfa ro gamma sigmak w gl ga psi sigmae ;
alfa=0.45;
ro=0.015;
gamma=0.33;
sigmak=0.1;
sigmae=0.01;
w=1.08;
psi=0.003;
gl=0.134;
ga=0.016;
model;
L*c^(-alfa)=((1/(1+ro))*c(+1)^(-alfa)*((gamma*A(+1)*k^(gamma-1))/(1+omega))-sigmak);
q=(A*k^gamma)/(1+omega);
k=((A)*k(-1)^(gamma))/(1+omega)-c-sigmak*k(-1);
L=L(-1)*(1+gl);
A=A(-1)*(1+ga);
omega=psi*exp(Tat)+psi*exp(Tat)^2;
Tat=w*Tat(-1)+e;
end;
initval;
L=6838;
k=19.74;
q=63.69;
c=44.34;
A=3.8;
omega=0.9;
end;
shocks;
var e=sigmae^2;
end;
steady;
stoch_simul(hp_filter=1600,order=1,irf=40);
system gives error as follow;
Residuals of the static equations:
Equation number 1 : 0.1
Equation number 2 : 63.69
Equation number 3 : 0
Equation number 4 : 0
Equation number 5 : 0
Equation number 6 : 0
Equation number 7 : 0
Error using print_info (line 74)
Impossible to find the steady state. Either the model doesn't have a steady
state, there are an infinity of steady states, or the guess values are too
far from the solution
Error in steady (line 92)
print_info(info,options_.noprint, options_);
Error in climate (line 157)
steady;
Error in dynare (line 180)
evalin('base',fname) ;