[MITgcm-support] Energy diagnostics?

[Martin Losch]

Hi Jody,

I am not aware of energy flux diagnostics [(E+P)u] in the code. I  
agree with you that for closing the budget they need to be computed  
correctly (at each time step and then averaged), horizontal averaging  
will also be an issue.

Martin

On Feb 15, 2009, at 10:25 PM, Klymak Jody wrote:

>
> Hi all,
>
> Before I spend a couple of days doing this, has anyone gone through  
> the effort of calculating energy diagnostics?
>
> I'm running barotropic flow over an obstacle (all in 2-D).  If I  
> take every timestep I can get a decent energy budget:
>
> (d/dt)\int E dV = \oint (E+P) u \cdot dA + \int \epsilon dV
>
> where \epsilon is specified and numerical dissipation.  However (E 
> +P)u changes so much over my output timesteps the residual when  
> compared w/ dE/dt is unrealistically large.      Averaging the  
> linear and non-linear flux over the dump interval should give a much  
> better estimate.
>
> Has anyone done this?  If not, I'll share my code.
>
> Thanks,  Jody
>
> --
> Jody Klymak
> http://web.uvic.ca/~jklymak/
>
>
>
>
>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support