[MITgcm-devel] seaice velocity mask

[Martin Losch]

Gael,

I just went back and checked, that I actually had seaiceMaskU/V
initialised in seaice_init_fixed.F (Revision 1.4), but Patrick removed
this part (Revision 1.5), probably, because they are re-initialised in
seaice_ini_varia.F (maybe because it's easier to taf to recognize
dependencies?).

I suggest, that we add a flag ALLOW_SEAICE_MASKRHS (undefined by
default) that is put ifdefs around the update of seaiceMaskU/V in
seaice_dynsolver.F and also in seaice_ini_varia.F and re-introduce the
initialisation of the fields into seaice_ini_fixed.F. Will that have
negative consequences for the adjoint?

After this works I guess that these fields need not be stored anymore
(except when ALLOW_SEAICE_MASKRHS is defined) and these directives can
be either removed (with a warning when both ALLOW_TAMC_AUTODIFF and
ALLOW_SEAICE_MASKRHS are defined) or also put between #ifdef
ALLOW_SEAICE_MASKRHS.

What do you think? and what does Patrick think about this?

Martin