   hipLaunchKernelGGL(PPMC_cuda, dim1dGrid, dim1dBlock, 0, 0, dev_conserved, Q_Lz, Q_Rz, nx, ny, nz, n_ghost, dz, dt, gama, 2, n_fields);
   CudaCheckError();
   #endif //PPMC
-
-
+  
+  
   // Step 2: Calculate the fluxes again
   #ifdef EXACT
   hipLaunchKernelGGL(Calculate_Exact_Fluxes_CUDA, dim1dGrid, dim1dBlock, 0, 0, Q_Lx, Q_Rx, F_x, nx, ny, nz, n_ghost, gama, 0, n_fields);

    // If memory is single allocated: memory_allocated becomes true and successive timesteps won't allocate memory.
    // If the memory is not single allocated: memory_allocated remains Null and memory is allocated every timestep.
    memory_allocated = true;

  }


  temp_potential = host_grav_potential;

  if ( !memory_allocated ){
    size_t global_free, global_total;
    CudaSafeCall( cudaMemGetInfo( &global_free, &global_total ) );
    
    // allocate memory on the GPU
    chprintf( " Allocating Hydro Memory: nfields: %d   n_cells: %d   nx: %d  ny: %d  nz: %d \n", n_fields, n_cells, nx, ny, nz );
    chprintf( " Memory needed: %f GB    Free: %f GB    Total:  %f GB  \n", n_fields*n_cells*sizeof(Real)/1e9, global_free/1e9, global_total/1e9  );
    dev_conserved = d_conserved;
    CudaSafeCall( cudaMalloc((void**)&Q_Lx,  n_fields*n_cells*sizeof(Real)) );
    CudaSafeCall( cudaMalloc((void**)&Q_Rx,  n_fields*n_cells*sizeof(Real)) );
