CTITLESAMPFKF  -- COSINE-BELL TAPER ON TRANSFORMED TRACES
C***********************************************************************
C                 COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991            *
C***********************************************************************
CA    AUTHOR      JAMES SUN/BRUCE VERWEST
CA    DESIGNER    JAMES SUN/BRUCE VERWEST
CA    LANGUAGE    VS FORTRAN
CA    SYSTEM      IBM / CRAY
CA    WRITTEN     07/01/87
C     REVISED     02-20-89 JJC FOR SPARC PRODUCTION.
CA                                                                      00080009
CA    SUBROUTINE SAMPFKF(P)                                             00080009
CA                                                                      00080009
CA    THIS SUBROUTINE PERFORMS A COSINE-BELL TAPER ON THE TEMPORAL      00080009
CA    FOURIER TRANSFORMED INPUT TRACE FROM THE HIGH-PASS FREQUENCY TO   00080009
CA    THE HIGH-CUT FREQUENCY.                                           00080009
CA                                                                      00080009
CA**********************************************************************
CA    ARGUMENTS:                                                        00080009
CA                                                                      00080009
CA    P            INPUT/         COMPLEX ARRAY WHICH CONTAINS THE      00080009
CA                 OUTPUT         TEMPORAL FOURIER TRANSFORMED TRACE    00080009
CA                                UPON INPUT AND THE COSINE-BELL TAPERED00080009
CA                                RESULT UPON EXIT                      00080009
CA
CA**********************************************************************
C                                                                       00080009
      SUBROUTINE SAMPFKF(P)
C
      COMMON /CMPFKC/ IF1,LNT,LW,LW2,IW1,ALPHA,SCALE,CMIN,VMUTE,
     +                AFFR,IKHHI
      COMMON /CMPFKT/ NT,DT,NW,DW,NWD2,NWD21,NWP2
      COMMON /CMPFKH/ NH,DH,NKH,DKH,NKHD2,NKHD21,NKH2,IHBEG
C
      COMPLEX P(LNT/2,1)
C
      PI  = 3.14159265
      MW  = LW-IW1+1
      JW1 = IW1+1
      AMW = 1./FLOAT(MW)
      AW  = 1.
      DO 100 IW = JW1,LW
        AW   = AW-AMW
        WGHT = .5*(1.-COS(AW*PI))
  100   CALL CSSCAL(NH,WGHT,P(IW,1),LNT/2)
C
      RETURN
      END