-
Notifications
You must be signed in to change notification settings - Fork 2
/
namelist.in
50 lines (42 loc) · 2.58 KB
/
namelist.in
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
39
40
41
42
43
44
45
46
47
48
49
! positions and velocities of the bodies in the solar system (in km, km/s)
! time below is 2014-02-02 20:18:00
&initial_state
x(1:10)=0.d0, 1.563021412664830d+07, -9.030189258080004d+07, -1.018974476358996d+08 ,
-2.443763125844157d+08, -2.35165468275322006d+08, -1.011712827283427d+09 ,
2.934840841770302d+09, 4.055112581124043d+09, 9.514009594170194d+08
y(1:10)=0.d0, 4.327888220902108d+07, 5.802615456116644d+07, 1.065689158175689d+08,
4.473211564076996d+07, 7.421837640432589d+08, -1.077496255617324d+09 ,
6.048399137411513d+08, -1.914578873112663d+09, -4.776029500570151d+09
z(1:10)=0.d0, 2.102123103174893d+06, 6.006513603716755d+06, -3.381951053601424d+03,
6.935657388967808d+06, 2.179850895804323d+06, 5.901251900068215d+07 ,
-3.576451387567792d+07, -5.400973716179796d+07, 2.358627841705075d+08
ux(1:10)=0.d0, -5.557001175482630d+01, -1.907374632532257d+01, -2.201749257051057d+01 ,
-3.456935754608896d+00, -1.262559929908801d+01, 6.507898648442419d+00 ,
-1.433852081777671d+00, 2.275119229131818d+00, 5.431808363374300d+00
uy(1:10)=0.d0, 1.840863017229157d+01, -2.963461693326599d+01, -2.071074857788741d+01,
-2.176307370133160d+01, -3.332552395475581d+00, -6.640809674126991d+00 ,
6.347897341634990d+00, 4.942356914027413d+00, -2.387056445508962d-02
uz(1:10)=0.d0 ,6.602621285552567d+00, 6.946391255404438d-01, 1.575245213712245d-03 ,
-3.711433859326417d-01 ,2.962741332356101d-01, -1.434198106014633d-01 ,
4.228261484335974d-02 ,-1.548950389954096d-01, -1.551877289694926d+00
! product of Gm for sun, then all other bodies (km^3 s-2)
gm=1.327124400d11,
22032.09d0, 324858.63d0, 398600.440d0,
42828.3d0, 126686511d0, 37931207.8d0 ,
5793966d0, 6835107d0, 872.4d0/
&run_vars
outputfile01='/tmp/output.nc',
run_forward_in_time=.true., ! if true run in the forward time direction
general_relativity=.false., ! add corrections for general relativity
dt=30.d0, ! days ! output interval
tfinal=1.d4, ! years ! run-time
interval_io=10., ! years for io
! which bodies a body interacts with:
! interact(1:10)=.true.,.true.,.true.,.true.,.true.,
! .true.,.true.,.true.,.true.,.true./
! this example is everything apart from mars, jupiter, saturn, and uranus, etc
! interact(1:10)=.true.,.true.,.true.,.true.,.false.,
! .false.,.false.,.false.,.false.,.true./
! this example is only the sun
interact(1:10)=.true.,.false.,.false.,.false.,.false.,
.false.,.false.,.false.,.false.,.false./