Ошибка в программа (Fortran)

[Далила]

Программа, в принципе, работает правильно. Ошибка не в алгоритме. Если запускать процесс 10 раз, то 8 из них будут верны, а 2 раза будет выдавать не те числа.

Задача состоит в том, чтобы проверить признак делимости на 9(сумма цифр числа, делящегося на 9, также делится на 9) для натуральных чисел, не превосходящих заданного N.

program test
include 'omp_lib.h'
integer j, k, l, s, num, p

n = 100
l = 0
k = 2
do while(k>1)
k = int(n / (10**l))
l = l + 1
end do
!$OMP PARALLEL NUM_THREADS(4) private(p,s,j)
!$OMP DO SCHEDULE (STATIC,1) 
do p = 1, n
j = l
s = 0
num = p
do while(j >= 0)
s = s + int(num / 10**j)
num = num - (10**j) * int(num / 10**j)
j = j - 1
enddo
if (mod(s,9) .eq. 0) then
print *, p, s, omp_get_thread_num()+1
end if
end do
!$OMP END DO
!$OMP END PARALLEL

end

[Vago]

Сообщение от Далила

Если запускать процесс 10 раз, то 8 из них будут верны, а 2 раза будет выдавать не те числа.

Распараллеливание в этих "10 разах" включено было? Или выключено?

[Далила]

Было включено

[Vago]

Без него - никак? (у меня сейчас нет этих библиотек).

1. Попробуйте отключить (в один поток у меня работает всегда правильно).
2. Если, всё-таки, никак - мне кажется, что num тоже в списке private() должно быть...

[Далила]

Без него, к сожалению, никак.

"...num тоже в списке private() должно быть..." и это пробовала, но тоже не помогает...