     if ( dev_conserved[            id] < density_floor ){
       if (dev_conserved[            id] > 0){  
         dens_0 = dev_conserved[            id];
-        printf("###Thread density change  %f -> %f \n", dens_0, density_floor );
+        // printf("###Thread density change  %f -> %f \n", dens_0, density_floor );
         dev_conserved[            id] = density_floor;
         // Scale the conserved values to the new density
         dev_conserved[1*n_cells + id] *= (density_floor / dens_0);

  gx = 0.0;
  gy = 0.0;
  gz = 0.0;
  
  #ifdef COUPLE_DELTA_E_KINETIC
  Real Ekin_0, Ekin_1;

      }
      else{
        dens_0 = dev_conserved[            id];
        // printf("###Thread Negative density change  %f \n", dens_0);
        dens_0 = Average_Cell_Single_Field( 0, xid, yid, zid, nx, ny, nz, n_cells, dev_conserved );    
        // printf("### New value:   %f \n", dens_0);
      }
    }
    // #ifdef COOLING_GRACKLE
