From bbbde8145635965d81c32959c7b9e61b506d66dd Mon Sep 17 00:00:00 2001 From: tradke Date: Mon, 7 Jul 2003 21:30:40 +0000 Subject: Get the input variable indices for the interpolator as INTEGER and then copy it into the CCTK_INT array. git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinAnalysis/Extract/trunk@81 5301f0c2-dbc4-4cee-b2f5-8d7afba4d129 --- src/D3_to_D2.F | 29 +++++++++++++++++++---------- 1 file changed, 19 insertions(+), 10 deletions(-) diff --git a/src/D3_to_D2.F b/src/D3_to_D2.F index 15218b7..6c7597c 100644 --- a/src/D3_to_D2.F +++ b/src/D3_to_D2.F @@ -77,7 +77,7 @@ c Local variables, passed on c Local variables, here only - INTEGER :: i, j, num_arrays, ierror + INTEGER :: i, j, num_arrays, vindex, ierror INTEGER :: interp_handle, param_table_handle, coord_system_handle CCTK_POINTER, dimension(3) :: interp_coords INTEGER, dimension(8) :: in_array_indices @@ -153,14 +153,22 @@ c ---------------------------------------------------- interp_coords(2) = CCTK_PointerTo(ys) interp_coords(3) = CCTK_PointerTo(zs) - call CCTK_VarIndex(in_array_indices(1), "extract::g00") - call CCTK_VarIndex(in_array_indices(2), "admbase::gxx") - call CCTK_VarIndex(in_array_indices(3), "admbase::gxy") - call CCTK_VarIndex(in_array_indices(4), "admbase::gxz") - call CCTK_VarIndex(in_array_indices(5), "admbase::gyy") - call CCTK_VarIndex(in_array_indices(6), "admbase::gyz") - call CCTK_VarIndex(in_array_indices(7), "admbase::gzz") - call CCTK_VarIndex(in_array_indices(8), "staticconformal::psi") + call CCTK_VarIndex(vindex, "extract::g00") + in_array_indices(1) = vindex + call CCTK_VarIndex(vindex, "admbase::gxx") + in_array_indices(2) = vindex + call CCTK_VarIndex(vindex, "admbase::gxy") + in_array_indices(3) = vindex + call CCTK_VarIndex(vindex, "admbase::gxz") + in_array_indices(4) = vindex + call CCTK_VarIndex(vindex, "admbase::gyy") + in_array_indices(5) = vindex + call CCTK_VarIndex(vindex, "admbase::gyz") + in_array_indices(6) = vindex + call CCTK_VarIndex(vindex, "admbase::gzz") + in_array_indices(7) = vindex + call CCTK_VarIndex(vindex, "staticconformal::psi") + in_array_indices(8) = vindex out_arrays(1) = CCTK_PointerTo(g00s) out_arrays(2) = CCTK_PointerTo(gxxs) @@ -200,7 +208,8 @@ c print*,'1 done' c Calculate radial derivatives and project onto sphere c ---------------------------------------------------- - call CCTK_VarIndex(in_array_indices(1), "extract::temp3d") + call CCTK_VarIndex(vindex, "extract::temp3d") + in_array_indices(1) = vindex if (conformal_state > 0) then c print*,'1 1/2' -- cgit v1.2.3