diff --git a/ElectromagneticScattering/scattererTmatrix.pro b/ElectromagneticScattering/scattererTmatrix.pro
index 3e62be26e38ebfd3abe401a19bf6fc824de60373..8d570c170109c4183fd7f9e86b629a80ea0318f8 100644
--- a/ElectromagneticScattering/scattererTmatrix.pro
+++ b/ElectromagneticScattering/scattererTmatrix.pro
@@ -199,6 +199,7 @@ PostProcessing {
   For pe In {1:p_max}
     { Name VPWM_postpro~{pe}; NameOfFormulation VPWMN_helmholtz_vector; NameOfSystem T;
       Quantity {
+        // { Name epsr_xx; Value { Local { [ CompXX[epsilonr[]] ]; In All_domains; Jacobian JVol; } } }
         { Name E_scat           ; Value { Local { [{u}]; In All_domains; Jacobian JVol; } } }
         { Name H_scat           ; Value { Local { [siwt*I[]/(mur[]*mu0*omega0)*{Curl u}]; In All_domains; Jacobian JVol; } } }
         { Name Mnm_source~{pe}  ; Value { Local { [   Mnm_source~{pe}[] ]; In All_domains; Jacobian JVol; } } }
@@ -227,6 +228,7 @@ PostOperation {
   For pe In {1:p_max}
     {Name VPWM_postop~{pe}; NameOfPostProcessing VPWM_postpro~{pe} ;
       Operation {
+        // Print [ epsr_xx , OnElementsOf All_domains, File "epsr_xx.pos"];
         If (flag_plotcuts==1)
           Print [ E_scat , OnGrid
                 {(r_pml_in-1*nm)*Sin[$B]*Cos[$C], (r_pml_in-1*nm)*Sin[$B]*Sin[$C], (r_pml_in-1*nm)*Cos[$B]} 
diff --git a/ElectromagneticScattering/scattererTmatrix_data.geo b/ElectromagneticScattering/scattererTmatrix_data.geo
index 159bb893c68661fb0f99e9738fd936cd1998cc9d..b6013410f73eae4b649512c2549a77552e976d9d 100644
--- a/ElectromagneticScattering/scattererTmatrix_data.geo
+++ b/ElectromagneticScattering/scattererTmatrix_data.geo
@@ -132,7 +132,7 @@ DefineConstant[
   space2pml    = {lambda_bg/2 , Name StrCat[pp3  , "0space around scatterer [nm]"] , Highlight Str[colorppNO]    , Closed 1,Min (lambda_bg/10.), Max (3*lambda_bg)},
   pml_size     = {lambda_bg , Name StrCat[pp3  , "1PML thickness [nm]"] , Highlight Str[colorppWA]  , Closed 0, Min (lambda_bg/10), Max (3*lambda_bg)},
   paramaille   = {5. , Name StrCat[pp3  , "2mesh size"], Highlight Str[colorppWA]  , Closed 0},
-  refine_scat  = {1. , Name StrCat[pp3  , "3scatterer mesh refinement"], Highlight Str[colorppWA]  , Closed 0}
+  refine_scat  = {2. , Name StrCat[pp3  , "3scatterer mesh refinement"], Highlight Str[colorppWA]  , Closed 0}
   is_FEM_o2    = {1 , Name StrCat[pp3   , "4Interpolation order "] , Choices {0="order 1",1="order 2"}, Closed 0}
 ];