/[MITgcm]/MITgcm/pkg/my82/my82_calc_visc.F
ViewVC logotype

Contents of /MITgcm/pkg/my82/my82_calc_visc.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.3 - (show annotations) (download)
Tue Jan 20 00:26:04 2009 UTC (16 years, 1 month ago) by jmc
Branch: MAIN
CVS Tags: checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q
Changes since 1.2: +5 -6 lines
add missing "_d 0"

1 C $Header: /u/gcmpack/MITgcm/pkg/my82/my82_calc_visc.F,v 1.2 2007/07/31 23:01:58 ce107 Exp $
2 C $Name: $
3 #include "MY82_OPTIONS.h"
4
5 SUBROUTINE MY82_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 MY82_CALC_VISC |
13 C | o Add contrubution to net viscosity from MY82 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 "GRID.h"
22 #include "MY82.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 CALC_COMMON_FACTORS
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_MY82
34
35 C == Local variables ==
36 C I, J, K - Loop counters
37 INTEGER i,j
38 CEOP
39 DO j=jMin,jMax
40 DO i=iMin,iMax
41 KappaRU(i,j,k) = max ( KappaRU(i,j,k),
42 & KappaRU(i,j,k) - viscAr + _maskW(i,j,k,bi,bj) *
43 & 0.5 _d 0*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i-1,j,k,bi,bj)) )
44 ENDDO
45 ENDDO
46
47 DO j=jMin,jMax
48 DO i=iMin,iMax
49 KappaRV(i,j,k) = max ( KappaRV(i,j,k),
50 & KappaRV(i,j,k) - viscAr + _maskS(i,j,k,bi,bj) *
51 & 0.5 _d 0*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i,j-1,k,bi,bj)) )
52 ENDDO
53 ENDDO
54
55 #endif /* ALLOW_MY82 */
56
57 RETURN
58 END

  ViewVC Help
Powered by ViewVC 1.1.22