ZH: display tx in period power grid

This commit is contained in:
Eric Teunis de Boone 2022-12-23 12:28:33 +01:00
parent 0d8bce9953
commit 82c10059cb

View file

@ -177,7 +177,7 @@ if __name__ == "__main__":
os.makedirs(fig_subdir, exist_ok=True) os.makedirs(fig_subdir, exist_ok=True)
# Read in antennas from file # Read in antennas from file
_, __, antennas = beacon.read_beacon_hdf5(antennas_fname) _, tx, antennas = beacon.read_beacon_hdf5(antennas_fname)
# Read original REvent # Read original REvent
ev = REvent(fname) ev = REvent(fname)
# .. patch in our antennas # .. patch in our antennas
@ -315,7 +315,12 @@ if __name__ == "__main__":
fig.colorbar(sc, ax=axs) fig.colorbar(sc, ax=axs)
if fig_dir: if fig_dir:
fig.tight_layout()
fig.savefig(path.join(fig_dir, __file__+f'.maxima.run{r}.pdf')) fig.savefig(path.join(fig_dir, __file__+f'.maxima.run{r}.pdf'))
if True:
axs.plot(tx.x/1e3, tx.y/1e3, marker='X', color='k')
fig.tight_layout()
fig.savefig(path.join(fig_dir, __file__+f'.maxima.run{r}.with_tx.pdf'))
## Save ks to file ## Save ks to file
best_idx = np.argmax(maxima_per_loc) best_idx = np.argmax(maxima_per_loc)