Added debug output

src/make_cut.f90

@@ -334,6 +334,10 @@ zmean=zmean/params%mpe
 select case (matrule)
 case (1)                                                                        ! Assign material properties of volumetric majority material
   matel=params%materialn(sum(maxloc(vol_lsf)))
+  write(*,*) 'zmean: ',zmean
+  write(*,*) 'ztrans: ',ztrans
+  write(*,*) 'matel: ',matel
+  write(*,*) 'transnum: ',mat(matel)%transnum
   if (zmean.lt.mat(matel)%ztrans) matel=mat(matel)%transnum
   matnum=matel
   viscosity=mat(matel)%viscosity