From f0428118d714b53c364fd42b5c8fd0ed102e99b0 Mon Sep 17 00:00:00 2001
From: schnetter <schnetter@40f6ab95-0e4f-0410-8daa-ee8d7420be1d>
Date: Mon, 21 Aug 2006 14:18:45 +0000
Subject: Correct error in setting surface radii when the parameters
 set_spherical or set_elliptic are used.

git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/SphericalSurface/trunk@24 40f6ab95-0e4f-0410-8daa-ee8d7420be1d
---
 src/radius.c | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/src/radius.c b/src/radius.c
index 102ae1d..0356259 100644
--- a/src/radius.c
+++ b/src/radius.c
@@ -69,7 +69,8 @@ void SphericalSurface_Set (CCTK_ARGUMENTS)
       
       for (j=0; j<nphi[n]; ++j) {
         for (i=0; i<ntheta[n]; ++i) {
-          sf_radius[n] = radius[n];
+          int const ind = i + maxntheta * (j + maxnphi * n);
+          sf_radius[ind] = radius[n];
         }
       }
       
@@ -131,12 +132,13 @@ void SphericalSurface_Set (CCTK_ARGUMENTS)
       
       for (j=0; j<nphi[n]; ++j) {
         for (i=0; i<ntheta[n]; ++i) {
+          int const ind = i + maxntheta * (j + maxnphi * n);
           CCTK_REAL const theta = sf_origin_theta[n] + i * sf_delta_theta[n];
           CCTK_REAL const phi = sf_origin_phi[n] + j * sf_delta_phi[n];
           CCTK_REAL const x2 = pow (sin(theta) * cos(phi), 2);
           CCTK_REAL const y2 = pow (sin(theta) * sin(phi), 2);
           CCTK_REAL const z2 = pow (cos(theta)           , 2);
-          sf_radius[n] = 1.0 / sqrt (x2 / rx2 + y2 / ry2 + z2 / rz2);
+          sf_radius[ind] = 1.0 / sqrt (x2 / rx2 + y2 / ry2 + z2 / rz2);
         }
       }
       
-- 
cgit v1.2.3