Changeset 157 in lmdz_wrf for branches/LMDZ_WRFmeas_develop/WRFV3/lmdz/diagphy_mod.F90

Timestamp:

Aug 8, 2014, 11:54:57 AM (10 years ago)

Author:

lfita

Message:

Adding 'threscheck_var3D', subroutine to check accordingly a threshold

File:

• branches/LMDZ_WRFmeas_develop/WRFV3/lmdz/diagphy_mod.F90 (modified) (1 diff)

branches/LMDZ_WRFmeas_develop/WRFV3/lmdz/diagphy_mod.F90

@@ -388 +388 @@
 END SUBROUTINE check_var3D
 
+SUBROUTINE threscheck_var3D(funcn, varn, var, sizev, zsize, thresvalue, kthres, 
+  stoprun)
+!  Subroutine to check the consistency of a 3D LMDSZ - variable (klon, klev) !
+!    * thresvalue: threshold for the variable
+
+  IMPLICIT NONE
+
+#include "dimensions.h"
+
+  INTEGER, INTENT(IN)                                    :: sizev, zsize
+  CHARACTER(LEN=50), INTENT(IN)                          :: funcn, varn
+  REAL, DIMENSION(sizev,zsize), INTENT(IN)               :: var
+  REAL, INTENT(IN)                                       :: thresvalue
+  LOGICAL, INTENT(IN)                                    :: stoprun
+
+! Local
+  INTEGER                                                :: i, k, wrongi, xpt, ypt
+  CHARACTER(LEN=50)                                      :: errmsg
+  LOGICAL                                                :: found
+  REAL, DIMENSION(sizev*zsize)                           :: wrongvalues
+  INTEGER, DIMENSION(sizev*zsize,2)                      :: wronggridpt
+
+!!!!!!! Variables
+! funcn: at which functino of part of the program variable is checked
+! varn: name of the variable
+! var: variable to check
+! sizev: size of the variable
+! zsize: vertical size of the variable
+! thresvalue: threshold attenaible value for the variable
+! kthres: kind of threshold= -1, below wrong, 1: above wrong
+! stoprun: Should the run stop if it founds a problem?
+
+  errmsg = 'ERROR -- error -- ERROR -- error'
+
+  found = .FALSE.
+  wrongi = 0
+  DO i=1,sizev
+    DO k=1,zsize
+      IF (kthres == -1 .AND. var(i,k) < thresvalue ) THEN
+        IF (wrongi == 0) found = .TRUE.
+        wrongi = wrongi + 1
+        wrongvalues(wrongi) = var(i,k)
+        wronggridpt(wrongi,1) = i
+        wronggridpt(wrongi,2) = k
+      ELSE IF (kthres == 1 .AND. var(i,k) > thresvalue ) THEN
+        IF (wrongi == 0) found = .TRUE.
+        wrongi = wrongi + 1
+        wrongvalues(wrongi) = var(i,k)
+        wronggridpt(wrongi,1) = i
+        wronggridpt(wrongi,2) = k
+      END IF
+    END DO
+  END DO
+
+  IF (found) THEN
+    PRINT *,TRIM(errmsg)
+    PRINT *,"  at '" // TRIM(funcn) // "' variable '" //TRIM(varn)//                 &
+      "' is wrong (threshold: ',thresvalue,') in Nvalues= ",wrongi,                  &
+      ' at i (x,y) k value___'
+    DO i=1,wrongi
+       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
+      STOP
+    END IF
+  END IF
+
+  RETURN
+
+END SUBROUTINE threscheck_var3D
 
 !C======================================================================