Changeset 62 in lmdz_wrf

Timestamp:

Jul 24, 2014, 11:37:02 AM (10 years ago)

Author:

lfita

Message:

Adding check_var3D to check variables after thermals

File:

• branches/LMDZ_WRFmeas/WRFV3/lmdz/physiq.F90 (modified) (4 diffs)

branches/LMDZ_WRFmeas/WRFV3/lmdz/physiq.F90

@@ -1301 +1301 @@
 ! Lluis
       INTEGER                                            :: llp
+      CHARACTER(LEN=50)                                  :: lvarname, lfunc
+      REAL                                               :: largest
+
       llp = 734
+      lfunc = 'physiq'
+      largest = 10.e5
 
 !c======================================================================
@@ -2951 +2956 @@
 
       endif
+      lvarname = 't_seri'
+      fname='after thermals'
+      CALL check_var3D(fname, varname, t_seri, klon, klev, largest, .FALSE.)
+      lvarname = 'q_seri'
+      CALL check_var3D(fname, varname, q_seri, klon, klev, largest, .FALSE.)
+      lvarname = 'd_t_ajsb'
+      CALL check_var3D(fname, varname, d_t_ajsb, klon, klev, largest, .FALSE.)
+      lvarname = 'd_q_ajsb'
+      CALL check_var3D(fname, varname, d_q_ajsb, klon, klev, largest, .FALSE.)
+
 !c
 !c===================================================================
@@ -3039 +3054 @@
         call writefield_phy('q_seri',q_seri,llm)
       endif
-
 !c
 !c-------------------------------------------------------------------
@@ -4350 +4364 @@
       RETURN
       END SUBROUTINE gr_fi_ecrit
+
+SUBROUTINE check_var3D(funcn, varn, var, sizev, zsize, bigvalue, stoprun)
+!  Subroutine to check the consistency of a 3D LMDSZ - variable (klon, klev) !
+!    * NaN value: by definition is variable /= variable
+!    * bigvalue: 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)                                       :: bigvalue
+  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,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
+! bigvalue: biggest attenaible value for the variable
+! 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 (var(i,k) /= var(i,k) .OR. ABS(var(i,k)) > bigvalue ) 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 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 check_var3D
+