statistical equilibrium, spin up and total kinetic energy

Discussion about analysis, visualization, and collaboration tools and techniques

Moderators: arango, robertson

Post Reply
Message
Author
jul10lara

statistical equilibrium, spin up and total kinetic energy

#1 Unread post by 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!!!!!

:P

User avatar
kate
Posts: 4091
Joined: Wed Jul 02, 2003 5:29 pm
Location: CFOS/UAF, USA

Re: statistical equilibrium, spin up and total kinetic energy

#2 Unread post by 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

Re: statistical equilibrium, spin up and total kinetic energy

#3 Unread post by 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!!!

Post Reply