potential3dg.f90

Ground state of the time-independent Gross-Pitaevskii equation

real(kind(1.d0)) function potentialV(x,y,z)
  !
  ! User supplied potential
  !
  ! This function takes on input the coordinates x,y,z and returns the value
  ! of the potential at this point 
  !
  ! all values in atomic units
  !
  use system3Dg
  implicit none
  integer, parameter :: dp=kind(1.d0)
  
  real(dp), intent(in) :: x,y,z

  real(dp) :: wx2,wy2,wz2
  
  wx2 = (1.36784d-14)**2
  wy2 = wx2
  wz2 = 8.d0*wx2
  
  potentialV = 0.5d0*mass*(wx2*x**2 + wy2*y**2 + wz2*z**2)
  
  return
end function potentialV