Orbital elements

analysis.rmd

@@ -27,6 +27,14 @@ ox <- baseg + geom_line(aes(y=oxidizer)) + ylab("Oxidizer (litres)") + ggtitle("
 at <- baseg + geom_line(aes(y=available_thrust)) + ylab("Available thrust (Newtons)") + ggtitle("Available thrust")
 thrust <- baseg + geom_line(aes(y=current_thrust)) + ylab("Current thrust (Newtons)") + ggtitle("Current thrust")
 g_force <- baseg + geom_line(aes(y=g_force)) + ylab("G-force (m/s²)") + ggtitle("G-force")
+ap <- baseg + geom_line(aes(y=apoapsis)) + ylab("Apoapsis height (m)") + ggtitle("Apoapsis")
+pa <- baseg + geom_line(aes(y=periapsis)) + ylab("Periapsis height (m)") + ggtitle("Periapsis")
+orad <- baseg + geom_line(aes(y=orbital_radius)) + ylab("Orbital radius (m)") + ggtitle("Orbital radius")
+os <- baseg + geom_line(aes(y=orbit_speed)) + ylab("Orbital speed (m/s)") + ggtitle("Orbital speed")
+op <- baseg + geom_line(aes(y=orbit_period)) + ylab("Orbital period (s)") + ggtitle("Orbital period")
+
+
+
 ```
 
 ```{r drawgraphs, echo=FALSE}
@@ -40,4 +48,9 @@ ox
 at
 thrust
 g_force
+ap
+pa
+orad
+os
+op
 ```

monitor.py

@@ -42,7 +42,7 @@ button.rect_transform.position = (0, 0)
 button_clicked = conn.add_stream(getattr, button, 'clicked')
 
 f = open(datafile, 'w')
-f.write('mission_time, current_altitude, v_speed, h_speed, mass, electric_charge, liquid_fuel, oxidizer, available_thrust, current_thrust, g_force\n')
+f.write('mission_time, current_altitude, v_speed, h_speed, mass, electric_charge, liquid_fuel, oxidizer, available_thrust, current_thrust, g_force, apoapsis, periapsis, orbital_radius, orbit_speed, orbit_period\n')
 
 vessel = conn.space_center.active_vessel
 body = vessel.orbit.body
@@ -61,7 +61,13 @@ while button_clicked() == False:
     ox = vessel.resources.amount('Oxidizer')
     g_force = vessel.flight(r_frame).g_force
     timept = vessel.met
-    f.write('{0}, {1}, {2}, {3}, {4}, {5}, {6}, {7}, {8}, {9}, {10}\n'.format(timept, current_altitude, v_speed, h_speed, vm, e_charge, lf, ox, at, c_thrust, g_force))
+    vo = vessel.orbit
+    ap = vo.apoapsis
+    pa = vo.periapsis
+    orad = vo.radius
+    os = vo.speed
+    op = vo.period
+    f.write('{0}, {1}, {2}, {3}, {4}, {5}, {6}, {7}, {8}, {9}, {10}, {11}, {12}, {13}, {14}, {15}\n'.format(timept, current_altitude, v_speed, h_speed, vm, e_charge, lf, ox, at, c_thrust, g_force, ap, pa, orad, os, op)) 
     time.sleep(0.1)
 
 f.close()