--- MITgcm/pkg/aim/aim_external_forcing.F	2001/01/26 16:23:18	1.1
+++ MITgcm/pkg/aim/aim_external_forcing.F	2001/02/02 21:36:29	1.2
@@ -0,0 +1,261 @@
+C     $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/aim/Attic/aim_external_forcing.F,v 1.2 2001/02/02 21:36:29 adcroft Exp $
+C     $Name:  $
+
+#include "AIM_OPTIONS.h"
+
+CStartOfInterface
+      SUBROUTINE AIM_EXTERNAL_FORCING_U(
+     I           iMin, iMax, jMin, jMax,bi,bj,kLev,
+     I           myCurrentTime,myThid)
+C     /==========================================================\
+C     | S/R AIM_EXTERNAL_FORCING_U                               |
+C     | o Add AIM tendency terms to U tendency.                  |
+C     \==========================================================/
+
+C     == Global data ==
+#include "SIZE.h"
+#include "EEPARAMS.h"
+#include "PARAMS.h"
+#include "GRID.h"
+#include "DYNVARS.h"
+
+#ifdef ALLOW_AIM
+#include "atparam0.h"
+#include "atparam1.h"
+      INTEGER NGP
+      INTEGER NLON
+      INTEGER NLAT
+      INTEGER NLEV
+      PARAMETER ( NLON=IX, NLAT=IL, NLEV=KX, NGP=NLON*NLAT )
+#include "com_physvar.h"
+#endif /* ALLOW_AIM */
+
+C     == Routine arguments ==
+C     iMin - Working range of tile for applying forcing.
+C     iMax
+C     jMin
+C     jMax
+C     kLev
+      INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
+      _RL myCurrentTime
+      INTEGER myThid
+CEndOfInterface
+
+#ifdef ALLOW_AIM
+C     == Local variables ==
+C     Loop counters
+      INTEGER I, J
+      INTEGER I2, I3, kAtm
+
+      DO J=1,sNy
+       DO I=1,sNx
+        I2 = sNx*(J-1)+I
+        I3 = sNx*(J-1)+mod(I+sNx-2,sNx)+1
+         Katm = _KD2KA( Klev )
+         gU(i,j,kLev,bi,bj) =
+     &   gU(i,j,kLev,bi,bj) + 0.5*(UT_PBL(I2,katm)+UT_PBL(I3,katm))
+       ENDDO
+      ENDDO
+
+      _EXCH_XYZ_R8( gU , myThid) 
+
+#endif /* ALLOW_AIM */
+
+      RETURN
+      END
+CStartOfInterface
+      SUBROUTINE AIM_EXTERNAL_FORCING_V(
+     I           iMin, iMax, jMin, jMax,bi,bj,kLev,
+     I           myCurrentTime,myThid)
+C     /==========================================================\
+C     | S/R EXTERNAL_FORCING_V                                   |
+C     | o Add AIM tendency to meridional velocity.               |
+C     \==========================================================/
+      IMPLICIT rEAL*8 (A-H,O-Z)
+
+C     == Global data ==
+#include "SIZE.h"
+#include "EEPARAMS.h"
+#include "PARAMS.h"
+#include "GRID.h"
+#include "DYNVARS.h"
+#include "FFIELDS.h"
+
+#ifdef ALLOW_AIM
+#include "atparam0.h"
+#include "atparam1.h"
+      INTEGER NGP
+      INTEGER NLON
+      INTEGER NLAT
+      INTEGER NLEV
+      PARAMETER ( NLON=IX, NLAT=IL, NLEV=KX, NGP=NLON*NLAT )
+#include "com_physvar.h"
+#endif /* ALLOW_AIM */
+
+C     == Routine arguments ==
+C     iMin - Working range of tile for applying forcing.
+C     iMax
+C     jMin
+C     jMax
+C     kLev
+      INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
+      _RL myCurrentTime
+      INTEGER myThid
+CEndOfInterface
+
+#ifdef ALLOW_AIM
+C     == Local variables ==
+C     Loop counters
+      INTEGER I, J
+      INTEGER I2, I3, JN, kAtm
+
+      DO J=1,sNy
+       DO I=1,sNx
+        I2 = sNx*(J-1)+I
+        JN = J
+        IF ( JN .EQ. sNy ) JN = sNy-1
+        I3 = sNx*(JN )+I
+         Katm = _KD2KA( Klev )
+         gV(i,j,kLev,bi,bj) =
+     &   gV(i,j,kLev,bi,bj) + 0.5*(VT_PBL(I2,katm)+VT_PBL(I3,katm))
+       ENDDO
+      ENDDO
+
+      _EXCH_XYZ_R8( gV , myThid) 
+
+#endif /* ALLOW_AIM */
+
+      RETURN
+      END
+CStartOfInterface
+      SUBROUTINE AIM_EXTERNAL_FORCING_T(
+     I           iMin, iMax, jMin, jMax,bi,bj,kLev,
+     I           myCurrentTime,myThid)
+C     /==========================================================\
+C     | S/R AIM_EXTERNAL_FORCING_T                               |
+C     | o Add AIM tendency to T                                  | 
+C     \==========================================================/
+
+C     == Global data ==
+#include "SIZE.h"
+#include "EEPARAMS.h"
+#include "PARAMS.h"
+#include "GRID.h"
+#include "DYNVARS.h"
+
+#ifdef ALLOW_AIM
+#include "atparam0.h"
+#include "atparam1.h"
+      INTEGER NGP
+      INTEGER NLON
+      INTEGER NLAT
+      INTEGER NLEV
+      PARAMETER ( NLON=IX, NLAT=IL, NLEV=KX, NGP=NLON*NLAT )
+#include "com_physvar.h"
+#endif
+
+C     == Routine arguments ==
+C     iMin - Working range of tile for applying forcing.
+C     iMax
+C     jMin
+C     jMax
+C     kLev
+      INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
+      _RL myCurrentTime
+      INTEGER myThid
+CEndOfInterface
+
+#ifdef ALLOW_AIM
+C     == Local variables ==
+C     Loop counters
+      INTEGER I, J
+
+C--   Forcing term
+      _RL thetaLim,kT,ka,ks,term1,term2,thetaEq,termP,rSurf
+      _RL pGround
+
+C--   Forcing term
+      pGround = 1. _d 5
+      RD      = 287. _d 0
+      CPAIR   = 1004. _d 0
+      DO J=1,sNy
+       DO I=1,sNx
+        I2 = sNx*(J-1)+I
+        katm = _KD2KA( Klev )
+        gT(i,j,kLev,bi,bj) = gT(i,j,kLev,bi,bj) +(
+     &                           TT_PBL(I2,katm)
+     &                          +TT_CNV(I2,katm)
+     &                          +TT_LSC(I2,katm)
+     &                          +TT_RSW(I2,katm)
+     &                          +TT_RLW(I2,katm)
+     &        )*((pGround/rC(kLev))**(RD/CPAIR))
+       ENDDO
+      ENDDO
+
+      _EXCH_XYZ_R8( gT , myThid) 
+
+#endif /* ALLOW_AIM */
+
+      RETURN
+      END
+CStartOfInterface
+      SUBROUTINE AIM_EXTERNAL_FORCING_S(
+     I           iMin, iMax, jMin, jMax,bi,bj,kLev,
+     I           myCurrentTime,myThid)
+C     /==========================================================\
+C     | S/R AIM_EXTERNAL_FORCING_S                               |
+C     | o Add AIM tendency to S.                                 |
+C     \==========================================================/
+
+C     == Global data ==
+#include "SIZE.h"
+#include "EEPARAMS.h"
+#include "PARAMS.h"
+#include "GRID.h"
+#include "DYNVARS.h"
+
+#ifdef ALLOW_AIM
+#include "atparam0.h"
+#include "atparam1.h"
+      INTEGER NGP
+      INTEGER NLON
+      INTEGER NLAT
+      INTEGER NLEV
+      PARAMETER ( NLON=IX, NLAT=IL, NLEV=KX, NGP=NLON*NLAT )
+#include "com_physvar.h"
+#endif
+
+
+C     == Routine arguments ==
+C     iMin - Working range of tile for applying forcing.
+C     iMax
+C     jMin
+C     jMax
+C     kLev
+      INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
+      _RL myCurrentTime
+      INTEGER myThid
+CEndOfInterface
+
+#ifdef ALLOW_AIM
+C     == Local variables ==
+C     Loop counters
+      INTEGER I, J
+
+      DO J=1,sNy
+       DO I=1,sNx
+        katm = _KD2KA( kLev )
+        I2 = sNx*(J-1)+I
+            gS(i,j,kLev,bi,bj) = gS(i,j,kLev,bi,bj)
+     &                       +QT_PBL(I2,katm)
+     &                       +QT_CNV(I2,katm)
+     &                       +QT_LSC(I2,katm)
+       ENDDO
+      ENDDO
+
+      _EXCH_XYZ_R8( gS , myThid) 
+
+#endif /* ALLOW_AIM */
+
+      RETURN
+      END