Changeset 50 in lmdz_wrf for branches/LMDZ_WRFmeas

Timestamp:

Jul 23, 2014, 4:33:35 PM (10 years ago)

Author:

lfita

Message:

Adding similar checks in 'diagtpq'

File:

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

branches/LMDZ_WRFmeas/WRFV3/lmdz/diagphy_mod.F90

@@ -441 +441 @@
 !$OMP THREADPRIVATE(h_vcol_pre, h_dair_pre, h_qw_pre, h_ql_pre)
 !$OMP THREADPRIVATE(h_qs_pre, qw_pre, ql_pre, qs_pre , ec_pre)
+
+! L. Fita, LMD July 2014
+      CHARACTER(LEN=50)                                  :: errmsg, fname, varname
+      LOGICAL                                            :: found
+      REAL                                               :: largest
+
+      fname = 'diagetpq'
+      errmsg = 'ERROR -- error -- ERROR -- error'
+      largest = 10.e3
 !c======================================================================
 !C
@@ -583 +592 @@
       ec_pre (idiag)    = ec_tot
 !C
+      IF (d_h_vcol .NE. d_h_vcol .OR. ABS(d_h_vcol) > largest) THEN
+        PRINT *,TRIM(errmsg)
+        PRINT *,'  ' // TRIM(fname) // ': Wrong d_h_vcol= ',d_h_vcol,' !!!'
+        PRINT *,'    d_h_vcol: Heat flux (W/m2) define as the Enthalpy change' //    &
+          ' (J/m2) during one time step (dtime) for the whole atmosphere (air,' //   &
+          ' watter vapour, liquid and solid)'
+        PRINT *,'    d_h_vcol = (h_vcol_tot - h_vcol_pre(idiag2) )'
+        PRINT *,'    h_vcol_tot= ',h_vcol_tot
+        IF (h_vcol_tot .NE. h_vcol_tot .OR. ABS(h_vcol_tot) > largest) THEN
+          PRINT *,'      h_vcol_tot = h_dair_tot+h_qw_tot+h_ql_tot+h_qs_tot'
+          PRINT *,'    airetot= ',airetot
+          IF (airetot .NE. airetot .OR. ABS(airetot) > largest) THEN
+            varname = 'airephy'
+            CALL check_var(fname, varname, airephy, klon, largest, .FALSE.)
+          END IF
+          PRINT *,'    h_dair_tot= ',h_dair_tot
+          IF (h_dair_tot .NE. h_dair_tot .OR. ABS(h_dair_tot) > largest) THEN
+            varname = 'zh_dair_col'
+            CALL check_var(fname, varname, zh_dair_col, klon, largest, .FALSE.)
+          END IF
+          PRINT *,'    h_qw_tot= ',h_qw_tot
+          IF (h_qw_tot .NE. h_qw_tot .OR. ABS(h_qw_tot) > largest) THEN
+            varname = 'zh_qw_col'
+            CALL check_var(fname, varname, zh_qw_col, klon, largest, .FALSE.)
+          END IF
+          PRINT *,'    h_ql_tot= ',h_ql_tot
+          IF (h_ql_tot .NE. h_ql_tot .OR. ABS(h_ql_tot) > largest) THEN
+            varname = 'zh_ql_col'
+            CALL check_var(fname, varname, zh_ql_col, klon, largest, .FALSE.)
+          END IF
+          PRINT *,'    h_qs_tot= ',h_qs_tot
+          IF (h_qs_tot .NE. h_qs_tot .OR. ABS(h_qs_tot) > largest) THEN
+            varname = 'zh_qs_col'
+            CALL check_var(fname, varname, zh_qs_col, klon, largest, .FALSE.)
+          END IF
+        END IF
+        STOP
+      END IF
+
       RETURN 
       END SUBROUTINE diagetpq