ncpus.F

Go to the documentation of this file.

C> @file
C> Set number of cpus.
C>
C> ### Program History Log
C> Date | Programmer | Comments
C> -----|------------|---------
C> 94-08-19 | Iredell | Initial.
C> 98-11-09 | Vuong | Add doc>block and remove cray references.
C> 1998-12-18 | Iredell | IBM SMP version.
C> 2010-11-16 | Slovacek | Linux must have different call.
C> 2012-11-01 | Mirvis | Multi-threading on LINUX-IBM/TIDE.
C>
C> @author Iredell @date 94-08-19
 
C> Set number of CPUs - the number of processors over which
C> to parallelize.
C>
C> @return Number of CPUs assigned.
C>
C> @author Iredell @date 94-08-19
      FUNCTION ncpus()
      INTEGER tid
C    Obtain thread number
#ifdef openmp
      INTEGER omp_get_num_threads,omp_get_thread_num
!$OMP PARALLEL PRIVATE(TID)
      tid = omp_get_thread_num()
!      PRINT *, '...............thread # ', TID
      if (tid. eq. 0) then
        ncpus=omp_get_num_threads()
!        PRINT *, 'totaly #------------------- of threads = ',NCPUS
      endif
!$OMP END PARALLEL
#else
      tid = 0
      ncpus = 1
#endif
      RETURN
      END