/[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.40 by dfer, Thu Sep 17 16:57:13 2009 UTC revision 1.41 by mlosch, Fri Nov 13 08:07:50 2009 UTC
# Line 345  C These are (powers of) length scales Line 345  C These are (powers of) length scales
345           L4rdt = L4rdt_D(i,j,bi,bj)           L4rdt = L4rdt_D(i,j,bi,bj)
346           L5 = (L2*L3)           L5 = (L2*L3)
347    
348    #ifndef AUTODIFF_DISABLE_REYNOLDS_SCALE
349  C Velocity Reynolds Scale  C Velocity Reynolds Scale
350           IF ( viscAhRe_max.GT.0. .AND. KE(i,j).GT.0. ) THEN           IF ( viscAhRe_max.GT.0. .AND. KE(i,j).GT.0. ) THEN
351             Uscl=SQRT(KE(i,j)*L2)*viscAhRe_max             Uscl=SQRT(KE(i,j)*L2)*viscAhRe_max
# Line 356  C Velocity Reynolds Scale Line 357  C Velocity Reynolds Scale
357           ELSE           ELSE
358             U4scl=0.             U4scl=0.
359           ENDIF           ENDIF
360    #endif /* ndef AUTODIFF_DISABLE_REYNOLDS_SCALE */
361    
362  cph-leith#ifndef ALLOW_AUTODIFF_TAMC  cph-leith#ifndef ALLOW_AUTODIFF_TAMC
363  #ifndef AUTODIFF_DISABLE_LEITH  #ifndef AUTODIFF_DISABLE_LEITH
# Line 471  C These are (powers of) length scales Line 473  C These are (powers of) length scales
473           L4rdt = L4rdt_Z(i,j,bi,bj)           L4rdt = L4rdt_Z(i,j,bi,bj)
474           L5 = (L2*L3)           L5 = (L2*L3)
475    
476    #ifndef AUTODIFF_DISABLE_REYNOLDS_SCALE
477  C Velocity Reynolds Scale (Pb here at CS-grid corners !)  C Velocity Reynolds Scale (Pb here at CS-grid corners !)
478           IF ( viscAhRe_max.GT.0. .OR. viscA4Re_max.GT.0. ) THEN           IF ( viscAhRe_max.GT.0. .OR. viscA4Re_max.GT.0. ) THEN
479             keZpt=0.25 _d 0*( (KE(i,j)+KE(i-1,j-1))             keZpt=0.25 _d 0*( (KE(i,j)+KE(i-1,j-1))
# Line 486  C Velocity Reynolds Scale (Pb here at CS Line 489  C Velocity Reynolds Scale (Pb here at CS
489             Uscl =0.             Uscl =0.
490             U4scl=0.             U4scl=0.
491           ENDIF           ENDIF
492    #endif /* ndef AUTODIFF_DISABLE_REYNOLDS_SCALE */
493    
494  #ifndef AUTODIFF_DISABLE_LEITH  #ifndef AUTODIFF_DISABLE_LEITH
495  C This is the vector magnitude of the vorticity gradient squared  C This is the vector magnitude of the vorticity gradient squared
Legend:
Removed from v.1.40  
changed lines
  Added in v.1.41
  ViewVC Help
Powered by ViewVC 1.1.22