[MITgcm-devel] [MITgcm-cvs] MITgcm/pkg/salt_plume CVS Commit

[Martin Losch]

what I forgot to say, is of course that the problem goes away, when I remove the masks for taux/y in seaice_get_dynforcing.F but that does not seem to be a good solution.

What also works is replacing
if (tmpscal4.GT.0) with
if (tmpscal3.LT.0. _d 0)
because then the discrimintant tmpscal4 is not only larger 0, but also larger than tmpscal2^2 so that the outer sqrt is also defined