Changeset 5578 for LMDZ6

Timestamp:

Mar 18, 2025, 11:32:47 AM (5 months ago)

Author:

Laurent Fairhead

Message:

Reindexing the do loop to prevent the matrix index reaching zero and crashing in debug mode.
Gives the same results as previously

File:

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

LMDZ6/trunk/libf/phylmd/tropopause_m.f90

@@ -20 +20 @@
   USE lmdz_reprobus_wrappers, ONLY: itroprep
   USE lmdz_cppkeys_wrapper, ONLY: CPPKEY_REPROBUS
+  USE print_control_mod, ONLY: lunout
 
 !-------------------------------------------------------------------------------
@@ -93 +94 @@
   DO i = 1, klon
     !--- UPPER TROPOPAUSE: |PV|=2PVU POINT STARTING FROM TOP
-    DO kt=klev-1,1,-1; savkt = kt; IF(ALL(ABS(pvor_cen(i,kt-nadj:kt))<=pv0)) EXIT; END DO
+!    DO kt=klev-1,1,-1
+!      savkt = kt
+!      IF (kt-nadj == 0) THEN
+!        WRITE(lunout,*)'ABORT_PHYSIC tropopause_m kt= ',kt
+!        call abort_physic("tropopause_m", " kt = nadj", 1)
+!      ENDIF
+!      IF(ALL(ABS(pvor_cen(i,kt-nadj:kt))<=pv0)) THEN 
+!        EXIT
+!      ENDIF
+!    END DO
+    DO kt=klev-1,nadj+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