Divergence of Horizontal Wind using Spherical Coordinates

I'm wanting to compute divergence of the horizontal wind field from U and V components of wind (meters per second), but using a spherical coordinate system based X and Y (latitude/longitude):

Am I able to do this using the divergence function in MATLAB or will I need to do this outside of the function? Thanks.