[MITgcm-support] energy diagnostics

[Jeff Scott]

Hi Martin,

Several years ago Alistair and I spent a few months checking the KE 
balance in the MITGCM, albeit offline (dumping everything every time step, 
running sychronously). It took a few months because it was rather a pain 
in the neck, mostly because of Adams-Bashforth.

We never did any inline calculations - but I did however document the 
calculations and I could send you this if you'd like.

Jeff






On Tue, 1 Apr 2008, Martin Losch wrote:

> Hi,
>
> I am sure that this has been discussed on this list before, but because the 
> search function on http://mitgcm.org/search.html still does not work for the 
> support mailing list I cannot find anything about energy diagnostics.
> I want to diagnose energy production (u*grad(p)), redistribution ((u*grad) 
> u), and dissipation (u*dissipation, that is horizontal friction and bottom 
> and side drag) terms in the horizontal momentum equations. As far as i can 
> see, there are no built-in diagnostics for this yet.
> Is there anybody who has done this and is willing to share her/his code?
> Otherwise my plan was to use gU/V and gU/Vdiss (and mulitply by U/V) and save 
> just before calling timestep (in dynamics.F), but is there an explicit 
> tendency due to grad(p)? Do I have to save this contribution in 
> momentum_correction_step.F to be consistent with the pressure correction 
> method?
>
> Martin
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support