var/convertor/kma_wave2grid/RFTF3M.inc

   1 C***********************************************************************
   2       SUBROUTINE RFTF3M (N,INC,LOT, C, WA,IFAC, CH)
   3       DIMENSION       C(INC,N)   ,CH(INC,N)  ,WA(N)    ,IFAC(*)
   4         CALL RFTF9M (N,INC,LOT, NA,C, WA, IFAC, CH)
   5 C
   6            CF = 1.0/FLOAT(N)
   7            N4=(N/4)*4
   8       IF (NA .EQ. 1) THEN
   9         IF (N4 .GE. 4) THEN
  10           DO 117 K=1,N4,4
  11           DO 117 L=1,LOT
  12             C(L,K  ) = CF*C(L,K  )
  13             C(L,K+1) = CF*C(L,K+1)
  14             C(L,K+2) = CF*C(L,K+2)
  15             C(L,K+3) = CF*C(L,K+3)
  16   117     CONTINUE
  17         ENDIF
  18         IF (N4 .NE. N) THEN
  19           DO 118 K=N4+1,N
  20           DO 118 L=1,LOT
  21             C(L,K) = CF*C(L,K)
  22   118     CONTINUE
  23         ENDIF
  24       ELSE
  25         IF (N4 .GE. 4) THEN
  26           DO 119 K=1,N4,4
  27           DO 119 L=1,LOT
  28             C(L,K  ) = CF*CH(L,K  )
  29             C(L,K+1) = CF*CH(L,K+1)
  30             C(L,K+2) = CF*CH(L,K+2)
  31             C(L,K+3) = CF*CH(L,K+3)
  32   119     CONTINUE
  33         ENDIF
  34         IF (N4 .NE. N) THEN
  35           DO 120 K=N4+1,N
  36           DO 120 L=1,LOT
  37             C(L,K) = CF*CH(L,K)
  38   120     CONTINUE
  39         ENDIF
  40       ENDIF
  41       RETURN
  42       END SUBROUTINE RFTF3M
  43