Divergence function

I am working with the divergence function. I have a meridional flux and zonal flux in g/ms. Which are the units of the value of the divergence results of the next command

d=divergence(Fm,Fz)?

I have a problem between the results found with the grads using the comand hdivg and with matlab using the comand divergence.

The pattern is the same but there is a factor of 1e6 between the values. I don't understand why?. Can help me anybody?

Thank you

[E]

I had a similar problem with curlz function. I found it tend to scale with mesh grid size