SUBROUTINE FPW86SUP ( RHOT, GRHOT, WK, WK0, WK1, WK2, WK3, WK4, WK5, WK6, WK7 )

Argument Definitions (+ indicates altered content)
REAL*8             GRHOT(*),   RHOT(*),   +WK(*),     +WK0(*)
REAL*8             WK1(*),     WK2(*),     WK3(*),     WK4(*)
REAL*8             WK5(*),     WK6(*),     WK7(*)
Description
                                                                        * 
  Copyright (C) 2010 Jacek Kobus                    * 
                                                                        * 
  This program is free software; you can redistribute it and/or modify  * 
  it under the terms of the GNU General Public License version 2 as     * 
  published by the Free Software Foundation.                            * 
                                                                        * 
 
### fpw86sup ### 
    Calculates exchange potential according to a formula of Perdew and 
    Wang Yue (PRB 33(12) (1986) 8800) and returns it in grhot array
Source file:fpw86sup.f
Intrinsic Functions Called
GENERIC*8          ABS,        SQRT
External Functions and Subroutines Called
REAL*8             FDFTPOT,    FF,         FFDP,       FFP
SUBROUTINE         COPY,       N2F,        NFNG,       SCAL
Parameter Variables Used
REAL*8             CONST13            (CONST13 = 1.D0/3.D0)
REAL*8             CONST43            (CONST43 = 4.D0/3.D0)
Local Variables (+ indicates altered content)
INTEGER           +I
REAL*8            +AKF,       +FDFTCOEFF, +S,         +T
REAL*8            +U,         +W1,        +W2
Referenced Common Block Variables (+ indicates altered content)
CONSFP             REAL*8             HALF,      PII
CONSFP             REAL*8             PRECIS,    THREE
CONSFP             REAL*8             TWO
CONSINT            INTEGER            IONE
GRIDI              INTEGER            MXSIZE
SCFR               REAL*8             ALPHAF