Tweaked wall texture to try to prevent sharp edges at bounary of ejecta

MintonGroup · Nov 2, 2021 · d9778c9 · d9778c9
1 parent 5456394
commit d9778c9
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/crater/crater_realistic_topography.f90 b/src/crater/crater_realistic_topography.f90
@@ -350,7 +350,7 @@ subroutine complex_wall_texture(user,surf,crater,domain,deltaMtot)
                                               xynoise * ybar + offset * rn(2))* znoise
          end do
          newdem = max(newdem + noise * areafrac,crater%melev - crater%floordepth)
-         if (r > 1.1_DP) newdem = max(newdem,crater%melev + crater%ejrim * r**(-3))
+         if (r > 1.1_DP) newdem = max(newdem,newdem + areafrac * crater%ejrim * r**(-3))
 
          elchange  = newdem - surf(xpi,ypi)%dem
          deltaMtot = deltaMtot + elchange