fields_f = fields_f + f_f_k * dt
!$omp parallel do default(shared) private(i,j,k) LOOP fields_f(i, j, k) = fields_f(i, j, k) + f_f_k(i, j, k) * dt ENDLOOP !$omp end parallel do