Joule

ElectromagneticScattering/scattering.pro

@@ -460,6 +460,7 @@ PostProcessing {
     { Name PW_postpro  ; NameOfFormulation PW_helmholtz_vector; NameOfSystem P;
       Quantity {
         { Name E_tot             ; Value { Local { [{u}+Einc[]]; In All_domains; Jacobian JVol; } } }
+        { Name normE_tot         ; Value { Local { [Norm[{u}+Einc[]]]; In All_domains; Jacobian JVol; } } }
         { Name E_scat            ; Value { Local { [{u}]; In All_domains; Jacobian JVol; } } }
         { Name E_scat_sph        ; Value { Local { [Vector[
                                     (X[]*CompX[{u}] + Y[]*CompY[{u}] + Z[]*CompZ[{u}])/r3D_sph[],
@@ -566,13 +567,13 @@ PostOperation {
   If (flag_study==RES_PW)
     {Name PW_postop; NameOfPostProcessing PW_postpro ;
       Operation {
-        Print [ E_tot   , OnElementsOf Domain     , File StrCat[myDir,"E_tot.pos"]];
-        Print [ E_scat  , OnElementsOf All_domains, File StrCat[myDir,"E_scat.pos"]];
-        // Print [ E_scat_sph    , OnElementsOf All_domains, File StrCat[myDir,"E_scat_sph.pos"]];
-        Print [ H_scat  , OnElementsOf Domain     , File StrCat[myDir,"H_scat.pos"]];
-        Print [ H_tot   , OnElementsOf All_domains, File StrCat[myDir,"H_tot.pos"]];
-        Print [ Poy_dif , OnElementsOf All_domains, File StrCat[myDir,"Poy_dif.pos"]];
-        Print [ Poy_tot , OnElementsOf Domain     , File StrCat[myDir,"Poy_tot.pos"]];
+        Print [ E_tot    , OnElementsOf Domain     , File StrCat[myDir,"E_tot.pos"]];
+        Print [ E_tot    , OnElementsOf Domain     , File StrCat[myDir,"E_tot.pos"]];
+        Print [ normE_tot, OnElementsOf Scat_In    , File StrCat[myDir,"normE_tot.pos"]];
+        Print [ H_scat   , OnElementsOf Domain     , File StrCat[myDir,"H_scat.pos"]];
+        Print [ H_tot    , OnElementsOf All_domains, File StrCat[myDir,"H_tot.pos"]];
+        Print [ Poy_dif  , OnElementsOf All_domains, File StrCat[myDir,"Poy_dif.pos"]];
+        Print [ Poy_tot  , OnElementsOf Domain     , File StrCat[myDir,"Poy_tot.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]}