Edited to only provide debug output from process 0

src/refine_surface.f90

@@ -350,7 +350,7 @@ if (iproc.eq.0 .and. params%debug>=1) write(*,'(a,i2,a,i4,a)') shift//'S.',is,':
 do i=1,size(ed%n1)
    i1=ed(i)%n1
    i2=ed(i)%n2
-   if (i1.eq.i2) then
+   if (i1 == i2 .and. iproc == 0) then
      print *,'*******************************************************************'
      print *,'Nodes for edge ',i,' have the same spatial coordinates'
      print *,'The edges nodes are are node i1: ',i1,' and i2: ',i2