[MITgcm-support] Fwd: Inter-tile discontinuities

[Pochini, Enrico]

Hi Martin,

thank you for your answer. I used the data.seaice and SEAICE_OPTIONS.h from
verification/lab_sea. I've tried to set the parameters as you said and I
also tried with SEAICEuseKrylov=.TRUE., which automatically sets
LSR_ERROR=1.E-12 and SEAICEnonLinIterMax=10. However I still have the issue.
I attach a link to the code and namelist which I used. In STDOUT it seems
that the solver for the sea ice dynamics does not converge.
Eventually the model crashes...

Enrico

 results72.nc
<https://drive.google.com/a/inogs.it/file/d/1IdUfelBlDrti9iK5RysnXMjYgSaH0-LW/view?usp=drive_web>

 results_seaice72.nc
<https://drive.google.com/a/inogs.it/file/d/1VrUQut_oD9yWhV9v4QHQG4mZ_-miETmc/view?usp=drive_web>

 sysdiag_48.nc
<https://drive.google.com/a/inogs.it/file/d/1DjQkpt4AAgPO8U_x_yDE1S3YOJqtCq76/view?usp=drive_web>

 SEAICE_OPTIONS.h
<https://drive.google.com/a/inogs.it/file/d/1HDxBg7PnZ1h-WcqSdInaDOBC47LLSUfJ/view?usp=drive_web>

 SEAICE_SIZE.h
<https://drive.google.com/a/inogs.it/file/d/1QkgiKhPcp7yJ0dI1Xk3rb61U_CrGNmY6/view?usp=drive_web>

 data.seaice
<https://drive.google.com/a/inogs.it/file/d/1Lw4bM0-T1IH-0ysFH7SgsH8RA7FQjnP7/view?usp=drive_web>

 data
<https://drive.google.com/a/inogs.it/file/d/16KvqeT5ZcD8otrsoyoyOg-vC6XxCECOF/view?usp=drive_web>

 STDOUT.0000
<https://drive.google.com/a/inogs.it/file/d/1Bu28s7ce58z4SKM2827WVPas9goOH6VH/view?usp=drive_web>

 job_err.err
<https://drive.google.com/a/inogs.it/file/d/1wkXF6ip0mPz1FFDMwxafHSctTHNPcJCa/view?usp=drive_web>

Il giorno mar 29 ott 2019 alle ore 17:17 Martin Losch <Martin.Losch at awi.de>
ha scritto:

> Hi Enrico,
>
> this is most likely linked to the dynamics solver. For historical reasons,
> the default parameters are not very fortunate. Without knowing your
> “data.seaice” and “SEAICE_OPTIONS.h”, I am guessing that you can fix this
> problem by reducing the linear tolerance (by setting LSR_ERROR to values of
> 1.e-5 or smaller).
> I also suggest to increase the number of nonlinear iterations to 10
> (SEAICEnonLinIterMax, default is 2). These unfortunate defaults are part of
> an issue (<https://github.com/MITgcm/MITgcm/issues/171>) on github, but I
> never got around to changing them to something more sensible (and nobody
> seems to care).
>
> Alternatively you could use the EVP code, but then make sure that you use
> at least mEVP (again, the defaults are not great), see the documentation (<
> https://mitgcm.readthedocs.io/en/latest/phys_pkgs/seaice.html#more-stable-variants-of-elastic-viscous-plastic-dynamics-evp-mevp-and-aevp
> >)
>
> Martin
>
> > On 29. Oct 2019, at 14:42, Pochini, Enrico <epochini at inogs.it> wrote:
> >
> > Dear all,
> >
> > (I've attached the file as gdrive link..)
> > I'm having issues with the SEAICE module i guess. After I activated it I
> started to observe discontinuities at tiles edge, where the marine ice is
> present. The discontinuities grow in time. I see discontinuities in surface
> fields and shallower levels of 3D fields.
> > What could cause this issue?
> >
> > Thanks,
> >
> > Enrico
> >  results24.nc
> > _______________________________________________
> > MITgcm-support mailing list
> > MITgcm-support at mitgcm.org
> > http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20191104/b672e9df/attachment.html>