def main():
'''
thicknesses in microns (um)
'''
eps = putil.getEps('al', 1500.)
print eps
exit()
# Set layer indices
ni = 1.
nf = 4.
nt = sqrt(-100 + 1j)
ns_film = np.array([nf])
ns_bare = np.array([ni])
# Set layer thicknesses
ds = np.array([1]) # um
# Set excitation propeties
wls = np.linspace(1., 2., 100) # um
# ths = np.linspace(0.0, 90, 1000) * pi/180.
ths = np.array([45. * pi / 180])
pol = 'p'
# Collect data
R_bare = pl.zeros((len(wls), len(ths)), dtype='float')
R_film = pl.zeros((len(wls), len(ths)), dtype='float')
for ith, th in enumerate(ths):
for iwl, wl in enumerate(wls):
R_film[iwl, ith] = tm.solvestack(ni, nt, ns_film, ds, wl, pol,
th)[0]
R_bare[iwl, ith] = tm.solvestack(ni, nt, ns_bare, ds, wl, pol,
th)[0]
R = R_film / R_bare
# Plot data
pl.figure(figsize=(10, 7.5))
pl.plot(wls, R[:, 0], 'r', lw=2, label='R')
#pl.xlim(ths[0], ths[-1])
#pl.ylim(0, 1)
pl.ylabel(r'$R$', fontsize=18)
pl.xlabel(r'$\theta$', fontsize=18)
pl.show()
return
def main():
nm = 1e-9
um = 1e-6
cm = 1e-3
wls = np.linspace(1500,6000,100)
# e_agst = putil.getEps('a-gete', wls)
# e_cgst = putil.getEps('c-gete', wls)
e_agst = putil.getEps('a-gst225', wls)
e_cgst = putil.getEps('c-gst225', wls)
n_agst = sqrt(e_agst)
n_agst = n_agst.real + 1j * n_agst.imag
n_cgst = sqrt(e_cgst)
n_cgst = n_cgst.real + 1j * n_cgst.imag
e_al = putil.getEps('al',wls)
n_al = sqrt(e_al)
n_al = n_al.real + 1j * n_al.imag
Nths = 10
ths = np.linspace(10,40,Nths) * pi/180.
pol = 'p'
ds = np.array([175]) # nm
R = np.zeros((len(wls),Nths,2), dtype='float')
T = np.zeros((len(wls),Nths,2), dtype='float')
A = np.zeros((len(wls),Nths,2), dtype='float')
plt.figure(figsize=(5,4))
for phase in ['a','c']:
for ip, p in enumerate(['p','s']):
for ith, th in enumerate(ths):
for iwl, wl in enumerate(wls):
ni = 1+0j
ns_a = np.array([n_agst[iwl]])
ns_c = np.array([n_cgst[iwl]])
nt = n_al[iwl]
# reflection from aluminum with no GST
ref = tm.solvestack(ni, nt, np.array([1]), ds, wl, p, th)[0]
if phase == 'a':
R[iwl,ith,ip], T[iwl,ith,ip], A[iwl,ith,ip] = tm.solvestack(ni, nt, ns_a, ds, wl, p, th)
else:
R[iwl,ith,ip], T[iwl,ith,ip], A[iwl,ith,ip] = tm.solvestack(ni, nt, ns_c, ds, wl, p, th)
R[iwl,ith,ip] /= ref
Rav = np.sum(np.sum(R,axis=1), axis=1) / (2*Nths)
c = 'b' if phase == 'a' else 'r'
plt.plot(wls/1e3, Rav, c, lw=2, label='R')
plt.ylabel(r'$R$', fontsize=18)
plt.xlabel(r'$\lambda$ ($\mu m$)', fontsize=18)
plt.legend(('amor.','crys.'), loc='best')
plt.tight_layout()
plt.show()
