diff --git a/dG3D/benchmarks/CrackPhaseField3D/pf.py b/dG3D/benchmarks/CrackPhaseField3D/pf.py
index 8cb587a4b99b427091a86502e5b45cc4e013ebe4..fda338eba7c6e618153c3d18452630253cdd0135 100644
--- a/dG3D/benchmarks/CrackPhaseField3D/pf.py
+++ b/dG3D/benchmarks/CrackPhaseField3D/pf.py
@@ -80,7 +80,7 @@ if Solving_mode == 1:
 	ftime =1.0   		# Final time
 	tol=1.e-5  		# Relative tolerance for NR scheme
 	tolAbs = 1.e-6		# Absolute tolerance for NR scheme
-	nstepArch= 1	# Number of step between 2 archiving
+	nstepArch= 10	# Number of step between 2 archiving
 	nstepArchEnergy = 1	# Number of step between 2 energy computation and archiving
 	nstepArchForce = 1	# Number of step between 2 force archiving
 	MaxIter = 200		# Maximum number of iterations
@@ -95,7 +95,7 @@ elif Solving_mode == 2:
 	ftime =1.0e-2   	# Final time
 	tol=1.e-4   		# Relative tolerance for NR scheme
 	tolAbs = 1.e-6		# Absolute tolerance for NR scheme
-	nstepArch= 1 	# Number of step between 2 archiving
+	nstepArch= 10 	# Number of step between 2 archiving
 	nstepArchEnergy = 1 # Number of step between 2 energy computation and archiving
 	nstepArchForce = 1	# Number of step between 2 force archiving
 	MaxIter = 200		# Maximum number of iterations
@@ -232,10 +232,11 @@ if LoadingMode == 0:
 		mysolver.setPathFollowingIncrementAdaptation(True,4,0.3)
 		mysolver.setPathFollowingControlType(0)
 		mysolver.setPathFollowingCorrectionMethod(0)
-		mysolver.setPathFollowingArcLengthStep(1e-1)
-		mysolver.setBoundsOfPathFollowingArcLengthSteps(0,2.);
-
+		mysolver.setPathFollowingArcLengthStep(3e-2)
+		mysolver.setBoundsOfPathFollowingArcLengthSteps(0,3e-2);
 
+	mysolver.archivingForceOnPhysicalGroup("Face", 400, 1, nstepArchForce)
+	mysolver.archivingForceOnPhysicalGroup("Face", 200, 1, nstepArchForce)
 	# mysolver.pathFollowing(True,1)
 	# mysolver.setPathFollowingIncrementAdaptation(False,4)
 	# mysolver.setPathFollowingLocalSteps(9.0e-1,1.0e-3)
@@ -321,7 +322,7 @@ mysolver.archivingNodeDisplacement(111,1, nstepArchForce)
 mysolver.solve()
 
 check = TestCheck()
-check.equal(-6.645770e+05,mysolver.getArchivedForceOnPhysicalGroup("face", 400, 1),1.e-4)
+check.equal(-2966.175372583341,mysolver.getArchivedForceOnPhysicalGroup("Face", 400, 1),1.e-4)