Changeset 4649 for LMDZ6/trunk/libf/phylmd

Timestamp:

Aug 28, 2023, 7:58:11 AM (15 months ago)

Author:

fhourdin

Message:

Supression de depassements de tableaux

Routines non utilisees dans les configurations standard

Location:

LMDZ6/trunk/libf/phylmd

Files:

• vdif_kcay.F90 (modified) (3 diffs)
• yamada.F90 (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/vdif_kcay.F90

@@ -221 +221 @@
   END IF
 
-  DO ilev = 1, nlev
+  DO ilev = 2, nlay
     DO igrid = 1, ngrid
       q2(igrid, ilev) = amax1(q2(igrid,ilev), q2min)
@@ -591 +591 @@
   IF (1==1) THEN
 
+    sss=0.
+    sssq=0.
+    ! WARNING : travail sur le point ig=1 ????
     DO ilev = 2, klev - 1
       sss = sss + plev(1, ilev-1) - plev(1, ilev+1)
@@ -624 +627 @@
 
 
-    DO ilev = 1, nlev
+    DO ilev = 2, klev-1
       DO igrid = 1, ngrid
         q2(igrid, ilev) = max(q2(igrid,ilev), q2min)

LMDZ6/trunk/libf/phylmd/yamada.F90

@@ -128 +128 @@
     l0(ig) = 100.
   END DO
-  DO k = 2, klev - 1
+  DO k = 2, klev
     DO ig = 1, ngrid
       l(ig, k) = l0(ig)*kap*zlev(ig, k)/(kap*zlev(ig,k)+l0(ig))
@@ -139 +139 @@
       sqz(ig) = 1.E-10
     END DO
-    DO k = 2, klev - 1
+    !DO k = 2, klev - 1
+    DO k = 2, klev
       DO ig = 1, ngrid
         q2(ig, k) = l(ig, k)**2*zz(ig, k)