diff --git a/DiffractionGratings/grating2D_conical.pro b/DiffractionGratings/grating2D_conical.pro
index d36f95e8104003732f425d135eee99551ae54203..f061331cdd4b7ab11b9ffbb945a3ead40f996e41 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"},
-  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}
+  phideg = {  30 , Name StrCat[pp2, "3incident plane wave angle (phi) [deg]"] , Highlight Str[colorpp2], Closed close_menu},
+  psideg = {  10 , Name StrCat[pp2, "4incident plane wave angle (psi) [deg]"] , Highlight Str[colorpp2], Closed close_menu}
 ];
 lambda0    = lambda0    * nm;
 lambda_min = lambda_min * nm;
diff --git a/DiffractionGratings/grating3D_data_retrieve_2D_lamellar.geo b/DiffractionGratings/grating3D_data_retrieve_2D_lamellar.geo
index ed3c4a729bdb24d5b302ed45622ce8c951700c7b..83f95d704f1c440ed906799040ed370010be504c 100644
--- a/DiffractionGratings/grating3D_data_retrieve_2D_lamellar.geo
+++ b/DiffractionGratings/grating3D_data_retrieve_2D_lamellar.geo
@@ -9,8 +9,8 @@ DefineConstant[
     FLAG_TOTAL    = {0     , Name StrCat[pp1,"/0Formulation"],Choices {0="scattered field",1="total field"}},
     lambda0       = {1000  , Name StrCat[pp1,"/1lambda0 [nm]"]},
     thetadeg      = {30    , Name StrCat[pp1,"/2theta0 [deg]"]},
-    phideg        = {90     , Name StrCat[pp1,"/3phi0 [deg]"]},
-    psideg        = {0     , Name StrCat[pp1,"/4psi0 [deg]"]},
+    phideg        = {30     , Name StrCat[pp1,"/3phi0 [deg]"]},
+    psideg        = {10    , Name StrCat[pp1,"/4psi0 [deg]"]},
     period_x      = {10    , Name StrCat[pp2,"/1X period [nm]"]},
     period_y      = {1000  , Name StrCat[pp2,"/2Y period [nm]"]},
     thick_L_1     = {200    , Name StrCat[pp2,"/3thickness layer 1 [nm] (superstrate)"]},