#2107 closed Bug (fixed)

non-zero ice concentration on land

Reported by: clem Owned by:
Priority: low Milestone:
Component: LIM3 Version: trunk
Severity: major Keywords:
Cc:

Description

Context

In some configurations, specifically when using BDY with sea ice, the model crashes due to a floating point exception in the ice thermodynamics.
It ends up that the crash results from the advection scheme (Ultimate-Macho). For some reason (I do not know why), the scheme sometime produces non-zero ice fields in a few points on land. Therefore, a non-zero ice concentration on land may occur. A freezing temperature is then calculated at 273.15K (since sss=0) and thus a division by 0 occurs in icethd_dh.F90
Note that the problem may also arise for Prather advection scheme but it is hidden since the ice fields are land-masked in this case at the end of the routine.

In theory this problem also concerns NEMO 3.6 but no such behavior has been observed so far except maybe in a recent ticket (#2106)

Fix

The advected fields resulting from the Ultimate-Macho scheme must be land-masked (as for Prather). I suppose this is not a clean fix but it does the job

Commit History (3)

ChangesetAuthorTimeChangeLog
9883clem2018-07-05T17:57:31+02:00

solve issue raised in ticket #2107 and maybe in ticket #2106

9882clem2018-07-05T17:51:31+02:00

finish solving the issue raised in ticket #2107

9866clem2018-07-02T18:19:04+02:00

solve problem described in ticket #2107

Change History (6)

comment:1 Changed 22 months ago by clem

In 9866:

solve problem described in ticket #2107

comment:2 Changed 22 months ago by clem

It turns out that NEMO 3.6 is not affected since Ultimate-Macho has been developed for NEMO4 only.
This issue is solved at r9866

comment:3 Changed 21 months ago by clem

Actually Prather scheme has the same issue. Therefore a debug is needed in nemo 3.6 as well.

comment:4 Changed 21 months ago by clem

In 9882:

finish solving the issue raised in ticket #2107

comment:5 Changed 21 months ago by clem

In 9883:

solve issue raised in ticket #2107 and maybe in ticket #2106

comment:6 Changed 15 months ago by clem

  • Resolution set to fixed
  • Status changed from new to closed
Note: See TracTickets for help on using tickets.