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

closed as off-topic by Kuba, bobthechemist, Oleksandr R., Yves Klett, m_goldberg Apr 24 at 15:23

  • The question does not concern the technical computing software Mathematica by Wolfram Research. Please see the help center to find out about the topics that can be asked here.
If this question can be reworded to fit the rules in the help center, please edit the question.

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
This question appears to be off-topic because OP has not responded to the comments to clarify the question. –  Kuba Apr 24 at 12:07
add comment