pkg/atm_ocn_coupler/coupler.F

C $Header: /u/gcmpack/MITgcm/pkg/aim_ocn_coupler/coupler.F,v 1.2 2004/05/21 19:31:50 jmc Exp $
C $Name:  $

      PROGRAM COUPLER
C     /==========================================================\
C     | PROGRAM COUPLER                                          |
C     | o Main routine for 'Coupler' component. 'Coupler'        |
C     |   component coordiantes the exchange of data between     |
C     |   component models in a coupled model experiment.        |
C     |==========================================================|
C     | This version uses the MIT Coupler "checkpoint1" library  |
C     | calls.                                                   |
C     \==========================================================/
      IMPLICIT NONE
#include "mpif.h"
#include "CPL_PARAMS.h"

C     == Local variables ==
C     I                - Loop counter
C     MPI_COMM_Coupler
C     rc               - MPI return code
      INTEGER I
      INTEGER rc

C     Initialise the coupler component
      CALL INITIALISE

C     Perform registration with other components
      CALL ACCEPT_COMPONENT_REGISTRATIONS

C     Coordinate the transfer configuration information 
C     between components
      CALL RECEIVE_COMPONENT_CONFIGS
      CALL SENDOUT_COMPONENT_CONFIGS

      DO I=1,nCouplingSteps

C      Receive updated state
       CALL CPL_RECV_OCN_FIELDS
       CALL CPL_RECV_ATM_FIELDS

C      Send out fields
       CALL CPL_SEND_ATM_FIELDS
       CALL CPL_SEND_OCN_FIELDS

      ENDDO

C     o Finalize MPI
C     First wait for everybody to finish. Nobody should call
C     MPI_Finalize before all the component modules are
C     ready to finish. On some systems once one participant
C     gets to MPI_Finalize then its unclear what will
C     happen after that. If everybody does on MPI_Barrier
C     on COMM_WORLD then we will be OK.
      CALL MPI_BARRIER( MPI_COMM_WORLD, rc )
      CALL MPI_FINALIZE(rc)

      STOP
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/aim_ocn_coupler/coupler.F,v 1.2 2004/05/21 19:31:50 jmc Exp $
2	C $Name: $
3
4	PROGRAM COUPLER
5	C /==========================================================\
6	C \| PROGRAM COUPLER \|
7	C \| o Main routine for 'Coupler' component. 'Coupler' \|
8	C \| component coordiantes the exchange of data between \|
9	C \| component models in a coupled model experiment. \|
10	C \|==========================================================\|
11	C \| This version uses the MIT Coupler "checkpoint1" library \|
12	C \| calls. \|
13	C \==========================================================/
14	IMPLICIT NONE
15	#include "mpif.h"
16	#include "CPL_PARAMS.h"
17
18	C == Local variables ==
19	C I - Loop counter
20	C MPI_COMM_Coupler
21	C rc - MPI return code
22	INTEGER I
23	INTEGER rc
24
25	C Initialise the coupler component
26	CALL INITIALISE
27
28	C Perform registration with other components
29	CALL ACCEPT_COMPONENT_REGISTRATIONS
30
31	C Coordinate the transfer configuration information
32	C between components
33	CALL RECEIVE_COMPONENT_CONFIGS
34	CALL SENDOUT_COMPONENT_CONFIGS
35
36	DO I=1,nCouplingSteps
37
38	C Receive updated state
39	CALL CPL_RECV_OCN_FIELDS
40	CALL CPL_RECV_ATM_FIELDS
41
42	C Send out fields
43	CALL CPL_SEND_ATM_FIELDS
44	CALL CPL_SEND_OCN_FIELDS
45
46	ENDDO
47
48	C o Finalize MPI
49	C First wait for everybody to finish. Nobody should call
50	C MPI_Finalize before all the component modules are
51	C ready to finish. On some systems once one participant
52	C gets to MPI_Finalize then its unclear what will
53	C happen after that. If everybody does on MPI_Barrier
54	C on COMM_WORLD then we will be OK.
55	CALL MPI_BARRIER( MPI_COMM_WORLD, rc )
56	CALL MPI_FINALIZE(rc)
57
58	STOP
59	END