--- MITgcm/pkg/diagnostics/diagnostics_fill_state.F	2005/01/28 01:05:02	1.15
+++ MITgcm/pkg/diagnostics/diagnostics_fill_state.F	2005/02/10 15:01:34	1.16
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/diagnostics/diagnostics_fill_state.F,v 1.15 2005/01/28 01:05:02 jmc Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/diagnostics/diagnostics_fill_state.F,v 1.16 2005/02/10 15:01:34 jmc Exp $
 C $Name:  $
 
 #include "DIAG_OPTIONS.h"
@@ -162,6 +162,46 @@
          CALL DIAGNOSTICS_FILL(tmpMk,'UV_VEL_Z',0,Nr,0,1,1,myThid)
         ENDIF
       
+        IF ( DIAGNOSTICS_IS_ON('WU_VEL  ',myThid) ) THEN
+         DO bj = myByLo(myThid), myByHi(myThid)
+          DO bi = myBxLo(myThid), myBxHi(myThid)
+           DO K=1,Nr
+            km1 = MAX(k-1,1)
+            DO j = 1,sNy
+             DO i = 1,sNx+1
+              tmpMk(i,j,K,bi,bj) = 0.25 _d 0
+     &           *(uVel(i,km1,K,bi,bj)+uVel(i,j,K,bi,bj))
+     &           *(wVel(i-1,j,K,bi,bj)*rA(i-1,j,bi,bj)
+     &            +wVel( i ,j,K,bi,bj)*rA( i ,j,bi,bj)
+     &            )*recip_rAw(i,j,bi,bj)
+             ENDDO
+            ENDDO
+           ENDDO
+          ENDDO
+         ENDDO
+         CALL DIAGNOSTICS_FILL(tmpMk,'WU_VEL  ',0,Nr,0,1,1,myThid)
+        ENDIF
+
+        IF ( DIAGNOSTICS_IS_ON('WV_VEL  ',myThid) ) THEN
+         DO bj = myByLo(myThid), myByHi(myThid)
+          DO bi = myBxLo(myThid), myBxHi(myThid)
+           DO K=1,Nr
+            km1 = MAX(k-1,1)
+            DO j = 1,sNy+1
+             DO i = 1,sNx
+              tmpMk(i,j,K,bi,bj) = 0.25 _d 0
+     &           *(vVel(i,km1,K,bi,bj)+vVel(i,j,K,bi,bj))
+     &           *(wVel(i,j-1,K,bi,bj)*rA(i,j-1,bi,bj)
+     &            +wVel(i, j ,K,bi,bj)*rA(i, j ,bi,bj)
+     &            )*recip_rAs(i,j,bi,bj)
+             ENDDO
+            ENDDO
+           ENDDO
+          ENDDO
+         ENDDO
+         CALL DIAGNOSTICS_FILL(tmpMk,'WV_VEL  ',0,Nr,0,1,1,myThid)
+        ENDIF
+
 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
 
         IF ( DIAGNOSTICS_IS_ON('UVELTH  ',myThid) ) THEN