Removed test code that causes CTEM to stop during an otherwise normal…

… run
MintonGroup · Nov 21, 2022 · 546615a · 546615a
1 parent 19e4cf2
commit 546615a
Show file tree

Hide file tree

Showing 3 changed files with 6 additions and 6 deletions.
diff --git 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
@@ -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
@@ -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)