From 9e617b9c86fbf6d9755efc9d2a739675dd865fe8 Mon Sep 17 00:00:00 2001
From: Geoffroy Lesur <geoffroy.lesur@univ-grenoble-alpes.fr>
Date: Sat, 2 Nov 2024 18:26:32 +0100
Subject: [PATCH] assume full sphere symetry when SPHERICAL coordinates with 1
 dimension but COMPONENTS>=2

---
 src/gridHost.cpp | 24 +++++++++++++++++-------
 1 file changed, 17 insertions(+), 7 deletions(-)

diff --git a/src/gridHost.cpp b/src/gridHost.cpp
index d31a232c..ad61f043 100644
--- a/src/gridHost.cpp
+++ b/src/gridHost.cpp
@@ -201,15 +201,25 @@ void GridHost::MakeGrid(Input &input) {
       // dir >= DIMENSIONS/ Init simple uniform grid
       for(int i = 0 ; i < np_tot[dir] ; i++) {
         // Initialize to default values
-        xstart[dir] = -0.5;
-        xend[dir] = 0.5;
-
+        #if GEOMETRY != SPHERICAL
+          xstart[dir] = -0.5;
+          xend[dir] = 0.5;
+        #elif GEOMETRY == SPHERICAL
+          if(dir == JDIR) {
+            xstart[dir] = 0;
+            xend[dir] = M_PI;
+          }
+          if(dir == KDIR) {
+            xstart[dir] = -0.5;
+            xend[dir] = 0.5;
+          }
+        #endif
         this->xbeg[dir] = xstart[dir];
         this->xend[dir] = xend[dir];
-        dx[dir](i) = 1.0;
-        x[dir](i)=0.0;
-        xl[dir](i)=-0.5;
-        xr[dir](i)=0.5;
+        dx[dir](i) = xend[dir] - xstart[dir];
+        x[dir](i)=0.5*(xstart[dir]+xend[dir]);
+        xl[dir](i)=xstart[dir];
+        xr[dir](i)=xend[dir];
       }
     }
   }