diff --git a/src/gray_tables.f90 b/src/gray_tables.f90 index ffe60bc..d8c024e 100644 --- a/src/gray_tables.f90 +++ b/src/gray_tables.f90 @@ -410,7 +410,8 @@ contains psi_n = sol(1) ! Handle spurious or no solutions - if (info /= 1 .or. psi_n < 0 .or. psi_n > 1) cycle + if (info /= 1) cycle + if (psi_n < 0 .or. psi_n > 1) cycle end block ! Compute and print the ψ_n(R,z) = ψ_n₀ contour