From 6f1bb35c25d7c3b6601d2efff132d29a743e5cd0 Mon Sep 17 00:00:00 2001 From: huixingjian Date: Thu, 11 Jan 2024 19:59:29 +0100 Subject: [PATCH] re-add z0*cos --- src/laser/MultiLaser.cpp | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/laser/MultiLaser.cpp b/src/laser/MultiLaser.cpp index 9da99b3ec8..b7df328cb8 100644 --- 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);