[MITgcm-support] Diagnose terms in vorticity equation

[Matthew Mazloff]

Yes, I always do it post-run using all the momentum term diagnostics. But I suppose you could save disk space by introducing these diagnostics into the code. 

Matt


> On Feb 1, 2023, at 9:01 AM, Yilang Xu <yxu at whoi.edu> wrote:
> 
> 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$>_______________________________________________
> 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!Glete6cuQ3CyNPRNM2B0oAoFh9fTRU6gQ3pAQipz6P3wxcFuNEWHUPl5S7ByUz_pQjWw3npOJ0KWXg$

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20230201/49142139/attachment-0001.html>