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

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

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

revision 1.4 by adcroft, Wed Jun 10 16:05:39 1998 UTC revision 1.5 by cnh, Tue Aug 18 16:32:41 1998 UTC
# Line 6  CStartOfInterFace Line 6  CStartOfInterFace
6        SUBROUTINE CALC_DIFFUSIVITY(        SUBROUTINE CALC_DIFFUSIVITY(
7       I        bi,bj,iMin,iMax,jMin,jMax,K,       I        bi,bj,iMin,iMax,jMin,jMax,K,
8       I        maskC,maskUp,KapGM,K33,       I        maskC,maskUp,KapGM,K33,
9       O        KappaZT,KappaZS,       O        KappaRT,KappaRS,
10       I        myThid)       I        myThid)
11    
12  C     /==========================================================\  C     /==========================================================\
13  C     | SUBROUTINE CALC_DIFFUSIVITY                              |  C     | SUBROUTINE CALC_DIFFUSIVITY                              |
14  C     | o Calculate net diffusivity                              |  C     | o Calculate net diffusivity                              |
 C     |==========================================================|  
15  C     \==========================================================/  C     \==========================================================/
16        IMPLICIT NONE        IMPLICIT NONE
17    
# Line 33  C Line 32  C
32        _RS maskUp(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RS maskUp(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
33        _RL KapGM(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL KapGM(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
34        _RL K33(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nz)        _RL K33(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nz)
35        _RL KappaZT(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nz)        _RL KappaRT(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nz)
36        _RL KappaZS(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nz)        _RL KappaRS(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nz)
37        INTEGER myThid        INTEGER myThid
38  CEndOfInterface  CEndOfInterface
39    
# Line 44  C     I, J, K - Loop counters Line 43  C     I, J, K - Loop counters
43    
44        DO j=jMin,jMax        DO j=jMin,jMax
45         DO i=iMin,iMax         DO i=iMin,iMax
46          KappaZT(i,j,k) = maskC(i,j)*maskUp(i,j)*          KappaRT(i,j,k) = maskC(i,j)*maskUp(i,j)*
47       &      ( diffKzT + KapGM(i,j)*K33(i,j,k) )       &      ( diffKrT + KapGM(i,j)*K33(i,j,k) )
48         ENDDO         ENDDO
49        ENDDO        ENDDO
50    
51        DO j=jMin,jMax        DO j=jMin,jMax
52         DO i=iMin,iMax         DO i=iMin,iMax
53          KappaZS(i,j,k) = maskC(i,j)*maskUp(i,j)*          KappaRS(i,j,k) = maskC(i,j)*maskUp(i,j)*
54       &      ( diffKzS + KapGM(i,j)*K33(i,j,k) )       &      ( diffKrS + KapGM(i,j)*K33(i,j,k) )
55         ENDDO         ENDDO
56        ENDDO        ENDDO
57    
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22