pragma omp parallel for before do_step() and grid_copy()