From cc75cf41aea4bf1f0748bb17d06128b787548d96 Mon Sep 17 00:00:00 2001
From: mhirabay <mhirabay@purdue.edu>
Date: Thu, 29 Dec 2016 14:11:40 +0000
Subject: [PATCH] Call init_porosity_stack here

---
 src/init/init_surf.f90 | 5 ++++-
 1 file changed, 4 insertions(+), 1 deletion(-)

diff --git a/src/init/init_surf.f90 b/src/init/init_surf.f90
index 8a0609f7..838f799f 100644
--- a/src/init/init_surf.f90
+++ b/src/init/init_surf.f90
@@ -38,6 +38,9 @@ subroutine init_surf(user,surf)
    !if (user%docrustal_thinning) surf%mantle = 0._DP
 
    if (user%doregotrack) call init_regolith_stack(user,surf)
+   
+   ! If doporosity, call init_porosity_stack to define the porolayer linked list. 
+ 	if (user%doporosity)  call init_porosity_stack(user,surf)
 
    return
-   end subroutine init_surf
+end subroutine init_surf