From b48afe77422e94346cd0413357a740e5353c47ec Mon Sep 17 00:00:00 2001
From: ryoji <ryoji@535fb057-194f-0410-b5a5-c63992f15602>
Date: Fri, 16 Feb 2001 19:44:11 +0000
Subject: For new interp.

git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/RotatingDBHIVP/trunk@22 535fb057-194f-0410-b5a5-c63992f15602
---
 src/RotatingDBHIVP.F | 25 ++++++++++++++++++++-----
 1 file changed, 20 insertions(+), 5 deletions(-)

(limited to 'src')

diff --git a/src/RotatingDBHIVP.F b/src/RotatingDBHIVP.F
index 1ff0c1e..28073c7 100644
--- a/src/RotatingDBHIVP.F
+++ b/src/RotatingDBHIVP.F
@@ -641,14 +641,29 @@ c
          enddo
       enddo
       
-      call CCTK_InterpHandle (handle, "simple_local")
-      
       npoints = nx*ny*nz
 
-      call CCTK_Interp (ierror,cctkGH,handle,npoints,3,10,10,
-     $     ne,nq,np,abseta,q,phi,
+!     Interpolator handle.
+
+      handle = -1
+
+      if (interpolation_order .eq. 1) then
+        call CCTK_InterpHandle (handle, "first-order uniform cartesian")
+      else if (interpolation_order .eq. 2) then
+        call CCTK_InterpHandle (handle, "second-order uniform cartesian")
+      else if (interpolation_order .eq. 3) then
+        call CCTK_InterpHandle (handle, "third-order uniform cartesian")
+      endif
+
+      if (handle .lt. 0) then
+        call CCTK_WARN (0, "Couldn't get handle for interpolation operator")
+      endif
+
+      call CCTK_InterpLocal (ierror,cctkGH,handle,npoints,3,10,10,
+     $     ne,nq,np,etagrd,qgrd,phigrd-pi,
+     $     CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL,
+     $     abseta,q,phi,
      $     CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL,
-     $     etagrd(1),qgrd(1),phigrd(1)-pi,deta,dq,dphi,
      $     psisph,detapsisph,dqpsisph,dphipsisph,detaetapsisph,
      $     detaqpsisph,detaphipsisph,dqqpsisph,dqphipsisph,dphiphipsisph,
      $     CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL,
-- 
cgit v1.2.3