[MITgcm-support] Coriolis terms with partially filled cells

[Martin Losch]

Dear Noriyuki,

I agree with you, the weighting with A_c*h_c is missing in the code (\Delta{r}_f/c is constant and drops out). I don’t know of updated code, but you can easily try it yourself like this:
         uCoriolisTerm(i,j) =
     &     0.5*( _fCori( i ,j,bi,bj)
     &          *0.5*( vFld( i ,j)+vFld( i ,j+1) )* _rA(i,j,bi,bj)* _hFacC(i,j,k,bi,bj)
     &         + _fCori(i-1,j,bi,bj)
     &          *0.5*( vFld(i-1,j)+vFld(i-1,j+1) )*_ rA(i-1,j,bi,bj)* _hFacC(i-1,j,k,bi,bj)
     &   )*recip_hfacw(i,j,k,bi,bj)*recip_rAw(i,j,bi,bj)

And similarly for the vCoriolisTerm (replace “w” with “s” and beware of the indices) and some of the non-hydrostatic terms.
It would be interesting to learn about the effects of this modification.

Martin

On Dec 26, 2014, at 3:59 PM, Noriyuki Yamamoto <nymmto at kugi.kyoto-u.ac.jp> wrote:

> Hi all.
> 
> I'm using partially filled cells to represent smooth topography.
> And I found that codes for energy conserving coriolis terms (in
> mom_u_coriolis.F, mom_v_coriolis.F)
> do not correspond to that in the online document
> (http://mitgcm.org/public/r2_manual/latest/online_documents/node54.html)
> especially when one uses partially filled cell or variable grid spacing.
> 
> Please tell me if there's already a modified version.
> Cheers, Noriyuki.
> 
> -- 
> Noriyuki Yamamoto
> PhD Student - Physical Oceanography Group 
> Division of Earth and Planetary Sciences,
> Graduate School of Science, Kyoto University.
> Mail:nymmto at kugi.kyoto-u.ac.jp
> Tel:+81-75-753-3924
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support