add threaded region check to seed dispersal call

ESCOMP · Oct 19, 2023 · 732d3b9 · 732d3b9
1 parent c811ae9
commit 732d3b9
Showing 1 changed file with 13 additions and 0 deletions.
diff --git a/src/utils/clmfates_interfaceMod.F90 b/src/utils/clmfates_interfaceMod.F90
@@ -2669,8 +2669,21 @@ subroutine WrapGlobalSeedDispersal(this,is_restart_flag)
 
    logical :: set_restart_flag ! local logical variable to pass to IsItDispersalTime
                                ! if optional is_restart_flag is true
+#ifdef _OPENMP
+   logical, external :: omp_in_parallel
+#endif
 
    type (neighbor_type),  pointer :: neighbor
+
+   ! Check to see if we are not in a threaded region.  Fail the run if this returns true.
+#ifdef _OPENMP
+   write(iulog,*) 'omp check: ', omp_in_parallel()
+   if (omp_in_parallel()) then
+      call endrun(msg='clmfates interface error: MPI routine called within threaded region'//&
+           errMsg(sourcefile, __LINE__))
+   end if
+#endif
+
    ! This should only be run once per day
    if(is_beg_curr_day()) then