/[MITgcm]/MITgcm/model/src/calc_viscosity.F
ViewVC logotype

Contents of /MITgcm/model/src/calc_viscosity.F

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

Revision 1.9 - (show annotations) (download)
Thu Oct 8 20:03:18 2009 UTC (14 years, 7 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint62v, checkpoint62u, checkpoint62t, checkpoint65a, 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.8: +3 -3 lines 
allow to specify a vertical profile of vertical viscosity
1 C $Header: /u/gcmpack/MITgcm/model/src/calc_viscosity.F,v 1.8 2008/10/20 23:51:39 jmc Exp $
2 C $Name: $
3
4 #include "PACKAGES_CONFIG.h"
5 #include "CPP_OPTIONS.h"
6
7 CBOP
8 C !ROUTINE: CALC_VISCOSITY
9 C !INTERFACE:
10 SUBROUTINE CALC_VISCOSITY(
11 I bi,bj, iMin,iMax,jMin,jMax,
12 O KappaRU, KappaRV,
13 I myThid )
14
15 C !DESCRIPTION: \bv
16 C *==========================================================*
17 C | SUBROUTINE CALC_VISCOSITY
18 C | o Calculate net vertical viscosity
19 C *==========================================================*
20 C \ev
21
22 C !USES:
23 IMPLICIT NONE
24 C == GLobal variables ==
25 #include "SIZE.h"
26 #include "EEPARAMS.h"
27 #include "PARAMS.h"
28 #include "DYNVARS.h"
29 #include "GRID.h"
30
31 C !INPUT/OUTPUT PARAMETERS:
32 C == Routine arguments ==
33 C iMin,iMax,jMin,jMax :: Range of points for which calculation
34 C bi,bj :: current tile indices
35 C KappaRU :: Total vertical viscosity for zonal flow.
36 C KappaRV :: Total vertical viscosity for meridional flow.
37 C myThid :: my Thread Id number
38 INTEGER iMin,iMax,jMin,jMax
39 INTEGER bi,bj
40 _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
41 _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
42 INTEGER myThid
43
44 C !LOCAL VARIABLES:
45 C == Local variables ==
46 C i, j, k :: Loop counters
47 INTEGER i,j,k
48 CEOP
49
50 DO k = 1,Nr
51
52 DO j = 1-Oly, sNy+Oly
53 DO i = 1-Olx, sNx+Olx
54 KappaRU(i,j,k) = viscArNr(k)
55 KappaRV(i,j,k) = viscArNr(k)
56 ENDDO
57 ENDDO
58
59 #ifdef ALLOW_KPP
60 IF (useKPP) CALL KPP_CALC_VISC(
61 I bi,bj, iMin,iMax,jMin,jMax, k,
62 O KappaRU,KappaRV,
63 I myThid)
64 #endif
65
66 #ifdef ALLOW_PP81
67 IF (usePP81) CALL PP81_CALC_VISC(
68 I bi,bj, iMin,iMax,jMin,jMax, k,
69 O KappaRU,KappaRV,
70 I myThid)
71 #endif
72
73 #ifdef ALLOW_MY82
74 IF (useMY82) CALL MY82_CALC_VISC(
75 I bi,bj, iMin,iMax,jMin,jMax, k,
76 O KappaRU,KappaRV,
77 I myThid)
78 #endif
79
80 #ifdef ALLOW_GGL90
81 IF (useGGL90) CALL GGL90_CALC_VISC(
82 I bi,bj, iMin,iMax,jMin,jMax, k,
83 O KappaRU,KappaRV,
84 I myThid)
85 #endif
86
87 C-- end of k loop
88 ENDDO
89
90 RETURN
91 END
  ViewVC Help
Powered by ViewVC 1.1.22