Changeset 2769 for LMDZ5/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Jan 24, 2017, 3:22:44 PM (7 years ago)

Author:

lguez

Message:

Added test on pressure levels. It is strange that this test could be
useful. paprs is input from the dynamical part of the code, and is
computed simply from surface pressure and hybrid coefficients. So the
test should always pass, unless there is a NaN in surface
pressure. Experience shows that this test, which is also made deep in
the ozone regridding procedures, sometimes fails, usually because of a
damaged restart file.

File:

• LMDZ5/trunk/libf/phylmd/physiq_mod.F90 (modified) (2 diffs)

LMDZ5/trunk/libf/phylmd/physiq_mod.F90

@@ -16 +16 @@
        d_u, d_v, d_t, d_qx, d_ps)
 
+    use assert_m, only: assert
     USE ioipsl, only: histbeg, histvert, histdef, histend, histsync, &
          histwrite, ju2ymds, ymds2ju, getin
@@ -1125 +1126 @@
     ENDIF
 
-    !======================================================================
+    ! Quick check on pressure levels:
+    do i = 1, klon
+       call assert(paprs(i, nbp_lev + 1) < paprs(i, nbp_lev), &
+            "physiq_mod paprs bad order")
+    end do
 
     IF (first) THEN