[cig-commits] [commit] devel, master: added a comment about Radau's approximation of Clairaut's equation for first-order terms of ellipticity. (fe1cc1c)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/specfem3d_globe

On branches: devel,master
Link       : https://github.com/geodynamics/specfem3d_globe/compare/bc58e579b3b0838a0968725a076f5904845437ca...be63f20cbb6f462104e949894dbe205d2398cd7f

>---------------------------------------------------------------

commit fe1cc1ce8635a6943d36461cdd1b583752f4d3dc
Author: Dimitri Komatitsch <komatitsch at lma.cnrs-mrs.fr>
Date:   Wed Apr 30 21:41:31 2014 +0200

    added a comment about Radau's approximation of Clairaut's equation for first-order terms of ellipticity.


>---------------------------------------------------------------

fe1cc1ce8635a6943d36461cdd1b583752f4d3dc
 ...ximation_of_Clairaut_ellipticity_1953_first_page.pdf | Bin 0 -> 86029 bytes
 src/shared/make_ellipticity.f90                         |   3 +++
 2 files changed, 3 insertions(+)

diff --git a/doc/Jeffreys_The_figures_of_rotating_planets_about_Radau_approximation_of_Clairaut_ellipticity_1953_first_page.pdf b/doc/Jeffreys_The_figures_of_rotating_planets_about_Radau_approximation_of_Clairaut_ellipticity_1953_first_page.pdf
new file mode 100644
index 0000000..7e076c6
Binary files /dev/null and b/doc/Jeffreys_The_figures_of_rotating_planets_about_Radau_approximation_of_Clairaut_ellipticity_1953_first_page.pdf differ
diff --git a/src/shared/make_ellipticity.f90 b/src/shared/make_ellipticity.f90
index bbe960c..a6f9f39 100644
--- a/src/shared/make_ellipticity.f90
+++ b/src/shared/make_ellipticity.f90
@@ -143,6 +143,9 @@
 
   do i=2,NR
     call intgrl(i_rho,r,1,i,rho,s1,s2,s3)
+!! DK DK Radau approximation of Clairaut's equation for first-order terms of ellipticity, see e.g. Jeffreys H.,
+!! DK DK The figures of rotating planets, Mon. Not. R. astr. Soc., vol. 113, p. 97-105 (1953).
+!! DK DK The Radau approximation is mentioned on page 97.
     call intgrl(i_radau,r,1,i,radau,s1,s2,s3)
     z=(2.0d0/3.0d0)*i_radau/(i_rho*r(i)*r(i))
     eta(i)=(25.0d0/4.0d0)*((1.0d0-(3.0d0/2.0d0)*z)**2.0d0)-1.0d0