/[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.30 by mlosch, Wed Nov 14 15:55:26 2007 UTC revision 1.31 by mlosch, Wed Nov 14 16:00:33 2007 UTC
# Line 315  CCCCCCCCCCCCCCC Divergence Point Calcula Line 315  CCCCCCCCCCCCCCC Divergence Point Calcula
315  C These are (powers of) length scales  C These are (powers of) length scales
316           IF (useAreaViscLength) THEN           IF (useAreaViscLength) THEN
317            L2=rA(i,j,bi,bj)            L2=rA(i,j,bi,bj)
318              L4rdt=0.03125 _d 0*recip_dt*L2**2
319           ELSE           ELSE
320            L2=2. _d 0/((recip_DXF(I,J,bi,bj)**2+recip_DYF(I,J,bi,bj)**2))            L2=2. _d 0/((recip_DXF(I,J,bi,bj)**2+recip_DYF(I,J,bi,bj)**2))
321              L4rdt=recip_dt/( 6. _d 0*(recip_DXF(I,J,bi,bj)**4
322         &                             +recip_DYF(I,J,bi,bj)**4)
323         &                   +8. _d 0*((recip_DXF(I,J,bi,bj)
324         &                             *recip_DYF(I,J,bi,bj))**2) )
325           ENDIF           ENDIF
326           L3=(L2**1.5)           L3=(L2**1.5)
327           L4=(L2**2)           L4=(L2**2)
328           L5=(L2**2.5)           L5=(L2*L3)
329    
330           L2rdt=0.25 _d 0*recip_dt*L2           L2rdt=0.25 _d 0*recip_dt*L2
331    
          IF (useAreaViscLength) THEN  
           L4rdt=0.03125 _d 0*recip_dt*L2**2  
          ELSE  
           L4rdt=recip_dt/( 6. _d 0*(recip_DXF(I,J,bi,bj)**4  
      &                            +recip_DYF(I,J,bi,bj)**4)  
      &                   +8. _d 0*((recip_DXF(I,J,bi,bj)  
      &                             *recip_DYF(I,J,bi,bj))**2) )  
          ENDIF  
   
332  C Velocity Reynolds Scale  C Velocity Reynolds Scale
333           IF ( viscAhRe_max.GT.0. .AND. KE(i,j).GT.0. ) THEN           IF ( viscAhRe_max.GT.0. .AND. KE(i,j).GT.0. ) THEN
334             Uscl=sqrt(KE(i,j)*L2)*viscAhRe_max             Uscl=sqrt(KE(i,j)*L2)*viscAhRe_max
# Line 446  CCCCCCCCCCCCC Vorticity Point Calculatio Line 442  CCCCCCCCCCCCC Vorticity Point Calculatio
442  C These are (powers of) length scales  C These are (powers of) length scales
443           IF (useAreaViscLength) THEN           IF (useAreaViscLength) THEN
444            L2=rAz(i,j,bi,bj)            L2=rAz(i,j,bi,bj)
445              L4rdt=0.125 _d 0*recip_dt*rAz(i,j,bi,bj)**2
446           ELSE           ELSE
447            L2=2. _d 0/((recip_DXV(I,J,bi,bj)**2+recip_DYU(I,J,bi,bj)**2))            L2=2. _d 0/((recip_DXV(I,J,bi,bj)**2+recip_DYU(I,J,bi,bj)**2))
448              L4rdt=recip_dt/
449         &     ( 6. _d 0*(recip_DXV(I,J,bi,bj)**4+recip_DYU(I,J,bi,bj)**4)
450         &      +8. _d 0*((recip_DXV(I,J,bi,bj)*recip_DYU(I,J,bi,bj))**2))
451           ENDIF           ENDIF
452    
453           L3=(L2**1.5)           L3=(L2**1.5)
454           L4=(L2**2)           L4=(L2**2)
455           L5=(L2**2.5)           L5=(L2*L3)
456    
457           L2rdt=0.25 _d 0*recip_dt*L2           L2rdt=0.25 _d 0*recip_dt*L2
          IF (useAreaViscLength) THEN  
           L4rdt=0.125 _d 0*recip_dt*rAz(i,j,bi,bj)**2  
          ELSE  
           L4rdt=recip_dt/  
      &     ( 6. _d 0*(recip_DXV(I,J,bi,bj)**4+recip_DYU(I,J,bi,bj)**4)  
      &      +8. _d 0*((recip_DXV(I,J,bi,bj)*recip_DYU(I,J,bi,bj))**2))  
          ENDIF  
458    
459  C Velocity Reynolds Scale (Pb here at CS-grid corners !)  C Velocity Reynolds Scale (Pb here at CS-grid corners !)
460           IF ( viscAhRe_max.GT.0. .OR. viscA4Re_max.GT.0. ) THEN           IF ( viscAhRe_max.GT.0. .OR. viscA4Re_max.GT.0. ) THEN
Legend:
Removed from v.1.30  
changed lines
  Added in v.1.31
  ViewVC Help
Powered by ViewVC 1.1.22