-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathfof2_median_plot.py
39 lines (30 loc) · 1.16 KB
/
fof2_median_plot.py
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
import matplotlib.pyplot as plt
import matplotlib.dates as mdates
import numpy as np
import pandas as pd
from pandas.plotting import register_matplotlib_converters
register_matplotlib_converters()
import datetime
#initialize variables
y1 = pd.read_excel('data/fof2_subplot.xlsx', delimiter=',', header=None)
y2 = pd.read_excel('data/fof2_median.xlsx', delimiter=',', header=None)
customdate = datetime.datetime(2015, 3, 15, 0) #date start
x = [customdate + datetime.timedelta(hours=i) for i in range(len(y1))]
#plot config
fig, host = plt.subplots(figsize=(9, 4))
plt.grid(linestyle='-')
p1, = host.plot(x, y1, "blue", label="foF2")
p2, = host.plot(x, y2, "lime", label="Median foF2")
host.set_xlabel("Maret 2015")
host.set_ylabel("foF2 $(Mhz)$")
host.yaxis.label.set_color('black')
host.yaxis.label.set_color('black')
tkw = dict(size=3, width=1.5)
host.tick_params(axis='y', colors='black', **tkw)
host.tick_params(axis='y', colors='black', **tkw)
host.tick_params(axis='x', **tkw)
lines = [p1, p2]
host.legend(lines, [l.get_label() for l in lines], loc="lower right")
formatter = mdates.DateFormatter('%d')
plt.gcf().axes[0].xaxis.set_major_formatter(formatter)
plt.show()