potential2ds.f90

Ground state of the time-independent Gross-Pitaevskii equation

real(kind(1.d0)) function potentialV0(x,z)
  !
  ! User supplied additional potential
  !
  ! This function takes on input the coordinates x,z and returns the value
  ! of the potential at this point *** exclusive of the harmonic trap ***
  !
  implicit none
  integer, parameter :: dp=kind(1.d0)
  
  real(dp), intent(in) :: x,z
  
  potentialV0 = 0.d0
  
  return
end function potentialV0