diff --git a/NonLinearEVP/NonLinearEVP_data.geo b/NonLinearEVP/NonLinearEVP_data.geo
index 2da8246aa6e202b0aabef864f19e1ecc697080b5..e63b9c56f361c5df87ef8a9877239a1ef3561c21 100644
--- a/NonLinearEVP/NonLinearEVP_data.geo
+++ b/NonLinearEVP/NonLinearEVP_data.geo
@@ -3,6 +3,8 @@
 ////   _data.geo               ////
 ///////////////////////////////////
 
+deg2rad  = Pi/180;
+
 pp0        = "1Geometry/0";
 pp1        = "2Polarization-Bloch/0";
 pp3        = "3Eigenvalue problem parameters/0";
@@ -30,9 +32,9 @@ DefineConstant[
   flag_Hparallel = {1     , Name StrCat[pp1 , "1polarization case"] , Choices {0="E //",1="H //"} },
   kx             = {0.75  , Name StrCat[pp1 , "2kx [Pi\a]"] , Highlight Str[colorpp1]  , Closed close_menu} , 
 
-  eps_oo_1       = {1     , Name StrCat[pp2 , "0eps_oo_1 [ - ]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
-  om_d_1         = {1.1   , Name StrCat[pp2 , "1om_d_1 [2Pic\a]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
-  gam_1          = {0.05  , Name StrCat[pp2 , "2gam_1 [2Pic\a]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
+  eps_oo_1       = {1     , Name StrCat[pp2 , "0epsilon_inf_Drude [ - ]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
+  om_d_1         = {1.1   , Name StrCat[pp2 , "1Omega_Drude [2Pic\a]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
+  gam_1          = {0.05  , Name StrCat[pp2 , "2Gamma_Drude [2Pic\a]"] , Highlight Str[colorpp2]  , Closed close_menu} , 
   
   neig          =  {1     , Name StrCat[pp3 , "0Number of eigenvalues [int]"] , Highlight Str[colorpp3]  , Closed close_menu} ,
   eig_target_re  = {0.0077, Name StrCat[pp3 , "1EV real part target [2Pic\a]"] , Highlight Str[colorpp3]  , Closed close_menu} , 
@@ -51,7 +53,7 @@ DefineConstant[
   flag_outEigvec = {1     , Name StrCat[pp4, "output eigenvector?"], Choices{0,1}}
 ];
 
-
+// normalized units so eigenvalues are close to unity
 cel      = a_lat/(2*Pi);
 epsf     = 8.854187817e-3;
 muf      = 400.*Pi;
@@ -60,9 +62,7 @@ epsilon0 = epsf*nm;
 mu0      = muf*nm;
 norm     = a_lat/(2.*Pi*cel);
 
-Printf("norm %f",norm);    
-
-// Normalization
+// Normalize
 d_sq           = d_sq          * a_lat;
 space2pml      = space2pml     * a_lat;
 pmlsize        = pmlsize       * a_lat;