/[MITgcm]/MITgcm/eesupp/src/eedie.F
ViewVC logotype

Diff of /MITgcm/eesupp/src/eedie.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.6 by heimbach, Tue Apr 10 22:35:24 2001 UTC revision 1.13 by jmc, Thu Oct 11 19:15:18 2012 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3    
4  #include "CPP_EEOPTIONS.h"  #include "CPP_EEOPTIONS.h"
5    #ifdef USE_LIBHPM
6    # include "f_hpm.h"
7    #endif
8    
9  CStartOfInterface  CBOP
10        SUBROUTINE EEDIE        SUBROUTINE EEDIE
11  C     /==========================================================\  C     *==========================================================*
12  C     | SUBROUTINE EEDIE                                         |  C     | SUBROUTINE EEDIE                                         |
13  C     | o Close execution "environment", particularly perform    |  C     | o Close execution "environment", particularly perform    |
14  C     |   steps to terminate parallel processing.                |  C     |   steps to terminate parallel processing.                |
15  C     |==========================================================|  C     *==========================================================*
16  C     | Note: This routine can also be compiled with CPP         |  C     | Note: This routine can also be compiled with CPP         |
17  C     | directives set so that no multi-processing is initialised|  C     | directives set so that no multi-processing is initialised|
18  C     | This is OK and should work fine.                         |  C     | This is OK and should work fine.                         |
19  C     \==========================================================/  C     *==========================================================*
20        IMPLICIT NONE        IMPLICIT NONE
21    
22  C     == Global variables ==  C     == Global variables ==
23  #include "SIZE.h"  #include "SIZE.h"
24  #include "EEPARAMS.h"  #include "EEPARAMS.h"
25  #include "EESUPPORT.h"  #include "EESUPPORT.h"
26  CEndOfInterface  CEOP
27    
28  C     == Local variables ==  C     == Local variables ==
29  C     msgBuf       - I/O Buffer  C     msgBuf       :: I/O Buffer
30  C     nThreadsDone - Used to count number of completed threads.  C     nThreadsDone :: Used to count number of completed threads.
31  C     I            - Loop counter.  C     I            :: Loop counter.
32        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
33        INTEGER nThreadsDone        INTEGER nThreadsDone
34        INTEGER I        INTEGER I
35  #ifdef ALLOW_USE_MPI  #ifdef ALLOW_USE_MPI
36  C     mpiRC            - Error code reporting variable used  C     mpiRC        :: Error code reporting variable used with MPI.
 C                        with MPI.  
37        INTEGER mpiRC        INTEGER mpiRC
38  #endif /* ALLOW_USE_MPI */  #endif /* ALLOW_USE_MPI */
39    
40          IF ( eeBootError ) THEN
41    C--   Skip ended threads counting if earlier error was found
42            WRITE(msgBuf,'(2A)')
43         &   'EEDIE: earlier error in multi-proc/thread setting'
44            CALL PRINT_ERROR( msgBuf, 1 )
45            fatalError = .TRUE.
46    
47          ELSE
48  C--   Check that all the threads have ended  C--   Check that all the threads have ended
49  C     No thread should reach this loop before all threads have set  C     No thread should reach this loop before all threads have set
50  C     threadIsComplete to TRUE. If they do then either there is a bug  C     threadIsComplete to TRUE. If they do then either there is a bug
51  C     in the code or the behaviour of the parallel compiler directives  C     in the code or the behaviour of the parallel compiler directives
52  C     are not right for this code. In the latter case different  C     are not right for this code. In the latter case different
53  C     directives may be available or the compiler itself may have a  C     directives may be available or the compiler itself may have a
54  C     bug or you may need a different parallel compiler for main.F  C     bug or you may need a different parallel compiler for main.F
55        nThreadsDone = 0          nThreadsDone = 0
56        DO I = 1, nThreads          DO I = 1, nThreads
57         IF ( threadIsComplete(I) ) nThreadsDone = nThreadsDone+1           IF ( threadIsComplete(I) ) nThreadsDone = nThreadsDone+1
58        ENDDO          ENDDO
59        IF ( nThreadsDone .LT. nThreads ) THEN          IF ( nThreadsDone .LT. nThreads ) THEN
60         WRITE(msgBuf,'(A,I5,A,I5,A)')           WRITE(msgBuf,'(A,I5,A)')
61       & 'S/R EEDIE: Only ',nThreadsDone,       &    'S/R EEDIE: Only',nThreadsDone,' threads have completed,'
62       & ' threads have completed, ',nThreads,           CALL PRINT_ERROR( msgBuf, 1 )
63       & ' are expected for this configuration!'           WRITE(msgBuf,'(A,I5,A)')
64         WRITE(*,*) msgBuf       &    'S/R EEDIE:',nThreads,' are expected for this config !'
65         WRITE(*,*)           CALL PRINT_ERROR( msgBuf, 1 )
66       & 'Possibly you have different setenv PARALLEL and nThreads?'           eeEndError = .TRUE.
67         eeEndError = .TRUE.           fatalError = .TRUE.
68         fatalError = .TRUE.          ENDIF
69    
70    C--   end if/else eebootError
71        ENDIF        ENDIF
72    
73    #ifdef USE_LIBHPM
74          CALL F_HPMTERMINATE(myProcId)
75    #endif
76    
77    C--   Flush IO-unit before MPI termination
78          CALL MDS_FLUSH( errorMessageUnit, 1 )
79    c#ifdef ALLOW_USE_MPI
80          CALL MDS_FLUSH( standardMessageUnit, 1 )
81    c#endif /* ALLOW_USE_MPI */
82    
83  #ifdef ALLOW_USE_MPI  #ifdef ALLOW_USE_MPI
84  C--  C- Note: since MPI_INIT is always called, better to also always terminate MPI
85    C        (even if usingMPI=F) --> comment out test on usingMPI
86    c     IF ( usingMPI ) THEN
87    
88  C--   MPI style multiple-process termination  C--   MPI style multiple-process termination
89  C--   ======================================  C--   ======================================
90  #ifndef ALWAYS_USE_MPI  #ifdef COMPONENT_MODULE
91        IF ( usingMPI ) THEN         IF ( useCoupler) CALL MPI_BARRIER( MPI_COMM_WORLD, mpiRC )
92    #endif
93    #ifdef ALLOW_OASIS
94           IF ( useOASIS ) CALL OASIS_FINALIZE
95  #endif  #endif
96         CALL MPI_FINALIZE  ( mpiRC )         CALL MPI_FINALIZE  ( mpiRC )
97         IF ( mpiRC .NE. MPI_SUCCESS ) THEN         IF ( mpiRC .NE. MPI_SUCCESS ) THEN
# Line 74  C--   ================================== Line 102  C--   ==================================
102       &       mpiRC       &       mpiRC
103          CALL PRINT_ERROR( msgBuf, 1 )          CALL PRINT_ERROR( msgBuf, 1 )
104         ENDIF         ENDIF
 C  
 #ifndef ALWAYS_USE_MPI  
       ENDIF  
 #endif  
 #endif /* ALLOW_USE_MPI */  
105    
106    c     ENDIF
107    #endif /* ALLOW_USE_MPI */
108    
109        RETURN        RETURN
110        END        END
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.22