[MITgcm-support] Diagnose terms in vorticity equation

[Matthew Mazloff]

Hello

Once you have closed the momentum budget you can close the vorticity budget
by taking the curl.
When we do this we use the stencil:
   vort3(I,J)=
     &      recip_rAz(I,J,bi,bj)*(
     &      ( vFld(I,J)*dyC(I,J,bi,bj)
     &       -vFld(I-1,J)*dyC(I-1,J,bi,bj) )
     &     -( uFld(I,J)*dxC(I,J,bi,bj)
     &       -uFld(I,J-1)*dxC(I,J-1,bi,bj) )
     &                           )

As this makes pxy - pyx = 0
and otherwise it doesn't always (e.g. dx is often a function of y).

-Matt



On Wed, Jan 25, 2023 at 6:36 AM Yilang Xu <yxu at whoi.edu> wrote:

> Dear MITgcm community,
>
>
>
> I am wondering if anyone has diagnosed the terms in vorticity equation in
> MITgcm. I see that there are diagnostics for momentum budget, and I am
> curious if there are similar things for vorticity equation.
>
> Thank you very much.
>
>
>
> Best,
>
> Yilang
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
>
> https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!Hn9HembHHugP5pPzonRkldVlrpJux6OKounalValIwp9nJMXH2HsglZDt7FRQ9sfBvCc2LrSEldRTw$
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20230125/a557a733/attachment.html>