Yokoya elliptic (reopened)

Replaced the existing look-up table method with an analytical calculation method. I reopened the merge request to change the value of qr_th from 0.9 to 0.8. 0.8 is corresponds to Migliorati's paper. (Please refer to FIG. 3 of the paper.)