OpenMP实现Fortran即时矩阵向量乘法并行出错求助

我最近在用Fortran实现一个即时计算模式的矩阵-向量乘法——简单说就是不提前存储矩阵的元素，而是在计算过程中实时生成对应位置的元素值。串行版本测试下来完全正常，但用OpenMP做并行化之后，结果就不对了，翻来覆去查了好久都没找到问题，想请各位大佬帮忙排查一下！

目前用到的模块代码如下：

Module module1
    integer:: N, M
    integer, allocatable:: Jeven_A(:), Jeven_B(:), Jodd_A(:), Jodd_B(:), s(:)
END module module1

对应的subroutine核心是基于上述模块里的数组参数，实时计算矩阵元素并完成与输入向量的乘法运算，串行逻辑已经反复验证过是正确的。

内容的提问来源于stack exchange，提问作者user45857