Version with debugging (commented out) of strain softening function

src/vrm.f90

@@ -304,6 +304,16 @@ select case (trim(plasticity_type))
                    straintot,e2d,e2dprev,flag_vrm_pb)
     endif
     friction_angle = phi*180.d0/pi
+
+    !if (friction_angle.lt.2.d0) then
+    !  write (*,*) 'plasticity_parameters: ',plasticity_parameters
+    !  write (*,*) 'friction_angle: ',friction_angle
+    !  write (*,*) 'phi: ',phi
+    !  write (*,*) 'plasticity_ss_type_phi: ',plasticity_ss_type_phi
+    !  write (*,*) 'e2d: ',e2d
+    !  write (*,*) 'e2dprev: ',e2dprev
+    !endif
+
     sin_phi = sin(phi)
     cos_phi = cos(phi)
     zeta    = cos_theta-sin_phi*sin_theta/sqrt3