From 1ee5804ac5e4d8fbbf677645b1cc0c2cff11033d Mon Sep 17 00:00:00 2001
From: mhirabay <mhirabay@purdue.edu>
Date: Thu, 29 Dec 2016 14:10:47 +0000
Subject: [PATCH] Initialize the porosity linked list

---
 src/init/init_porosity_stack.f90 | 46 ++++++++++++++++++++++++++++++++
 1 file changed, 46 insertions(+)
 create mode 100644 src/init/init_porosity_stack.f90

diff --git a/src/init/init_porosity_stack.f90 b/src/init/init_porosity_stack.f90
new file mode 100644
index 00000000..702e76c3
--- /dev/null
+++ b/src/init/init_porosity_stack.f90
@@ -0,0 +1,46 @@
+!**********************************************************************************************************************************
+!
+!  Unit Name   : init_porosity_stack
+!  Unit Type   : subroutine
+!  Project     : CTEM
+!  Language    : Fortran 2003
+!
+!  Description : Initializes the porosity stack to null
+!  
+!
+!  Input
+!    Arguments : porolayer
+!
+!  Output
+!    Arguments : porolayer
+!           
+! 
+!  Notes       :  
+!
+!**********************************************************************************************************************************
+subroutine init_porosity_stack(user,surf)
+   use module_globals
+   use module_util
+   use module_init, EXCEPT_THIS_ONE => init_porosity_stack
+   implicit none
+
+   ! Arguments
+   type(usertype),intent(in) :: user
+   type(surftype),dimension(:,:),intent(inout) :: surf
+   type(regodatatype) :: bedrock    
+   integer(I4B) :: xp, yp
+
+   ! Internal variables
+   logical :: initstat
+
+   do yp = 1, user%gridsize
+      do xp = 1, user%gridsize
+
+      	call util_init_list(surf(xp,yp)%porolayer, initstat)
+         
+      end do
+   end do
+
+   return
+end subroutine init_porosity_stack
+