k0 intop1(elw.k) Wave number of incident field k0x k0*sin(theta0) Wave vector x-component of incident field k0y -k0*cos(theta0) Wave vector y-component of incident field Pin 1[Pa]*exp(-i*(k0x*x+k0y*y)) Incident pressure field R (elw.p_t-Pin)/Pin Reflection coefficient absR aveop1(abs(R)) Absolute value of reflection coefficient alpha 1-absR^2 Absorption coefficient