diff --git a/DiffractionGratings/grating3D_postplot_Mmatrix.py b/DiffractionGratings/grating3D_postplot_Mmatrix.py
index f209b638a953679cc6911fe862db91e61d1e190a..ed8739d8823879993ec78cbd62509102b2623055 100644
--- a/DiffractionGratings/grating3D_postplot_Mmatrix.py
+++ b/DiffractionGratings/grating3D_postplot_Mmatrix.py
@@ -86,19 +86,19 @@ TOT2_sin = T2_sin+R2_sin+Qscat_sin
 
 M[1-1,1-1,:,:] = 0.5*(np.abs(rpp)**2+np.abs(rsp)**2+np.abs(rps)**2+np.abs(rss)**2)
 M[1-1,2-1,:,:] = 0.5*(np.abs(rpp)**2+np.abs(rsp)**2-np.abs(rps)**2-np.abs(rss)**2)
-M[1-1,3-1,:,:] = np.real(rpp*np.conj(rps)+rsp*np.conj(rss)) #*-1
+M[1-1,3-1,:,:] = np.real(rpp*np.conj(rps)+rsp*np.conj(rss))
 M[1-1,4-1,:,:] = np.imag(rpp*np.conj(rps)+rsp*np.conj(rss))
 M[2-1,1-1,:,:] = 0.5*(np.abs(rpp)**2-np.abs(rsp)**2+np.abs(rps)**2-np.abs(rss)**2)
 M[2-1,2-1,:,:] = 0.5*(np.abs(rpp)**2-np.abs(rsp)**2-np.abs(rps)**2+np.abs(rss)**2)
-M[2-1,3-1,:,:] = np.real(rpp*np.conj(rps)-rsp*np.conj(rss)) #*-1
+M[2-1,3-1,:,:] = np.real(rpp*np.conj(rps)-rsp*np.conj(rss))
 M[2-1,4-1,:,:] = np.imag(rpp*np.conj(rps)-rsp*np.conj(rss))
-M[3-1,1-1,:,:] = np.real(rpp*np.conj(rsp)+rps*np.conj(rss)) #*-1
-M[3-1,2-1,:,:] = np.real(rpp*np.conj(rsp)-rps*np.conj(rss)) #*-1
+M[3-1,1-1,:,:] = np.real(rpp*np.conj(rsp)+rps*np.conj(rss))
+M[3-1,2-1,:,:] = np.real(rpp*np.conj(rsp)-rps*np.conj(rss))
 M[3-1,3-1,:,:] = np.real(rpp*np.conj(rss)+rps*np.conj(rsp))
-M[3-1,4-1,:,:] = np.imag(rpp*np.conj(rss)-rps*np.conj(rsp)) #*-1
+M[3-1,4-1,:,:] = np.imag(rpp*np.conj(rss)-rps*np.conj(rsp))
 M[4-1,1-1,:,:] =-np.imag(rpp*np.conj(rsp)+rps*np.conj(rss))
 M[4-1,2-1,:,:] =-np.imag(rpp*np.conj(rsp)-rps*np.conj(rss))
-M[4-1,3-1,:,:] =-np.imag(rpp*np.conj(rss)+rps*np.conj(rsp)) #*-1
+M[4-1,3-1,:,:] =-np.imag(rpp*np.conj(rss)+rps*np.conj(rsp))
 M[4-1,4-1,:,:] = np.real(rpp*np.conj(rss)-rps*np.conj(rsp))
 
 
@@ -147,7 +147,6 @@ for i in range(4):
             norm = mpl.colors.Normalize(vmin=-1,vmax=1)
             sm = mpl.cm.ScalarMappable(cmap=plt.cm.bwr, norm=norm)
             sm.set_array([])
-            # plt.colorbar(sm, ticks=np.linspace(-1,1,3))
         ax.set_rorigin(which_orig)
 plt.subplots_adjust(top=0.92, bottom=0.08, left=0.10, right=0.95, hspace=0.25,wspace=0.35)
 plt.savefig('Mmatrix.jpg')