diff --git a/src/Makefile.am b/src/Makefile.am
index abc481ca..1bdb5932 100755
--- a/src/Makefile.am
+++ b/src/Makefile.am
@@ -7,7 +7,7 @@ OPTREPORT = -qopt-report=5
 IPRODUCTION = -g -traceback -no-wrap-margin -assume byterecl -O3 -qopt-prefetch=0 -sox $(PAR) $(SIMDVEC) $(HEAPARR)
 #IDEBUG = -O0 -g -traceback -debug all -nogen-interfaces -assume byterecl -m64 -heap-arrays -FR -no-pie -no-ftz -fpe-all=0 -mp1 -fp-model strict -fpe0 -align all -pad -ip -prec-div -prec-sqrt -assume protect-parens -CB -no-wrap-margin -init=snan,arrays 
 IDEBUG = -O0 -g -traceback -debug all -nogen-interfaces -assume byterecl -m64 -heap-arrays -FR -no-pie -no-ftz -fpe-all=0 -mp1 -fp-model strict -fpe0 -align all -pad -ip -prec-sqrt -assume protect-parens -CB -no-wrap-margin -init=snan,arrays 
-AM_FCFLAGS =  $(IDEBUG)
+AM_FCFLAGS =  $(IPRODUCTION)
 #ifort debug flags
 
 #gfortran optimized flags
diff --git a/src/ejecta/ejecta_ray_func.f90 b/src/ejecta/ejecta_ray_func.f90
index 693f6fd9..459fefde 100644
--- a/src/ejecta/ejecta_ray_func.f90
+++ b/src/ejecta/ejecta_ray_func.f90
@@ -38,7 +38,7 @@ function ejecta_ray_func(theta,thetar,r,n,w) result(ans)
     c = w / r
     b = thetar 
     dtheta = min(2*pi - abs(theta - b),abs(theta - b))
-    a = sqrt(2 * pi) / (n * c * erf(pi / (2 *sqrt(2._DP) * c)))
+    a = sqrt(2 * pi) / (n * c * erf(pi / (2 *sqrt(2._DP) * c))) !this is the intensity function
     ans = a * exp(-dtheta**2 / (2 * c**2))
 
 !return
diff --git a/src/regolith/regolith_streamtube.f90 b/src/regolith/regolith_streamtube.f90
index c5689b43..0ae14b88 100644
--- a/src/regolith/regolith_streamtube.f90
+++ b/src/regolith/regolith_streamtube.f90
@@ -240,10 +240,10 @@ subroutine regolith_streamtube(user,surf,crater,domain,ejb,ejtble,xp,yp,xpi,ypi,
    tots    = 0.0_DP
    depthb = crater%imp / 2.0_DP
 
-   if (eradc<=user%testimp) then
-      write(*,*) lrad/crater%frad, user%testimp, crater%frad, rm, deltar, eradc, eradi, erado, ebh, newlayer%meltfrac
-      stop
-   end if
+   ! if (eradc<=user%testimp) then
+   !    write(*,*) lrad/crater%frad, user%testimp, crater%frad, rm, deltar, eradc, eradi, erado, ebh, newlayer%meltfrac
+   !    stop
+   ! end if
 
    call regolith_shock_damage_zone(crater,rm,eradi,depthb,xsfints)