fix interp_free issue if there is only one handle

Nek5000 · Feb 14, 2021 · 19906d4 · 19906d4
1 parent 002ed02
commit 19906d4
Showing 1 changed file with 6 additions and 5 deletions.
diff --git a/core/interp.f b/core/interp.f
@@ -197,11 +197,12 @@ subroutine interp_free(ih)
 
       common /intp_h/ ih_intp(2,INTP_HMAX)
 
-      ih_intp1 = ih_intp(1,ih)
-      ih_intp2 = ih_intp(2,ih)
-
-      call fgslib_findpts_free(ih_intp1)
-      call fgslib_findpts_free(ih_intp2)
+      if(ih_intp(2,ih) .ne. ih_intp(1,ih)) then
+        call fgslib_findpts_free(ih_intp(2,ih))
+        ih_intp(2,ih) = -1
+      endif
+      call fgslib_findpts_free(ih_intp(1,ih))
+      ih_intp(1,ih) = -1
 
       return
       end