Changeset 170 in lmdz_wrf

Timestamp:

Aug 12, 2014, 11:22:27 AM (10 years ago)

Author:

lfita

Message:

Adding real x,y coordinates of the wrong value

File:

• branches/LMDZ_WRFmeas_develop/WRFV3/lmdz/abort_gcm.F90 (modified) (6 diffs)

branches/LMDZ_WRFmeas_develop/WRFV3/lmdz/abort_gcm.F90

@@ -57 +57 @@
 #endif
 #include "iniprint.h"
+#include "dimensions.h"
  
 !C
@@ -77 +78 @@
       CHARACTER(LEN=1), INTENT(IN)                       :: if1c
 
+! Local
+      INTEGER                                            :: xpt, ypt
+
       write(lunout,*) 'in abort_gcm'
 #ifdef CPP_IOIPSL
@@ -87 +91 @@
 !c     call histclo(4)
 !c     call histclo(5)
+      ypt = INT(wronggridpt(i,1)/wiim) + 1
+      xpt = wronggridpt(i,1) - (ypt-1)*wiim
+
       write(lunout,*) 'Stopping in ', modname
-      write(lunout,*) 'Stopping at point ', point,' wrong: ',val, if1c, thres
+      write(lunout,*) 'Stopping at point ', point,'( ',xpt,', ',ypt,' ) wrong: ',    &
+        val, if1c, thres
       write(lunout,*) 'Reason = ',message
       if (ierr .eq. 0) then
@@ -109 +117 @@
 #endif
 #include "iniprint.h"
+#include "dimensions.h"
  
 !C
@@ -129 +138 @@
       CHARACTER(LEN=1), INTENT(IN)                       :: if1c
 
+! Local
+      INTEGER                                            :: xpt, ypt
+
       write(lunout,*) 'in abort_gcm'
 #ifdef CPP_IOIPSL
@@ -139 +151 @@
 !c     call histclo(4)
 !c     call histclo(5)
+      ypt = INT(wronggridpt(i,1)/wiim) + 1
+      xpt = wronggridpt(i,1) - (ypt-1)*wiim
+
       write(lunout,*) 'Stopping in ', modname
-      write(lunout,*) 'Stopping at point (klon, klev):', pk, ',', pl,' wrong: ',val, &
-        if1c, thres
+      write(lunout,*) 'Stopping at point (klon, klev):', pk, '( ',xpt,', ',ypt,' ) ',&
+        pl,' wrong: ',val, if1c, thres
       write(lunout,*) 'Reason = ',message
       if (ierr .eq. 0) then