[MITgcm-support] Diagnose terms in vorticity equation
Yilang Xu
yxu at whoi.edu
Wed Feb 1 12:01:51 EST 2023
Hi Matt,
Thanks for the explanation. I track down the code you list to mom_calc_relvort3.F. I guess the terms in the vorticity equation can be coded for diagnostic purpose in a similar way. Or, the vorticity budget can be calculated post-run, since the velocities, viscosities, momHDiv, and momVort3 can be diagnosed.
Cheers,
Yilang
From: MITgcm-support <mitgcm-support-bounces at mitgcm.org> on behalf of Matthew Mazloff <mmazloff at ucsd.edu>
Reply-To: "mitgcm-support at mitgcm.org" <mitgcm-support at mitgcm.org>
Date: Wednesday, January 25, 2023 at 14:27
To: "mitgcm-support at mitgcm.org" <mitgcm-support at mitgcm.org>
Cc: Yilang Xu <yxu at whoi.edu>
Subject: Re: [MITgcm-support] Diagnose terms in vorticity equation
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<mailto: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<mailto:MITgcm-support at mitgcm.org>
https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!Mih3wA!Hn9HembHHugP5pPzonRkldVlrpJux6OKounalValIwp9nJMXH2HsglZDt7FRQ9sfBvCc2LrSEldRTw$<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/20230201/da2a83db/attachment.html>
More information about the MITgcm-support
mailing list