diff --git a/example/porous_media/fromPoro2Perm.py b/example/porous_media/fromPoro2Perm.py
deleted file mode 100644
index 2688b2e9..00000000
--- a/example/porous_media/fromPoro2Perm.py
+++ /dev/null
@@ -1,74 +0,0 @@
-import numpy as np
-import matplotlib.pyplot as plt
-import cv2
-from .utils import make_plot, save_img
-
-
-# %% 
-path= '/home/grzegorz/GITHUB/LBM/TCLB/example/porous_media/'
-image = cv2.imread(path+'porous_media_contrast_1024x256.png') 
-# image = cv2.imread('japan_1024x640.png') 
-
-gray_image = cv2.cvtColor(image, cv2.COLOR_BGR2GRAY) # convert to grayscale
-# gray_image = image
-# np_gray_image = np.array(gray_image)
-# Convert to float
-float_image = gray_image.astype(np.float32)
-
-# Rescale to 0-1 range
-normalized_image = float_image / 255.0
-
-shape = gray_image.shape
-print(f"input shape: {shape}")
-
-
-# %% 
-# Parameters for the Kozeny-Carman equation
-K = 1.  # Kozeny constant, arbitrary choice for demonstration
-S = 1.  # Specific surface area, arbitrary choice for demonstration
-
-# Epsilon values to prevent division by zero in the denominator
-
-
-# Porosity range from 0 to 1 (not including 1 to avoid division by zero in the original equation)
-phi = normalized_image #np.linspace(0.0, 0.9, 200)
-
-# Prepare the plot
-
-
-# Calculate and plot permeability for each epsilon value
-
-
-
-
-
-    
-    
-# Define the smooth limiter function based on the given C code
-
-
-
-# Adjusted Kozeny-Carman equation to include epsilon in the denominator
-def Kozeny_Carman(phi,K,S, epsilon = 1E-3):
-    return (phi**3) / (K * ((1 - phi)**2 + epsilon)) * (1/S**2)
-
-rephi = rescale(phi, 0.01, 0.99)
-k = Kozeny_Carman(rephi,K,S)
-
-make_plot(rephi, 'Input Porosity', 'gist_gray')
-make_plot(k, 'Permability', 'gist_gray')
-
-# ks = rescale(k, 0.01, 1000) 
-# make_plot(ks, 'RescaledPermability')
-# h, w, c = k.shape
-h, w = k.shape
-save_img(k, h, w, 'gist_gray')
-
-# %%
-
-# x = np.linspace(0,10,1000)
-# y = smooth_limiter(x, max=10, steepness= 1)
-
-# plt.plot(x,y)
-# plt.plot(x,x)
-# plt.grid()