diff --git a/DiffractionGratings/grating2D.pro b/DiffractionGratings/grating2D.pro
index 84d0a5fdb2021ed47436d2aa83edecf85de27467..4531aeadbdbce3644cd427ac7cec69b583bf96fb 100644
--- a/DiffractionGratings/grating2D.pro
+++ b/DiffractionGratings/grating2D.pro
@@ -486,9 +486,9 @@ PostOperation {
         Print[ Q_layer_dep[layer_dep] , OnGlobal, Format FrequencyTable, File > StrCat[myDir, "absorption-Q_layer_dep.txt"]];
         Print[ Q_layer_cov[layer_cov] , OnGlobal, Format FrequencyTable, File > StrCat[myDir, "absorption-Q_layer_cov.txt"]];
         Print[ Hz_tot    , OnElementsOf Plot_domain, File StrCat[myDir, Sprintf("Hz_tot_lambda%.2fnm_1.pos", lambda0/nm)] , Name Sprintf("Hz_tot_%.2fnm.pos", lambda0/nm)];
-        // debug
-        Print[ u1    , OnElementsOf Plot_domain, File StrCat[myDir, Sprintf("u1%.2fnm_1.pos", lambda0/nm)] , Name Sprintf("u1%.2fnm.pos", lambda0/nm)];
-        Print[ E1    , OnElementsOf Plot_domain, File StrCat[myDir, Sprintf("E1%.2fnm_1.pos", lambda0/nm)] , Name Sprintf("E1%.2fnm.pos", lambda0/nm)];
+        // // debug
+        // Print[ u1    , OnElementsOf Plot_domain, File StrCat[myDir, Sprintf("u1%.2fnm_1.pos", lambda0/nm)] , Name Sprintf("u1%.2fnm.pos", lambda0/nm)];
+        // Print[ E1    , OnElementsOf Plot_domain, File StrCat[myDir, Sprintf("E1%.2fnm_1.pos", lambda0/nm)] , Name Sprintf("E1%.2fnm.pos", lambda0/nm)];
         Echo[ Str["For i In {0:2}",
                   "  View[i].LineWidth = 4;View[i].ColormapNumber = 15;View[0].AxesFormatX = '%.1f';",
                   "EndFor"], File StrCat[myDir,"tmp0.geo"]] ;
diff --git a/ElectromagneticScattering/scattering_data.geo b/ElectromagneticScattering/scattering_data.geo
index 66db18ef76b2d7c5bcba9bc6cfb64ea3a79f04ff..8cea03ed354067b62804b7b297c237710caf08fd 100644
--- a/ElectromagneticScattering/scattering_data.geo
+++ b/ElectromagneticScattering/scattering_data.geo
@@ -10,7 +10,7 @@ cel      = 1.0/(Sqrt[epsilon0 * mu0]);
 deg2rad  = Pi/180.;
 
 pp0        = "1Geometry/0";
-pp1        = "2Study Type/0";
+pp1        = "0Study Type/0";
 pp2        = "3Electromagnetic parameters/0";
 pp3        = "4Mesh size and PMLs parameters/0";
 pp4        = "5Postpro options/0";
@@ -34,7 +34,7 @@ RES_QNM   = 3;
 
 DefineConstant[
   flag_shape = {ELL , Name StrCat[pp0, "0Scatterer shape"], 
-    Choices {ELL="ellispoid",PARALL="parallelepiped",CYL="cylinder",CONE="cone",TOR="split torus"},Closed 0,GmshOption "Reset", Autocheck 0}];
+    Choices {ELL="ellispoid",PARALL="parallelepiped",CYL="cylinder",CONE="cone",TOR="split torus"},Closed 0}];
 
 If (flag_shape==ELL)
   DefineConstant[
@@ -103,7 +103,7 @@ DefineConstant[
         Choices {RES_PW="Plane Wave Response",
                  RES_TMAT="T-Matrix", 
                  RES_GREEN="Green's Tensor",
-                 RES_QNM="Quasi-Normal Modes"},Closed 0}];
+                 RES_QNM="Quasi-Normal Modes"},Closed 0, GmshOption "Reset", Autocheck 0}];
 
 DefineConstant[
   epsr_In_re  = { 9., Name StrCat[pp2  , "0scatterer permittivity (real) []"] , Highlight Str[colorppOK] , Closed 0},