statistical equilibrium, spin up and total kinetic energy

[jul10lara]

Hi every body

I have run a 10 perpetual years simulation because
I want to get statistical equilibrium for the model,
but I don't know if the model has get it or not.
So the question is: how do I know if the simulation has reach statistical equilibrium?
I heard the steady state is reached when the change in total kinetic energy
becomes minimum, and this can be plotted.
So the other question is: How do I know the total (or mean) kinetic energy of the
simulation, and how can I plot this?

Thank you!!!!!

[kate]

When you run ROMS, there is much written to standard out, including TKE every NINFO timesteps. I set NINFO to 1 and get vast quantities of output. How to plot it? Well, pick a standard tool. I don't know if the guys have developed this plot in Matlab. I'd do it in NCL or ask my programmer to do it in Python, but you could use gnuplot if that was your thing.

[jul10lara]

I forgot to say I use ROMS_tools xD
My problem is solved.
I just have to run the matlab scripts roms_diags and plot_diags. With this scrypts I can see plotted the kinetic energy. If it reaches the asymptote, or fluctuates around its average, it means the model has reached the statistical equilibrium.

Thank you!!!