[MITgcm-support] vorticity balance

[Matthew Mazloff]

Hi Marco

Taking curl of the momentum equation is the proper method. Due to discretization the vorticity budget may not close to numerical precision, but being sure you have closed the momentum budget first will give you confidence you are doing things correctly. 

And note that if you are using spherical coordinates the discretization will cause the pressure term to not vanish, so you need to account for that

Matt


On Jan 29, 2015, at 6:34 AM, marco reale <reale.marco82 at gmail.com> wrote:

> hi everyone,
> 
> I’m just curios to know whether my approach is correct or not : I want to do a vorticity balance of a basin. In order to compute each term of vorticity equation my idea is to apply the curl definition as mom_cal_relVort.F to each term of momentum buget equation , namely UM_advec+UmHydroPX+Um_diss+GradientEta+Um_corios+UM_ext.After I get each term  in order to have a vorticity not a variation I multiply each term for frequency of shapshot set in data.diagnostics.Have you ever done this kind of computation using MITGCM ? Is it correct this approach or not?
> 
> cheers
> 
> Marco
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support