/[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.6 by m_bates, Thu Aug 22 17:12:18 2013 UTC revision 1.7 by jmc, Thu Aug 22 23:01:47 2013 UTC
# Line 376  C      The bottom of the grid cell is sh Line 376  C      The bottom of the grid cell is sh
376  C      integration level, so, advance the depth.  C      integration level, so, advance the depth.
377         IF (-rF(k+1).LE. GM_K3D_EadyMinDepth) CYCLE         IF (-rF(k+1).LE. GM_K3D_EadyMinDepth) CYCLE
378    
379  C      Don't bother going any deeper since the top of the  C      Do not bother going any deeper since the top of the
380  C      cell is deeper than the bottom integration level  C      cell is deeper than the bottom integration level
381         IF (-rF(k).GE.GM_K3D_EadyMaxDepth) EXIT         IF (-rF(k).GE.GM_K3D_EadyMaxDepth) EXIT
382    
# Line 679  C     This is necessary to find the eige Line 679  C     This is necessary to find the eige
679          ENDDO          ENDDO
680         ENDDO         ENDDO
681  C     This fudge is necessary to avoid division by zero in gmredi_calc_eigs.  C     This fudge is necessary to avoid division by zero in gmredi_calc_eigs.
682  C     It does not affect the end result since it's in the overlap region.  C     It does not affect the end result since it is in the overlap region.
683         j=1-Oly         j=1-Oly
684         DO i=1-Olx,sNx+Olx         DO i=1-Olx,sNx+Olx
685          N2W(i,j,k) = GM_K3D_minN2          N2W(i,j,k) = GM_K3D_minN2
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22