refactor

ElectromagneticScattering/scattering.pro

@@ -479,17 +479,6 @@ PostProcessing {
     }
   EndIf
   If (flag_study==RES_TMAT)
-    { Name VPWMN_postpro_test; NameOfFormulation VPWMN_helmholtz_vector; NameOfSystem A;
-            Quantity {
-          { Name E_scat            ; Value { Local { [$PE]; In All_domains; Jacobian JVol; } } } 
-          { Name E_scat_sph        ; Value { Local { [Vector[
-                            (X[]*CompX[{u}] + Y[]*CompY[{u}] + Z[]*CompZ[{u}])/r3D_sph[],
-                            (X[]*Z[]*CompX[{u}] + Y[]*Z[]*CompY[{u}] - (X[]^2+Y[]^2)*CompZ[{u}])/(r3D_sph[]*r2D_sph[]),
-                            (-Y[]*CompX[{u}] + X[]*CompY[{u}])/r2D_sph[]
-                            ]];
-                        In All_domains; Jacobian JVol; } } }
-    }
-    }
     For pe In {1:p_max}
       // { Name VPWM_postpro~{pe}; NameOfFormulation VPWM_helmholtz_vector~{pe};NameOfSystem M~{pe};
       { Name VPWM_postpro~{pe}; NameOfFormulation VPWMN_helmholtz_vector; NameOfSystem A;
@@ -597,18 +586,6 @@ PostOperation {
     }
   EndIf
   If (flag_study==RES_TMAT)
-    {Name VPWMN_postop_test ; NameOfPostProcessing VPWMN_postpro_test;
-      Operation {
-        // Print [ E_scat   , OnElementsOf Domain     , File StrCat[myDir,"Escat.pos"]];
-          For kr In {0:nb_cuts-1}
-            Print [ E_scat_sph , OnGrid 
-                  {r_sph~{kr}*Sin[$B]*Cos[$C], r_sph~{kr}*Sin[$B]*Sin[$C], r_sph~{kr}*Cos[$B]} 
-                  {r_sph~{kr}, {sph_scan : Pi-sph_scan+(Pi-2.0*sph_scan)/(10*(npts_theta-1.0)) : (Pi-2.0*sph_scan)/(npts_theta-1.0)}, 
-                  {sph_scan : 2.0*Pi-sph_scan+(2.0*Pi-2.0*sph_scan)/(10*(npts_phi-1.0)) : (2.0*Pi-2.0*sph_scan)/(npts_phi-1.0)} }, 
-                  File StrCat[myDir,StrCat[StrCat["E_scat_onsphere_sph_M",Sprintf["%g",pe]],Sprintf["_r%g.dat",kr]]], Format Table];
-          EndFor
-      }
-    }
     For pe In {1:p_max}
       {Name VPWM_postop~{pe}; NameOfPostProcessing VPWM_postpro~{pe} ;
         Operation {