/[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.7 - (show annotations) (download)
Thu Sep 16 09:35:11 2004 UTC (19 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint57m_post, checkpoint57g_pre, checkpoint57s_post, checkpoint58b_post, checkpoint57b_post, checkpoint57g_post, checkpoint56b_post, checkpoint57y_post, checkpoint57r_post, checkpoint57d_post, checkpoint57i_post, checkpoint59, checkpoint58, checkpoint55, checkpoint57, checkpoint56, checkpoint58f_post, checkpoint57n_post, checkpoint58d_post, checkpoint58a_post, checkpoint57z_post, checkpoint54f_post, checkpoint58y_post, checkpoint58t_post, checkpoint55i_post, checkpoint58m_post, checkpoint57l_post, checkpoint57t_post, checkpoint55c_post, checkpoint57v_post, checkpoint57f_post, checkpoint60, checkpoint61, checkpoint57a_post, checkpoint57h_pre, checkpoint58w_post, checkpoint57h_post, checkpoint57y_pre, checkpoint55g_post, checkpoint58o_post, checkpoint57c_post, checkpoint58p_post, checkpoint58q_post, checkpoint55d_post, checkpoint58e_post, mitgcm_mapl_00, checkpoint55d_pre, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, checkpoint55h_post, checkpoint58n_post, checkpoint57e_post, checkpoint55b_post, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint55f_post, checkpoint59c, checkpoint59b, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59h, checkpoint59k, checkpoint57p_post, checkpint57u_post, checkpoint57q_post, eckpoint57e_pre, checkpoint58k_post, checkpoint58v_post, checkpoint56a_post, checkpoint58l_post, checkpoint57h_done, checkpoint57j_post, checkpoint57f_pre, checkpoint58g_post, checkpoint58x_post, checkpoint59j, checkpoint58h_post, checkpoint56c_post, checkpoint58j_post, checkpoint57a_pre, checkpoint55a_post, checkpoint57o_post, checkpoint57k_post, checkpoint57w_post, checkpoint58i_post, checkpoint57x_post, checkpoint58c_post, checkpoint58u_post, checkpoint58s_post, checkpoint55e_post, checkpoint61d, checkpoint61b, checkpoint61c, checkpoint61a
Changes since 1.6: +8 -1 lines 
o prepare addition of new packages GGL90 and OPPS
  - new parameters useGGL90 and useOPPS
  - include call to the new routines at the (hopefully) right places
1 C $Header: /u/gcmpack/MITgcm/model/src/calc_viscosity.F,v 1.6 2004/09/02 09:13:49 mlosch 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,K,
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 bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
34 C myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
35 C KappaRU :: Total vertical viscosity for zonal flow.
36 C KappaRV :: Total vertical viscosity for meridional flow.
37 INTEGER bi,bj,iMin,iMax,jMin,jMax,K
38 _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
39 _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
40 INTEGER myThid
41
42 C !LOCAL VARIABLES:
43 C == Local variables ==
44 C I, J, K - Loop counters
45 INTEGER i,j
46 CEOP
47
48 DO j = 1-Oly, sNy+Oly
49 DO i = 1-Olx, sNx+Olx
50 KappaRU(i,j,k) = viscAr
51 ENDDO
52 ENDDO
53
54 DO j = 1-Oly, sNy+Oly
55 DO i = 1-Olx, sNx+Olx
56 KappaRV(i,j,k) = viscAr
57 ENDDO
58 ENDDO
59
60 #ifdef ALLOW_KPP
61 IF (useKPP) CALL KPP_CALC_VISC(
62 I bi,bj,iMin+1,iMax,jMin+1,jMax,K,
63 O KappaRU,KappaRV,
64 I myThid)
65 #endif
66
67 #ifdef ALLOW_PP81
68 IF (usePP81) CALL PP81_CALC_VISC(
69 I bi,bj,iMin+1,iMax,jMin+1,jMax,K,
70 O KappaRU,KappaRV,
71 I myThid)
72 #endif
73
74 #ifdef ALLOW_MY82
75 IF (useMY82) CALL MY82_CALC_VISC(
76 I bi,bj,iMin+1,iMax,jMin+1,jMax,K,
77 O KappaRU,KappaRV,
78 I myThid)
79 #endif
80
81 #ifdef ALLOW_GGL90
82 IF (useGGL90) CALL GGL90_CALC_VISC(
83 I bi,bj,iMin+1,iMax,jMin+1,jMax,K,
84 O KappaRU,KappaRV,
85 I myThid)
86 #endif
87
88
89 RETURN
90 END
  ViewVC Help
Powered by ViewVC 1.1.22