Changeset 4180 for trunk/LMDZ.COMMON/libf/evolution/orbit.F90

Timestamp:

Apr 10, 2026, 7:17:55 PM (6 hours ago)

Author:

jbclement

Message:

PEM:

• Rework layering-related logic, especially clarify interactions between surface and subsurface water tendencies and disable CO2 lag resistance (inconsistent without updating pressure and mass balance + PCM).
• Prevent simultaneous activation of layering and ice flows.
• Add warning when flux_geo /= 0 while soil is disabled.
• Add new utility function "is_lvl_enabled" for displaying.
• Replace deprecated 'minieps' with 'eps'/'tol'.

JBC

File:

• trunk/LMDZ.COMMON/libf/evolution/orbit.F90 (modified) (2 diffs)

trunk/LMDZ.COMMON/libf/evolution/orbit.F90

@@ -15 +15 @@
 ! DEPENDENCIES
 ! ------------
-use numerics, only: dp, di, k4, largest_nb, minieps
+use numerics, only: dp, di, k4, largest_nb, eps
 
 ! DECLARATION
@@ -477 +477 @@
 call print_msg('Obl. (current|accepted min|accepted max): '//real2str(obliquity)//' | '//real2str(min_obl)//real2str(max_obl),LVL_NFO)
 
-max_ecc = min(eccentricity + max_change_ecc,1. - minieps) ! ecc < 1.
+max_ecc = min(eccentricity + max_change_ecc,1. - eps) ! ecc < 1.
 min_ecc = max(eccentricity - max_change_ecc,0._dp) ! ecc >= 0.
 call print_msg('Ecc. (current|accepted min|accepted max): '//real2str(eccentricity)//' | '//real2str(min_ecc)//' | '//real2str(max_ecc),LVL_NFO)