#############################################################################
# This script combines the results from the K-X-Gamma-L-U sweeps
# and plots the bandstructure.
# Properties:     
# a:            period used to normalize the frequency (fs_all_norm=fs_all*a/c)
#               in this case the woodpile 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 parameter sweeps
a=getnamed("::model","pile_a");  # get structure period, used to normalize frequency
fs_all = matrix(50, 80);

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

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

sweepname="Gamma-L"; 
resonance=getsweepresult(sweepname,"fs");
fs_all(1:50,41:60)=resonance.fs;

sweepname="L-U"; 
resonance=getsweepresult(sweepname,"fs");
fs_all(1:50,61:80)=resonance.fs;

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