I am new to Dynare and am trying to simulate an extension of the Gertler-Kiyotaki model (so many of the equations are the same as or similar to ones in the original paper). All of the computations go through and the steady state is found but many of the moments appear to have the wrong sign and no plots are being displayed (in fact, in the panel on the right side all of the variables say NaN in all periods). Also note that none of the variables are in exponential form (so the equations in the code are precisely those derived). I’ve attached my .mod file.

Thanks, the order=1 suggestion worked. When you say “shock size of 100%,” you’re referring to the fact that I set var e =1? If I change it to .01 then it becomes 1%? Also, suppose I wanted to plot some function of variables I defined (let’s say S/N); then can I simply type something like “SoverN = S/N; rplot S/N;” ?