Added ability to alter ejecta thickness based on the relative slope o…

…f the impact site to the downrange location.

src/crater/crater_populate.f90

@@ -173,7 +173,7 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt
          if (user%doseismic) call seismic_shake(user,surf,crater,domain)
 
          ! Generate dynamic diffusion
-         if (user%dosoftening) call crater_soften(user,surf,crater,domain)
+         !if (user%dosoftening) call crater_soften(user,surf,crater,domain)
 
          ! find the average height and slope at crater location
          call crater_averages(user,surf,crater)
@@ -194,12 +194,11 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt
                call ejecta_table_define(user,crater,domain,ejb,ejtble)
                call ejecta_interpolate(crater,domain,crater%frad,ejb(1:ejtble),ejtble,crater%ejrim)
             end if
-            !call ejecta_emplace(user,surf,crater,domain,ejb(1:ejtble),ejtble,ejbmass)
+            call ejecta_emplace(user,surf,crater,domain,ejb(1:ejtble),ejtble,ejbmass)
          else
             ejtble = 0
          end if
 
-
          ! Place crater onto the surface
          if (crater%fcrat > domain%smallest_crater) then
             call crater_emplace(user,surf,crater,domain,ejbmass)
@@ -215,6 +214,7 @@ subroutine crater_populate(user,surf,crater,domain,prod,production_list,vdist,nt
          end if
 
 
+
          ! Collapse any remaining unstable slopes
          if (user%docollapse) call crater_slope_collapse(user,surf,crater,domain)