k0 intop1(acpr.k) "Wave number"
kx_e sin(theta0) "Plane wave direction, x-component"
ky_e -cos(theta0) "Plane wave direction, y-component"
kx k0*kx_e "Incident wave number, x-component"
ky k0*ky_e "Incident wave number, y-component"
Pin intop2(-acpr.I_by) "Projected incident power"
Pout intop2(acpr.I_sy) "Projected reflected power"
alpha 1-Pout/Pin "Absorption coefficient"
Zn aveop1(acpr.p_t/(-acpr.vy))/intop1(acpr.Z) "Surface impedance (normalized)"