diff --git a/src/gray_project.f90 b/src/gray_project.f90 index 74e51c0..4cb4b19 100644 --- a/src/gray_project.f90 +++ b/src/gray_project.f90 @@ -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