#############################################################################
# This script combines the results from the Gamma-M-K-Gamma sweep
# and plots the bandstructure.
# Properties:     
# a:            period used to normalize the frequency (fs_all_norm=fs_all*a/c)
#               in this case the square lattice constant
# fs_all:       Frequencies of bands in units of Hz
# fs_all_norm:  Frequencies of bands in units of Hz * a / c
#############################################################################

runsweep; # run all sweeps

# get ax from model
a = getnamed("::model","a");
fs_all = matrix(50, 30);

# get fs data from the sweeps
sweepname="Gamma-M";
resonance=getsweepresult(sweepname,"fs");
fs_all(1:50,1:10)=resonance.fs;

sweepname="M-K";
resonance=getsweepresult(sweepname,"fs");
fs_all(1:50,11:20)=resonance.fs;

sweepname="K-Gamma";
resonance=getsweepresult(sweepname,"fs");
fs_all(1:50,21:30)=resonance.fs;

fs_all_norm = fs_all*a/c;

plot(1:30,transpose(fs_all_norm),"k (Gamma-M-K-Gamma)","f (Hz*a/c)","Bandstructure","plot points, color=blue");