Cacher is the code snippet organizer for pro developers

import seaborn as sns

colors = ['blue', 'red']  
sns.set_palette(sns.xkcd_palette(colors))

fig, axarr = plt.subplots(16,1,sharex=True, sharey=True, figsize=(16,32))
fig.suptitle('TITEL', fontsize=18)

fig.add_subplot(111, frameon=False)
plt.tick_params(labelcolor='none', top=False, bottom=False, left=False, right=False)

for k in range(0,16):             
    axarr[k].plot(NH3[k*250:(k+1)*250], label='Messung')     
    axarr[k].plot(yreg[k*250:(k+1)*250], label='Simulation') 
    axarr[k].grid()
    axarr[k].legend()

plt.xlabel('Zeit [min]')
plt.ylabel('$NH_3 [mg/Nm^3]$')

fig.subplots_adjust(top=.95)
plt.show()