/[MITgcm]/MITgcm/pkg/ggl90/ggl90_calc_visc.F
ViewVC logotype

Diff of /MITgcm/pkg/ggl90/ggl90_calc_visc.F

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

revision 1.5 by jmc, Thu Oct 8 20:07:18 2009 UTC revision 1.6 by gforget, Fri Aug 6 18:37:05 2010 UTC
# Line 46  C     I, J, K - Loop counters Line 46  C     I, J, K - Loop counters
46    
47        DO j=jMin,jMax        DO j=jMin,jMax
48         DO i=iMin,iMax         DO i=iMin,iMax
 #ifdef ALLOW_GGL90_SMOOTH  
49          KappaRU(i,j,k) = KappaRU(i,j,k) +          KappaRU(i,j,k) = KappaRU(i,j,k) +
50       & ((       &                   (GGL90viscArU(i,j,k,bi,bj) - viscArNr(k) )
      &   p4 *(GGL90viscAr(i  ,j  ,k,bi,bj) * mskCor(i  ,j  ,bi,bj)  
      &       +GGL90viscAr(i-1,j  ,k,bi,bj) * mskCor(i-1,j  ,bi,bj))  
      &  +p8 *(GGL90viscAr(i-1,j-1,k,bi,bj) * mskCor(i-1,j-1,bi,bj)  
      &       +GGL90viscAr(i-1,j+1,k,bi,bj) * mskCor(i-1,j+1,bi,bj)  
      &       +GGL90viscAr(i  ,j-1,k,bi,bj) * mskCor(i  ,j-1,bi,bj)  
      &       +GGL90viscAr(i  ,j+1,k,bi,bj) * mskCor(i  ,j+1,bi,bj))  
      &  )  
      & /(p4 * 2. _d 0  
      &  +p8 *(      maskC(i-1,j-1,k,bi,bj) * mskCor(i-1,j-1,bi,bj)  
      &       +      maskC(i-1,j+1,k,bi,bj) * mskCor(i-1,j+1,bi,bj)  
      &       +      maskC(i  ,j-1,k,bi,bj) * mskCor(i  ,j-1,bi,bj)  
      &       +      maskC(i  ,j+1,k,bi,bj) * mskCor(i  ,j+1,bi,bj))  
      &  ) - viscArNr(k) )  
      &  *maskC(i  ,j,k,bi,bj)*mskCor(i  ,j,bi,bj)  
      &  *maskC(i-1,j,k,bi,bj)*mskCor(i-1,j,bi,bj)  
 #else  
         KappaRU(i,j,k) = KappaRU(i,j,k) + _maskW(i,j,k,bi,bj) *  
      &                   (.5 _d 0*(GGL90viscAr(i,j,k,bi,bj)  
      &                            +GGL90viscAr(i-1,j,k,bi,bj))  
      &                    - viscArNr(k) )  
 #endif  
51         ENDDO         ENDDO
52        ENDDO        ENDDO
53    
54        DO j=jMin,jMax        DO j=jMin,jMax
55         DO i=iMin,iMax         DO i=iMin,iMax
 #ifdef ALLOW_GGL90_SMOOTH  
         KappaRV(i,j,k) = KappaRV(i,j,k) +  
      & ((  
      &   p4 *(GGL90viscAr(i  ,j  ,k,bi,bj) * mskCor(i  ,j  ,bi,bj)  
      &       +GGL90viscAr(i  ,j-1,k,bi,bj) * mskCor(i  ,j-1,bi,bj))  
      &  +p8 *(GGL90viscAr(i-1,j  ,k,bi,bj) * mskCor(i-1,j  ,bi,bj)  
      &       +GGL90viscAr(i-1,j-1,k,bi,bj) * mskCor(i-1,j-1,bi,bj)  
      &       +GGL90viscAr(i+1,j  ,k,bi,bj) * mskCor(i+1,j  ,bi,bj)  
      &       +GGL90viscAr(i+1,j-1,k,bi,bj) * mskCor(i+1,j-1,bi,bj))  
      &  )  
      & /(p4 * 2. _d 0  
      &  +p8 *(      maskC(i-1,j  ,k,bi,bj) * mskCor(i-1,j  ,bi,bj)  
      &       +      maskC(i-1,j-1,k,bi,bj) * mskCor(i-1,j-1,bi,bj)  
      &       +      maskC(i+1,j  ,k,bi,bj) * mskCor(i+1,j  ,bi,bj)  
      &       +      maskC(i+1,j-1,k,bi,bj) * mskCor(i+1,j-1,bi,bj))  
      &  ) - viscArNr(k) )  
      &   *maskC(i,j  ,k,bi,bj)*mskCor(i,j  ,bi,bj)  
      &   *maskC(i,j-1,k,bi,bj)*mskCor(i,j-1,bi,bj)  
 #else  
56          KappaRV(i,j,k) = KappaRV(i,j,k) + _maskS(i,j,k,bi,bj) *          KappaRV(i,j,k) = KappaRV(i,j,k) + _maskS(i,j,k,bi,bj) *
57       &                   (.5 _d 0*(GGL90viscAr(i,j,k,bi,bj)       &                   (GGL90viscArV(i,j,k,bi,bj) - viscArNr(k) )
      &                            +GGL90viscAr(i,j-1,k,bi,bj))  
      &                    - viscArNr(k) )  
 #endif  
58         ENDDO         ENDDO
59        ENDDO        ENDDO
60    
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22