diff --git a/j2plast.c b/j2plast.c
index 1714213b921cd7b7b13083bd1896287ae5339e5a..9193b1e425c87532af34e02e35f0e3e5eed5c83a 100644
--- a/j2plast.c
+++ b/j2plast.c
@@ -1115,7 +1115,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
 	//POWER LAW HARDENING
 	case H_POW:
 		if(p<DP_MIN){
-			printf("Should not be here: p: %.10e\n",p);
+			//printf("Should not be here: p: %.10e\n",p);
 			*R=0.;
 			if(hexp<1.0){
 				*dR = 1e20;
@@ -1136,7 +1136,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
 	//POWER LAW HARDENING FOR DRUCKER PRAGER
 	case H_POW_DP:
 		if(p<DP_MIN){
-			printf("Should not be here: p: %.10e\n",p);
+			//printf("Should not be here: p: %.10e\n",p);
 			*R=0.;
 			if(hexp<1.0){
 				*dR = 1e20;
@@ -1159,7 +1159,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
 	case H_EXP:
 	
 		if(p<DP_MIN){
-			printf("Should not be here: p: %.10e\n",p);
+			//printf("Should not be here: p: %.10e\n",p);
 			*R=0.;
 			*dR = hmod1*hexp;
 			*ddR = -hexp*(*dR);
@@ -1178,7 +1178,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
 	case H_EXP_DP:
 	
 		if(p<DP_MIN){
-			printf("Should not be here: p: %.10e\n",p);
+			//printf("Should not be here: p: %.10e\n",p);
 			*R=0.;
 			*dR = hmod1*hexp;
 			*ddR = -hexp*(*dR);
@@ -1195,7 +1195,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
 	//SWIFT LAW
 	case H_SWIFT:
 		if(p<DP_MIN){
-			printf("Should not be here: p: %.10e\n",p);
+			//printf("Should not be here: p: %.10e\n",p);
 			*R=0.;
 			*dR = sy0*hmod1*hexp;
 			*ddR = (*dR)*(hexp-1.)*hmod1;
@@ -1212,7 +1212,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
 	//SWIFT LAW FOR DRUCKER PRAGER
 	case H_SWIFT_DP:
 		if(p<DP_MIN){
-			printf("Should not be here: p: %.10e\n",p);
+			//printf("Should not be here: p: %.10e\n",p);
 			*R=0.;
 			*dR = sy0*hmod1*hexp;
 			*ddR = (*dR)*(hexp-1.)*hmod1;
@@ -1230,7 +1230,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
 	case H_LINEXP:
 	
 		if(p<DP_MIN){
-			printf("Should not be here: p: %.10e\n",p);
+			//printf("Should not be here: p: %.10e\n",p);
 			*R = 0;
 			*dR = hmod1 + hmod2*hexp;
 			*ddR = -hmod2*hexp*hexp;
@@ -1247,7 +1247,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
 	case H_LINEXP_DP:
 	
 		if(p<DP_MIN){
-			printf("Should not be here: p: %.10e\n",p);
+			//printf("Should not be here: p: %.10e\n",p);
 			*R = 0;
 			*dR = hmod1 + hmod2*hexp;
 			*ddR = -hmod2*hexp*hexp;
@@ -1263,7 +1263,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
 	//POWER LAW HARDENING EXTRAPOLATED AFTER 16% DEFO TO MIMIC DIGIMAT TO ABAQUS
 	case H_POW_EXTRAPOL:
 		if(p<DP_MIN){
-			printf("Should not be here: p: %.10e\n",p);
+			//printf("Should not be here: p: %.10e\n",p);
 			*R=0.;
 			if(hmod1<1){
 				*dR = 1e20;
@@ -1299,7 +1299,7 @@ void j2hard (double p, double sy0, int htype, double hmod1, double hmod2, double
 	//POWER LAW HARDENING EXTRAPOLATED AFTER 16% DEFO TO MIMIC DIGIMAT TO ABAQUS DRUCKER PRAGER
 	case H_POW_EXTRAPOL_DP:
 		if(p<DP_MIN){
-			printf("Should not be here: p: %.10e\n",p);
+			//printf("Should not be here: p: %.10e\n",p);
 			*R=0.;
 			if(hmod1<1){
 				*dR = 1e20;