Take the 2-minute tour ×
Mathematica Stack Exchange is a question and answer site for users of Mathematica. It's 100% free, no registration required.

How do you change Curl in Mathematica so that it uses theta as latitude rather than as colatittude? Changing theta to Pi/2 - theta doesn't seem to get all the signs right in the final result.

Basically, I want to use the curl in a right handed spherical coordinate system with unit vectors going eastward, northward, and vertical. If I take the Curl[R[r, theta, phi]], where R is a vector{R, Theta, Phi} and change theta to Pi/2 - theta so that the angle used in the result from Mathematica's Curl is measured from the equator instead of the pole (i.e., it is now the Latitude instead of coLatitude), the result does not quite agree with what I find documented in the literature. Two of the three vector components agree in terms and signs but the third Mathematica component gives

Tan[theta]*Theta[r, theta, phi]/r Sec[theta]*D[Theta[r, theta, phi], phi]/r + 
   Phi[r, theta, phi]/r

whereas the literature gives

Tan[theta]*Theta[r, theta, phi]/r + Sec[theta]*D[Theta[r, theta, phi], phi]/r - 
   Phi[r, theta, phi]/r

which agrees in terms but not in sign. Multiplying by -1 still leaves one term in disagreement.

So, I guess a better form of the question is how do you change the Mathematica curl from its current Eastward, Southward, Vertical coordinate system to an Eastward, Northward, Vertical coordinate system that uses Latitude rather than colatitude for theta.?

share|improve this question
6  
I think I know what you mean, but I don't want to waste time guessing wrong. Could you be more specific about what you are trying to do? Please include the code that you tried and the result you think is incorrect. –  Jens Feb 25 '13 at 1:47
add comment

Know someone who can answer? Share a link to this question via email, Google+, Twitter, or Facebook.

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.