Replace pause statement to write message

geodynamics · Oct 18, 2024 · b88001f · b88001f
1 parent 05087f7
commit b88001f
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/src/Fortran_libraries/SOLVER_src/direct_solver/m_ludcmp_band.f90 b/src/Fortran_libraries/SOLVER_src/direct_solver/m_ludcmp_band.f90
@@ -115,7 +115,8 @@ subroutine ludcmp_band(n, nband, band_a, band_lu ,indx, d)
           aamax = max(aamax, abs(band_a(k,j)))
 !          aamax = max(aamax, abs(a(i,j)))
         end do
-        if (aamax.eq.0.0d0) pause 'singular matrix in ludcmp_band'
+        if (aamax.eq.0.0d0) write(*,*)                                  &
+     &                    '!!!!!!  singular matrix in ludcmp  !!!!!!'
         vv(i)=1.0d0 / aamax
         idx_org(i) = i
       end do