diff --git a/src/make_cut.f90 b/src/make_cut.f90 index 19445a1ead2ba43c704a2f350379ea3ed64938e8..38004d2b472503d3a0b87147bb52f19a606dad54 100644 --- a/src/make_cut.f90 +++ b/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