src/gray_tables.f90: flux-averages do not go past ρ=1

This commit is contained in:
Michele Guerini Rocco 2025-02-10 10:30:41 +01:00
parent 5043dd843d
commit 9efaf045e2
Signed by: rnhmjoj
GPG Key ID: BFBAF4C975F76450

View File

@ -178,7 +178,7 @@ contains
! Δρ_p for the uniform ρ_p grid
drho_p = 1.0_wp_ / (n_grid - 1)
do i = 0, n_grid
do i = 0, n_grid-1
rho_p = i * drho_p
rho_t = equil%pol2tor(rho_p)
call equil%flux_average(rho_p, &