[MITgcm-devel] obcs_seaice_sponge.F

[Menemenlis, Dimitris (3248)]

Jean-Michel and Patrick, I know that it is no longer actively supported but, for the record, there
seems to be a problem with obcs_sponge.F, and it dates back to very first check in 10+ years ago.

According to comments, the relaxation time scale is supposed to be linear
and controlled by U/Vrelaxobcsinner/bound:

C U/Vrelaxobcsinner/bound :: relaxation time scale (in seconds) on the boundary
C                            (bound) and at the innermost grid point of the sponge
C                            layer (inner); relaxation time scales in-between
C                            are linearly interpolated from these values

But as coded, the relaxation is quadratic because there is an additional attenuation term due
to linearly decreasing u/v/t/srelax as it goes away from the boundary.

Additionally, as coded, relaxation will be zero if spongeThickness = 1.

I will fix these two problems in obcs_seaice_sponge.F before checking in,

Let me know if you want me to do anything about these two problem in obcs_sponge.F

Cheers

Dimitris Menemenlis