[MITgcm-support] application of Gauss (or Stokes) theorem to vorticity eon on C-grid

[钱钰坤]

Hi Paolo,


I am also going to do this but a little worried about the staggerred grids.  I guess 'no_slip_sides=.TRUE.' (in parm01 of data file) is also important to ensure zero circulation at boundaries?



  Best regards 
 
-----------------------------------------
Yu-Kun Qian (钱钰坤)South China Sea Institute of Oceanology
Sun Yat-sen University
Xingang West Road, Haizhu District
Guangzhou, P.R. China
Email: qianyukun at scsio.ac.cn, qianyk at mail3.sysu.edu.cn
   
 

 
 
 
------------------ Original ------------------
From:  "Paola Cessi"<pcessi at ucsd.edu>;
Date:  Tue, Sep 5, 2023 08:04 AM
To:  "MITgcm-support"<mitgcm-support at mitgcm.org>; 

Subject:  [MITgcm-support] application of Gauss (or Stokes) theorem to vorticity eon on C-grid

 

Dear colleagues,

I am trying to calculate the circulation balance over a semi-enclosed area, by taking the vorticity balance and integrating over said area. One of the terms in the vorticity balance is div (fu, fv), which, when integrated over the area, should give the integral over the boundary of (fu,fv) \cdot n, where n is the normal to the boundary. Over the closed portion of the boundary, i.e. on the solid walls, (fu,fv) \cdot n should be zero. However, because the vorticity points are not on the velocity boundary faces, these terms do not vanish. The residual left from this contribution dominates the balance, so it is not a small error.


I have tried different versions of the vorticity definitions, i.e. the volume flux form as in  (equation 13) of Adcroft et al https://journals.ametsoc.org/view/journals/mwre/132/12/mwr2823.1.xml, or the simple derivatives and neither of them work. 


Does anyone have any words of wisdom or references that could be helpful?


Many thanks,
 Paola Cessi

--------------------------------------------------------------------------------------------
Paola Cessi                                                         Tel: +1 858 534 0622 
Scripps Institution of Oceanography                  Fax: +1 858 534 8045
9500 Gilman Drive #0213                             e-mail: pcessi at ucsd.edu 
La Jolla, CA 92093-0213          Web: http://pordlabs.ucsd.edu/pcessi
USA              
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20230905/3d04dcbd/attachment-0001.html>