From 09b04ede1e44a309077b5e94001f71c902f3adde Mon Sep 17 00:00:00 2001
From: Guillaume Demesy <guillaume.demesy@fresnel.fr>
Date: Thu, 9 Jan 2020 23:42:58 +0100
Subject: [PATCH] Joule
---
ElectromagneticScattering/scattering.pro | 15 ++++++++-------
1 file changed, 8 insertions(+), 7 deletions(-)
diff --git a/ElectromagneticScattering/scattering.pro b/ElectromagneticScattering/scattering.pro
index 8768a0c..b4dde40 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]}
--
GitLab