Changeset 5699

Timestamp:

Jun 13, 2025, 9:01:23 PM (4 weeks ago)

Author:

yann meurdesoif

Message:

Convection GPU porting : remove zilch function (set to 0 flattened array) or similar...
YM

Location:

LMDZ6/trunk/libf/phylmd

Files:

• cv3_routines.f90 (modified) (3 diffs)
• cv3p_mixing.f90 (modified) (1 diff)
• cv_routines.f90 (modified) (6 diffs)
• cva_driver.f90 (modified) (3 diffs)

LMDZ6/trunk/libf/phylmd/cv3_routines.f90

@@ -2520 +2520 @@
 ! ---  TO REPRESENT EQUAL PROBABILITIES OF MIXING
 ! =====================================================================
-
-  CALL zilch(asum, nloc*nd)
-  CALL zilch(csum, nloc*nd)
-  CALL zilch(csum, nloc*nd)
+  asum(1:nloc,1:nd) = 0.
+  csum(1:nloc,1:nd) = 0.
 
   DO il = 1, ncum
@@ -2688 +2686 @@
 
 ! MAF: renormalisation de MENT
-  CALL zilch(zm, nloc*na)
+  zm(1:nloc,1:na) = 0.
+  
   DO jm = 1, nl
     DO im = 1, nl
@@ -2919 +2918 @@
     IF (num1<=0) GO TO 400
 
-    CALL zilch(wdtrain, ncum)
+    wdtrain(1:ncum) = 0.0
 
 

LMDZ6/trunk/libf/phylmd/cv3p_mixing.f90

@@ -294 +294 @@
 ! =====================================================================
 
-  CALL zilch(csum, nloc*nd)
-
+  csum(:,:) = 0.
+  
   DO il = 1, ncum
     lwork(il) = .FALSE.

LMDZ6/trunk/libf/phylmd/cv_routines.f90

@@ -719 +719 @@
   ! J Teixeira fix
 
-  CALL zilch(byp, ncum)
+  byp(1:ncum) = 0
+
   DO i = 1, ncum
     lcape(i) = .TRUE.
@@ -885 +886 @@
   ! =====================================================================
 
-  DO i = 1, ncum*nlp
-    nent(i, 1) = 0
-    m(i, 1) = 0.0
-  END DO
+!ym very bad
+!ym  DO i = 1, ncum*nlp
+!ym    nent(i, 1) = 0
+!ym    m(i, 1) = 0.0
+!ym  END DO
+  DO k = 1, nlp
+    DO i = 1, ncum 
+      nent(i, k) = 0
+      m(i, k) = 0.0
+    ENDDO
+  ENDDO
 
   DO k = 1, nlp
@@ -907 +915 @@
   ! -------------------------------------------------------------------
 
-  CALL zilch(work, ncum)
-
+  work(1:ncum) = 0.
+ 
   DO j = minorig + 1, nl
     DO i = 1, ncum
@@ -1004 +1012 @@
   ! =====================================================================
 
-  CALL zilch(bsum, ncum*nlp)
+  bsum(1:ncum,1:nlp) = 0.
   DO ij = 1, ncum
     lwork(ij) = .FALSE.
@@ -1180 +1188 @@
 
 
-  CALL zilch(wdtrain, ncum)
+  wdtrain(1:ncum) = 0.
   DO i = nl + 1, 1, -1
 
@@ -1440 +1448 @@
     IF (num1<=0) GO TO 1500
 
-    CALL zilch(amp1, ncum)
-    CALL zilch(ad, ncum)
+    amp1(1:ncum)=0.
+    ad(1:ncum)=0.
 
     DO k = i + 1, nl + 1

LMDZ6/trunk/libf/phylmd/cva_driver.f90

@@ -1004 +1004 @@
 !      END IF
       IF (iflag_mix>=1) THEN
-        CALL zilch(supmax, nloc*nd)
+        supmax(:,:)=0.
         if (prt_level >= 9) &
              PRINT *, 'cva_driver -> cv3p_mixing'
@@ -1016 +1016 @@
 
       ELSE
-        CALL zilch(supmax, nloc*nd)
+        supmax(:,:)=0.
       END IF
     END IF
@@ -1093 +1093 @@
                         unk, vnk, hp, tv, tvp, ep, clw, m, sig, &
                         ment, qent, uent, vent, nent, sigij, elij, ments, qents, traent)
-        CALL zilch(hent, nloc*nd*nd)
+        hent(1:nloc,1:nd,1:nd) = 0.
       ELSE
 !!jyg:  Essais absurde pour voir