C $Header: /u/gcmpack/MITgcm/pkg/thsice/thsice_get_velocity.F,v 1.3 2012/03/06 15:37:24 jmc Exp $
C $Name:  $

#include "THSICE_OPTIONS.h"
#ifdef ALLOW_SEAICE
# include "SEAICE_OPTIONS.h"
#endif /* ALLOW_SEAICE */


CBOP
C !ROUTINE: THSICE_GET_VELOCITY

C !INTERFACE: ==========================================================
      SUBROUTINE THSICE_GET_VELOCITY(
     O                  uLoc, vLoc,
     I                  bi, bj, myTime, myIter, myThid )

C !DESCRIPTION: \bv
C     *===========================================================*
C     | SUBROUTINE THSICE_GET_VELOCITY
C     | o load seaice velocity from pkg/seaice common block
C     *===========================================================*
C \ev

C !USES: ===============================================================
      IMPLICIT NONE

C     === Global variables ===

#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
c#include "THSICE_SIZE.h"
c#include "THSICE_PARAMS.h"
#ifdef ALLOW_SEAICE
# include "SEAICE_SIZE.h"
# include "SEAICE.h"
#endif /* ALLOW_SEAICE */

C !INPUT PARAMETERS: ===================================================
C     === Routine arguments ===
C     uLoc/vLoc :: current ice velocity on C-grid [m/s]
C     bi,bj     :: Tile indices
C     myTime    :: Current time in simulation (s)
C     myIter    :: Current iteration number
C     myThid    :: My Thread Id number
      _RL     uLoc(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      _RL     vLoc(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
      INTEGER bi,bj
      _RL     myTime
      INTEGER myIter
      INTEGER myThid

#ifdef ALLOW_THSICE
C !LOCAL VARIABLES: ====================================================
C     === Local variables ===
C     i,j,      :: Loop counters
      INTEGER i, j
CEOP

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|


#ifdef ALLOW_SEAICE
      IF ( useSEAICE ) THEN
        DO j=1-OLy,sNy+OLy
         DO i=1-OLx,sNx+OLx
          uLoc(i,j) = UICE(i,j,bi,bj)
          vLoc(i,j) = VICE(i,j,bi,bj)
         ENDDO
        ENDDO
      ELSE
#else /* ALLOW_SEAICE */
      IF ( .TRUE. ) THEN
#endif /* ALLOW_SEAICE */
C-    set ice velocity to zero
        DO j=1-OLy,sNy+OLy
         DO i=1-OLx,sNx+OLx
          uLoc(i,j) = 0.
          vLoc(i,j) = 0.
         ENDDO
        ENDDO
      ENDIF

#endif /* ALLOW_THSICE */

      RETURN
      END