From 19b366ca8a5998b2b0798265ad31b2d4d0734717 Mon Sep 17 00:00:00 2001
From: daminton <daminton@purdue.edu>
Date: Thu, 9 Feb 2017 18:58:19 +0000
Subject: [PATCH] Added new user input that allows you to truncate the ejecta
 blanket at a specific radius

---
 src/io/io_input.f90 | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/src/io/io_input.f90 b/src/io/io_input.f90
index a2c5e1db..6d5b5cd5 100644
--- a/src/io/io_input.f90
+++ b/src/io/io_input.f90
@@ -86,7 +86,7 @@ subroutine io_input(infile,user)
    user%doporosity = .false.
    user%soften_factor = 0.43_DP
    user%soften_slope = 1.8_DP
-   user%discontinuous = .true.
+   user%ejecta_truncation = 10.0_DP
    
    open(unit=LUN,file=infile,status="old",iostat=ierr)
    if (ierr /= 0) then
@@ -303,11 +303,11 @@ subroutine io_input(infile,user)
             call io_get_token(line, ilength, ifirst, ilast, ierr)
             token = line(ifirst:ilast)
             read(token, *) user%testtally
-         case ("DISCONTINUOUS") ! Do discontinuous ejecta blanket
+         case ("EJECTA_TRUNCATION") ! Size of ejecta blankets to model in terms of crater diameter 
             ifirst = ilast + 1
             call io_get_token(line, ilength, ifirst, ilast, ierr)
             token = line(ifirst:ilast)
-            read(token, *) user%discontinuous
+            read(token, *) user%ejecta_truncation
          ! Porosity model
          case ("POROSITYFLG")
             ifirst = ilast + 1