re-add z0*cos

Hi-PACE · Jan 11, 2024 · 6f1bb35 · 6f1bb35
1 parent 615061b
commit 6f1bb35
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/laser/MultiLaser.cpp b/src/laser/MultiLaser.cpp
@@ -1030,7 +1030,7 @@ MultiLaser::InitLaserSlice (const amrex::Geometry& geom, const int islice, const
                         arr(i, j, k, comp + 1 ) = 0._rt;
                     }
                     // Compute envelope for time step 0
-                    Complex diffract_factor = 1._rt + I * (z-zfoc*std::cos(propagation_angle_yz)+z0) * 2._rt/( k0 * w0 * w0 );
+                    Complex diffract_factor = 1._rt + I * (z-zfoc*std::cos(propagation_angle_yz)+z0*std::cos(propagation_angle_yz)) * 2._rt/( k0 * w0 * w0 );
                     Complex inv_complex_waist_2 = 1._rt /( w0 * w0 * diffract_factor );
                     Complex prefactor = a0/diffract_factor;
                     Complex time_exponent = z*z/(L0*L0);