diff --git a/DiffractionGratings/grating3D.pro b/DiffractionGratings/grating3D.pro
index dd56ea54ac78a458f48aadde72dd00d6cae400b4..1af89f2d01326b65c322bc30ff60d7755c436d65 100644
--- a/DiffractionGratings/grating3D.pro
+++ b/DiffractionGratings/grating3D.pro
@@ -181,35 +181,6 @@ Function{
   H1d[Omega_super] = Hr[];
   H1d[Omega_subs]  = Ht[];
 
-  // ////////////////////////:
-  // k1s[]  = (Abs[theta0]>1e-12) ?      k1[]/\zhat[] / Norm[    k1[]/\zhat[]] : Vector[-Sin[phi0], Cos[phi0], 0.] ;
-  // k2s[]  = (Abs[theta0]>1e-12) ?  Re[k2[]]/\zhat[] / Norm[ Re[k2[]]/\zhat[]] : Vector[-Sin[phi0], Cos[phi0], 0.] ;
-  // k1rs[] = (Abs[theta0]>1e-12) ?     k1r[]/\zhat[] / Norm[    k1r[]/\zhat[]] : Vector[-Sin[phi0], Cos[phi0], 0.] ;
-  // k1p[]  = (Abs[theta0]>1e-12) ?      k1[]/\k1s[]  / Norm[     k1[]/\k1s[]]  : Vector[ Cos[phi0], Sin[phi0], 0.] ;
-  // k2p[]  = (Abs[theta0]>1e-12) ?  Re[k2[]]/\k2s[]  / Norm[ Re[k2[]]/\k2s[]]  : Vector[ Cos[phi0], Sin[phi0], 0.] ;
-  // k1rp[] = (Abs[theta0]>1e-12) ?     k1r[]/\k1rs[] / Norm[    k1r[]/\k1rs[]] : Vector[-Cos[phi0],-Sin[phi0], 0.] ;
-  // rsa[] = (CompZ[k1[]]-CompZ[k2[]])/(CompZ[k1[]]+CompZ[k2[]]);
-  // tsa[] =           2.*CompZ[k1[]] /(CompZ[k1[]]+CompZ[k2[]]);
-  // rpa[] = (CompZ[k1[]]*epsr2[]-CompZ[k2[]]*epsr1[])/(CompZ[k1[]]*epsr2[]+CompZ[k2[]]*epsr1[]);
-  // tpa[] =                  (2.*CompZ[k1[]]*epsr2[])/(CompZ[k1[]]*epsr2[]+CompZ[k2[]]*epsr1[]) * Sqrt[epsr1[]/epsr2[]];
-  // AmplEi[] =       Cos[psi0]* k1p[] +       Sin[psi0]* k1s[];
-  // AmplEr[] = rpa[]*Cos[psi0]*k1rp[] + rsa[]*Sin[psi0]*k1rs[];
-  // AmplEt[] = tpa[]*Cos[psi0]* k2p[] + tsa[]*Sin[psi0]* k2s[];
-  // Eia[]     = Ae*AmplEi[] * Exp[I[]*k1[]*XYZ[]];
-  // Eta[]     = Ae*AmplEt[] * Exp[I[]*k2[]*XYZ[]];
-  // Era[]     = Ae*AmplEr[] * Exp[I[]*k1r[]*XYZ[]];
-  // Hia[]     =  1/(om0*mu0*mur[]) * k1[]  /\ Ei[];
-  // Hra[]     =  1/(om0*mu0*mur[]) * k1r[] /\ Er[];
-  // Hta[]     =  1/(om0*mu0*mur[]) * k2[]  /\ Et[];
-  // E1a[Omega_super]  = Eia[]+Era[];
-  // E1a[Omega_subs]   = Eta[];
-  // E1da[Omega_super] = Era[];
-  // E1da[Omega_subs]  = Eta[];
-  // H1a[Omega_super]  = Hia[]+Hra[];
-  // H1a[Omega_subs]   = Hta[];
-  // H1da[Omega_super] = Hra[];
-  // H1da[Omega_subs]  = Hta[];
-  // /////////////////////
   source[] = (om0/cel)^2*(epsr[]-epsr_annex[])*E1[];
   
   // Bloch phase shifts