Changeset 59 in lmdz_wrf for branches/LMDZ_WRFmeas/WRFV3

Timestamp:

Jul 23, 2014, 5:58:14 PM (10 years ago)

Author:

lfita

Message:

Adding real x,y coordinates of the wrong values

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/diagphy_mod.F90 (modified) (6 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/diagphy_mod.F90

@@ -269 +269 @@
   IMPLICIT NONE
 
+#include "dimensions.h"
+
   INTEGER, INTENT(IN)                                    :: sizev
   CHARACTER(LEN=50), INTENT(IN)                          :: funcn, varn
@@ -276 +278 @@
 
 ! Local
-  INTEGER                                                :: i, wrongi
+  INTEGER                                                :: i, wrongi, xpt, ypt
   CHARACTER(LEN=50)                                      :: errmsg
   LOGICAL                                                :: found
@@ -306 +308 @@
     PRINT *,TRIM(errmsg)
     PRINT *,"  at '" // TRIM(funcn) // "' variable '" //TRIM(varn)//                 &
-      "' is wrong in Nvalues= ",wrongi,' at i value___'
+      "' is wrong in Nvalues= ",wrongi,' at i (x, y) value___'
     DO i=1,wrongi
-      PRINT *,wronggridpt(i), wrongvalues(i)
+       ypt = INT(wronggridpt(i,1)/wiim) + 1
+       xpt = wronggridpt(i,1) - (ypt-1)*wiim
+      PRINT *,wronggridpt(i), '(',xpt,', ',ypt,')', wrongvalues(i)
     END DO
     IF (stoprun) THEN
@@ -326 +330 @@
   IMPLICIT NONE
 
+#include "dimensions.h"
+
   INTEGER, INTENT(IN)                                    :: sizev, zsize
   CHARACTER(LEN=50), INTENT(IN)                          :: funcn, varn
@@ -333 +339 @@
 
 ! Local
-  INTEGER                                                :: i, k, wrongi
+  INTEGER                                                :: i, k, wrongi, xpt, ypt
   CHARACTER(LEN=50)                                      :: errmsg
   LOGICAL                                                :: found
@@ -367 +373 @@
     PRINT *,TRIM(errmsg)
     PRINT *,"  at '" // TRIM(funcn) // "' variable '" //TRIM(varn)//                 &
-      "' is wrong in Nvalues= ",wrongi,' at i k value___'
+      "' is wrong in Nvalues= ",wrongi,' at i (x,y) k value___'
     DO i=1,wrongi
-      PRINT *,wronggridpt(i,1), wronggridpt(i,2), wrongvalues(i)
+       ypt = INT(wronggridpt(i,1)/wiim) + 1
+       xpt = wronggridpt(i,1) - (ypt-1)*wiim
+      PRINT *,wronggridpt(i,1), '(',xpt,', ',ypt,')', wronggridpt(i,2), wrongvalues(i)
     END DO
     IF (stoprun) THEN