Changeset 5706 for LMDZ6/trunk/libf/phylmd/cv3_routines.f90

Timestamp:

Jun 16, 2025, 4:18:19 PM (3 weeks ago)

Author:

yann meurdesoif

Message:

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

YM

File:

• LMDZ6/trunk/libf/phylmd/cv3_routines.f90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/cv3_routines.f90

@@ -614 +614 @@
 
 ! Compute icbmax.
-
+  
+  !ym do not do that, independance between column !
   icbmax = 2
   DO i = 1, len
@@ -700 +701 @@
 
 ! Re-compute icbsmax (icbsmax2):                          !convect3
-!                                                         !convect3
-  icbsmax2 = 2                                            !convect3
-  DO i = 1, len                                           !convect3
-    icbsmax2 = max(icbsmax2, icbs(i))                     !convect3
-  END DO                                                  !convect3
+!
+!ym column independant, do not use reduction                                                         !convect3
+!ym  icbsmax2 = 2                                            !convect3
+!ym  DO i = 1, len                                           !convect3
+!ym    icbsmax2 = max(icbsmax2, icbs(i))                     !convect3
+!ym  END DO                                                  !convect3
 
 ! initialization outputs:
 
-  DO k = 1, icbsmax2                                      ! convect3
+!ym  DO k = 1, icbsmax2                                      ! convect3
+  DO k = 1, nd                                      ! convect3
     DO i = 1, len                                         ! convect3
-      tp(i, k) = 0.0                                      ! convect3
-      tvp(i, k) = 0.0                                     ! convect3
-      clw(i, k) = 0.0                                     ! convect3
+      IF (k<=MAX(2,icbs(i))) THEN
+        tp(i, k) = 0.0                                      ! convect3
+        tvp(i, k) = 0.0                                     ! convect3
+        clw(i, k) = 0.0                                     ! convect3
+      ENDIF
     END DO                                                ! convect3
   END DO                                                  ! convect3
@@ -718 +723 @@
 ! tp and tvp below cloud base:
 
-  DO k = minorig, icbsmax2 - 1
+!ym  DO k = minorig, icbsmax2 - 1
+  DO k = minorig, nd
     DO i = 1, len
-      tp(i, k) = tnk(i) - (gz(i,k)-gznk(i))*cpinv(i)
-      tvp(i, k) = tp(i, k)*(1.+qnk(i)/eps-qnk(i))        !whole thing (convect3)
+      IF (k<=MAX(2,icbs(i))-1) THEN
+        tp(i, k) = tnk(i) - (gz(i,k)-gznk(i))*cpinv(i)
+        tvp(i, k) = tp(i, k)*(1.+qnk(i)/eps-qnk(i))        !whole thing (convect3)
+      ENDIF
     END DO
   END DO