Removed debug output

src/make_cut.f90

@@ -331,11 +331,8 @@ enddo
 select case (matrule)
 case (1)                                                                        ! Assign material properties of volumetric majority material
   matel=params%materialn(sum(maxloc(vol_lsf)))
-  write (*,*) 'Before c1 test 1'
   if (zmean.lt.mat(matel)%ztrans) matel=mat(matel)%transnum
-  write (*,*) 'After c1 test 1'
   matnum=matel
-  write (*,*) 'After c1 assignment'
   viscosity=mat(matel)%viscosity
   density=mat(matel)%density
   penal=mat(matel)%penalty
@@ -349,11 +346,8 @@ case (1)
 
 case(2)                                                                         ! Assign material properties of volumetric minority material
   matel=params%materialn(sum(minloc(vol_lsf)))
-  write (*,*) 'Before c2 test 1'  
   if (zmean.lt.mat(matel)%ztrans) matel=mat(matel)%transnum
-  write (*,*) 'After c2 test 1'
   matnum=matel
-  write (*,*) 'After c2 assignment'
   viscosity=mat(matel)%viscosity
   density=mat(matel)%density
   penal=mat(matel)%penalty
@@ -378,9 +372,7 @@ case default
   do i=1,nlsf
     matel=params%materialn(i)
     if (zmean.lt.mat(matel)%ztrans) matel=mat(matel)%transnum
-    write (*,*) 'Before c0 test 1'
     if (vol_lsf(i).ge.0.015625-eps) matnum=matnum+matel
-    write (*,*) 'After c0 test 1'
     viscosity=viscosity+vol_lsf(i)*mat(matel)%viscosity
     density=density+vol_lsf(i)*mat(matel)%density
     penal=penal+vol_lsf(i)*mat(matel)%penalty
@@ -391,9 +383,7 @@ case default
     heat=heat+vol_lsf(i)*mat(matel)%heat
   enddo
   matel=params%materialn(0)
-  write (*,*) 'Before c0 test 2'
   if (vol_lsf0.ge.0.015625-eps) matnum=matnum+100
-  write (*,*) 'After c0 test 2'
   viscosity=viscosity+vol_lsf0*mat(matel)%viscosity
   density=density+vol_lsf0*mat(matel)%density
   penal=penal+vol_lsf0*mat(matel)%penalty