diff --git a/makedata.py b/makedata.py index 1a46ba5028a74070bd671729f2f0f74546b61472..aea1fe8a1fe4056c567db4826de462afe9931657 100755 --- a/makedata.py +++ b/makedata.py @@ -7,10 +7,12 @@ import numpy as np a = 1 -with open("../sagittarius_a/results/power_1420.ecsv", "r") as f: +with open("doppler.txt", "r") as f: for line in f.readlines(): + a += 1 if line[0] == '#' or line[0] == 't': print(line, end='', flush=True) else: print(line, end='', flush=True) - sleep(1) + if a > 1000: + sleep(1) diff --git a/power_plot.py b/power_plot.py index 859ef0dfa2770bb4926970ee0bbaeab97e164f1c..200315ba4e73ce5baccde3be56e31c654c86b33d 100755 --- a/power_plot.py +++ b/power_plot.py @@ -1,6 +1,6 @@ #!/usr/bin/env python3 -from datetime import datetime +from datetime import datetime, timedelta import sys import os @@ -51,7 +51,7 @@ class PowerPlot: ) fig.suptitle(f"Maximum frequency", fontsize=16) fig.tight_layout() - (self.totalplot,) = self.ax.plot([], [], "-", color="blue") + (self.totalplot,) = self.ax.plot([], [], ".", color="blue") (self.totalplothighlight,) = self.ax.plot( [], [], ".", color="red", markersize=10 ) @@ -117,7 +117,8 @@ class PowerPlotMainWindow(QMainWindow): line = sys.stdin.readline() values = line.split(",") - self.time_dt = np.array(date2num(datetime.utcfromtimestamp(float(values[0])))) + first_time = datetime.utcfromtimestamp(float(values[0])) + self.time_dt = np.array(date2num(first_time)) self.last_x_buffer = np.empty(last_seconds + 1) self.last_x_buffer[:] = np.nan @@ -126,8 +127,8 @@ class PowerPlotMainWindow(QMainWindow): self.data = np.array([float(values[1])]) - self.start_time = datetime(2023, 12, 13, 9, 15) - self.stop_time = datetime(2023, 12, 13, 11, 30) + self.start_time = first_time + self.stop_time = first_time + timedelta(hours=2) plt.ion()