diff --git a/lyceanem/electromagnetics/beamforming.py b/lyceanem/electromagnetics/beamforming.py
index 2063b36..5041222 100644
--- 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
index 3627217..5924175 100644
--- 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(