Call init_porosity_stack here

MintonGroup · Dec 29, 2016 · cc75cf4 · cc75cf4
1 parent 1ee5804
commit cc75cf4
Showing 1 changed file with 4 additions and 1 deletion.
diff --git 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