Changeset 5705 for LMDZ6/trunk/libf/phylmd/cv_routines.f90

Timestamp:

Jun 13, 2025, 9:47:39 PM (4 weeks ago)

Author:

yann meurdesoif

Message:

Convection GPU porting : suppress potential dependency between columns, may change results (cv_closure)

YM

File:

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

LMDZ6/trunk/libf/phylmd/cv_routines.f90

@@ -808 +808 @@
   ! Compute icbmax.
   ! -------------------------------------------------------------------
-
-  icbmax = 2
-  DO i = 1, ncum
-    icbmax = max(icbmax, icb(i))
-  END DO
+  
+!ym independance betwwen column
+!ym  icbmax = 2
+!ym  DO i = 1, ncum
+!ym    icbmax = max(icbmax, icb(i))
+!ym  END DO
 
   ! =====================================================================
@@ -836 +837 @@
 
   ! dtpbl = average of tvp-tv in the PBL (k=nk to icb-1).
-
-  DO k = minorig, icbmax
+!ym independance betwwen column
+!ym  DO k = minorig, icbmax
+  DO k = minorig, nd
     DO i = 1, ncum
-      IF ((k>=nk(i)) .AND. (k<=(icb(i)-1))) THEN
-        dtpbl(i) = dtpbl(i) + (tvp(i,k)-tv(i,k))*dph(i, k)
-      END IF
+      IF (k<=MAX(2,icb(i))) THEN
+        IF ((k>=nk(i)) .AND. (k<=(icb(i)-1))) THEN
+          dtpbl(i) = dtpbl(i) + (tvp(i,k)-tv(i,k))*dph(i, k)
+        END IF
+      ENDIF
     END DO
   END DO