src/equilibrium.f90: set q=0 for ψ>1 in numeric case

2024-02-02 17:13:34 +01:00 · 2024-02-02 17:13:34 +01:00 · c99fc0b343
commit c99fc0b343
parent 85c84c9a48
1 changed files with 7 additions and 1 deletions
--- a/src/equilibrium.f90
+++ b/src/equilibrium.f90
@ -1195,7 +1195,13 @@ contains

      case (EQ_EQDSK_FULL, EQ_EQDSK_PARTIAL)
        ! Numerical data
+        if (psin < 1) then
          fq = q_spline%eval(psin)
+        else
+          ! q is undefined or infinite
+          fq = 0
+        end if
+
    end select
  end function fq