diff --git a/ElectromagneticScattering/scattering.pro b/ElectromagneticScattering/scattering.pro index 8768a0c1d6ba0932935c3c89b4f333b9a693ae0d..b4dde40106c78d868bf25f980451d950f0d1da2c 100644 --- a/ElectromagneticScattering/scattering.pro +++ b/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]}