Corrected label to new format in line with vtk mesh requirements, and…

… added area normalisation component to time domain scattering.
LyceanEM · Dec 17, 2024 · eb7089c · eb7089c
1 parent 240236b
commit eb7089c
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 1 deletion.
diff --git a/lyceanem/electromagnetics/beamforming.py b/lyceanem/electromagnetics/beamforming.py
@@ -1353,7 +1353,7 @@ def GPUBeamformingMap(Etheta, Ephi, DirectivityMap, az_range, el_range, waveleng
 def create_display_mesh(
     field_data,
     field_radius=1.0,
-    label="Poynting Vector (Magnitude)",
+    label='Poynting_Vector_(Magnitude_(W/m2))',
     log_type="amplitude",
     plot_max=None,
 ):

diff --git a/lyceanem/models/time_domain.py b/lyceanem/models/time_domain.py
@@ -112,6 +112,8 @@ def calculate_scattering(
             else:
                 conformal_E_vectors = desired_E_axis.reshape(num_sources, 3)
 
+    # convert from field strengths to diffraction integral components
+    conformal_E_vectors = conformal_E_vectors * aperture_coords.point_data['Area'].reshape(-1, 1)
     if scattering == 0:
         # only use the aperture point cloud, no scattering required.
         scatter_points = meshio.Mesh(