Changeset 52 in lmdz_wrf for branches/LMDZ_WRFmeas

Timestamp:

Jul 23, 2014, 5:21:24 PM (10 years ago)

Author:

lfita

Message:

Adding 'check_vars3D' for LMDZ 3D variables (klon,klev) !

File:

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

branches/LMDZ_WRFmeas/WRFV3/lmdz/diagphy_mod.F90

@@ -306 +306 @@
     PRINT *,TRIM(errmsg)
     PRINT *,"  at '" // TRIM(funcn) // "' variable '" //TRIM(varn)//                 &
-      "' is wrong at i value___"
+      "' is wrong in Nvalues= ",wrongi,' at i value___"
     DO i=1,wrongi
       PRINT *,wronggridpt(i), wrongvalues(i)
@@ -318 +318 @@
 
 END SUBROUTINE check_var
+
+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
+
+  INTEGER, INTENT(IN)                                    :: sizev, zsize
+  CHARACTER(LEN=50), INTENT(IN)                          :: funcn, varn
+  REAL, DIMENSION(sizev), INTENT(IN)                     :: var
+  REAL, INTENT(IN)                                       :: bigvalue
+  LOGICAL, INTENT(IN)                                    :: stoprun
+
+! Local
+  INTEGER                                                :: i, wrongi
+  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 k value___"
+    DO i=1,wrongi
+      PRINT *,wronggridpt(i,:), wrongvalues(i)
+    END DO
+    IF (stoprun) THEN
+      STOP
+    END IF
+  END IF
+
+  RETURN
+
+END SUBROUTINE check_var3D
+
 
 !C======================================================================
@@ -611 +673 @@
             varname = 'zh_dair_col'
             CALL check_var(fname, varname, zh_dair_col, klon, largest*10.e6, .FALSE.)
+            PRINT *,'      zh_dair_col = func(q,ql,qs,paprs)'
+            varname = 'q'
+            CALL check_var3D(fname, varname, q, klon, klev, largest, .FALSE.)
+            varname = 'ql'
+            CALL check_var3D(fname, varname, ql, klon, klev, largest, .FALSE.)
+            varname = 'qs'
+            CALL check_var3D(fname, varname, qs, klon, klev, largest, .FALSE.)
+            varname = 'paprs'
+            CALL check_var3D(fname, varname, paprs, klon, klev, largest*10.e6, .FALSE.)
           END IF
           PRINT *,'    h_qw_tot= ',h_qw_tot
@@ -631 +702 @@
       END IF
 
+      IF (qs_tot .NE. qs_tot .OR. ABS(qs_tot) > largest) THEN
+        PRINT *,TRIM(errmsg)
+        PRINT *,'  ' // TRIM(fname) // ': Wrong qs_tot= ',qs_tot,' !!!'
+        PRINT *,'    qs: total mass of solid watter (kg/m2)'
+        PRINT *,'    qs = zqs_col(i)*airephy(i)'
+        varname = 'airephy'
+        CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
+        varname = 'zqs_col'
+        CALL check_var(fname, varname, zqs_col, klon, largest*10.e5, .FALSE.)
+        STOP
+      END IF
       RETURN 
       END SUBROUTINE diagetpq