/[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.4 - (hide annotations) (download)
Thu May 30 02:27:00 2002 UTC (22 years ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint47e_post, checkpoint46l_post, checkpoint46g_pre, checkpoint47c_post, checkpoint50c_post, checkpoint46f_post, checkpoint48e_post, checkpoint50c_pre, checkpoint46b_post, checkpoint48i_post, checkpoint46l_pre, checkpoint51, checkpoint50, checkpoint50d_post, checkpoint50b_pre, checkpoint51f_post, checkpoint48b_post, checkpoint51d_post, checkpoint48c_pre, checkpoint47d_pre, checkpoint47a_post, checkpoint48d_pre, checkpoint51j_post, checkpoint47i_post, checkpoint47d_post, checkpoint46d_pre, checkpoint48d_post, checkpoint48f_post, checkpoint45d_post, checkpoint46j_pre, checkpoint48h_post, checkpoint51b_pre, checkpoint46a_post, checkpoint47g_post, checkpoint46j_post, checkpoint51h_pre, checkpoint46k_post, checkpoint48a_post, checkpoint50f_post, checkpoint50a_post, checkpoint50f_pre, checkpoint47j_post, branch-exfmods-tag, branchpoint-genmake2, checkpoint46e_pre, checkpoint48c_post, checkpoint45b_post, checkpoint46b_pre, checkpoint51b_post, checkpoint51c_post, checkpoint46c_pre, checkpoint46, checkpoint47b_post, checkpoint46h_pre, checkpoint46m_post, checkpoint46a_pre, checkpoint50g_post, checkpoint45c_post, checkpoint46g_post, checkpoint50h_post, checkpoint50e_pre, checkpoint50i_post, checkpoint51i_pre, checkpoint47f_post, checkpoint50e_post, checkpoint46i_post, checkpoint46c_post, checkpoint50d_pre, checkpoint46e_post, checkpoint51e_post, checkpoint47, checkpoint48, checkpoint49, checkpoint46h_post, checkpoint51f_pre, checkpoint48g_post, checkpoint47h_post, checkpoint51g_post, checkpoint46d_post, checkpoint50b_post, checkpoint51a_post
Branch point for: branch-exfmods-curt, branch-genmake2
Changes since 1.3: +6 -6 lines 
Extended calculation to full overlap
(needed to avoid partial recomputation in adjoint mode).
CVS ----------------------------------------------------------------------
1 heimbach 1.4 C $Header: /u/gcmpack/MITgcm/model/src/calc_viscosity.F,v 1.3 2001/09/26 18:09:14 cnh Exp $
2     C $Name: $
3 heimbach 1.1
4     #include "CPP_OPTIONS.h"
5    
6 cnh 1.3 CBOP
7     C !ROUTINE: CALC_VISCOSITY
8     C !INTERFACE:
9 heimbach 1.1 SUBROUTINE CALC_VISCOSITY(
10     I bi,bj,iMin,iMax,jMin,jMax,K,
11     O KappaRU,KappaRV,
12     I myThid)
13    
14 cnh 1.3 C !DESCRIPTION: \bv
15     C *==========================================================*
16     C | SUBROUTINE CALC_VISCOSITY
17     C | o Calculate net vertical viscosity
18     C *==========================================================*
19     C \ev
20    
21     C !USES:
22 heimbach 1.1 IMPLICIT NONE
23     C == GLobal variables ==
24     #include "SIZE.h"
25     #include "EEPARAMS.h"
26     #include "PARAMS.h"
27     #include "DYNVARS.h"
28     #include "GRID.h"
29    
30 cnh 1.3 C !INPUT/OUTPUT PARAMETERS:
31 heimbach 1.1 C == Routine arguments ==
32     C bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
33     C myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
34 cnh 1.3 C KappaRU :: Total vertical viscosity for zonal flow.
35     C KappaRV :: Total vertical viscosity for meridional flow.
36 heimbach 1.1 INTEGER bi,bj,iMin,iMax,jMin,jMax,K
37     _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
38     _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
39     INTEGER myThid
40    
41 cnh 1.3 C !LOCAL VARIABLES:
42 heimbach 1.1 C == Local variables ==
43     C I, J, K - Loop counters
44     INTEGER i,j
45 cnh 1.3 CEOP
46 heimbach 1.1
47 heimbach 1.4 DO j = 1-Oly, sNy+Oly
48     DO i = 1-Olx, sNx+Olx
49 heimbach 1.1 KappaRU(i,j,k) = viscAr
50     ENDDO
51     ENDDO
52    
53 heimbach 1.4 DO j = 1-Oly, sNy+Oly
54     DO i = 1-Olx, sNx+Olx
55 heimbach 1.1 KappaRV(i,j,k) = viscAr
56     ENDDO
57     ENDDO
58    
59     #ifdef ALLOW_KPP
60     IF (useKPP) CALL KPP_CALC_VISC(
61     I bi,bj,iMin+1,iMax,jMin+1,jMax,K,
62     O KappaRU,KappaRV,
63     I myThid)
64     #endif
65    
66    
67     RETURN
68     END
  ViewVC Help
Powered by ViewVC 1.1.22