/[MITgcm]/MITgcm/pkg/gmredi/gmredi_k3d.F
ViewVC logotype

Diff of /MITgcm/pkg/gmredi/gmredi_k3d.F

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

revision 1.5 by m_bates, Thu Aug 8 22:39:36 2013 UTC revision 1.6 by m_bates, Thu Aug 22 17:12:18 2013 UTC
# Line 88  C     small  :: a small number (to avoid Line 88  C     small  :: a small number (to avoid
88    
89  C     dfdy    :: gradient of the Coriolis paramter, df/dy, 1/(m*s)  C     dfdy    :: gradient of the Coriolis paramter, df/dy, 1/(m*s)
90  C     dfdx    :: gradient of the Coriolis paramter, df/dx, 1/(m*s)  C     dfdx    :: gradient of the Coriolis paramter, df/dx, 1/(m*s)
91  C     gradf   :: total gradient of the Coriolis paramter, SQRT(df/dx**2+df/dy**2), 1/(m*s)  C     gradf   :: gradient of the Coriolis paramter at a cell centre, 1/(m*s)
92  C     RRhines :: The Rhines scale (m)  C     RRhines :: The Rhines scale (m)
93  C     Rmix    :: Mixing length  C     Rmix    :: Mixing length
94  C     N2      :: Square of the buoyancy frequency (1/s**2)  C     N2      :: Square of the buoyancy frequency (1/s**2)
# Line 250  C       Not limited Line 250  C       Not limited
250        ENDDO        ENDDO
251        DO j=1-Oly,sNy+Oly        DO j=1-Oly,sNy+Oly
252         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
253          gradf(i,j) = SQRT( dfdx(i,j)**2 + dfdy(i,j)**2 )          gradf(i,j) = recip_rSphere*fCoriCos(i,j,bi,bj)
254         ENDDO         ENDDO
255        ENDDO        ENDDO
256    
# Line 413  C       If the minimum depth for the int Line 413  C       If the minimum depth for the int
413  C       bottom then give the eady growth rate a dummy, non-zero value  C       bottom then give the eady growth rate a dummy, non-zero value
414  C       to avoid floating point exceptions.  These points are taken care  C       to avoid floating point exceptions.  These points are taken care
415  C       of by setting K3D=GM_K3D_smallK later.  C       of by setting K3D=GM_K3D_smallK later.
416          IF (kLow_C(i,j).NE.0          IF (-r_Low(i,j,bi,bj).LE.GM_K3D_EadyMinDepth) THEN
      &       .AND. -r_Low(i,j,bi,bj).LT.GM_K3D_EadyMinDepth) THEN  
417            eady(i,j) = small            eady(i,j) = small
418    
419  C       Otherwise, multiply eady by the various constants to get the  C       Otherwise, multiply eady by the various constants to get the
# Line 469  C     Calculate the diffusivity (on the Line 468  C     Calculate the diffusivity (on the
468         DO j=1-Oly,sNy+Oly         DO j=1-Oly,sNy+Oly
469          DO i=1-Olx,sNx+Olx          DO i=1-Olx,sNx+Olx
470           IF (k.LE.kLow_C(i,j)) THEN           IF (k.LE.kLow_C(i,j)) THEN
471             IF (-r_Low(i,j,bi,bj).LT.GM_K3D_EadyMinDepth) THEN             IF (-r_Low(i,j,bi,bj).LE.GM_K3D_EadyMinDepth) THEN
472               K3D(i,j,k,bi,bj) = GM_K3D_smallK               K3D(i,j,k,bi,bj) = GM_K3D_smallK
473             ELSE             ELSE
474               IF (urms(i,j,k).EQ.0.0) THEN               IF (urms(i,j,k).EQ.0.0) THEN
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.22