FIxed issue that happens when a crater isn't big enough to affect the…

… surface
MintonGroup · Dec 14, 2022 · 88dd83b · 88dd83b
1 parent a6249f2
commit 88dd83b
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 0 deletions.
diff --git a/src/crater/crater_populate.f90 b/src/crater/crater_populate.f90
@@ -247,6 +247,7 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt
 
          ! Place crater onto the surface
          call crater_emplace(user,surf,crater,domain,ejbmass)
+         if (abs(ejbmass) < 2*tiny(1.0_DP)) cycle
 
          call ejecta_distance_estimate(user,crater,domain,crater%ejdis) ! Fast but imprecise estimate of the total ejecta distance
                                                                         ! For very steep size distributions, only a fraction of the

diff --git a/src/ejecta/ejecta_emplace.f90 b/src/ejecta/ejecta_emplace.f90
@@ -307,6 +307,7 @@ subroutine ejecta_emplace(user,surf,crater,domain,ejb,ejtble,deltaMtot,age,age_r
    fmasscons = (-deltaMtot)/ ejbmass
    cumulative_elchange = cumulative_elchange * fmasscons
    crater%ejrim = crater%ejrim * fmasscons
+   if (abs(fmasscons) < tiny(1.0_DP)) return
    ejb(:)%thick = ejb(:)%thick + log(fmasscons)
    maxhits = 1
    ! Create box for soften calculation (will be no bigger than the grid itself)