Changeset 5708

Timestamp:

Jun 16, 2025, 4:26:48 PM (4 weeks ago)

Author:

yann meurdesoif

Message:

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

YM

File:

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

LMDZ6/trunk/libf/phylmd/cv3p1_closure.f90

@@ -154 +154 @@
   ! compute icbmax:
 
-  icbmax = 2
-  DO il = 1, ncum
-    icbmax = max(icbmax, icb(il))
-  END DO
+!ym break the column independance
+!ym  icbmax = 2
+!ym  DO il = 1, ncum
+!ym    icbmax = max(icbmax, icb(il))
+!ym  END DO
+
   ! if(prt.level.GE.20) print*,'cv3p1_param apres 200'
 
   ! update sig and w0 below cloud base:
-
-  DO k = 1, icbmax
-    DO il = 1, ncum
-      IF (k<=icb(il)) THEN
-        sig(il, k) = beta*sig(il, k) - 2.*alpha*buoy(il, icb(il))*buoy(il, &
-          icb(il))
-        sig(il, k) = amax1(sig(il,k), 0.0)
-        w0(il, k) = beta*w0(il, k)
-      END IF
+!ym column independance
+!ym  DO k = 1, icbmax
+  DO k = 1, nd
+    DO il = 1, ncum
+      IF (k<=MAX(2,icb(il))) THEN
+        IF (k<=icb(il)) THEN
+          sig(il, k) = beta*sig(il, k) - 2.*alpha*buoy(il, icb(il))*buoy(il, &
+            icb(il))
+          sig(il, k) = amax1(sig(il,k), 0.0)
+          w0(il, k) = beta*w0(il, k)
+        END IF
+      ENDIF
     END DO
   END DO