aboutsummaryrefslogtreecommitdiff
path: root/src/AHFinder_gau.F
diff options
context:
space:
mode:
Diffstat (limited to 'src/AHFinder_gau.F')
-rw-r--r--src/AHFinder_gau.F7
1 files changed, 5 insertions, 2 deletions
diff --git a/src/AHFinder_gau.F b/src/AHFinder_gau.F
index ad79044..42a8acb 100644
--- a/src/AHFinder_gau.F
+++ b/src/AHFinder_gau.F
@@ -34,6 +34,8 @@
integer npoints
integer param_table_handle,interp_handle,coord_system_handle,sum_handle
integer ierror
+ character(len=128) :: operator
+ CCTK_INT nchars
CCTK_INT rerror,error1,error2
@@ -67,7 +69,7 @@
CCTK_REAL :: m1p, m2p, x1p, x2p, y1p, y2p, z1p, z2p
CCTK_REAL :: r12, r22, sigma
CCTK_REAL, dimension(nx,ny,nz) :: fac
- CCTK_INT :: use_att, apower, nchars
+ CCTK_INT :: use_att, apower
logical :: use_rot_att
logical :: str_comp
logical :: gaussf_exists
@@ -685,7 +687,8 @@
call CCTK_WARN(0,"Cannot create parameter table for interpolator")
endif
- call CCTK_InterpHandle (interp_handle,"Lagrange polynomial interpolation")
+ call CCTK_FortranString (nchars, interpolation_operator, operator)
+ call CCTK_InterpHandle (interp_handle, operator)
if (interp_handle .lt. 0) then
call CCTK_WARN(0,"Cannot get handle for interpolation ! Forgot to activate an implementation providing interpolation operators ??")
endif