From 2cf4bd64ddb984fa28939f45b778eb0ba6e3b9ac Mon Sep 17 00:00:00 2001 From: daminton Date: Thu, 9 Feb 2017 18:56:27 +0000 Subject: [PATCH] Overhauled crater formation routine. All raised rims are now part of the ejecta blanket --- src/ejecta/ejecta_rootfind.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/ejecta/ejecta_rootfind.f90 b/src/ejecta/ejecta_rootfind.f90 index 9e7f952b..716da6cf 100644 --- a/src/ejecta/ejecta_rootfind.f90 +++ b/src/ejecta/ejecta_rootfind.f90 @@ -153,7 +153,7 @@ subroutine ejecta_rootfind(user,crater,domain,erad,lrad,vejsq,ejang,firstrun) end if niter = i erad = resultat - if ((erad <= crater%frad).and.(lrad >= crater%frad).and.(erad >= 0._DP)) exit + if ((erad <= crater%rad).and.(lrad >= crater%rad).and.(erad >= 0._DP)) exit factor = 0.5_DP * (factor + 1.0_DP) ! Failed. Try again with a new factor end do everything