Fix normalize intensity (#338)

LASY-org · Dec 19, 2024 · a26c7c5 · a26c7c5
1 parent acca4f6
commit a26c7c5
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/lasy/utils/laser_utils.py b/lasy/utils/laser_utils.py
@@ -130,7 +130,8 @@ def normalize_peak_intensity(peak_intensity, grid):
         if input_peak_intensity == 0.0:
             print("Field is zero everywhere, normalization will be skipped")
         else:
-            grid.set_temporal_field(np.sqrt(peak_intensity / input_peak_intensity))
+            field *= np.sqrt(peak_intensity / input_peak_intensity)
+            grid.set_temporal_field(field)
 
 
 def get_full_field(laser, theta=0, slice=0, slice_axis="x", Nt=None):