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	cnh	1.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