/[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.11 - (show annotations) (download)
Wed Dec 24 19:09:33 2014 UTC (9 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint65z, checkpoint65x, checkpoint65y, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65i, checkpoint65n, checkpoint65l, checkpoint65m, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint65o
Changes since 1.10: +42 -19 lines 
- add one more level to vertical viscosity local arrays (Nr+1, previously Nr)
  since no-slip bottom BC uses viscosity @ k+1 to update velocity @ level k
- for now and until vertical mixing scheme are updated to fill up Nr+1 level,
  just copy Nr value to Nr+1
1 C $Header: /u/gcmpack/MITgcm/model/src/calc_viscosity.F,v 1.10 2014/07/30 03:31:35 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+1)
41 _RL kappaRV(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr+1)
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 INTEGER ki
49 CEOP
50
51 DO k = 1,Nr+1
52 ki = MIN(k,Nr)
53
54 DO j = 1-OLy, sNy+OLy
55 DO i = 1-OLx, sNx+OLx
56 kappaRU(i,j,k) = viscArNr(ki)
57 kappaRV(i,j,k) = viscArNr(ki)
58 ENDDO
59 ENDDO
60
61 #ifdef ALLOW_KPP
62 IF ( useKPP .AND. k.LE.Nr ) THEN
63 CALL KPP_CALC_VISC(
64 I bi,bj, iMin,iMax,jMin,jMax, k,
65 O kappaRU, kappaRV,
66 I myThid)
67 ENDIF
68 #endif
69
70 #ifdef ALLOW_PP81
71 IF ( usePP81 .AND. k.LE.Nr ) THEN
72 CALL PP81_CALC_VISC(
73 I bi,bj, iMin,iMax,jMin,jMax, k,
74 O kappaRU, kappaRV,
75 I myThid)
76 ENDIF
77 #endif
78
79 #ifdef ALLOW_KL10
80 IF ( useKL10 .AND. k.LE.Nr ) THEN
81 CALL KL10_CALC_VISC(
82 I bi,bj, iMin,iMax,jMin,jMax, k,
83 O kappaRU, kappaRV,
84 I myThid)
85 ENDIF
86 #endif
87
88 #ifdef ALLOW_MY82
89 IF ( useMY82 .AND. k.LE.Nr ) THEN
90 CALL MY82_CALC_VISC(
91 I bi,bj, iMin,iMax,jMin,jMax, k,
92 O kappaRU, kappaRV,
93 I myThid)
94 ENDIF
95 #endif
96
97 #ifdef ALLOW_GGL90
98 IF ( useGGL90 .AND. k.LE.Nr ) THEN
99 CALL GGL90_CALC_VISC(
100 I bi,bj, iMin,iMax,jMin,jMax, k,
101 O kappaRU, kappaRV,
102 I myThid)
103 ENDIF
104 #endif
105
106 IF ( k.EQ.Nr+1 .AND.
107 & ( usePP81 .OR. useKL10 .OR. useMY82 .OR. useGGL90 )
108 & ) THEN
109 DO j = 1-OLy, sNy+OLy
110 DO i = 1-OLx, sNx+OLx
111 kappaRU(i,j,k) = kappaRU(i,j,ki)
112 kappaRV(i,j,k) = kappaRV(i,j,ki)
113 ENDDO
114 ENDDO
115 ENDIF
116
117 C-- end of k loop
118 ENDDO
119
120 RETURN
121 END
  ViewVC Help
Powered by ViewVC 1.1.22