--- MITgcm_contrib/plumes/external_forcing.F 2004/05/13 22:21:45 1.1 +++ MITgcm_contrib/plumes/external_forcing.F 2004/05/25 18:11:58 1.2 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/plumes/external_forcing.F,v 1.1 2004/05/13 22:21:45 molod Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/plumes/external_forcing.F,v 1.2 2004/05/25 18:11:58 molod Exp $ C $Name: $ #include "PACKAGES_CONFIG.h" @@ -74,6 +74,12 @@ #endif /* ALLOW_FIZHI */ C AMM +C AMM + CALL PLUME_TENDENCY_APPLY_U( + & iMin,iMax, jMin,jMax, bi,bj, kLev, + & myCurrentTime, myThid ) +C AMM + C Add windstress momentum impulse into the top-layer IF ( kLev .EQ. kSurface ) THEN DO j=jMin,jMax @@ -162,6 +168,13 @@ & myCurrentTime, myThid ) #endif /* ALLOW_FIZHI */ C AMM + +C AMM + CALL PLUME_TENDENCY_APPLY_V( + & iMin,iMax, jMin,jMax, bi,bj, kLev, + & myCurrentTime, myThid ) +C AMM + C Add windstress momentum impulse into the top-layer IF ( kLev .EQ. kSurface ) THEN DO j=jMin,jMax @@ -258,6 +271,13 @@ & myCurrentTime, myThid ) #endif /* ALLOW_FIZHI */ C AMM + +C AMM + CALL PLUME_TENDENCY_APPLY_T( + & iMin,iMax, jMin,jMax, bi,bj, kLev, + & myCurrentTime, myThid ) +C AMM + C Add heat in top-layer IF ( kLev .EQ. kSurface ) THEN DO j=jMin,jMax