deleted omp branch, omp code now in model_omp.cpp