Changeset 5707

Timestamp:

Jun 16, 2025, 4:20:42 PM (7 days ago)

Author:

yann meurdesoif

Message:

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

YM

File:

• LMDZ6/trunk/libf/phylmd/cv3_routines.f90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/cv3_routines.f90

@@ -2134 +2134 @@
 ! compute icbmax:
 
-  icbmax = 2
-  DO i = 1, ncum
-    icbmax = max(icbmax, icb(i))
-  END DO
+!ym  icbmax = 2
+!ym  DO i = 1, ncum
+!ym    icbmax = max(icbmax, icb(i))
+!ym  END DO
 
 ! update sig and w0 below cloud base:
 
-  DO k = 1, icbmax
+!ym  DO k = 1, icbmax
+  DO k = 1, nd
     DO i = 1, ncum
-      IF (k<=icb(i)) THEN
-        sig(i, k) = beta*sig(i, k) - &
-                    2.*alpha*buoy(i, icb(i))*buoy(i, icb(i))
-        sig(i, k) = max(sig(i,k), 0.0)
-        w0(i, k) = beta*w0(i, k)
-      END IF
+      IF (k<=MAX(2,icb(i))) THEN
+        IF (k<=icb(i)) THEN
+          sig(i, k) = beta*sig(i, k) - &
+                      2.*alpha*buoy(i, icb(i))*buoy(i, icb(i))
+          sig(i, k) = max(sig(i,k), 0.0)
+          w0(i, k) = beta*w0(i, k)
+        END IF
+      ENDIF
     END DO
   END DO