diff --git a/src/main/CTEM.f90 b/src/main/CTEM.f90
index d08e77b9..4f64f5c4 100644
--- a/src/main/CTEM.f90
+++ b/src/main/CTEM.f90
@@ -155,6 +155,15 @@ program CTEM
    end do
 end if
 
+! If doporosity is true, then destroy the linked list for porosity
+if (user%doporosity) then
+   do yp = 1, user%gridsize
+      do xp = 1, user%gridsize
+         call util_destroy_list(surf(xp,yp)%porolayer)
+      end do
+   end do
+end if
+
 ! Deallocate all the allocatables
 deallocate(seedarr)
 deallocate(surf,prod,vdist,pdist,crtscl,truedist,truelist,obsdist,obslist,nflux,production_list)