/[MITgcm]/MITgcm/pkg/kpp/kpp_calc_visc.F
ViewVC logotype

Annotation of /MITgcm/pkg/kpp/kpp_calc_visc.F

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

Revision 1.4 - (hide annotations) (download)
Thu Oct 8 20:06:22 2009 UTC (14 years, 8 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint62v, checkpoint62u, checkpoint62t, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint62c, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62w, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint63g, checkpoint64, checkpoint65, checkpoint62, checkpoint63, checkpoint63p, checkpoint63q, checkpoint63r, checkpoint63s, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint62b, checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint64i, checkpoint64h, checkpoint64k, checkpoint64j, checkpoint64m, checkpoint64l, checkpoint64o, checkpoint64n, checkpoint64a, checkpoint64c, checkpoint64b, checkpoint64e, checkpoint64d, checkpoint64g, checkpoint64f, checkpoint61w, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.3: +6 -6 lines 
modif for vertical profile of (background) viscosity:
 for now, just uses 1rst value: viscArNr(1)
1 jmc 1.4 C $Header: /u/gcmpack/MITgcm/pkg/kpp/kpp_calc_visc.F,v 1.3 2007/05/03 21:31:35 jmc Exp $
2 jmc 1.3 C $Name: $
3 heimbach 1.1
4     #include "KPP_OPTIONS.h"
5    
6 jmc 1.4 SUBROUTINE KPP_CALC_VISC(
7 heimbach 1.1 I bi,bj,iMin,iMax,jMin,jMax,K,
8     U KappaRU,KappaRV,
9     I myThid)
10    
11 jmc 1.4 C *==========================================================*
12 heimbach 1.1 C | SUBROUTINE KPP_CALC_VISC |
13     C | o Add contrubution to net viscosity from KPP mixing |
14 jmc 1.4 C *==========================================================*
15 heimbach 1.1 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 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_KPP
34     #include "KPP.h"
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) = max ( KappaRU(i,j,k),
43 jmc 1.4 & KappaRU(i,j,k) - viscArNr(1) + _maskW(i,j,k,bi,bj) *
44 heimbach 1.1 & 0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i-1,j,k,bi,bj)) )
45     ENDDO
46     ENDDO
47    
48     DO j=jMin,jMax
49     DO i=iMin,iMax
50     KappaRV(i,j,k) = max ( KappaRV(i,j,k),
51 jmc 1.4 & KappaRV(i,j,k) - viscArNr(1) + _maskS(i,j,k,bi,bj) *
52 heimbach 1.1 & 0.5*(KPPviscAz(i,j,k,bi,bj)+KPPviscAz(i,j-1,k,bi,bj)) )
53     ENDDO
54     ENDDO
55    
56     #endif /* ALLOW_KPP */
57    
58     RETURN
59     END
  ViewVC Help
Powered by ViewVC 1.1.22