ACCIS/constellation.py at main · SIOSlab/ACCIS · GitHub

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
50
51
52
53
54
55
56
import math

import build.accis as accis

from accis_io import *
from accis_run import *

sat1 = accis_get_sat('sat')

p1 = sat1.get_param()

a = p1["Orbit Semi-Major Axis (km)"]

mu = 3.986E5;

mean_mot = math.degrees(math.sqrt(mu / a**3))

dt = 300

df = -dt * mean_mot
dW =  dt * 360 / 86400

n_sat = 4

sats = []

for j in range(n_sat):

    sat = accis_get_sat('sat')

    p = sat.get_param()

    p["Orbit True Anomaly at Epoch (deg)"] = j * df

    p["Orbit RAAN (deg)"] = j * dW

    p["Satellite ID"] = j + 1

    p["Random Number Generator Seed"] = 11 * j + 14

    sat.set_param(p)

    sats.append(sat)

t_mins = 100

t_save = 10

steps_per_min = 120

steps = t_mins * steps_per_min
save_cadence = t_save * steps_per_min

accis_run(sats, steps, save_cadence)

accis_save_results(sats)