dgoldberg/streamice/streamice_taub.F

C $Header: /u/gcmpack/MITgcm_contrib/dgoldberg/streamice/streamice_visc_beta.F,v 1.3 2012/09/04 21:11:44 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
          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	C $Header: /u/gcmpack/MITgcm_contrib/dgoldberg/streamice/streamice_visc_beta.F,v 1.3 2012/09/04 21:11:44 dgoldberg Exp $
2	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	IF (STREAMICE_hmask(i,j,bi,bj).eq.1) THEN
44
45	streamice_taubx (i,j,bi,bj) =
46	& tau_beta_eff_streamice(i,j,bi,bj) * U_streamice(i,j,bi,bj)
47	streamice_tauby (i,j,bi,bj) =
48	& tau_beta_eff_streamice(i,j,bi,bj) * V_streamice(i,j,bi,bj)
49
50	ENDIF
51	ENDDO
52	ENDDO
53	ENDDO
54	ENDDO
55
56
57	#endif
58	RETURN
59	END