From 0c7d692477cce86beefdf88659f254db9d85f176 Mon Sep 17 00:00:00 2001 From: daminton Date: Wed, 8 Feb 2017 14:58:25 +0000 Subject: [PATCH] Updated code to conserve mass by doing collapse stage before ejecta --- src/crater/crater_populate.f90 | 7 +++---- 1 file changed, 3 insertions(+), 4 deletions(-) diff --git a/src/crater/crater_populate.f90 b/src/crater/crater_populate.f90 index 18d6e77a..5d4c2662 100644 --- a/src/crater/crater_populate.f90 +++ b/src/crater/crater_populate.f90 @@ -180,6 +180,9 @@ 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) + ! Collapse any remaining unstable slopes + if (user%docollapse) call crater_slope_collapse(user,surf,crater,domain,ejbmass) + 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 ! craters are retained. The full ejecta_table_define function @@ -208,10 +211,6 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt nsincetally = nsincetally + 1 if (.not.user%testflag) call io_updatePbar("") - - ! Collapse any remaining unstable slopes - if (user%docollapse) call crater_slope_collapse(user,surf,crater,domain) - !if (user%docrustal_thinning) call crust_thin(user,surf,crater,domain,mdepth) ! Find out if the current crater is the largest or smallest and if so record it