Changeset 2110 for LMDZ5/trunk/libf/dyn3dpar

Timestamp:

Aug 27, 2014, 5:54:44 PM (10 years ago)

Author:

lguez

Message:

Abort if surface pressure becomes negative. The call to abort_gcm
already was in the sequential version but not in dyn3dpar nor in
dyn3dmem. In dyn3dmem, there is a variable checksum_all, which is
always true. The call to MPI_ALLREDUCE which should update
checksum_all is commented out. I do not know why (performance?).

Non-ASCII characters in comments are not always rendered properly and
they risk being lost. See revision [1740].

Bug fix in procedure convect2: the dimension len must be declared
before the array idcum which has this dimension. Bug fix in procedure
icefrac: the dimensions nl and len must be declared before the arrays
which have these dimensions.

File:

• LMDZ5/trunk/libf/dyn3dpar/integrd_p.F (modified) (2 diffs)

LMDZ5/trunk/libf/dyn3dpar/integrd_p.F

@@ -137 +137 @@
        
         IF( .NOT. checksum ) THEN
-         write(lunout,*) "integrd: negative surface pressure ",
-     &                                                ps(stop_it)
+           write(lunout,*) "integrd: ps = ", ps(stop_it)
          write(lunout,*) " at node ij =", stop_it
          ! since ij=j+(i-1)*jjp1 , we have
@@ -145 +144 @@
          write(lunout,*) " lon = ",rlonv(i)*180./pi, " deg",
      &                   " lat = ",rlatu(j)*180./pi, " deg"
+         call abort_gcm("integrd_p", "negative surface pressure", 1)
         ENDIF