Hi,
Can anyone help me to short it out ?
In the Northern edge gradient BC shows some Wrong-boundary index (possibly).
Its doesn’t match with Southern BC. South is OK.
I think it will:
cff2=0.5*(h(i ,Jend+1)+zeta(i ,Jend+1,knew)
& +h(i-1,Jstr-1)+zeta(i-1,Jend+1,knew))
to
cff2=0.5*(h(i ,Jend+1)+zeta(i ,Jend+1,knew)
& +h(i-1,Jend+1)+zeta(i-1,Jend+1,knew))
The Northern boundary depth correction reads h(i-1, Jstr-1)
the Southern boundary row instead of h(i-1, Jend+1),
so cff2 (the water column depth used to rescale momentum) is
computed from the wrong location, corrupting
the barotropic correction at the Northern OBC.
Thank you for your kind response.
Best regards