diff --git a/check_delaunay.f90 b/check_delaunay.f90
index 87e4098de80b5944c137850c531fa21b559aa6c1..88234c544e31695be904c98193411a3d409d7a35 100644
--- a/check_delaunay.f90
+++ b/check_delaunay.f90
@@ -172,15 +172,20 @@ do while (nflip > 0)
          j2=ed(iedge)%m2
 
          dn1n2=dist(surface%x (i1),surface%x (i2),surface%y (i1),surface%y (i2),surface%z (i1),surface%z (i2), &
-                    surface%xn(i1),surface%xn(i2),surface%yn(i1),surface%yn(i2),surface%zn(i1),surface%zn(i2), params%distance_exponent)
+                    surface%xn(i1),surface%xn(i2),surface%yn(i1),surface%yn(i2),surface%zn(i1),surface%zn(i2), &
+                    params%distance_exponent)
          dn1m1=dist(surface%x (i1),surface%x (j1),surface%y (i1),surface%y (j1),surface%z (i1),surface%z (j1), &
-                    surface%xn(i1),surface%xn(j1),surface%yn(i1),surface%yn(j1),surface%zn(i1),surface%zn(j1), params%distance_exponent)
+                    surface%xn(i1),surface%xn(j1),surface%yn(i1),surface%yn(j1),surface%zn(i1),surface%zn(j1), &
+                    params%distance_exponent)
          dn1m2=dist(surface%x (i1),surface%x (j2),surface%y (i1),surface%y (j2),surface%z (i1),surface%z (j2), &
-                    surface%xn(i1),surface%xn(j2),surface%yn(i1),surface%yn(j2),surface%zn(i1),surface%zn(j2), params%distance_exponent)
+                    surface%xn(i1),surface%xn(j2),surface%yn(i1),surface%yn(j2),surface%zn(i1),surface%zn(j2), &
+                    params%distance_exponent)
          dn2m1=dist(surface%x (i2),surface%x (j1),surface%y (i2),surface%y (j1),surface%z (i2),surface%z (j1), &
-                    surface%xn(i2),surface%xn(j1),surface%yn(i2),surface%yn(j1),surface%zn(i2),surface%zn(j1), params%distance_exponent)
+                    surface%xn(i2),surface%xn(j1),surface%yn(i2),surface%yn(j1),surface%zn(i2),surface%zn(j1), &
+                    params%distance_exponent)
          dn2m2=dist(surface%x (i2),surface%x (j2),surface%y (i2),surface%y (j2),surface%z (i2),surface%z (j2), &
-                    surface%xn(i2),surface%xn(j2),surface%yn(i2),surface%yn(j2),surface%zn(i2),surface%zn(j2), params%distance_exponent)
+                    surface%xn(i2),surface%xn(j2),surface%yn(i2),surface%yn(j2),surface%zn(i2),surface%zn(j2), &
+                    params%distance_exponent)
 
          dist1=(dn1m1**2+dn2m1**2-dn1n2**2)/(2.d0*dn1m1*dn2m1)
          dist2=(dn1m2**2+dn2m2**2-dn1n2**2)/(2.d0*dn1m2*dn2m2)