- function [rh,rv] = ref_coef (psi, epsp, epspp)
- eps = epsp - i .* epspp;
- psirad = psi.*(pi./180.);
- arg1 = eps - (cos(psirad).^2);
- arg2 = sqrt(arg1);
- arg3 = sin(psirad);
- arg4 = eps.*arg3;
- rv = (arg4-arg2)./(arg4+arg2);
- rh = (arg3-arg2)./(arg3+arg2);
- return