Phi_E = pi/2 # 1st Euler angle
Theta_E_list = [0, pi/4, pi/2] # 2nd Eulet angle
# Incidence angles
Phi_i_list = numpy.linspace(0, 0.999*pi/2, 37) # array of Phi_i values
Phi_i_deg = Phi_i_list*180/pi
Kx_list = front.get_Kx_from_Phi(Phi_i_list) # array of Kx values
data = []
for Theta_E in Theta_E_list:
R = Berreman4x4.rotation_Euler((Phi_E, Theta_E, 0))
back.setMaterial(uniaxialMaterialRef.rotated(R))
l = []
for Kx in Kx_list:
Jr = s.getJones(Kx)[0]
l.append(Berreman4x4.extractEllipsoParam(Jr))
data.append(l)
data = numpy.array(data)
fig = pyplot.figure()
ax = fig.add_subplot(111)
ax.plot(Phi_i_deg, data[:,:,0,0].T)
ax.set_xlabel(r"$\theta_i$")
ax.set_ylabel(r"$\Psi_{pp}$")
pyplot.tight_layout()
##############################################################################
# We reproduce figure 6.17
print("We reproduce figure 6.17, p. 239...")
