Added new user flag that controls soften size

src/globals/module_globals.f90

@@ -52,7 +52,7 @@ module module_globals
 integer(I4B),parameter :: MAXLAYER=20          ! Maximum number of layers (you need roughly 1-2 layers per order of magnitude of 
                                                ! resolution
 real(DP),parameter :: TALLYCOVERAGE = 0.01_DP   ! The total area coverage to reach before a tally step is executed
-real(DP),parameter :: SUBPIXELCOVERAGE = 0.05_DP ! The total area coverage to reach before a subpixel evaluate step is executed: 0.05_DP
+real(DP),parameter :: SUBPIXELCOVERAGE = 0.025_DP ! The total area coverage to reach before a subpixel evaluate step is executed: 0.05_DP
 real(DP),parameter :: COOKIESIZE = 3.0_DP      ! Relative size of old crater to new crater that cookie cutting is applied
                                                ! Only craters smaller than COOKIESIZE times the new crater are cookie cut
 real(DP),parameter :: ALPHA = 0.125_DP
@@ -175,8 +175,9 @@ module module_globals
    real(DP) :: regcoh  ! target surface regolith layer cohesion
 
    ! Crater diffusion input parameters
-   real(DP) :: soften_factor ! K_{c,0} from Minton et al. (2017)
-   real(DP) :: soften_slope  ! power law index of diffusion model K_{c,i}
+   real(DP) :: soften_factor ! Kbar,0 from Minton et al. (2017)
+   real(DP) :: soften_slope  ! power law index of diffusion model Kbar
+   real(DP) :: soften_size  ! size of diffusion area in crater radii
 
    ! Ejecta softening variables
    logical           :: dosoftening  ! Set T to use the extra crater softening model