Changeset 1550 for LMDZ5/trunk/libf

Timestamp:

Jul 5, 2011, 11:44:55 AM (13 years ago)

Author:

lguez

Message:

Bug fix in "bilan_dyn_p". The index was out of bounds in the removed
assignment . Also, the removed assignment was useless.

Bug fix in "coefkzmin". The size of a dummy array cannot exceed the
size of the associated actual array. ("coefkzmin" is called by
"coef_diff_turb".) "km(:, klev+1)" and "kn(:, klev+1)" were not
defined in "coefkzmin" so this was maybe an innocuous bug.

Location:

LMDZ5/trunk/libf

Files:

• dyn3d/integrd.F (modified) (1 diff)
• dyn3d/vlsplt.F (modified) (2 diffs)
• dyn3dpar/bilan_dyn_p.F (modified) (1 diff)
• dyn3dpar/integrd_p.F (modified) (1 diff)
• phylmd/coefkzmin.F (modified) (1 diff)
• phylmd/hgardfou.F (modified) (1 diff)

LMDZ5/trunk/libf/dyn3d/integrd.F

@@ -89 +89 @@
         IF( ps(ij).LT.0. ) THEN
          PRINT*,' Au point ij = ',ij, ' , pression sol neg. ', ps(ij)
-         STOP' dans integrd'
+         print *, ' dans integrd'
+         stop 1
         ENDIF
       ENDDO

LMDZ5/trunk/libf/dyn3d/vlsplt.F

@@ -129 +129 @@
       IMPLICIT NONE
 c
-#include "dimensions.h"
-#include "paramet.h"
-#include "logic.h"
-#include "comvert.h"
-#include "comconst.h"
+      include "dimensions.h"
+      include "paramet.h"
+      include "logic.h"
+      include "comvert.h"
+      include "comconst.h"
+      include "iniprint.h"
 c
 c
@@ -351 +352 @@
 
       IF(n0.gt.0) THEN
-      PRINT*,'Nombre de points pour lesquels on advect plus que le'
+      if (prt_level > 2) PRINT *,
+     $        'Nombre de points pour lesquels on advect plus que le'
      &       ,'contenu de la maille : ',n0
 

LMDZ5/trunk/libf/dyn3dpar/bilan_dyn_p.F

@@ -389 +389 @@
          Q_cum(:,jjb:jje,:,:)=0.
          flux_uQ_cum(:,jjb:jje,:,:)=0.
-         flux_v_cum(:,jjb:jje,:)=0.
          if (pole_sud) jje=jj_end-1
          flux_v_cum(:,jjb:jje,:)=0.

LMDZ5/trunk/libf/dyn3dpar/integrd_p.F

@@ -128 +128 @@
          PRINT*,' Au point ij = ',stop_it, ' , pression sol neg. '
      &         , ps(stop_it)
-         STOP' dans integrd'
+         print *, ' dans integrd'
+         stop 1
         ENDIF
 

LMDZ5/trunk/libf/phylmd/coefkzmin.F

@@ -50 +50 @@
       REAL teta(klon,klev)
 
-      REAL km(klon,klev+1)
-      REAL kn(klon,klev+1)
+      REAL km(klon,klev)
+      REAL kn(klon,klev)
       integer knon
 

LMDZ5/trunk/libf/phylmd/hgardfou.F

@@ -45 +45 @@
          jbad = 0
          DO i = 1, klon
-         IF (zt(i).GT.370.0) THEN
+         IF (zt(i) > 370.) THEN
             jbad = jbad + 1
             jadrs(jbad) = i