From beccd00b300ff0fcac9180e1796b270eab674a21 Mon Sep 17 00:00:00 2001 From: Michele Guerini Rocco Date: Tue, 3 Dec 2024 16:04:46 +0100 Subject: [PATCH] src/gray_core.f90: fix power at step 0 MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit The fraction of power in ray jk at the step zero is P₀(jk), not necessarily 1 with more than one ray. --- src/gray_core.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/gray_core.f90 b/src/gray_core.f90 index db73a26..5a2add2 100644 --- a/src/gray_core.f90 +++ b/src/gray_core.f90 @@ -237,7 +237,7 @@ contains do jk=1,params%raytracing%nray ! Save the step "zero" data call store_ray_data(params, equil, results%tables, & - i=0, jk=jk, s=stv(jk), P0=one, pos=yw(1:3,jk), & + i=0, jk=jk, s=stv(jk), P0=p0jk(jk), pos=yw(1:3,jk), & psi_n=-one, B=zero, b_n=[zero,zero,zero], k0=ak0, & N_pl=zero, N_pr=zero, N=yw(4:6,jk), N_pr_im=zero, & n_e=zero, T_e=zero, alpha=zero, tau=zero, dIds=zero, &