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

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

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

Revision 1.7 - (hide 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 mlosch 1.7 C $Header: /u/gcmpack/MITgcm/model/src/calc_viscosity.F,v 1.6 2004/09/02 09:13:49 mlosch Exp $
2 heimbach 1.4 C $Name: $
3 heimbach 1.1
4 edhill 1.5 #include "PACKAGES_CONFIG.h"
5 heimbach 1.1 #include "CPP_OPTIONS.h"
6    
7 cnh 1.3 CBOP
8     C !ROUTINE: CALC_VISCOSITY
9     C !INTERFACE:
10 heimbach 1.1 SUBROUTINE CALC_VISCOSITY(
11     I bi,bj,iMin,iMax,jMin,jMax,K,
12     O KappaRU,KappaRV,
13     I myThid)
14    
15 cnh 1.3 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 heimbach 1.1 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 cnh 1.3 C !INPUT/OUTPUT PARAMETERS:
32 heimbach 1.1 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 cnh 1.3 C KappaRU :: Total vertical viscosity for zonal flow.
36     C KappaRV :: Total vertical viscosity for meridional flow.
37 heimbach 1.1 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 cnh 1.3 C !LOCAL VARIABLES:
43 heimbach 1.1 C == Local variables ==
44     C I, J, K - Loop counters
45     INTEGER i,j
46 cnh 1.3 CEOP
47 heimbach 1.1
48 heimbach 1.4 DO j = 1-Oly, sNy+Oly
49     DO i = 1-Olx, sNx+Olx
50 heimbach 1.1 KappaRU(i,j,k) = viscAr
51     ENDDO
52     ENDDO
53    
54 heimbach 1.4 DO j = 1-Oly, sNy+Oly
55     DO i = 1-Olx, sNx+Olx
56 heimbach 1.1 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 mlosch 1.6 #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 mlosch 1.7 #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 heimbach 1.1
89     RETURN
90     END
  ViewVC Help
Powered by ViewVC 1.1.22