Added porosity plot function