[MITgcm-support] Non-Zero Meridional Velocity Kelvin Wave

[Martin Losch]

Hi Bertrand,

I didn’t read the thread very carefully but I found this

> I also checked my IC variables and those are perfectly symmetric about the equator (assuming that I am defining my initial U and V velocity at cell center, right?).

Just to avoid confusion: the initial conditions are actually specified at the corresponding grid points. I.e. dynamic topograph Eta is a the center of the C-cell, Uini and Vini are at the U- and V- points of the grid cell. If you assumed that they are all co-located, that may be the problem?

Martin


> On 17. Nov 2017, at 03:58, Bertrand Louis Rene Delorme <bdelorme at stanford.edu> wrote:
> 
> A little update on my problem (note also that the plot of Vm_dPHdy_Surf in my previous message is wrong):
> -          The full baroclinic pressure gradient force (ie: Vm_dPHdy-gravity*(ETAN(i) - ETAN(i-1) )/DYC)  and
> the Coriolis force (in Vm_Advec or Vm_Cori -- same) are the only forces occurring in the v-momentum equation 
> (AB_gV is non-zero but much smaller).
> -          I plotted those at the first time step (see enclosed). We can see that they do seem to have opposite spatial structures.
> Yet, looking at the difference of their amplitude ( ABS(f*u) - ABS(dPHI/dy) ), we can see that they do not balance
> each other perfectly. The induced acceleration is consistent with the plot of V.
>  
> So, I managed to find where the nonzero values in V come from, but still struggle to know why.
> If anyone has an idea, please let me know.
>  
> Thanks,
> Bertrand
>