diff options
| author | tradke <tradke@5301f0c2-dbc4-4cee-b2f5-8d7afba4d129> | 2003-07-07 21:30:40 +0000 |
|---|---|---|
| committer | tradke <tradke@5301f0c2-dbc4-4cee-b2f5-8d7afba4d129> | 2003-07-07 21:30:40 +0000 |
| commit | bbbde8145635965d81c32959c7b9e61b506d66dd (patch) | |
| tree | 9b45dba180dfe5b8fe6ff695cf6dbc56a89d170f | |
| parent | 9ba88a4e7c228c29823c9ec36b1fecc4edfa0642 (diff) | |
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
| -rw-r--r-- | src/D3_to_D2.F | 29 |
1 files 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' |