BUG: need to pass all densities

Author: aaronjridley

src/exchange_messages.cpp

@@ -18,8 +18,8 @@ bool Neutrals::exchange_old(Grid &grid) {
   int64_t nGCs = grid.get_nGCs();
 
   for (int iSpecies = 0; iSpecies < nSpecies; iSpecies++) {
-    if (species[iSpecies].DoAdvect)
-      DidWork = exchange_one_var(grid, species[iSpecies].density_scgc, false);
+    //if (species[iSpecies].DoAdvect)
+    DidWork = exchange_one_var(grid, species[iSpecies].density_scgc, false);
   }
 
   DidWork = exchange_one_var(grid, temperature_scgc, false);
@@ -1033,14 +1033,14 @@ bool exchange_one_var(Grid &grid,
   MPI_Barrier(aether_comm);
 
   // If this is a cubesphere grid, interpolate ghostcells to their proper location
-  if (grid.IsCubeSphereGrid & grid.gcInterpolationSet) {
-    report.print(3, "Interpolating Ghostcells to Proper Location");
-    var_scgc = interpolate_ghostcells(var_to_pass, grid);
-    var_to_pass = var_scgc;
-  }
+  //if (grid.IsCubeSphereGrid & grid.gcInterpolationSet) {
+  //  report.print(3, "Interpolating Ghostcells to Proper Location");
+  //  var_scgc = interpolate_ghostcells(var_to_pass, grid);
+  //  var_to_pass = var_scgc;
+  //}
 
   // Now we fill in the corners so that we don't have zero values there:
-  fill_corners(var_to_pass, nG);
+  //fill_corners(var_to_pass, nG);
 
   report.exit(function);
   return DidWork;