[cig-commits] [commit] devel: added a comment about Radau's approximation of Clairaut's equation for first-order terms of ellipticity. (fe1cc1c)
cig_noreply at geodynamics.org
cig_noreply at geodynamics.org
Wed Apr 30 12:44:09 PDT 2014
Repository : ssh://geoshell/specfem3d_globe
On branch : devel
Link : https://github.com/geodynamics/specfem3d_globe/compare/124baaa8969bcccf5b5896f1a595b2d49dd09e6f...60794754721e795e0a3cb34e89d65d86c81ecf10
>---------------------------------------------------------------
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
More information about the CIG-COMMITS
mailing list