/[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.11 by m_bates, Sat Sep 28 17:59:30 2013 UTC revision 1.12 by m_bates, Wed Oct 9 19:13:21 2013 UTC
# Line 416  C      We are in the integration depth r Line 416  C      We are in the integration depth r
416           IF ( (kLow_C(i,j).GE.k) .AND.           IF ( (kLow_C(i,j).GE.k) .AND.
417       &        (-hMixLayer(i,j,bi,bj).LE.-rC(k)) ) THEN       &        (-hMixLayer(i,j,bi,bj).LE.-rC(k)) ) THEN
418    
419             slopeC(i,j,k) = SQRT(SQRT(M4loc(i,j,k))/N2loc(i,j,k))             slopeC(i,j,k) = SQRT(M4loc(i,j,k))/N2loc(i,j,k)
420  C          Limit the slope.  Note, this is not all the Eady calculations.  C          Limit the slope.  Note, this is not all the Eady calculations.
421             IF (slopeC(i,j,k).LE.GM_K3D_maxSlope) THEN             IF (slopeC(i,j,k).LE.GM_K3D_maxSlope) THEN
422               eady(i,j) = eady(i,j)               eady(i,j) = eady(i,j)
# Line 425  C          Limit the slope.  Note, this Line 425  C          Limit the slope.  Note, this
425               slopeC(i,j,k) = GM_K3D_maxSlope               slopeC(i,j,k) = GM_K3D_maxSlope
426               eady(i,j) = eady(i,j)               eady(i,j) = eady(i,j)
427       &            + hfacC(i,j,k,bi,bj)*drF(k)*SQRT(M4loc(i,j,k))       &            + hfacC(i,j,k,bi,bj)*drF(k)*SQRT(M4loc(i,j,k))
428       &              *GM_K3D_maxSlope*GM_K3D_maxSlope       &              *GM_K3D_maxSlope
429             ENDIF             ENDIF
430             deltaH(i,j) = deltaH(i,j) + drF(k)             deltaH(i,j) = deltaH(i,j) + drF(k)
431           ENDIF           ENDIF
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.12
  ViewVC Help
Powered by ViewVC 1.1.22