1- For omp parallelisation double all omp_barrier 2- change comment from !$ to !! (because !$ is not compatible with omp parallelisation)