 
   }
 
-  
   #ifdef DYNAMIC_GPU_ALLOC
   // If memory is not single allocated then free the memory every timestep.
   Free_Memory_VL_3D();
   #endif
-
-  #ifdef TIME
-  printf("cpto: %6.2f  cpfr: %6.2f\n", cpto, cpfr);
-  printf("ppmx: %6.2f  ppmy: %6.2f  ppmz: %6.2f\n", ppmx, ppmy, ppmz);
-  printf("r1x:  %6.2f  r1y:  %6.2f  r1z:  %6.2f\n", r1x, r1y, r1z);
-  printf("r2x:  %6.2f  r2y:  %6.2f  r2z:  %6.2f\n", r2x, r2y, r2z);
-  printf("pcm:  %6.2f  ie:   %6.2f  cvu:  %6.2f\n", pcm, ie, cvu);
-  printf("buff: %6.2f  dti:  %6.2f\n", buff, dti);
-  printf("cool: %6.2f\n", cool);
-  printf("arrs: %6.2f  othr: %6.2f\n", arrays, other);
-  #endif
-  #ifdef TIME
-  cudaEventRecord(stop_VL, 0);
-  cudaEventSynchronize(stop_VL);
-  cudaEventElapsedTime(&elapsedTime, start_VL, stop_VL);
-  printf("Total time for VL step: %5.3f ms\n", elapsedTime);
-  #endif //TIME
 
   // return the maximum inverse timestep
   return max_dti;

#include"subgrid_routines_3D.h"
#include"io.h"

#define TIME

__global__ void Update_Conserved_Variables_3D_half(Real *dev_conserved, Real *dev_conserved_half, Real *dev_F_x, Real *dev_F_y,  Real *dev_F_z, int nx, int ny, int nz, int n_ghost, Real dx, Real dy, Real dz, Real dt, Real gamma, int n_fields);


    memory_allocated = true;
    #endif 
  }  

  // counter for which block we're on
  int block = 0;
