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

src/gray_project.f90: avoid division by zero

Browse Source
This commit is contained in:
Michele Guerini Rocco 2025-03-20 20:03:28 +01:00
parent fadf1ac7c7
commit 65b9d129ea
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_project.f90
View File

@ -269,7 +269,7 @@ contains
end if
! Signal that J_φ has alternating sign
if ((I_tot - abs(sum(J_phi*self%dA)))/ I_tot > 0.1) then
if (I_tot - abs(sum(J_phi*self%dA)) > 0.1_wp_*I_tot) then
rho_avg_J = -rho_avg_J
drho_avg_J = -drho_avg_J
rho_max_J = -rho_max_J