diff --git a/DiffractionGratings/grating2D_conical.pro b/DiffractionGratings/grating2D_conical.pro
index 51d0c97a6f38794c388e4ba5c045400b3f723cfe..d36f95e8104003732f425d135eee99551ae54203 100644
--- a/DiffractionGratings/grating2D_conical.pro
+++ b/DiffractionGratings/grating2D_conical.pro
@@ -18,8 +18,8 @@
 myDir = "res2D_conical/";
 DefineConstant[
   lambda0   = {lambda_min , Min lambda_min, Max lambda_max, Step (lambda_max-lambda_min)/(nb_lambdas-1), Name StrCat[pp2, "1wavelength [nm]"] , Loop 1, Highlight Str[colorpp2],Graph "200000200020", ServerAction "Reset GetDP/T0, GetDP/R0, GetDP/Lambda_step, GetDP/Omegaal absorption"},
-  phi_deg = {  45 , Name StrCat[pp2, "3incident plane wave angle (phi) [deg]"] , Highlight Str[colorpp2], Closed close_menu},
-  psi_deg = {  45 , Name StrCat[pp2, "4incident plane wave angle (psi) [deg]"] , Highlight Str[colorpp2], Closed close_menu}
+  phideg = {  45 , Name StrCat[pp2, "3incident plane wave angle (phi) [deg]"] , Highlight Str[colorpp2], Closed close_menu},
+  psideg = {  45 , Name StrCat[pp2, "4incident plane wave angle (psi) [deg]"] , Highlight Str[colorpp2], Closed close_menu}
 ];
 lambda0    = lambda0    * nm;
 lambda_min = lambda_min * nm;
@@ -140,9 +140,9 @@ Function{
   k1norm[] = k0*n1[];
   k2norm[] = k0*n2[];
 
-  phi0    =  phi_deg   * deg2rad;
+  phi0    =  phideg   * deg2rad;
   theta0  =  theta_deg * deg2rad;
-  psi0    =  psi_deg   * deg2rad;
+  psi0    =  psideg   * deg2rad;
   // Ae      =  1/Sqrt[ep0/mu0];
   Ae      = 1;