Changeset 3893

Timestamp:

Jan 25, 2016, 10:37:45 AM (9 years ago)

Author:

ymipsl

Message:

Bug fix for LMDZ : some operation was perform on whole range of an compressed array, making floating point exception when computing on invalid part of arrays.

YM

File:

• dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/yamada4.F90 (modified) (3 diffs)

dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/yamada4.F90

@@ -247 +247 @@
 
     DO k = 2, klev
-      q2(:, k) = l(:, k)**2*zz(:, k)
+      q2(:ngrid, k) = l(:ngrid, k)**2*zz(:ngrid, k)
     END DO
 
@@ -340 +340 @@
     ! print*,'Longueur ',l(:,:)
     DO k = 2, klev - 1
-      l(:, k) = max(l(:,k), 1.)
-      km(:, k) = l(:, k)*sqrt(q2(:,k))*sm(:, k)
-      q2(:, k) = q2(:, k) + dt*km(:, k)*m2(:, k)*(1.-rif(:,k))
-      q2(:, k) = min(max(q2(:,k),1.E-10), 1.E4)
-      q2(:, k) = 1./(1./sqrt(q2(:,k))+dt/(2*l(:,k)*b1))
-      q2(:, k) = q2(:, k)*q2(:, k)
+      DO ig = 1, ngrid
+        l(ig, k) = max(l(ig,k), 1.)
+        km(ig, k) = l(ig, k)*sqrt(q2(ig,k))*sm(ig, k)
+        q2(ig, k) = q2(ig, k) + dt*km(ig, k)*m2(ig, k)*(1.-rif(ig,k))
+        q2(ig, k) = min(max(q2(ig,k),1.E-10), 1.E4)
+        q2(ig, k) = 1./(1./sqrt(q2(ig,k))+dt/(2*l(ig,k)*b1))
+        q2(ig, k) = q2(ig, k)*q2(ig, k)
+      END DO
     END DO
 
@@ -377 +379 @@
   IF (iflag_pbl>=12) THEN
     ! print*,'YAMADA VDIF'
-    q2(:, 1) = q2(:, 2)
+    q2(1:ngrid, 1) = q2(1:ngrid, 2)
     CALL vdif_q2(dt, g, rconst, ngrid, plev, temp, kq, q2)
   END IF