#############################################################################
# This script combines the results from the Gamma-X sweep 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","ax");
fs_all = matrix(50, 20);

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

fs_all_norm = fs_all*a/c;


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