/[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.4 - (show 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 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
4 #include "CPP_OPTIONS.h"
5
6 CBOP
7 C !ROUTINE: CALC_VISCOSITY
8 C !INTERFACE:
9 SUBROUTINE CALC_VISCOSITY(
10 I bi,bj,iMin,iMax,jMin,jMax,K,
11 O KappaRU,KappaRV,
12 I myThid)
13
14 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 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 C !INPUT/OUTPUT PARAMETERS:
31 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 C KappaRU :: Total vertical viscosity for zonal flow.
35 C KappaRV :: Total vertical viscosity for meridional flow.
36 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 C !LOCAL VARIABLES:
42 C == Local variables ==
43 C I, J, K - Loop counters
44 INTEGER i,j
45 CEOP
46
47 DO j = 1-Oly, sNy+Oly
48 DO i = 1-Olx, sNx+Olx
49 KappaRU(i,j,k) = viscAr
50 ENDDO
51 ENDDO
52
53 DO j = 1-Oly, sNy+Oly
54 DO i = 1-Olx, sNx+Olx
55 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