Changeset 57 in lmdz_wrf

Timestamp:

Jul 23, 2014, 5:47:38 PM (10 years ago)

Author:

lfita

Message:

Adding for checking qw_tot, ql_tot, qs_tot, ec_tot

File:

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

branches/LMDZ_WRFmeas/WRFV3/lmdz/diagphy_mod.F90

@@ -704 +704 @@
       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)'
+      IF (qw_tot .NE. qw_tot .OR. ABS(qw_tot) > largest) THEN
+        PRINT *,TRIM(errmsg)
+        PRINT *,'  ' // TRIM(fname) // ': Wrong qw_tot= ',qw_tot,' !!!'
+        PRINT *,'    qw_tot: total mass of watter vapour (kg/m2)'
+        PRINT *,'    qw_tot = f(q,ql,qs,airephy,t,paprs)'
         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
+        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 = 't'
+        CALL check_var3D(fname, varname, t, klon, klev, largest, .FALSE.)
+        varname = 'paprs'
+        CALL check_var3D(fname, varname, paprs, klon, klev, largest, .FALSE.)
+        STOP
+      END IF
+
+      IF (ql_tot .NE. ql_tot .OR. ABS(ql_tot) > largest) THEN
+        PRINT *,TRIM(errmsg)
+        PRINT *,'  ' // TRIM(fname) // ': Wrong ql_tot= ',ql_tot,' !!!'
+        PRINT *,'    ql_tot: total mass of liquid watter (kg/m2)'
+        PRINT *,'    ql_tot = f(ql,airephy,t,paprs)'
+        varname = 'airephy'
+        CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
+        varname = 'ql'
+        CALL check_var3D(fname, varname, ql, klon, klev, largest, .FALSE.)
+        varname = 't'
+        CALL check_var3D(fname, varname, t, klon, klev, largest, .FALSE.)
+        varname = 'paprs'
+        CALL check_var3D(fname, varname, paprs, klon, klev, largest, .FALSE.)
+        STOP
+      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_tot: total mass of solid watter (kg/m2)'
+        PRINT *,'    qs_tot = f(qs,airephy,t,paprs)'
+        varname = 'airephy'
+        CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
+        varname = 'qs'
+        CALL check_var3D(fname, varname, qs, klon, klev, largest, .FALSE.)
+        varname = 't'
+        CALL check_var3D(fname, varname, t, klon, klev, largest, .FALSE.)
+        varname = 'paprs'
+        CALL check_var3D(fname, varname, paprs, klon, klev, largest, .FALSE.)
+        STOP
+      END IF
+
+      IF (ec_tot .NE. ec_tot .OR. ABS(ec_tot) > largest) THEN
+        PRINT *,TRIM(errmsg)
+        PRINT *,'  ' // TRIM(fname) // ': Wrong ec_tot= ',ec_tot,' !!!'
+        PRINT *,'    ec_tot: total cinetic energy (kg/m2)'
+        PRINT *,'    ec_tot = f(u,v,airephy,t,paprs)'
+        varname = 'airephy'
+        CALL check_var(fname, varname, airephy, klon, largest*10.e5, .FALSE.)
+        varname = 'u'
+        CALL check_var3D(fname, varname, u, klon, klev, largest, .FALSE.)
+        varname = 'v'
+        CALL check_var3D(fname, varname, v, klon, klev, largest, .FALSE.)
+        varname = 't'
+        CALL check_var3D(fname, varname, t, klon, klev, largest, .FALSE.)
+        varname = 'paprs'
+        CALL check_var3D(fname, varname, paprs, klon, klev, largest, .FALSE.)
+        STOP
+      END IF
+
       RETURN 
+
       END SUBROUTINE diagetpq