c     f77 -o fullfcc fullfcc.f; fullfcc > fcckpt.full.dat
      PROGRAM FULLfcc
      IMPLICIT DOUBLE PRECISION (A-H,O-Z)
      PARAMETER (NX=5,SX=0.5D0)
      DIMENSION b(3,3)
      
      b(1,1) = 1.d0
      b(2,1) = 1.d0
      b(3,1) = -1.d0
      
      b(1,2) = 1.
      b(2,2) = -1.
      b(3,2) = 1.d0
      
      b(1,3) = -1.d0
      b(2,3) = 1.d0
      b(3,3) = 1.d0
      
      write(6,*) nx*nx*nx
      do i=1,nx
      do j=1,nx 
      do k=1,nx
      c1 = (i-1+sx)/nx
      c2 = (j-1+sx)/nx
      c3 = (k-1+sx)/nx
      fkx = 2 * (c1*b(1,1) + c2*b(1,2) + c3*b(1,3))
      fky = 2 * (c1*b(2,1) + c2*b(2,2) + c3*b(2,3))
      fkz = 2 * (c1*b(3,1) + c2*b(3,2) + c3*b(3,3))
      write (6,10) fkx, fky, fkz, 1.d0/nx/nx/nx
10    format(4f18.14)
      enddo
      enddo
      enddo
      stop
      end