From a7a651d4dd602efcb37b05c04664a9f46f03cfae Mon Sep 17 00:00:00 2001
From: daminton <daminton@purdue.edu>
Date: Tue, 29 Nov 2016 15:48:29 +0000
Subject: [PATCH] added new function to compute critical slope angles inside
 new craters

---
 src/crater/crater_critical_slope.f90 | 36 ++++++++++++++++++++++++++++
 1 file changed, 36 insertions(+)
 create mode 100644 src/crater/crater_critical_slope.f90

diff --git a/src/crater/crater_critical_slope.f90 b/src/crater/crater_critical_slope.f90
new file mode 100644
index 00000000..9d94cde1
--- /dev/null
+++ b/src/crater/crater_critical_slope.f90
@@ -0,0 +1,36 @@
+!  Project     : CTEM
+!  Language    : Fortran 2003
+!
+!  Description : Forms the maximum slope for the initial crater morphology
+!  
+!
+!  Input
+!    Arguments : 
+!
+!  Output
+!    Arguments :
+!           
+! 
+!  Notes       :  
+!
+!**********************************************************************************************************************************
+
+
+function crater_critical_slope(user,crater,iradsq) result(critical)
+   use module_globals
+   use module_util
+   use module_crater, EXCEPT_THIS_ONE => crater_form_exterior_func
+   implicit none
+
+  ! Arguments
+   type(usertype),intent(in) :: user
+   type(cratertype),intent(in) :: crater
+   integer(I4B),intent(in) :: iradsq
+   real(DP) :: critical
+
+   ! Internal variables
+
+   critical = CRITSLP
+
+end function crater_critical_slope
+