istp/gray
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases 1 Wiki Activity

src/gray_core.f90: fix power at step 0

Browse Source 
The fraction of power in ray jk at the step zero is P₀(jk), not
necessarily 1 with more than one ray.
This commit is contained in:
Michele Guerini Rocco 2024-12-03 16:04:46 +01:00
parent 4a5c858ae7
commit beccd00b30
Signed by: rnhmjoj
GPG Key ID: BFBAF4C975F76450
1 changed files with 1 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
src/gray_core.f90
View File

@ -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, &