Overhauled crater formation routine. All raised rims are now part of …

…the ejecta blanket
MintonGroup · Feb 9, 2017 · 2cf4bd6 · 2cf4bd6
1 parent e650686
commit 2cf4bd6
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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