diff options
| author | ryoji <ryoji@535fb057-194f-0410-b5a5-c63992f15602> | 2001-02-16 19:44:11 +0000 |
|---|---|---|
| committer | ryoji <ryoji@535fb057-194f-0410-b5a5-c63992f15602> | 2001-02-16 19:44:11 +0000 |
| commit | b48afe77422e94346cd0413357a740e5353c47ec (patch) | |
| tree | e1ed08d65689bdea05d438409f6c6194e7ed423a /src | |
| parent | 28ac93e767b7ecb0616a7bb4f73c2d738a9075b9 (diff) | |
For new interp.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/RotatingDBHIVP/trunk@22 535fb057-194f-0410-b5a5-c63992f15602
Diffstat (limited to 'src')
| -rw-r--r-- | src/RotatingDBHIVP.F | 25 |
1 files changed, 20 insertions, 5 deletions
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, |