[MITgcm-support] funky ice dynamics in doubly periodic domain

[Jean Mensa]

Hi all,
I have a doubly periodic high resolution simulation forced with ECMWF
atmospheric forcing (forcings.png) running the seaice pkg with VP rheology.
Thermodynamics variables look good (ice.png) but ice velocity does
something funky and I am not sure why. During summer season ice thickness
is minimum (H_I) and concentration low (A_I) and ice dynamics seems to be
driven by wind stress. This is all good and fine except that during winter
ice velocity (S_I) seems to converge to some constant value decoupling from
wind stress. This does not seem right to me. Is that an artifact of the
doubly periodic configuration? Notice that the simulation runs fine and
gives no failed convergence warnings.
Thanks,
Jean


data.seaice,

 &SEAICE_PARM01
  SEAICEuseFluxForm = .TRUE.,
  SEAICEuseFlooding  = .TRUE.,
  SEAICE_EPS         = 1.E-8,
  SEAICE_drag =  0.00114,
  SEAICEuseDYNAMICS  = .TRUE.,
  SEAICEadvScheme    = 33,
  LSR_ERROR = 1.e-6,
  NPSEUDOTIMESTEPS = 20,
 &
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160222/7fcf6729/attachment-0001.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: forcings.png
Type: image/png
Size: 124520 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160222/7fcf6729/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ice.png
Type: image/png
Size: 77463 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20160222/7fcf6729/attachment-0003.png>