dgoldberg/streamice/streamice_taub.F

C $Header: /u/gcmpack/MITgcm_contrib/dgoldberg/streamice/streamice_taub.F,v 1.1 2012/09/18 17:04:40 dgoldberg Exp $
C $Name:  $

#include "STREAMICE_OPTIONS.h"

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

CBOP
      SUBROUTINE STREAMICE_TAUB ( myThid )

C     /============================================================\
C     | SUBROUTINE                                                 |   
C     | o                                                          |
C     |============================================================|
C     |                                                            |
C     \============================================================/
      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "GRID.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "STREAMICE.h"
#include "STREAMICE_CG.h"

C     !INPUT/OUTPUT ARGUMENTS
      INTEGER myThid
!       _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
!       _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)

#ifdef ALLOW_STREAMICE


C     LOCAL VARIABLES
      INTEGER i, j, bi, bj, k, l, umid, vmid
      _RL ux, uy, vx, vy, exx, eyy, exy, second_inv, unorm

      DO bj=myByLo(myThid),myByHi(myThid)
       DO bi=myBxLo(myThid),myBxHi(myThid)
        DO j=1,sNy
         DO i=1,sNx
cph(
          streamice_taubx (i,j,bi,bj) = 0.
          streamice_tauby (i,j,bi,bj) = 0.
cph)
          IF (STREAMICE_hmask(i,j,bi,bj).eq.1) THEN
  
           streamice_taubx (i,j,bi,bj) = 
     &      tau_beta_eff_streamice(i,j,bi,bj) * U_streamice(i,j,bi,bj)
           streamice_tauby (i,j,bi,bj) = 
     &      tau_beta_eff_streamice(i,j,bi,bj) * V_streamice(i,j,bi,bj)
           
          ENDIF
         ENDDO
        ENDDO
       ENDDO
      ENDDO


#endif
      RETURN
      END
1	heimbach	1.2	C $Header: /u/gcmpack/MITgcm_contrib/dgoldberg/streamice/streamice_taub.F,v 1.1 2012/09/18 17:04:40 dgoldberg Exp $
2	dgoldberg	1.1	C $Name: $
3
4			#include "STREAMICE_OPTIONS.h"
5
6			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
7
8			CBOP
9			SUBROUTINE STREAMICE_TAUB ( myThid )
10
11			C /============================================================\
12			C \| SUBROUTINE \|
13			C \| o \|
14			C \|============================================================\|
15			C \| \|
16			C \============================================================/
17			IMPLICIT NONE
18
19			C === Global variables ===
20			#include "SIZE.h"
21			#include "GRID.h"
22			#include "EEPARAMS.h"
23			#include "PARAMS.h"
24			#include "STREAMICE.h"
25			#include "STREAMICE_CG.h"
26
27			C !INPUT/OUTPUT ARGUMENTS
28			INTEGER myThid
29			! _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
30			! _RL taudx (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
31
32			#ifdef ALLOW_STREAMICE
33
34
35			C LOCAL VARIABLES
36			INTEGER i, j, bi, bj, k, l, umid, vmid
37			_RL ux, uy, vx, vy, exx, eyy, exy, second_inv, unorm
38
39			DO bj=myByLo(myThid),myByHi(myThid)
40			DO bi=myBxLo(myThid),myBxHi(myThid)
41			DO j=1,sNy
42			DO i=1,sNx
43	heimbach	1.2	cph(
44			streamice_taubx (i,j,bi,bj) = 0.
45			streamice_tauby (i,j,bi,bj) = 0.
46			cph)
47	dgoldberg	1.1	IF (STREAMICE_hmask(i,j,bi,bj).eq.1) THEN
48
49			streamice_taubx (i,j,bi,bj) =
50			& tau_beta_eff_streamice(i,j,bi,bj) * U_streamice(i,j,bi,bj)
51			streamice_tauby (i,j,bi,bj) =
52			& tau_beta_eff_streamice(i,j,bi,bj) * V_streamice(i,j,bi,bj)
53
54			ENDIF
55			ENDDO
56			ENDDO
57			ENDDO
58			ENDDO
59
60
61			#endif
62			RETURN
63			END