 
 void lbm_simulation_destroy(lbm_simulation* lbm_sim)
 {
-    lbm_vars_dealloc(&lbm_sim->h_vars);
     lbm_vars_cuda_dealloc(&lbm_sim->d_vars);
     HANDLE_ERROR(cudaFree(lbm_sim->d_pal_lat));
     free(lbm_sim);

    lbm_sim->nl = nx * ny * nz;

    lbm_sim->omega = omega;
   
    lbm_vars_alloc(&lbm_sim->h_vars, lbm_sim->nl);
    
    lbm_sim->switch_f0_f1 = false;


}

struct lbm_simulation{
    lbm_vars d_vars;
    double* d_pal_lat;
    dim3 dimComputationGrid, dimComputationBlock;
    size_t shared_mem_size;
