Maybe check your equations carefully again. For example, under
WS = (1-alpha)*(mc*alpha^alpha/Rk^alpha)^(1/(1-alpha)) = 0.245, But in main model,
WS=epsw/(epsw-1)*H1/H2 = 0.271;
That is why you have, for example,
Equation number 3 : -0.026 (i.e., residual not zero) after running your mod file. Thus, 0.245 - 0.271 = -0.026.
Also, for example, in steady-state (initval block)
W = WS = 0.24579. But if you use steady values for model equations, then
W=(1-alpha)*mc*A*K^alpha*N^(-alpha) = 0.27384; and
WS=epsw/(epsw-1)*H1/H2 = 0.27188;, so W \neq WS in
model; block, but W = WS in
These are some things you can check. Maybe this is helpful.