Changeset 2110 for LMDZ5/trunk/libf/dyn3dmem/integrd_loc.F

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/dyn3dmem/integrd_loc.F (modified) (2 diffs)

LMDZ5/trunk/libf/dyn3dmem/integrd_loc.F

@@ -147 +147 @@
 !     &                   MPI_LOGICAL,MPI_LOR,COMM_LMDZ,ierr)
       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
@@ -155 +154 @@
 !         write(lunout,*) " lon = ",rlonv(i)*180./pi, " deg",
 !     &                   " lat = ",rlatu(j)*180./pi, " deg"
+         call abort_gcm("integrd_loc", "negative surface pressure", 1)
       ENDIF