--- MITgcm/model/src/calc_gw.F	2000/03/24 18:31:31	1.3
+++ MITgcm/model/src/calc_gw.F	2001/03/06 17:10:29	1.7
@@ -1,3 +1,5 @@
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/calc_gw.F,v 1.7 2001/03/06 17:10:29 jmc Exp $
+C $Name:  $
 #include "CPP_OPTIONS.h"
 
       SUBROUTINE CALC_GW(     
@@ -14,7 +16,6 @@
 #include "EEPARAMS.h"
 #include "PARAMS.h"
 #include "GRID.h"
-#include "CG2D.h"
 #include "GW.h"
 #include "CG3D.h"
 
@@ -40,7 +41,7 @@
       _RL    flx_Dn(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
       _RL    flx_Up(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
 C     I,J,K - Loop counters
-      INTEGER i,j,k, kP1
+      INTEGER i,j,k, kP1, kUp
       _RL  wOverride
       _RS  hFacROpen
       _RS  hFacRClosed
@@ -177,6 +178,21 @@
        ENDDO
       ENDDO
 
+#ifdef ALLOW_OBCS
+      IF (useOBCS) THEN
+C--   This call is aesthetic: it makes the W field
+C     consistent with the OBs but this has no algorithmic
+C     impact. This is purely for diagnostic purposes.
+       DO bj=myByLo(myThid),myByHi(myThid)
+        DO bi=myBxLo(myThid),myBxHi(myThid)
+         DO K=1,Nr
+          CALL OBCS_APPLY_W( bi, bj, K, wVel, myThid )
+         ENDDO
+        ENDDO
+       ENDDO
+      ENDIF
+#endif /* ALLOW_OBCS */
+
 #endif /* ALLOW_NONHYDROSTATIC */
 
       RETURN