model/src/calc_viscosity.F


#include "CPP_OPTIONS.h"

CStartOfInterFace
      SUBROUTINE CALC_VISCOSITY( 
     I        bi,bj,iMin,iMax,jMin,jMax,K,
     I        maskUp,
     O        KappaRU,KappaRV,
     I        myThid)

C     /==========================================================\
C     | SUBROUTINE CALC_VISCOSITY                                |
C     | o Calculate net viscosity                                |
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     maskUp  - land/water mask for Wvel points (above tracer level)
C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
C
      INTEGER bi,bj,iMin,iMax,jMin,jMax,K
      _RS maskUp(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
      _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
      _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
      INTEGER myThid
CEndOfInterface

C     == Local variables ==
C     I, J, K - Loop counters
      INTEGER i,j

      DO j=jMin,jMax
       DO i=iMin,iMax
        KappaRU(i,j,k) = viscAr
       ENDDO
      ENDDO

      DO j=jMin,jMax
       DO i=iMin,iMax
        KappaRV(i,j,k) = viscAr
       ENDDO
      ENDDO

#ifdef ALLOW_KPP
      IF (useKPP) CALL KPP_CALC_VISC(
     I        bi,bj,iMin+1,iMax,jMin+1,jMax,K,
     I        maskUp,
     O        KappaRU,KappaRV,
     I        myThid)
#endif


      RETURN
      END
1	heimbach	1.1
2			#include "CPP_OPTIONS.h"
3
4			CStartOfInterFace
5			SUBROUTINE CALC_VISCOSITY(
6			I bi,bj,iMin,iMax,jMin,jMax,K,
7			I maskUp,
8			O KappaRU,KappaRV,
9			I myThid)
10
11			C /==========================================================\
12			C \| SUBROUTINE CALC_VISCOSITY \|
13			C \| o Calculate net viscosity \|
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 maskUp - land/water mask for Wvel points (above tracer level)
27			C myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
28			C
29			INTEGER bi,bj,iMin,iMax,jMin,jMax,K
30			_RS maskUp(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
31			_RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
32			_RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
33			INTEGER myThid
34			CEndOfInterface
35
36			C == Local variables ==
37			C I, J, K - Loop counters
38			INTEGER i,j
39
40			DO j=jMin,jMax
41			DO i=iMin,iMax
42			KappaRU(i,j,k) = viscAr
43			ENDDO
44			ENDDO
45
46			DO j=jMin,jMax
47			DO i=iMin,iMax
48			KappaRV(i,j,k) = viscAr
49			ENDDO
50			ENDDO
51
52			#ifdef ALLOW_KPP
53			IF (useKPP) CALL KPP_CALC_VISC(
54			I bi,bj,iMin+1,iMax,jMin+1,jMax,K,
55			I maskUp,
56			O KappaRU,KappaRV,
57			I myThid)
58			#endif
59
60
61			RETURN
62			END