c     -*-Fortran-*-
c     $Header$
      
#include "cctk.h"
#include "cctk_Arguments.h"
#include "cctk_Functions.h"
#include "cctk_Parameters.h"
      
      subroutine IDSWTE_uniform_charge_exact (CCTK_ARGUMENTS)
      implicit none
      DECLARE_CCTK_ARGUMENTS
      DECLARE_CCTK_FUNCTIONS
      DECLARE_CCTK_PARAMETERS
      
      integer i, j, k
      
      do k=1,cctk_lsh(3)
         do j=1,cctk_lsh(2)
            do i=1,cctk_lsh(1)
               
               if (r(i,j,k) .ge. radius) then
                  phi(i,j,k) = charge / r(i,j,k)
               else
                  phi(i,j,k) =
     $                 charge / (2*radius**3) * (3*radius**2 - r(i,j,k)**2)
               end if
               
c     Create time-symmetric initial data
               phi_p(i,j,k) = phi(i,j,k)
               
            end do
         end do
      end do
      
      end