/[MITgcm]/MITgcm/pkg/mom_common/mom_hdissip.F
ViewVC logotype

Contents of /MITgcm/pkg/mom_common/mom_hdissip.F

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

Revision 1.4 - (show annotations) (download)
Fri Sep 16 19:32:20 2005 UTC (18 years, 8 months ago) by baylor
Branch: MAIN
CVS Tags: checkpoint57s_post
Changes since 1.3: +62 -101 lines 
Move calculation of viscosities to separate file mom_common/mom_calc_visc.F.  This allows Leith, LeithD,
and Smagorinsky to be used simultaneously, and soon will allow them to be used in strain-tension,
fluxform, and nonhydrostatic modes.  Also, introduce viscC4Smag, a biharmonic Smagorinsky
viscosity as in Griffies and Hallberg.
1 C $Header: /u/gcmpack/MITgcm/pkg/mom_common/mom_hdissip.F,v 1.3 2005/04/14 14:22:51 baylor Exp $
2 C $Name: $
3
4 #include "MOM_COMMON_OPTIONS.h"
5
6 SUBROUTINE MOM_HDISSIP(
7 I bi,bj,k,
8 I hDiv,vort3,tension,strain,KE,
9 I hFacZ,
10 O uDissip,vDissip,
11 I myThid)
12 IMPLICIT NONE
13 C
14 C Calculate horizontal dissipation terms in terms of tension and strain
15 C
16 C Du = d/dx At Tension + d/dy As Strain
17 C Dv = d/dx As Strain - d/dy At Tension
18
19 C == Global variables ==
20 #include "SIZE.h"
21 #include "EEPARAMS.h"
22 #include "GRID.h"
23 #include "PARAMS.h"
24
25 C == Routine arguments ==
26 INTEGER bi,bj,k
27 _RL hDiv(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
28 _RL vort3(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
29 _RL tension(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
30 _RL strain(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
31 _RL KE(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
32 _RS hFacZ(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
33 _RL uDissip(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
34 _RL vDissip(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
35 INTEGER myThid
36
37 C == Local variables ==
38 INTEGER I,J
39 _RL viscAh_t(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
40 _RL viscAh_s(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
41 _RL viscA4_t(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
42 _RL viscA4_s(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
43 LOGICAL harmonic, biharmonic, useVariableViscosity
44
45 CALL MOM_CALC_VISC(
46 I bi,bj,k,
47 O viscAh_s,viscAh_t,viscA4_s,viscA4_t,
48 O harmonic,biharmonic,useVariableViscosity,
49 I hDiv,vort3,tension,strain,KE,
50 I myThid)
51
52 C - Laplacian and bi-harmonic terms
53 IF (harmonic) THEN
54 DO j=2-Oly,sNy+Oly-1
55 DO i=2-Olx,sNx+Olx-1
56
57 uDissip(i,j) =
58 & recip_dyg(i,j,bi,bj)*recip_dyg(i,j,bi,bj)
59 & *recip_dxc(i,j,bi,bj)
60 & *(
61 & dyf( i , j ,bi,bj)*dyf( i , j ,bi,bj)
62 & *viscAh_t( i , j )*tension( i , j )
63 & -dyf(i-1, j ,bi,bj)*dyf(i-1, j ,bi,bj)
64 & *viscAh_t(i-1, j )*tension(i-1, j )
65 & )
66 & +recip_dxc(i,j,bi,bj)*recip_dxc(i,j,bi,bj)
67 & *recip_dyg(i,j,bi,bj)
68 & *(
69 & dxv( i ,j+1,bi,bj)*dxv( i ,j+1,bi,bj)
70 & *viscAh_s(i,j+1)*strain( i ,j+1)
71 & -dxv( i , j ,bi,bj)*dxv( i , j ,bi,bj)
72 & *viscAh_s(i, j )*strain( i , j )
73 & )
74
75 vDissip(i,j) =
76 & recip_dyc(i,j,bi,bj)*recip_dyc(i,j,bi,bj)
77 & *recip_dxg(i,j,bi,bj)
78 & *(
79 & dyu(i+1, j ,bi,bj)*dyu(i+1, j ,bi,bj)
80 & *viscAh_s(i+1,j)*strain(i+1,j)
81 & -dyu( i , j ,bi,bj)*dyu( i , j ,bi,bj)
82 & *viscAh_s( i ,j)*strain( i ,j)
83 & )
84 & -recip_dxg(i,j,bi,bj)*recip_dxg(i,j,bi,bj)
85 & *recip_dyc(i,j,bi,bj)
86 & *(
87 & dxf( i , j ,bi,bj)*dxf( i , j ,bi,bj)
88 & *viscAh_t(i, j )*tension(i, j )
89 & -dxf( i ,j-1,bi,bj)*dxf( i ,j-1,bi,bj)
90 & *viscAh_t(i,j-1)*tension(i,j-1)
91 & )
92
93 ENDDO
94 ENDDO
95 ENDIF
96 IF (biharmonic) THEN
97 STOP 'MOM_HDISSIP: BIHARMONIC NOT ALLOWED WITH STRAIN-TENSION'
98 ENDIF
99
100 RETURN
101 END
  ViewVC Help
Powered by ViewVC 1.1.22