Hi grads's users, I have a question concerning the way hdivg function is done. I looked at the script code and it seems that the finite difference method does not exactly correspond to the horizontal divergence expression in spherical coordinates. More specifically the term in tan(latitude) is ignored. Is this term negligible ? Thanks in advance RM