[MITgcm-devel] lsr/obcs

[Martin Losch]

Hi Jean-Michel,

just so I don't misunderstand the code for obcs in seaice_lsr.F:
#ifndef OBCS_UVICE_OLD
[…]
            BU(I,J,bi,bj) = ONE
[…]
            rhsU(I,J,bi,bj) = uIce(I,J,bi,bj)
#endif

This is absolutely correct, but would be more consistent (in terms of units) if it were
BU =  seaiceMassU(I,J,bi,bj)/SEAICE_deltaTdyn
rhsU = seaiceMassU(I,J,bi,bj)/SEAICE_deltaTdyn * uIce (I,J,bi,bj)

Is that correct? I guess scaling these values by mass/dt only has an effect on the residual computations. Do you agree?

Martin