BUG: ammend prev commit: IC's and BCs were missing some important cells.

Author: abukowski21

src/neutrals_bcs.cpp

@@ -196,8 +196,8 @@ bool Neutrals::set_lower_bcs(Grid grid,
   // - This kind-of assumes nGCs=2, so may need to be updated.
   // - If the first_lower_gc is at iAlt = 1, this may cause issues.
   // - The equator-most (j-hat) grid cell will be entirely below min_alt!
-  for (int iLon = nGCs; iLon < nLons - nGCs; iLon++) {
-    for (int iLat = nGCs; iLat < nLats - nGCs; iLat++) {
+  for (int iLon = 0; iLon < nLons; iLon++) {
+    for (int iLat = 0; iLat < nLats; iLat++) {
       for (iSpecies = 0; iSpecies < nSpecies; iSpecies++) {
 
         // k-index of 1st lower ghost cell is not constant on the dipole grid.

src/neutrals_ics.cpp

@@ -175,7 +175,7 @@ bool Neutrals::initial_conditions(Grid grid,
         for (iLat = 0; iLat < nLats; iLat++) {
           for (int iSpecies = 0; iSpecies < nSpecies; iSpecies++) {
             species[iSpecies].density_scgc.subcube(
-              iLon, iLat, 0, iLon, iLat, grid.first_lower_gc(iLon, iLat)).fill(
+              iLon, iLat, 0, iLon, iLat, grid.first_lower_gc(iLon, iLat)+1).fill(
                 species[iSpecies].lower_bc_density);
           }
         }