pkg/ggl90/ggl90_calc_visc.F

C $Header: /u/gcmpack/MITgcm/pkg/ggl90/ggl90_calc_visc.F,v 1.2 2007/07/31 23:01:58 ce107 Exp $
C $Name:  $
#include "GGL90_OPTIONS.h"

      SUBROUTINE GGL90_CALC_VISC( 
     I        bi,bj,iMin,iMax,jMin,jMax,K,
     U        KappaRU,KappaRV,
     I        myThid)

CBOP
C     /==========================================================\
C     | SUBROUTINE GGL90_CALC_VISC                               |
C     | o Add contrubution to net viscosity from GGL90 mixing    |
C     \==========================================================/
      IMPLICIT NONE

C     == GLobal variables ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "GRID.h"

C     == Routine arguments ==
C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
C     myThid - Instance number for this innvocation of GGL90_CALC_VISC
C
      INTEGER bi,bj,iMin,iMax,jMin,jMax,K
      _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
      _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
      INTEGER myThid

#ifdef ALLOW_GGL90
#include "GGL90.h"

C     == Local variables ==
C     I, J, K - Loop counters
      INTEGER i,j
CEOP
      DO j=jMin,jMax
       DO i=iMin,iMax
        KappaRU(i,j,k) = KappaRU(i,j,k) + _maskW(i,j,k,bi,bj) *
     &                   (.5 _d 0*(GGL90viscAr(i,j,k,bi,bj)
     &                            +GGL90viscAr(i-1,j,k,bi,bj))
     &                    -viscAr)
       ENDDO
      ENDDO

      DO j=jMin,jMax
       DO i=iMin,iMax
        KappaRV(i,j,k) = KappaRV(i,j,k) + _maskS(i,j,k,bi,bj) *
     &                   (.5 _d 0*(GGL90viscAr(i,j,k,bi,bj)
     &                            +GGL90viscAr(i,j-1,k,bi,bj))
     &                    -viscAr)
       ENDDO
      ENDDO

#endif /* ALLOW_GGL90 */

      RETURN
      END
1	dfer	1.3	C $Header: /u/gcmpack/MITgcm/pkg/ggl90/ggl90_calc_visc.F,v 1.2 2007/07/31 23:01:58 ce107 Exp $
2	ce107	1.2	C $Name: $
3	mlosch	1.1	#include "GGL90_OPTIONS.h"
4
5			SUBROUTINE GGL90_CALC_VISC(
6			I bi,bj,iMin,iMax,jMin,jMax,K,
7			U KappaRU,KappaRV,
8			I myThid)
9
10			CBOP
11			C /==========================================================\
12			C \| SUBROUTINE GGL90_CALC_VISC \|
13			C \| o Add contrubution to net viscosity from GGL90 mixing \|
14			C \==========================================================/
15			IMPLICIT NONE
16
17			C == GLobal variables ==
18			#include "SIZE.h"
19			#include "EEPARAMS.h"
20			#include "PARAMS.h"
21			#include "DYNVARS.h"
22			#include "GRID.h"
23
24			C == Routine arguments ==
25			C bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
26			C myThid - Instance number for this innvocation of GGL90_CALC_VISC
27			C
28			INTEGER bi,bj,iMin,iMax,jMin,jMax,K
29			_RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
30			_RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
31			INTEGER myThid
32
33			#ifdef ALLOW_GGL90
34			#include "GGL90.h"
35
36			C == Local variables ==
37			C I, J, K - Loop counters
38			INTEGER i,j
39			CEOP
40			DO j=jMin,jMax
41			DO i=iMin,iMax
42	dfer	1.3	KappaRU(i,j,k) = KappaRU(i,j,k) + _maskW(i,j,k,bi,bj) *
43			& (.5 _d 0*(GGL90viscAr(i,j,k,bi,bj)
44			& +GGL90viscAr(i-1,j,k,bi,bj))
45			& -viscAr)
46	mlosch	1.1	ENDDO
47			ENDDO
48
49			DO j=jMin,jMax
50			DO i=iMin,iMax
51	dfer	1.3	KappaRV(i,j,k) = KappaRV(i,j,k) + _maskS(i,j,k,bi,bj) *
52			& (.5 _d 0*(GGL90viscAr(i,j,k,bi,bj)
53			& +GGL90viscAr(i,j-1,k,bi,bj))
54			& -viscAr)
55	mlosch	1.1	ENDDO
56			ENDDO
57
58			#endif /* ALLOW_GGL90 */
59
60			RETURN
61			END