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

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

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

Revision 1.1 - (hide annotations) (download)
Thu Sep 16 11:27:18 2004 UTC (19 years, 7 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint57t_post, checkpoint57o_post, checkpoint58e_post, checkpoint57v_post, checkpoint58u_post, checkpoint58w_post, checkpoint57m_post, checkpoint55c_post, checkpoint57s_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, checkpoint57e_post, checkpoint55h_post, checkpoint58n_post, checkpoint58x_post, checkpoint57g_pre, checkpoint55b_post, checkpoint58t_post, checkpoint58h_post, checkpoint56c_post, checkpoint57y_pre, checkpoint55, checkpoint57f_pre, checkpoint57a_post, checkpoint58q_post, checkpoint54f_post, checkpoint55g_post, checkpoint58j_post, checkpoint59e, checkpoint59d, checkpoint59a, checkpoint55f_post, checkpoint59c, checkpoint59b, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, eckpoint57e_pre, checkpoint57h_done, checkpoint58f_post, checkpoint57x_post, checkpoint57n_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint57f_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, checkpoint58g_post, checkpoint58o_post, checkpoint57z_post, checkpoint57c_post, checkpoint58y_post, checkpoint55e_post, checkpoint58k_post, checkpoint58v_post, checkpoint55a_post, checkpoint58s_post, checkpoint58p_post, checkpoint57j_post, checkpoint58b_post, checkpoint57h_pre, checkpoint58m_post, checkpoint57l_post, checkpoint57h_post, checkpoint56a_post, checkpoint55d_post 
o initial check-in of TKE-model of Gaspar et al. 1990
1 mlosch 1.1 C$Header: $
2     C$Name: $
3     #include "GGL90_OPTIONS.h"
4    
5     SUBROUTINE GGL90_CALC_VISC(
6     I bi,bj,iMin,iMax,jMin,jMax,K,
7     U KappaRU,KappaRV,
8     I myThid)
9    
10     CBOP
11     C /==========================================================\
12     C | SUBROUTINE GGL90_CALC_VISC |
13     C | o Add contrubution to net viscosity from GGL90 mixing |
14     C \==========================================================/
15     IMPLICIT NONE
16    
17     C == GLobal variables ==
18     #include "SIZE.h"
19     #include "EEPARAMS.h"
20     #include "PARAMS.h"
21     #include "DYNVARS.h"
22     #include "GRID.h"
23    
24     C == Routine arguments ==
25     C bi, bj, iMin, iMax, jMin, jMax - Range of points for which calculation
26     C myThid - Instance number for this innvocation of GGL90_CALC_VISC
27     C
28     INTEGER bi,bj,iMin,iMax,jMin,jMax,K
29     _RL KappaRU(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
30     _RL KappaRV(1-Olx:sNx+Olx,1-Oly:sNy+Oly,Nr)
31     INTEGER myThid
32    
33     #ifdef ALLOW_GGL90
34     #include "GGL90.h"
35    
36     C == Local variables ==
37     C I, J, K - Loop counters
38     INTEGER i,j
39     CEOP
40     DO j=jMin,jMax
41     DO i=iMin,iMax
42     KappaRU(i,j,k) = max ( KappaRU(i,j,k),
43     & KappaRU(i,j,k) - viscAr + _maskW(i,j,k,bi,bj) *
44     & 0.5*(GGL90viscAr(i,j,k,bi,bj)+GGL90viscAr(i-1,j,k,bi,bj)) )
45     ENDDO
46     ENDDO
47    
48     DO j=jMin,jMax
49     DO i=iMin,iMax
50     KappaRV(i,j,k) = max ( KappaRV(i,j,k),
51     & KappaRV(i,j,k) - viscAr + _maskS(i,j,k,bi,bj) *
52     & 0.5*(GGL90viscAr(i,j,k,bi,bj)+GGL90viscAr(i,j-1,k,bi,bj)) )
53     ENDDO
54     ENDDO
55    
56     #endif /* ALLOW_GGL90 */
57    
58     RETURN
59     END
  ViewVC Help
Powered by ViewVC 1.1.22