/[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.29 by jmc, Mon Oct 22 13:24:50 2007 UTC revision 1.30 by mlosch, Wed Nov 14 15:55:26 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)
           L4rdt=0.03125 _d 0*recip_dt*L2**2  
318           ELSE           ELSE
319            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))
           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) )  
320           ENDIF           ENDIF
321           L3=(L2**1.5)           L3=(L2**1.5)
322           L4=(L2**2)           L4=(L2**2)
323           L5=(L2*L3)           L5=(L2**2.5)
324    
325           L2rdt=0.25 _d 0*recip_dt*L2           L2rdt=0.25 _d 0*recip_dt*L2
326    
327             IF (useAreaViscLength) THEN
328              L4rdt=0.03125 _d 0*recip_dt*L2**2
329             ELSE
330              L4rdt=recip_dt/( 6. _d 0*(recip_DXF(I,J,bi,bj)**4
331         &                            +recip_DYF(I,J,bi,bj)**4)
332         &                   +8. _d 0*((recip_DXF(I,J,bi,bj)
333         &                             *recip_DYF(I,J,bi,bj))**2) )
334             ENDIF
335    
336  C Velocity Reynolds Scale  C Velocity Reynolds Scale
337           IF ( viscAhRe_max.GT.0. .AND. KE(i,j).GT.0. ) THEN           IF ( viscAhRe_max.GT.0. .AND. KE(i,j).GT.0. ) THEN
338             Uscl=sqrt(KE(i,j)*L2)*viscAhRe_max             Uscl=sqrt(KE(i,j)*L2)*viscAhRe_max
# Line 442  CCCCCCCCCCCCC Vorticity Point Calculatio Line 446  CCCCCCCCCCCCC Vorticity Point Calculatio
446  C These are (powers of) length scales  C These are (powers of) length scales
447           IF (useAreaViscLength) THEN           IF (useAreaViscLength) THEN
448            L2=rAz(i,j,bi,bj)            L2=rAz(i,j,bi,bj)
           L4rdt=0.125 _d 0*recip_dt*rAz(i,j,bi,bj)**2  
449           ELSE           ELSE
450            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))
           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))  
451           ENDIF           ENDIF
452    
453           L3=(L2**1.5)           L3=(L2**1.5)
454           L4=(L2**2)           L4=(L2**2)
455           L5=(L2*L3)           L5=(L2**2.5)
456    
457           L2rdt=0.25 _d 0*recip_dt*L2           L2rdt=0.25 _d 0*recip_dt*L2
458             IF (useAreaViscLength) THEN
459              L4rdt=0.125 _d 0*recip_dt*rAz(i,j,bi,bj)**2
460             ELSE
461              L4rdt=recip_dt/
462         &     ( 6. _d 0*(recip_DXV(I,J,bi,bj)**4+recip_DYU(I,J,bi,bj)**4)
463         &      +8. _d 0*((recip_DXV(I,J,bi,bj)*recip_DYU(I,J,bi,bj))**2))
464             ENDIF
465    
466  C Velocity Reynolds Scale (Pb here at CS-grid corners !)  C Velocity Reynolds Scale (Pb here at CS-grid corners !)
467           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.29  
changed lines
  Added in v.1.30
  ViewVC Help
Powered by ViewVC 1.1.22