diff options
| author | ryoji <ryoji@b6f3ac56-194f-0410-8878-cdf6079d7f1b> | 2001-02-16 19:35:33 +0000 |
|---|---|---|
| committer | ryoji <ryoji@b6f3ac56-194f-0410-8878-cdf6079d7f1b> | 2001-02-16 19:35:33 +0000 |
| commit | 2da338fe56ad62fab28ea2dd3a19212ae0cb97fa (patch) | |
| tree | 28ab73e197ac7517b56b6546bdc56b37dcc801d0 /src/IDAxiOddBrillBH.F | |
| parent | d07503e7826b1af47b8907cd2032fea75a03c318 (diff) | |
For new interp.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiOddBrillBH/trunk@28 b6f3ac56-194f-0410-8878-cdf6079d7f1b
Diffstat (limited to 'src/IDAxiOddBrillBH.F')
| -rw-r--r-- | src/IDAxiOddBrillBH.F | 24 |
1 files changed, 20 insertions, 4 deletions
diff --git a/src/IDAxiOddBrillBH.F b/src/IDAxiOddBrillBH.F index 2c8dd75..141e904 100644 --- a/src/IDAxiOddBrillBH.F +++ b/src/IDAxiOddBrillBH.F @@ -382,17 +382,33 @@ c Compute eta,q,phi at the each points of cartesian grid enddo enddo - call CCTK_InterpHandle (handle, "simple_local") c Find the local origin of the spatial coordinates c ------------------------------------------------ npoints = nx*ny*nz + +c 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_Interp (ierror,cctkGH,handle,npoints,2,6,6, - $ ne,nq,abseta,q, + call CCTK_InterpLocal (ierror,cctkGH,handle,npoints,2,6,6, + $ ne,nq,etagrd,qgrd, + $ CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL, + $ abseta,q, $ CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL, - $ etagrd(1),qgrd(1),deta,dq, $ psi2dv,detapsisph,dqpsisph,detaetapsisph,detaqpsisph, $ dqqpsisph, $ CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL,CCTK_VARIABLE_REAL, |