[MITgcm-support] Coriolis terms with partially filled cells

Noriyuki Yamamoto nymmto at kugi.kyoto-u.ac.jp
Tue Jan 6 06:25:45 EST 2015


Dear Martin

Thank you very much for your help!
I'll try it soon.

I have been trying to close the zonally integrated momentum budget and needed volume weighted scheme for Coriolis terms.
If the modification caused significant changes in the model result,  I'd report them later.

Noriyuki

2015/01/06 18:48、Martin Losch <Martin.Losch at awi.de> のメール:

> 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
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
> 




More information about the MITgcm-support mailing list