OpenMP Fortran代码线程不遵守OMP同步规则引发运行报错问题

• 修复Work变量的OpenMP共享属性声明
  你代码中Work是带SAVE属性的可分配数组，这类静态存储变量如果未在并行区域显式声明为SHARED，部分编译器会默认给每个线程创建独立私有副本。此时仅执行!$OMP SINGLE块的线程的私有Work副本被分配，其他线程访问的自身私有副本均为未分配状态，直接触发你遇到的报错。
  修复方案：可以在子程序开头添加OpenMP指令显式声明共享属性：

  !$OMP SHARED(Work)
  

  如果外层调用的并行块用了DEFAULT(NONE)，需要在外层并行块的属性列表中加入Work的SHARED声明。
• 移除不必要的SAVE属性
  你的Work数组每次进入子程序分配、退出前释放，完全不需要SAVE静态存储属性，删除SAVE关键字后，可避免编译器对静态变量的默认线程私有处理。如果删除SAVE后需要保持共享属性，仍需补充!$OMP SHARED(Work)声明。
• 验证同步屏障生效
  虽然!$OMP END SINGLE默认带有隐式同步屏障，你可以在分配Work的SINGLE块后手动加显式屏障进一步验证：

  !$OMP SINGLE
  ALLOCATE (Work(nr,nth) )
  !$OMP END SINGLE
  !$OMP BARRIER
  

  确认是否因为屏障失效导致其他线程提前进入循环访问未分配的Work。
• 调整Work数组的分配逻辑
  每次调用子程序都执行分配/释放操作会额外引入开销，且在多线程环境下容易出同步问题，建议将Work数组的分配逻辑移到外层并行区域执行之前，作为全局共享数组传入calcT子程序，完全规避并行区域内动态分配的同步问题。
• 新增调试打印确认分配状态
  在第一个!$OMP DO指令前添加如下调试代码，确认各线程的Work数组分配状态：

  !$OMP CRITICAL
  PRINT *, "Thread ", OMP_GET_THREAD_NUM(), " Work allocated: ", ALLOCATED(Work)
  !$OMP END CRITICAL
  

  如果打印结果显示仅线程0的Work为已分配状态，即可确认是多线程私有副本导致的问题。
• 排查并行区域嵌套问题
  确认calcT子程序的调用位置是否位于另一层OpenMP并行区域内部，如果存在嵌套并行场景，需要确认编译和运行时都开启了嵌套并行支持，否则内层并行指令的执行逻辑会出现异常。

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