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

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

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

revision 1.2 by heimbach, Mon Aug 13 23:27:53 2001 UTC revision 1.4.4.1 by edhill, Thu Oct 2 18:10:45 2003 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4    #include "PACKAGES_CONFIG.h"
5  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
6    
7  CStartOfInterFace  CBOP
8    C     !ROUTINE: CALC_VISCOSITY
9    C     !INTERFACE:
10        SUBROUTINE CALC_VISCOSITY(        SUBROUTINE CALC_VISCOSITY(
11       I        bi,bj,iMin,iMax,jMin,jMax,K,       I        bi,bj,iMin,iMax,jMin,jMax,K,
12       O        KappaRU,KappaRV,       O        KappaRU,KappaRV,
13       I        myThid)       I        myThid)
14    
15  C     /==========================================================\  C     !DESCRIPTION: \bv
16  C     | SUBROUTINE CALC_VISCOSITY                                |  C     *==========================================================*
17  C     | o Calculate net viscosity                                |  C     | SUBROUTINE CALC_VISCOSITY                                
18  C     \==========================================================/  C     | o Calculate net vertical viscosity                        
19        IMPLICIT NONE  C     *==========================================================*
20    C     \ev
21    
22    C     !USES:
23          IMPLICIT NONE
24  C     == GLobal variables ==  C     == GLobal variables ==
25  #include "SIZE.h"  #include "SIZE.h"
26  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 20  C     == GLobal variables == Line 28  C     == GLobal variables ==
28  #include "DYNVARS.h"  #include "DYNVARS.h"
29  #include "GRID.h"  #include "GRID.h"
30    
31    C     !INPUT/OUTPUT PARAMETERS:
32  C     == Routine arguments ==  C     == Routine arguments ==
33  C     bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation  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  C     myThid - Instance number for this innvocation of CALC_COMMON_FACTORS
35  C  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        INTEGER bi,bj,iMin,iMax,jMin,jMax,K
38        _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
39        _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)        _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
40        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
41    
42    C     !LOCAL VARIABLES:
43  C     == Local variables ==  C     == Local variables ==
44  C     I, J, K - Loop counters  C     I, J, K - Loop counters
45        INTEGER i,j        INTEGER i,j
46    CEOP
47    
48        DO j=jMin,jMax        DO j = 1-Oly, sNy+Oly
49         DO i=iMin,iMax         DO i = 1-Olx, sNx+Olx
50          KappaRU(i,j,k) = viscAr          KappaRU(i,j,k) = viscAr
51         ENDDO         ENDDO
52        ENDDO        ENDDO
53    
54        DO j=jMin,jMax        DO j = 1-Oly, sNy+Oly
55         DO i=iMin,iMax         DO i = 1-Olx, sNx+Olx
56          KappaRV(i,j,k) = viscAr          KappaRV(i,j,k) = viscAr
57         ENDDO         ENDDO
58        ENDDO        ENDDO
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.4.4.1
  ViewVC Help
Powered by ViewVC 1.1.22