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

Diff of /MITgcm/pkg/mom_common/mom_calc_visc.F

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

revision 1.5 by baylor, Tue Sep 20 20:41:43 2005 UTC revision 1.6 by baylor, Tue Sep 20 21:01:30 2005 UTC
# Line 91  C     == Routine arguments == Line 91  C     == Routine arguments ==
91  C     == Local variables ==  C     == Local variables ==
92        INTEGER I,J        INTEGER I,J
93        _RL smag2fac, smag4fac        _RL smag2fac, smag4fac
94        _RL viscAhRemax, viscA4Remax        _RL viscAhRe_max, viscA4Re_max
95        _RL Alin,Alinmin,grdVrt,grdDiv        _RL Alin,Alinmin,grdVrt,grdDiv
96        _RL recip_dt,L2,L3,L4,L5,L2rdt,L4rdt        _RL recip_dt,L2,L3,L4,L5,L2rdt,L4rdt
97        _RL Uscl,U4scl        _RL Uscl,U4scl
# Line 136  C     == Local variables == Line 136  C     == Local variables ==
136       &  .OR.(viscC2leithD.NE.0.)       &  .OR.(viscC2leithD.NE.0.)
137       &  .OR.(viscC2smag.NE.0.)       &  .OR.(viscC2smag.NE.0.)
138    
139        IF (harmonic) viscAhremax=50        IF (harmonic) viscAhre_max=viscAhremax
140    
141        biharmonic=        biharmonic=
142       &      (viscA4.NE.0.)       &      (viscA4.NE.0.)
# Line 147  C     == Local variables == Line 147  C     == Local variables ==
147       &  .OR.(viscC4leithD.NE.0.)       &  .OR.(viscC4leithD.NE.0.)
148       &  .OR.(viscC4smag.NE.0.)       &  .OR.(viscC4smag.NE.0.)
149    
150        IF (biharmonic) viscA4remax=50        IF (biharmonic) viscA4re_max=viscA4remax
151    
152        calcleith=        calcleith=
153       &      (viscC2leith.NE.0.)       &      (viscC2leith.NE.0.)
# Line 190  C These are (powers of) length scales Line 190  C These are (powers of) length scales
190       &                        *recip_DYF(I,J,bi,bj))**2) )       &                        *recip_DYF(I,J,bi,bj))**2) )
191    
192  C Velocity Reynolds Scale  C Velocity Reynolds Scale
193           Uscl=sqrt(KE(i,j)*L2*0.5)/viscAhRemax           Uscl=sqrt(KE(i,j)*L2*0.5)/viscAhRe_max
194           U4scl=0.125*L2*Uscl/viscA4Remax           U4scl=0.125*L2*Uscl/viscA4Re_max
195    
196           IF (useFullLeith.and.calcleith) THEN           IF (useFullLeith.and.calcleith) THEN
197  C This is the vector magnitude of the vorticity gradient squared  C This is the vector magnitude of the vorticity gradient squared
# Line 294  C These are (powers of) length scales Line 294  C These are (powers of) length scales
294    
295  C Velocity Reynolds Scale  C Velocity Reynolds Scale
296           Uscl=sqrt((KE(i,j)+KE(i,j+1)+KE(i+1,j)+KE(i+1,j+1))*L2*0.125)/           Uscl=sqrt((KE(i,j)+KE(i,j+1)+KE(i+1,j)+KE(i+1,j+1))*L2*0.125)/
297       &         viscAhRemax       &         viscAhRe_max
298           U4scl=0.125*L2*Uscl/viscA4Remax           U4scl=0.125*L2*Uscl/viscA4Re_max
299    
300  C This is the vector magnitude of the vorticity gradient squared  C This is the vector magnitude of the vorticity gradient squared
301           IF (useFullLeith.and.calcleith) THEN           IF (useFullLeith.and.calcleith) THEN
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22