Returns the gravitational acceleration at the Earth's surface as function of latitude, based on Clairaut's formula.
TODO Reference
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real(kind=rk), | intent(in) | :: | latitude | Latitude [rad] |
function gravityClairaut(latitude) result(grav)
!! Returns the gravitational acceleration at the Earth's surface as function
!! of latitude, based on Clairaut's formula.
!!
!! TODO Reference
real(kind=rk),intent(in) :: latitude !! Latitude [rad]
grav = 9.780327_rk*(1.0026454_rk &
-0.0026512_rk*cos(2*latitude)&
+0.0000058_rk*cos(2*latitude)**2)
endfunction gravityClairaut