diff options
Diffstat (limited to 'src/D3_to_D2.F')
-rw-r--r-- | src/D3_to_D2.F | 13 |
1 files changed, 10 insertions, 3 deletions
diff --git a/src/D3_to_D2.F b/src/D3_to_D2.F index 7bc7586..da2294e 100644 --- a/src/D3_to_D2.F +++ b/src/D3_to_D2.F @@ -4,7 +4,8 @@ c ======================================================================== SUBROUTINE D3_to_D2(cctkGH,conformal_state,do_ADMmass,do_momentum,do_spin, - & Psi_power,origin,myproc,interpolation_order,Dx,Dy,Dz,Psi, + & Psi_power,origin,myproc,interpolation_operator,interpolation_order, + & Dx,Dy,Dz,Psi, & g00,gxx,gxy,gxz,gyy,gyz,gzz,hxx,hxy,hxz,hyy,hyz,hzz, & x,y,z,eta,Nt,Np,theta,phi, & Psis,g00s,gxxs,gxys,gxzs,gyys,gyzs,gzzs,dPsis,dgxxs,dgxys,dgxzs, @@ -46,6 +47,8 @@ c Input variables & Extract_temp3d LOGICAL :: & do_ADMmass(2) + CCTK_STRING,INTENT(IN) :: + & interpolation_operator c Output variables @@ -77,7 +80,9 @@ c Local variables, here only CCTK_POINTER, dimension(8) :: out_arrays CCTK_INT, dimension(8) :: out_array_type_codes - character(30) options_string + character(128) options_string + character(128) operator + CCTK_INT nchars c ------------------------------------------------------------------------ @@ -120,7 +125,9 @@ c -------------------------------------------------------------------- param_table_handle = -1 coord_system_handle = -1 - 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 |