pkg/thsice/thsice_get_velocity.F

a4eca6e929 Jean*0001 #include "THSICE_OPTIONS.h"
                0002 #ifdef ALLOW_SEAICE
                0003 # include "SEAICE_OPTIONS.h"
                0004 #endif /* ALLOW_SEAICE */
                0005
                0006
                0007 CBOP
                0008 C !ROUTINE: THSICE_GET_VELOCITY
                0009
                0010 C !INTERFACE: ==========================================================
                0011       SUBROUTINE THSICE_GET_VELOCITY(
                0012      O                  uLoc, vLoc,
                0013      I                  bi, bj, myTime, myIter, myThid )
                0014
                0015 C !DESCRIPTION: \bv
                0016 C     *===========================================================*
                0017 C     | SUBROUTINE THSICE_GET_VELOCITY
                0018 C     | o load seaice velocity from pkg/seaice common block
                0019 C     *===========================================================*
                0020 C \ev
                0021
                0022 C !USES: ===============================================================
                0023       IMPLICIT NONE
                0024
                0025 C     === Global variables ===
                0026
                0027 #include "SIZE.h"
                0028 #include "EEPARAMS.h"
                0029 #include "PARAMS.h"
                0030 c#include "THSICE_SIZE.h"
                0031 c#include "THSICE_PARAMS.h"
                0032 #ifdef ALLOW_SEAICE
a34cef4f76 Jean*0033 # include "SEAICE_SIZE.h"
a4eca6e929 Jean*0034 # include "SEAICE.h"
                0035 #endif /* ALLOW_SEAICE */
                0036
                0037 C !INPUT PARAMETERS: ===================================================
                0038 C     === Routine arguments ===
                0039 C     uLoc/vLoc :: current ice velocity on C-grid [m/s]
                0040 C     bi,bj     :: Tile indices
                0041 C     myTime    :: Current time in simulation (s)
                0042 C     myIter    :: Current iteration number
                0043 C     myThid    :: My Thread Id number
                0044       _RL     uLoc(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0045       _RL     vLoc(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
                0046       INTEGER bi,bj
                0047       _RL     myTime
                0048       INTEGER myIter
                0049       INTEGER myThid
                0050
                0051 #ifdef ALLOW_THSICE
                0052 C !LOCAL VARIABLES: ====================================================
                0053 C     === Local variables ===
                0054 C     i,j,      :: Loop counters
                0055       INTEGER i, j
                0056 CEOP
                0057
                0058 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0059
                0060
                0061 #ifdef ALLOW_SEAICE
                0062       IF ( useSEAICE ) THEN
                0063         DO j=1-OLy,sNy+OLy
                0064          DO i=1-OLx,sNx+OLx
f16ded4b41 Mart*0065           uLoc(i,j) = UICE(i,j,bi,bj)
                0066           vLoc(i,j) = VICE(i,j,bi,bj)
a4eca6e929 Jean*0067          ENDDO
                0068         ENDDO
                0069       ELSE
                0070 #else /* ALLOW_SEAICE */
                0071       IF ( .TRUE. ) THEN
                0072 #endif /* ALLOW_SEAICE */
                0073 C-    set ice velocity to zero
                0074         DO j=1-OLy,sNy+OLy
                0075          DO i=1-OLx,sNx+OLx
                0076           uLoc(i,j) = 0.
                0077           vLoc(i,j) = 0.
                0078          ENDDO
                0079         ENDDO
                0080       ENDIF
                0081
                0082 #endif /* ALLOW_THSICE */
                0083
                0084       RETURN
                0085       END