cptime.F90

Go to the documentation of this file.

!RJ: make interfaceable; generalization
SUBROUTINE cptime(PVCP,PTCP)
USE parkind1, ONLY : jprd, jpim
IMPLICIT NONE
REAL(KIND=JPRD) :: PVCP
REAL(KIND=JPRD) :: PTCP
!
#if defined (NEWTIMER)
! this routine should work better with OpenMP
! But doesn't work on Cray - and in any case does not return 
! CPU time for all threads combined 
INTEGER(KIND=JPIM),SAVE :: IFIRST=1
INTEGER(KIND=JPIM),SAVE :: KFIRST,KTPS
INTEGER(KIND=JPIM) :: KTICK
!     Usage of Fortran95 intrinsic function system_clock for ELAPSED time,
!     thus taking into account the parallelism if inside an open-mp region. REK
IF(ifirst.EQ.1) THEN
  ifirst=0
  CALL system_clock(kfirst,ktps)
  pvcp=0.0_jprd
  ptcp=pvcp
ELSE
  CALL system_clock(ktick)
  pvcp=REAL(ktick-kfirst,kind=jprd)/REAL(ktps,kind=jprd)
  ptcp=pvcp
ENDIF
#else
INTEGER(KIND=JPIM),SAVE :: ifirst=0
REAL(KIND=JPRD),SAVE :: ZFIRST
REAL(KIND=JPRD) :: ZSEC
!     Usage of Fortran95 intrinsic function for CPU timing.
IF(ifirst.EQ.0) THEN
  ifirst=1
  CALL cpu_time(zfirst)
  pvcp=0.0_jprd
  ptcp=pvcp
ELSE
  CALL cpu_time(zsec)
  pvcp=zsec-zfirst
  ptcp=pvcp
ENDIF
#endif
!
RETURN
END SUBROUTINE cptime