def example1():
#generate cap window
lmax = 20
nwins = 20
theta = 25.
cap = SHWindow.from_cap(lmax,nwins,theta)
cap.info()
cap.plot(20,show=False,fname='cap_tapers.png')
cap.plot_couplingmatrix(30,5,show=False,fname='cap_coupling.png')
def example2():
#generate cap window
lmax = 15
nwins = 15
topo = np.loadtxt('topo.dat.gz')
dh_mask = topo > 0.
print(dh_mask.shape)
region = SHWindow.from_mask(lmax, nwins, dh_mask, sampling=2)
region.info()
region.plot(nwins,show=False,fname='continent_tapers.png')
region.plot_couplingmatrix(30,5,show=False,fname='continent_coupling.png')
