BUG: Dipole missing factor of 2 in number of GCs

Author: abukowski21

edu/examples/Dipole/MoreDipoleInfo.ipynb

@@ -53,7 +53,7 @@ def main(alt_minRE, alt_maxRE, lat_min, lat_max, origins, extent, nLatsPerBlock,
             origin  = -1*origin - extent
 
         lat0 = (2*(lat_max - lat_min))*origin    
-        dlat = extent * ((lat_max - lat_min)*2) / (nLatsPerBlock - nGCs)
+        dlat = extent * (2*(lat_max - lat_min)) / (nLatsPerBlock - nGCs*2)
 
         # Put latitudes down evenly (centers & corners)
         # - This forms the invariant latitudes which field lines must pass thru

edu/examples/Dipole/dipole.py

@@ -335,7 +335,7 @@ bool Grid::init_dipole_grid(Quadtree quadtree_ion, Planets planet) {
 
   precision_t lat0 = 2.0 * (max_lat - min_lat) * lat_origin;
   precision_t dlat = 2.0 * size_up_norm(1) * (max_lat -  min_lat) /
-                     (nLats - nGCs);
+                     (nLats - nGCs * 2);
 
   arma_vec lat1d(nLats);
   arma_vec lat1dDown(nLats + 1);