/[MITgcm]/MITgcm/pkg/my82/my82_calc.F
ViewVC logotype

Diff of /MITgcm/pkg/my82/my82_calc.F

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

revision 1.5 by jmc, Tue Jan 20 00:26:04 2009 UTC revision 1.6 by jmc, Thu Oct 8 20:07:53 2009 UTC
# Line 11  C !INTERFACE: ========================== Line 11  C !INTERFACE: ==========================
11       I     bi, bj, myTime, myThid )       I     bi, bj, myTime, myThid )
12    
13  C !DESCRIPTION: \bv  C !DESCRIPTION: \bv
14  C     /==========================================================\  C     *==========================================================*
15  C     | SUBROUTINE MY82_CALC                                     |  C     | SUBROUTINE MY82_CALC                                     |
16  C     | o Compute all MY82 fields defined in MY82.h              |  C     | o Compute all MY82 fields defined in MY82.h              |
17  C     |==========================================================|  C     *==========================================================*
18  C     | This subroutine is based on SPEM code                    |  C     | This subroutine is based on SPEM code                    |
19  C     \==========================================================/  C     *==========================================================*
20        IMPLICIT NONE        IMPLICIT NONE
21  C  C
22  C--------------------------------------------------------------------  C--------------------------------------------------------------------
# Line 160  C     M. Satoh: Eq. (11.3.43) Line 160  C     M. Satoh: Eq. (11.3.43)
160           MYviscAr(I,J,K,bi,bj) = MYhbl(I,J,bi,bj)*tkel*SM(I,J,K)           MYviscAr(I,J,K,bi,bj) = MYhbl(I,J,bi,bj)*tkel*SM(I,J,K)
161           MYdiffKr(I,J,K,bi,bj) = MYhbl(I,J,bi,bj)*tkel*SH(I,J,K)           MYdiffKr(I,J,K,bi,bj) = MYhbl(I,J,bi,bj)*tkel*SH(I,J,K)
162  C     Set a minium (= background) value  C     Set a minium (= background) value
163           MYviscAr(I,J,K,bi,bj) = MAX(MYviscAr(I,J,K,bi,bj),viscAr)           MYviscAr(I,J,K,bi,bj) = MAX(MYviscAr(I,J,K,bi,bj),
164         &                               viscArnr(k) )
165           MYdiffKr(I,J,K,bi,bj) = MAX(MYdiffKr(I,J,K,bi,bj),           MYdiffKr(I,J,K,bi,bj) = MAX(MYdiffKr(I,J,K,bi,bj),
166       &                               diffKrNrT(k) )       &                               diffKrNrT(k) )
167  C     Set a maximum and mask land point  C     Set a maximum and mask land point
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22