pkg/ggl90/ggl90_calc_visc.F

89474f9a5c Mart*0001 #include "GGL90_OPTIONS.h"
                0002
310f4f2d26 Jean*0003 CBOP
                0004 C     !ROUTINE: GGL90_CALC_VISC
                0005
                0006 C     !INTERFACE:
5e48dccc42 Jean*0007       SUBROUTINE GGL90_CALC_VISC(
310f4f2d26 Jean*0008      I        bi, bj, iMin, iMax, jMin, jMax, k,
                0009      U        KappaRU, KappaRV,
                0010      I        myThid )
89474f9a5c Mart*0011
310f4f2d26 Jean*0012 C     !DESCRIPTION:
5e48dccc42 Jean*0013 C     *==========================================================*
89474f9a5c Mart*0014 C     | SUBROUTINE GGL90_CALC_VISC                               |
310f4f2d26 Jean*0015 C     | o Add contribution to net viscosity from GGL90 mixing    |
5e48dccc42 Jean*0016 C     *==========================================================*
89474f9a5c Mart*0017
310f4f2d26 Jean*0018 C     !USES:
                0019       IMPLICIT NONE
89474f9a5c Mart*0020 C     == GLobal variables ==
                0021 #include "SIZE.h"
                0022 #include "EEPARAMS.h"
                0023 #include "PARAMS.h"
310f4f2d26 Jean*0024 c#include "DYNVARS.h"
                0025 c#include "GRID.h"
5e48dccc42 Jean*0026 #include "GGL90.h"
89474f9a5c Mart*0027
310f4f2d26 Jean*0028 C     !INPUT/OUTPUT PARAMETERS:
                0029 C     bi, bj    :: current tile indices
                0030 C     iMin,iMax :: Range of points for which calculation is done
                0031 C     jMin,jMax :: Range of points for which calculation is done
                0032 C     k         :: current level index
                0033 C     KappaRU   :: vertical viscosity array for U-component
                0034 C     KappaRV   :: vertical viscosity array for V-component
                0035 C     myThid    :: my Thread Id number
                0036       INTEGER bi, bj, iMin, iMax, jMin, jMax, k
                0037       _RL KappaRU(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
                0038       _RL KappaRV(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
89474f9a5c Mart*0039       INTEGER myThid
                0040
                0041 #ifdef ALLOW_GGL90
310f4f2d26 Jean*0042 C     !LOCAL VARIABLES:
                0043 C     i, j      :: Loop counters
                0044       INTEGER i, j
                0045 CEOP
004d5ee949 Davi*0046
89474f9a5c Mart*0047       DO j=jMin,jMax
                0048        DO i=iMin,iMax
310f4f2d26 Jean*0049         KappaRU(i,j,k) = KappaRU(i,j,k)
                0050      &                 + ( GGL90viscArU(i,j,k,bi,bj) - viscArNr(k) )
89474f9a5c Mart*0051        ENDDO
                0052       ENDDO
                0053
                0054       DO j=jMin,jMax
                0055        DO i=iMin,iMax
310f4f2d26 Jean*0056         KappaRV(i,j,k) = KappaRV(i,j,k)
                0057      &                 + ( GGL90viscArV(i,j,k,bi,bj) - viscArNr(k) )
89474f9a5c Mart*0058        ENDDO
                0059       ENDDO
                0060
                0061 #endif /* ALLOW_GGL90 */
                0062
                0063       RETURN
                0064       END