From dc5f52c40f6cb17fbefeb2f8aae3d437df885d95 Mon Sep 17 00:00:00 2001
From: mhirabay <mhirabay@purdue.edu>
Date: Thu, 29 Dec 2016 14:15:55 +0000
Subject: [PATCH] Added the porosity layer.

---
 src/util/util_init_list.f90 | 13 +++++++------
 1 file changed, 7 insertions(+), 6 deletions(-)

diff --git a/src/util/util_init_list.f90 b/src/util/util_init_list.f90
index 635f1b7c..2159497d 100644
--- a/src/util/util_init_list.f90
+++ b/src/util/util_init_list.f90
@@ -54,7 +54,7 @@
 !                layer. 
 !
 !**********************************************************************************************************************************
-subroutine util_init_list(regolayer,initstat)
+subroutine util_init_list(regolayer, initstat)
    use module_globals
    use module_util, EXCEPT_THIS_ONE => util_init_list
    implicit none
@@ -74,9 +74,11 @@ subroutine util_init_list(regolayer,initstat)
          initstat = .true.
          nullify(regolayer%next)
          regolayer%regodata%thickness = VBIG ! This generates a buffer layer that the model should never reach if the run is structured properly
-         regolayer%regodata%comp = 0.0_DP
-         regolayer%regodata%meltfrac = 0.0_DP
-         regolayer%regodata%porosity = 0.0_DP
+         regolayer%regodata%comp      = 0.0_DP
+         regolayer%regodata%meltfrac  = 0.0_DP
+         regolayer%regodata%porosity  = 0.07_DP      ! This is also the initial bedrock porosity, Should be in the input file. 
+         regolayer%regodata%depth     = - 30000.0_DP ! The depth at the boundary of the crust. 
+                                                     ! This must be defined in the input file. But not in this version. 
       else
          write(*,*) 'util_init_list: Initialization failed. Exhausted memory.'
       end if
@@ -85,5 +87,4 @@ subroutine util_init_list(regolayer,initstat)
    end if
 
    return
-end subroutine util_init_list
-
+end subroutine util_init_list
\ No newline at end of file