Skip to content

Commit

Permalink
Check if a crater does not generate topography. If not, skip it.
Browse files Browse the repository at this point in the history
  • Loading branch information
daminton committed Jul 15, 2022
1 parent 7b427f6 commit 3df361b
Showing 1 changed file with 6 additions and 0 deletions.
6 changes: 6 additions & 0 deletions src/crater/crater_populate.f90
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down

0 comments on commit 3df361b

Please sign in to comment.