Skip to content
Commits on Source (2)
...@@ -174,8 +174,8 @@ contains ...@@ -174,8 +174,8 @@ contains
deallocate (work, iwork, stat=status) deallocate (work, iwork, stat=status)
call error_check (status, 'sh_diag: deallocation error') call error_check (status, 'sh_diag: deallocation error')
write (fo,'(/,a,i0,a,i0,a,5g14.6)') ' Mesh ',imesh,', Sector ',sect,': first five eigenvalues = ', evals(1:5) write (fo,'(/,a,i0,a,i0,a,5f14.4)') ' Mesh ',imesh,', Sector ',sect,': first five eigenvalues = ', evals(1:5)
write (fo,'(a,i0,a,i0,a,5g14.6,/)') ' Mesh ',imesh,', Sector ',sect,': final five eigenvalues = ', evals(ny-4:ny) write (fo,'(a,i0,a,i0,a,5f14.4,/)') ' Mesh ',imesh,', Sector ',sect,': final five eigenvalues = ', evals(ny-4:ny)
z = a z = a
evecs: do ii = 1, ny ! normalize eigenvectors evecs: do ii = 1, ny ! normalize eigenvectors
......
...@@ -174,8 +174,8 @@ contains ...@@ -174,8 +174,8 @@ contains
deallocate (work, iwork, stat=status) deallocate (work, iwork, stat=status)
call error_check (status, 'sh_diag: deallocation error') call error_check (status, 'sh_diag: deallocation error')
write (fo,'(/,a,i0,a,i0,a,5g14.6)') ' Mesh ',imesh,', Sector ',sect,': first five eigenvalues = ', evals(1:5) write (fo,'(/,a,i0,a,i0,a,5f14.4)') ' Mesh ',imesh,', Sector ',sect,': first five eigenvalues = ', evals(1:5)
write (fo,'(a,i0,a,i0,a,5g14.6,/)') ' Mesh ',imesh,', Sector ',sect,': final five eigenvalues = ', evals(ny-4:ny) write (fo,'(a,i0,a,i0,a,5f14.4,/)') ' Mesh ',imesh,', Sector ',sect,': final five eigenvalues = ', evals(ny-4:ny)
z = a z = a
evecs: do ii = 1, ny ! normalize eigenvectors evecs: do ii = 1, ny ! normalize eigenvectors
......