Adjusted peak ring radius scaling so it maxes out at 0.5 frad

MintonGroup · Nov 3, 2021 · 031f3b9 · 031f3b9
1 parent 481d2cd
commit 031f3b9
Showing 1 changed file with 1 addition and 0 deletions.
diff --git a/src/crater/crater_realistic_topography.f90 b/src/crater/crater_realistic_topography.f90
@@ -172,6 +172,7 @@ subroutine complex_peak(user,surf,crater,deltaMtot)
    !FWHM = 0.3_DP !Lansberg 
    a = crater%peakheight
    b = 0.003_DP * ((1e-3_DP * crater%fcrat)**(1.75_DP)) / (1e-3_DP * crater%fcrat) ! Make peak rings for sufficiently large craters
+   b = min(b, 0.5_DP)
    c = FWHM / (2 * sqrt(2 * log(2._DP)))
    !*********************