Opened 6 years ago
Closed 5 years ago
#2107 closed Bug (fixed)
non-zero ice concentration on land
| Reported by: | clem | Owned by: | systeam |
|---|---|---|---|
| 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)
| Changeset | Author | Time | ChangeLog |
|---|---|---|---|
| 9883 | clem | 2018-07-05T17:57:31+02:00 | solve issue raised in ticket #2107 and maybe in ticket #2106 |
| 9882 | clem | 2018-07-05T17:51:31+02:00 | finish solving the issue raised in ticket #2107 |
| 9866 | clem | 2018-07-02T18:19:04+02:00 | solve problem described in ticket #2107 |
Change History (6)
comment:1 Changed 6 years ago by clem
comment:2 Changed 6 years 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 6 years ago by clem
Actually Prather scheme has the same issue. Therefore a debug is needed in nemo 3.6 as well.
comment:4 Changed 6 years ago by clem
In 9882:
comment:5 Changed 6 years ago by clem
In 9883:
comment:6 Changed 5 years ago by clem
- Resolution set to fixed
- Status changed from new to closed
In 9866: