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

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

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

Revision 1.1 - (show annotations) (download)
Thu Sep 2 09:11:54 2004 UTC (19 years, 8 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, checkpoint54e_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 add two new packages
  - pp81 (Packanowski and Philander, 1981), Richardson number and
    stratification dependent mixing
  - my82 (Mellor and Yamada, 1982) level 2 turbulence closure scheme
1 C$Header: $
2 C$Name: $
3 #include "MY82_OPTIONS.h"
4
5 SUBROUTINE MY82_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 MY82_CALC_VISC |
13 C | o Add contrubution to net viscosity from MY82 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 CALC_COMMON_FACTORS
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_MY82
34 #include "MY82.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*(MYviscAr(i,j,k,bi,bj)+MYviscAr(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*(MYviscAr(i,j,k,bi,bj)+MYviscAr(i,j-1,k,bi,bj)) )
53 ENDDO
54 ENDDO
55
56 #endif /* ALLOW_MY82 */
57
58 RETURN
59 END
  ViewVC Help
Powered by ViewVC 1.1.22