Changeset 5344 for LMDZ6/trunk

Timestamp:

Nov 28, 2024, 5:16:07 PM (2 weeks ago)

Author:

Laurent Fairhead

Message:

Temporary fix on value of index so that the code does not crash in debug mode.
We still have to understand how it can be 0
TL, JG, DC, LF

File:

• LMDZ6/trunk/libf/phylmd/tropopause_m.f90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/tropopause_m.f90

@@ -44 +44 @@
   REAL,    PARAMETER :: w(5)=[0.1,0.25,0.3,0.25,0.1] !--- Vertical smoothing
   INTEGER, SAVE :: k0
+  INTEGER :: savkt
   LOGICAL, SAVE :: first=.TRUE.
 !$OMP THREADPRIVATE(k0,first)
@@ -92 +93 @@
   DO i = 1, klon
     !--- UPPER TROPOPAUSE: |PV|=2PVU POINT STARTING FROM TOP
-    DO kt=klev-1,1,-1; IF(ALL(ABS(pvor_cen(i,kt-nadj:kt))<=pv0)) EXIT; END DO
+    DO kt=klev-1,1,-1; savkt = kt; IF(ALL(ABS(pvor_cen(i,kt-nadj:kt))<=pv0)) EXIT; END DO
+    kt = savkt
+    IF (kt == 0 ) THEN
+       call abort_physic("dyn_tropopause", " kt = 1", 1)
+    ENDIF
     !--- LOWER TROPOPAUSE: |PV|=2PVU POINT STARTING FROM BOTTOM
     DO kb=k0,klev-1;   IF(ALL(ABS(pvor_cen(i,kb:kb+nadj))> pv0)) EXIT; END DO; kb=kb-1