pkg/aim/aim_aim2dyn_exchanges.F

C     $Header: /u/gcmpack/models/MITgcmUV/pkg/aim/aim_external_forcing.F,v 1.2 2001/02/02 21:36:29 adcroft Exp $
C     $Name:  $

#include "AIM_OPTIONS.h"

CStartOfInterface
      SUBROUTINE AIM_AIM2DYN_EXCHANGES(
     I           myCurrentTime,myThid)
C     /==========================================================\
C     | S/R AIM_AIM2DYN_EXCHANGES                                |
C     | o Do overlap exchanges for AIM related quantities.       |
C     |   AIM does not compute into overlaps so exchanges        |
C     |   are needed for AIM quantities.                         |
C     \==========================================================/
      IMPLICIT NONE

C     == Global data ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "DYNVARS.h"

#ifdef ALLOW_AIM
#include "AIM2DYN.h"
#endif /* ALLOW_AIM */

C     == Routine arguments ==
      _RL myCurrentTime
      INTEGER myThid
CEndOfInterface

C     == Local variables ==
C     I,J,K - Loop counters
C     bi,bj
      INTEGER I,J,K, bi,bj 

#ifdef ALLOW_AIM

      _EXCH_XY_R8( aim_drag, myThid )
      _EXCH_XYZ_R8( gT, myThid )
      _EXCH_XYZ_R8( gS, myThid )
      _EXCH_XYZ_R8( gTNM1, myThid )
      _EXCH_XYZ_R8( gSNM1, myThid )
      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
        DO K=1,Nr
         DO J=1-OLy,sNy+OLy
          DO I=1-OLx,sNx+OLx
           gT(i,j,k,bi,bj) = 
     &      maskC(i,j,k,bi,bj)*gT(i,j,k,bi,bj)
           gTNM1(i,j,k,bi,bj) = 
     &      maskC(i,j,k,bi,bj)*gTNM1(i,j,k,bi,bj)
           gS(i,j,k,bi,bj) = 
     &      maskC(i,j,k,bi,bj)*gS(i,j,k,bi,bj)
           gSNM1(i,j,k,bi,bj) = 
     &      maskC(i,j,k,bi,bj)*gSNM1(i,j,k,bi,bj)
          ENDDO
         ENDDO
        ENDDO
       ENDDO
      ENDDO

#endif /* ALLOW_AIM */

      RETURN
      END
1	C $Header: /u/gcmpack/models/MITgcmUV/pkg/aim/aim_external_forcing.F,v 1.2 2001/02/02 21:36:29 adcroft Exp $
2	C $Name: $
3
4	#include "AIM_OPTIONS.h"
5
6	CStartOfInterface
7	SUBROUTINE AIM_AIM2DYN_EXCHANGES(
8	I myCurrentTime,myThid)
9	C /==========================================================\
10	C \| S/R AIM_AIM2DYN_EXCHANGES \|
11	C \| o Do overlap exchanges for AIM related quantities. \|
12	C \| AIM does not compute into overlaps so exchanges \|
13	C \| are needed for AIM quantities. \|
14	C \==========================================================/
15	IMPLICIT NONE
16
17	C == Global data ==
18	#include "SIZE.h"
19	#include "EEPARAMS.h"
20	#include "PARAMS.h"
21	#include "GRID.h"
22	#include "DYNVARS.h"
23
24	#ifdef ALLOW_AIM
25	#include "AIM2DYN.h"
26	#endif /* ALLOW_AIM */
27
28	C == Routine arguments ==
29	_RL myCurrentTime
30	INTEGER myThid
31	CEndOfInterface
32
33	C == Local variables ==
34	C I,J,K - Loop counters
35	C bi,bj
36	INTEGER I,J,K, bi,bj
37
38	#ifdef ALLOW_AIM
39
40	_EXCH_XY_R8( aim_drag, myThid )
41	_EXCH_XYZ_R8( gT, myThid )
42	_EXCH_XYZ_R8( gS, myThid )
43	_EXCH_XYZ_R8( gTNM1, myThid )
44	_EXCH_XYZ_R8( gSNM1, myThid )
45	DO bj=myByLo(myThid),myByHi(myThid)
46	DO bi=myBxLo(myThid),myBxHi(myThid)
47	DO K=1,Nr
48	DO J=1-OLy,sNy+OLy
49	DO I=1-OLx,sNx+OLx
50	gT(i,j,k,bi,bj) =
51	& maskC(i,j,k,bi,bj)*gT(i,j,k,bi,bj)
52	gTNM1(i,j,k,bi,bj) =
53	& maskC(i,j,k,bi,bj)*gTNM1(i,j,k,bi,bj)
54	gS(i,j,k,bi,bj) =
55	& maskC(i,j,k,bi,bj)*gS(i,j,k,bi,bj)
56	gSNM1(i,j,k,bi,bj) =
57	& maskC(i,j,k,bi,bj)*gSNM1(i,j,k,bi,bj)
58	ENDDO
59	ENDDO
60	ENDDO
61	ENDDO
62	ENDDO
63
64	#endif /* ALLOW_AIM */
65
66	RETURN
67	END