potential2dg.f90

Ground state of the time-independent Gross-Pitaevskii equation

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

  real(dp) :: wx,wz
  
  wz = 1.51982d-14
  wx = 0.25d0*wz
  
  potentialV = 0.5d0*mass*(wx**2*x**2 + wz**2*z**2)
  
  return
end function potentialV