diff options
Diffstat (limited to 'src/D3_to_D2.F')
-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' |