diff --git a/src/vrm.f90 b/src/vrm.f90
index 32a4adef8c11f51b3eaddd62610268322d09248d..eb3217b7d6e097d34c4bf2d0812f0633589be83a 100644
--- a/src/vrm.f90
+++ b/src/vrm.f90
@@ -264,6 +264,18 @@ select case (trim(plasticity_type))
     fail  = 2.d0*viscosity*e2d-yield
     if (fail.gt.0.d0) then
       viscosity=0.5d0*yield/e2d
+      if (viscosity.le.params%viscositymin) then
+        write(*,*) 'VISCOSITY BELOW MINIMUM'
+        write(*,*) 'viscosity: ',viscosity
+        write(*,*) 'yield: ',yield
+        write(*,*) 'fail: ',fail
+        write(*,*) '-pressure: ',-pressure
+        write(*,*) 'sin_phi: ',sin_phi
+        write(*,*) 'c: ',c
+        write(*,*) 'cos_phi: ',cos_phi
+        write(*,*) 'zeta: ',zeta
+        write(*,*) 'e2d: ',e2d
+      endif
       is_plastic=.true.
     else
       is_plastic=.false.