forked from MichaelMurry/interplanetary
-
Notifications
You must be signed in to change notification settings - Fork 0
/
eqradius.m
39 lines (36 loc) · 957 Bytes
/
eqradius.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
function [ r_0 ] = eqradius( Planet )
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Description: This function defines the equatorial radius of the planet
% of interest
% Inputs:
% planet = string variable of the target planet
% Outputs
% r_0 = equatorial radius [km]
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Select equatorial radius
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
switch Planet
case 'Sun'
r_0 = 696000;
case 'Mercury'
r_0 = 2439.7;
case 'Venus'
r_0 = 6051.8;
case 'Earth'
r_0 = 6373.1363;
case 'Moon'
r_0 = 1737.4;
case 'Mars'
r_0 = 3396.19;
case 'Jupiter'
r_0 = 71492;
case 'Saturn'
r_0 = 60268.0;
case 'Uranus'
r_0 = 25559;
case 'Neptune'
r_0 = 24764;
case 'Pluto'
r_0 = 1195;
end
end