diff --git a/src/io/io_read_surf.f90 b/src/io/io_read_surf.f90
index 85fb3e32..e47da102 100644
--- a/src/io/io_read_surf.f90
+++ b/src/io/io_read_surf.f90
@@ -79,6 +79,9 @@ subroutine io_read_surf(user,surf)
    close(LUN)
 
    if (user%doregotrack) call io_read_regotrack(user,surf)
+   
+	! if doporosity, call io_read_porotrack to define the porosity linked list from the read files. 
+   if (user%doporosity)  call io_read_porotrack(user,surf)
 
    !if (user%docrustal_thinning) then
    !   recsize=sizeof(itmp)*user%gridsize*user%gridsize