From 75752fccdfedba8860e4a1170e9c6be15126e317 Mon Sep 17 00:00:00 2001 From: daminton Date: Fri, 2 Dec 2016 21:35:43 +0000 Subject: [PATCH] Added ability to alter ejecta thickness based on the relative slope of the impact site to the downrange location. --- src/crater/crater_populate.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/crater/crater_populate.f90 b/src/crater/crater_populate.f90 index 8606bdbd..885f3f6b 100644 --- a/src/crater/crater_populate.f90 +++ b/src/crater/crater_populate.f90 @@ -173,7 +173,7 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt if (user%doseismic) call seismic_shake(user,surf,crater,domain) ! Generate dynamic diffusion - if (user%dosoftening) call crater_soften(user,surf,crater,domain) + !if (user%dosoftening) call crater_soften(user,surf,crater,domain) ! find the average height and slope at crater location call crater_averages(user,surf,crater) @@ -194,12 +194,11 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt call ejecta_table_define(user,crater,domain,ejb,ejtble) call ejecta_interpolate(crater,domain,crater%frad,ejb(1:ejtble),ejtble,crater%ejrim) end if - !call ejecta_emplace(user,surf,crater,domain,ejb(1:ejtble),ejtble,ejbmass) + call ejecta_emplace(user,surf,crater,domain,ejb(1:ejtble),ejtble,ejbmass) else ejtble = 0 end if - ! Place crater onto the surface if (crater%fcrat > domain%smallest_crater) then call crater_emplace(user,surf,crater,domain,ejbmass) @@ -215,6 +214,7 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt end if + ! Collapse any remaining unstable slopes if (user%docollapse) call crater_slope_collapse(user,surf,crater,domain)