Parametrized some important values like Omega. Simplyfied the code a little. Wrote some useful comments. Readded rotation along tangents - i.e. using vecotor field for computation.