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

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

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

revision 1.2 by ce107, Tue Jul 31 23:01:58 2007 UTC revision 1.3 by jmc, Tue Jan 20 00:26:04 2009 UTC
# Line 2  C $Header$ Line 2  C $Header$
2  C $Name$  C $Name$
3  #include "MY82_OPTIONS.h"  #include "MY82_OPTIONS.h"
4    
5        SUBROUTINE MY82_CALC_VISC(        SUBROUTINE MY82_CALC_VISC(
6       I        bi,bj,iMin,iMax,jMin,jMax,K,       I        bi,bj,iMin,iMax,jMin,jMax,K,
7       U        KappaRU,KappaRV,       U        KappaRU,KappaRV,
8       I        myThid)       I        myThid)
# Line 18  C     == GLobal variables == Line 18  C     == GLobal variables ==
18  #include "SIZE.h"  #include "SIZE.h"
19  #include "EEPARAMS.h"  #include "EEPARAMS.h"
20  #include "PARAMS.h"  #include "PARAMS.h"
 #include "DYNVARS.h"  
21  #include "GRID.h"  #include "GRID.h"
22    #include "MY82.h"
23    
24  C     == Routine arguments ==  C     == Routine arguments ==
25  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
# Line 31  C Line 31  C
31        INTEGER myThid        INTEGER myThid
32    
33  #ifdef ALLOW_MY82  #ifdef ALLOW_MY82
 #include "MY82.h"  
34    
35  C     == Local variables ==  C     == Local variables ==
36  C     I, J, K - Loop counters  C     I, J, K - Loop counters
# Line 41  CEOP Line 40  CEOP
40         DO i=iMin,iMax         DO i=iMin,iMax
41          KappaRU(i,j,k) = max ( KappaRU(i,j,k),          KappaRU(i,j,k) = max ( KappaRU(i,j,k),
42       &         KappaRU(i,j,k) - viscAr + _maskW(i,j,k,bi,bj) *       &         KappaRU(i,j,k) - viscAr + _maskW(i,j,k,bi,bj) *
43       &         0.5*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i-1,j,k,bi,bj)) )       &    0.5 _d 0*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i-1,j,k,bi,bj)) )
44         ENDDO         ENDDO
45        ENDDO        ENDDO
46    
# Line 49  CEOP Line 48  CEOP
48         DO i=iMin,iMax         DO i=iMin,iMax
49          KappaRV(i,j,k) = max ( KappaRV(i,j,k),          KappaRV(i,j,k) = max ( KappaRV(i,j,k),
50       &         KappaRV(i,j,k) - viscAr + _maskS(i,j,k,bi,bj) *       &         KappaRV(i,j,k) - viscAr + _maskS(i,j,k,bi,bj) *
51       &         0.5*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i,j-1,k,bi,bj)) )       &    0.5 _d 0*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i,j-1,k,bi,bj)) )
52         ENDDO         ENDDO
53        ENDDO        ENDDO
54    

Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3

  ViewVC Help
Powered by ViewVC 1.1.22