From c88588d226e6e629c4ac9a8a3804a25df8d9a9c3 Mon Sep 17 00:00:00 2001
From: Dave Whipp <dwhipp@dal.ca>
Date: Wed, 26 Oct 2011 20:20:34 +0000
Subject: [PATCH] Added debug output
---
src/vrm.f90 | 12 ++++++++++++
1 file changed, 12 insertions(+)
diff --git a/src/vrm.f90 b/src/vrm.f90
index 32a4adef..eb3217b7 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.
--
GitLab