k0 intop1(acpr.k) "Wave number of incident field"
kx_e sin(theta0) "Plane wave direction, x"
ky_e -cos(theta0) "Plane wave direction, y"
k0x k0*kx_e "Wave vector x-component of incident field"
k0y k0*ky_e "Wave vector y-component of incident field"
p_in acpr.p_b "Incident pressure field"
p_re acpr.p_s "Reflected pressure field"
R p_re/p_in "Reflection coefficient"
absR abs(aveop1(R)) "Absolute and averaged value of the reflection coefficient"
alpha_R 1-absR^2 "Absorption coefficient (based on R)"
Pin intop2(-acpr.I_by) "Incident power"
Pre intop2(acpr.I_sy) "Reflected power"
alpha 1-Pre/Pin "Absorption coefficient"