var/convertor/kma_wave2grid/RFTB1M.inc

   1 C***********************************************************************
   2       SUBROUTINE RFTB1M (N,INC,LOT, C, WA,IFAC, CH)
   3       DIMENSION       C(INC,N)   ,CH(INC,N)  ,WA(N)    ,IFAC(*)
   4         CALL RFTB9M (N,INC,LOT, NA,C, WA, IFAC, CH)
   5       IF (NA .EQ. 1) RETURN
   6            N4 = (N/4)*4
   7       IF (N4 .GE. 4) THEN
   8          DO 117 K=1,N4,4
   9          DO 117 L=1,LOT
  10            C(L,K  ) = CH(L,K  )
  11            C(L,K+1) = CH(L,K+1)
  12            C(L,K+2) = CH(L,K+2)
  13            C(L,K+3) = CH(L,K+3)
  14   117    CONTINUE
  15       ENDIF
  16       IF (N4 .NE. N) THEN
  17          DO 118 K=N4+1,N
  18          DO 118 L=1,LOT
  19            C(L,K) = CH(L,K)
  20   118    CONTINUE
  21       ENDIF
  22       RETURN
  23       END SUBROUTINE RFTB1M
  24