From 3df361b3ab07a7ac62b778b8b3bb746a52cb8ec8 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Fri, 15 Jul 2022 10:57:20 -0400 Subject: [PATCH] Check if a crater does not generate topography. If not, skip it. --- src/crater/crater_populate.f90 | 6 ++++++ 1 file changed, 6 insertions(+) diff --git a/src/crater/crater_populate.f90 b/src/crater/crater_populate.f90 index e60b68ce..c3c779c3 100644 --- a/src/crater/crater_populate.f90 +++ b/src/crater/crater_populate.f90 @@ -240,6 +240,12 @@ 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) < tiny(ejbmass)) then ! Crater has no topography. Discard and move on. + ntrue = ntrue -1 + mass = mass - crater%impmass + cycle + end if + 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