Changeset 1550 for LMDZ5/trunk/libf/dyn3d

Timestamp:

Jul 5, 2011, 11:44:55 AM (14 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/dyn3d

Files:

• integrd.F (modified) (1 diff)
• vlsplt.F (modified) (2 diffs)

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