diff --git a/Plot-Scripts/PerforatedBilayerPlot.py b/Plot-Scripts/PerforatedBilayerPlot.py
new file mode 100644
index 0000000000000000000000000000000000000000..99cf4f4f97a24eece842e28a8caeff980df77608
--- /dev/null
+++ b/Plot-Scripts/PerforatedBilayerPlot.py
@@ -0,0 +1,205 @@
+import numpy as np
+import matplotlib.pyplot as plt
+import math
+import os
+import subprocess
+import fileinput
+import re
+import sys
+import matplotlib as mpl
+from mpl_toolkits.mplot3d import Axes3D
+import matplotlib.cm as cm
+import matplotlib.ticker as ticker
+from matplotlib.ticker import MultipleLocator,FormatStrFormatter,MaxNLocator
+import seaborn as sns
+import matplotlib.colors as mcolors
+import codecs
+import json
+# ----------------------------------------------------
+# --- Define Plot style:
+plt.style.use("seaborn")
+
+mpl.rcParams['text.usetex'] = True
+mpl.rcParams["font.family"] = "serif"
+mpl.rcParams["font.size"] = "8"
+mpl.rcParams['xtick.bottom'] = True
+mpl.rcParams['xtick.major.size'] = 3
+mpl.rcParams['xtick.minor.size'] = 1.5
+mpl.rcParams['xtick.major.width'] = 0.75
+mpl.rcParams['xtick.labelsize'] = 8
+mpl.rcParams['xtick.major.pad'] = 1
+
+mpl.rcParams['ytick.left'] = True
+mpl.rcParams['ytick.major.size'] = 3
+mpl.rcParams['ytick.minor.size'] = 1.5
+mpl.rcParams['ytick.major.width'] = 0.75
+mpl.rcParams['ytick.labelsize'] = 8
+mpl.rcParams['ytick.major.pad'] = 1
+
+mpl.rcParams['axes.titlesize'] = 8
+mpl.rcParams['axes.titlepad'] = 1
+mpl.rcParams['axes.labelsize'] = 8
+
+#Adjust Legend:
+mpl.rcParams['legend.frameon'] = True # Use frame for legend
+# mpl.rcParams['legend.framealpha'] = 0.5
+mpl.rcParams['legend.fontsize'] = 8 # fontsize of legend
+
+
+#Adjust grid:
+mpl.rcParams.update({"axes.grid" : True}) # Add grid
+mpl.rcParams['axes.labelpad'] = 3
+mpl.rcParams['grid.linewidth'] = 0.25
+mpl.rcParams['grid.alpha'] = 0.9 # 0.75
+mpl.rcParams['grid.linestyle'] = '-'
+mpl.rcParams['grid.color'] = 'gray'#'black'
+mpl.rcParams['text.latex.preamble'] = r'\usepackage{amsfonts}' # Makes Use of \mathbb possible.
+# ----------------------------------------------------------------------------------------
+textwidth = 6.26894 #textwidth in inch
+# width = textwidth * 0.5
+width = textwidth
+height = width / 1.618 # The golden ratio.
+
+
+dataset_numbers = [0, 1, 2, 3, 4 ,5]
+
+kappa_sim_upper = []
+
+#--- Choose wether to perforate upper (passive) or lower (active) layer
+# perforatedLayer = 'upper'
+perforatedLayer = 'lower'
+
+for dataset_number in dataset_numbers:
+
+ fig, ax = plt.subplots(figsize=(width,height))
+
+
+ Path_KappaUpper = './experiment/perforated-bilayer/results_' + perforatedLayer + '_' + str(dataset_number) + '/kappa_simulation.txt'
+ # Path_KappaLower = './experiment/perforated-bilayer/results_lower_' + str(dataset_number) + '/kappa_simulation.txt'
+
+ tmp = open(Path_KappaUpper).read().split(",")
+ # kappa_sim_lower = open(Path_KappaLower).read().split(",")
+ curvature = [float(i) for i in tmp]
+ # kappa_sim_upper[dataset_number] = np.array(curvature)
+ kappa_sim_upper.append(np.array(curvature))
+ # kappa_sim_lower = [float(i) for i in kappa_sim_lower]
+
+
+
+
+input = [0 , 0.05, 0.10, 0.15, 0.20, 0.25, 0.30] # Design parameter (volume ratio)
+
+
+
+# --------------- Plot Lines + Scatter -----------------------
+line_1 = ax.plot(input, kappa_sim_upper[0], # data
+ # color='forestgreen', # linecolor
+ marker='D', # each marker will be rendered as a circle
+ markersize=5, # marker size
+ # markerfacecolor='darkorange', # marker facecolor
+ markeredgecolor='black', # marker edgecolor
+ markeredgewidth=0.75, # marker edge width
+ # linestyle='dashdot', # line style will be dash line
+ linewidth=1.5, # line width
+ zorder=3,
+ label = r"$r = 0.12$ ")
+
+line_2 = ax.plot(input, kappa_sim_upper[1], # data
+ # color='forestgreen', # linecolor
+ marker='D', # each marker will be rendered as a circle
+ markersize=5, # marker size
+ # markerfacecolor='darkorange', # marker facecolor
+ markeredgecolor='black', # marker edgecolor
+ markeredgewidth=0.75, # marker edge width
+ # linestyle='dashdot', # line style will be dash line
+ linewidth=1.5, # line width
+ zorder=3,
+ label = r"$r = 0.17$ ")
+
+line_3 = ax.plot(input, kappa_sim_upper[2], # data
+ # color='forestgreen', # linecolor
+ marker='D', # each marker will be rendered as a circle
+ markersize=5, # marker size
+ # markerfacecolor='darkorange', # marker facecolor
+ markeredgecolor='black', # marker edgecolor
+ markeredgewidth=0.75, # marker edge width
+ # linestyle='dashdot', # line style will be dash line
+ linewidth=1.5, # line width
+ zorder=3,
+ label = r"$r = 0.22$ ")
+
+line_4 = ax.plot(input, kappa_sim_upper[3], # data
+ # color='forestgreen', # linecolor
+ marker='D', # each marker will be rendered as a circle
+ markersize=5, # marker size
+ # markerfacecolor='darkorange', # marker facecolor
+ markeredgecolor='black', # marker edgecolor
+ markeredgewidth=0.75, # marker edge width
+ # linestyle='dashdot', # line style will be dash line
+ linewidth=1.5, # line width
+ zorder=3,
+ label = r"$r = 0.34$ ")
+
+line_5 = ax.plot(input, kappa_sim_upper[4], # data
+ # color='orangered', # linecolor
+ marker='o', # each marker will be rendered as a circle
+ markersize=5, # marker size
+ # markerfacecolor='cornflowerblue', # marker facecolor
+ markeredgecolor='black', # marker edgecolor
+ markeredgewidth=0.75, # marker edge width
+ # linestyle='--', # line style will be dash line
+ linewidth=1.5, # line width
+ zorder=3,
+ alpha=0.8, # Change opacity
+ label = r"$r=0.43$ ")
+
+line_6 = ax.plot(input, kappa_sim_upper[5], # data
+ # color='forestgreen', # linecolor
+ marker='D', # each marker will be rendered as a circle
+ markersize=5, # marker size
+ # markerfacecolor='darkorange', # marker facecolor
+ markeredgecolor='black', # marker edgecolor
+ markeredgewidth=0.75, # marker edge width
+ # linestyle='dashdot', # line style will be dash line
+ linewidth=1.5, # line width
+ zorder=3,
+ label = r"$r = 0.49$ ")
+
+# --------------- Set Axes -----------------------
+# Plot - Title
+title = perforatedLayer + ' layer perforation'
+ax.set_title(title)
+# match dataset_number:
+# case 0:
+# ax.set_title(r"ratio $r = 0.12$")
+# case 1:
+# ax.set_title(r"ratio $r = 0.17$")
+# case 2:
+# ax.set_title(r"ratio $r = 0.22$")
+# case 3:
+# ax.set_title(r"ratio $r = 0.34$")
+# case 4:
+# ax.set_title(r"ratio $r = 0.43$")
+# case 5:
+# ax.set_title(r"ratio $r = 0.49$")
+
+ax.set_xlabel(r"Volume ratio ", labelpad=4)
+ax.set_ylabel(r"Curvature $\kappa$", labelpad=4)
+plt.tight_layout()
+
+
+# --- Set Legend
+legend = ax.legend()
+
+frame = legend.get_frame()
+frame.set_edgecolor('black')
+
+
+# --- Adjust left/right spacing:
+# plt.subplots_adjust(right=0.81)
+# plt.subplots_adjust(left=0.11)
+
+# ---------- Output Figure as pdf:
+fig.set_size_inches(width, height)
+fig.savefig('PerforatedBilayer_' + perforatedLayer + '.pdf')
+# plt.show()
\ No newline at end of file