!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
!                                                                       
!   FFTPACK 5.0                                                         
!   Copyright (C) 1995-2004, Scientific Computing Division,             
!   University Corporation for Atmospheric Research                     
!   Licensed under the GNU General Public License (GPL)                 
!                                                                       
!   Authors:  Paul N. Swarztrauber and Richard A. Valent                
!                                                                       
!   $Id: c1fm1f.f,v 1.2 2004/06/15 21:08:32 rodney Exp $                
!                                                                       
!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
                                                                        
      SUBROUTINE C1FM1F (N,INC,C,CH,WA,FNF,FAC) 
      COMPLEX       C(*) 
      REAL       CH(*),     WA(*),      FAC(*) 
!                                                                       
! FFTPACK 5.0 auxiliary routine                                         
!                                                                       
      INC2 = INC+INC 
      NF = FNF 
      NA = 0 
      L1 = 1 
      IW = 1 
      DO 125 K1=1,NF 
         IP = FAC(K1) 
         L2 = IP*L1 
         IDO = N/L2 
         LID = L1*IDO 
         NBR = 1+NA+2*MIN(IP-2,4) 
         write(*,*) wa(iw),wa(iw+1) 
         GO TO (52,62,53,63,54,64,55,65,56,66),NBR 
   52    CALL C1F2KF (IDO,L1,NA,C,INC2,CH,2,WA(IW)) 
         GO TO 120 
   62    CALL C1F2KF (IDO,L1,NA,CH,2,C,INC2,WA(IW)) 
         GO TO 120 
   53    CALL C1F3KF (IDO,L1,NA,C,INC2,CH,2,WA(IW)) 
         GO TO 120 
   63    CALL C1F3KF (IDO,L1,NA,CH,2,C,INC2,WA(IW)) 
         GO TO 120 
   54    CALL C1F4KF (IDO,L1,NA,C,INC2,CH,2,WA(IW)) 
         GO TO 120 
   64    CALL C1F4KF (IDO,L1,NA,CH,2,C,INC2,WA(IW)) 
         GO TO 120 
   55    CALL C1F5KF (IDO,L1,NA,C,INC2,CH,2,WA(IW)) 
         GO TO 120 
   65    CALL C1F5KF (IDO,L1,NA,CH,2,C,INC2,WA(IW)) 
         GO TO 120 
   56    CALL C1FGKF (IDO,IP,L1,LID,NA,C,C,INC2,CH,CH,                  &
     &     1,WA(IW))                                                    
         GO TO 120 
   66    CALL C1FGKF (IDO,IP,L1,LID,NA,CH,CH,2,C,C,                     &
     &     INC2,WA(IW))                                                 
  120    L1 = L2 
         IW = IW+(IP-1)*(IDO+IDO) 
         IF(IP .LE. 5) NA = 1-NA 
  125 END DO 
      RETURN 
      END