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

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

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

revision 1.4 by jmc, Thu Aug 22 23:01:47 2013 UTC revision 1.5 by m_bates, Sun May 18 02:49:52 2014 UTC
# Line 50  C     bi, bj    :: tile indices Line 50  C     bi, bj    :: tile indices
50  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
51  C     == Local variables ==  C     == Local variables ==
52        INTEGER i,j,k,kk,m        INTEGER i,j,k,kk,m
53  # ifdef use_lapack  # ifdef HAVE_LAPACK
54  C     info     :: error code from LAPACK  C     info     :: error code from LAPACK
55  C     idx      :: index used for sorting the eigenvalues  C     idx      :: index used for sorting the eigenvalues
56  C     a3d      :: lower diagonal of eigenvalue problem  C     a3d      :: lower diagonal of eigenvalue problem
# Line 114  C     Square of the Coriolis parameter Line 114  C     Square of the Coriolis parameter
114         ENDDO         ENDDO
115        ENDDO        ENDDO
116    
117  # ifdef use_lapack  # ifdef HAVE_LAPACK
118  C     Calculate the tridiagonal operator matrix for  C     Calculate the tridiagonal operator matrix for
119  C     f^2 d/dz 1/N^2 d/dz  C     f^2 d/dz 1/N^2 d/dz
120  C     a3d is the lower off-diagonal, b3d is the diagonal  C     a3d is the lower off-diagonal, b3d is the diagonal
# Line 393  C     Normalise the eigenvectors Line 393  C     Normalise the eigenvectors
393  # ifdef ALLOW_DIAGNOSTICS  # ifdef ALLOW_DIAGNOSTICS
394  C     Diagnostics  C     Diagnostics
395        IF ( useDiagnostics.AND.writediag ) THEN        IF ( useDiagnostics.AND.writediag ) THEN
396  #  ifdef use_lapack          #  ifdef HAVE_LAPACK        
397          CALL DIAGNOSTICS_FILL(a3d, 'GM_A3D  ',0,Nr,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(a3d, 'GM_A3D  ',0,Nr,0,1,1,myThid)
398          CALL DIAGNOSTICS_FILL(b3d, 'GM_B3D  ',0,Nr,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(b3d, 'GM_B3D  ',0,Nr,0,1,1,myThid)
399          CALL DIAGNOSTICS_FILL(c3d, 'GM_C3D  ',0,Nr,0,1,1,myThid)          CALL DIAGNOSTICS_FILL(c3d, 'GM_C3D  ',0,Nr,0,1,1,myThid)
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22