var/convertor/kma_wave2grid/RFFTFM.inc

   1 C***********************************************************************
   2       SUBROUTINE RFFTFM (N,INC,JUMP,LOT, R, WA,IFAC, WSAVE)
   3       DIMENSION       R(*)       ,WSAVE(*)   ,WA(N)    ,IFAC(*)
   4 cch?     REAL*8 CTIME1,CTIME2
   5 cch?     CALL XCLOCK(CTIME1,5)
   6 C
   7       IF (N .EQ. 1) RETURN
   8 C
   9       IF (JUMP .NE. 1) THEN
  10         CALL RFTF2M (N,INC,JUMP,LOT, R, WA,IFAC, WSAVE)
  11       ELSE
  12         CALL RFTF3M (N,INC,     LOT, R, WA,IFAC, WSAVE)
  13       ENDIF
  14 cch?     CALL XCLOCK(CTIME2,5)
  15 cch?     CTIME1=CTIME2-CTIME1
  16 cch?     CTIME2=1000.0*CTIME1/FLOAT(N*LOT)
  17 cch?     WRITE(6,500) INC,JUMP,N,LOT, CTIME1,CTIME2
  18 cch? 500 FORMAT(1H0,' RFFTBM: INC,JUMP,N,LOT =',4I6,'  TOTAL CPU =',F13.6
  19 cch?    *      ,' (SEC) ;  UNIT CPU =',F13.6,' (MSEC)')
  20       RETURN
  21       END SUBROUTINE RFFTFM
  22