1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
| var HALF_PI = Math.PI/2;
| module.exports = function(eccent, ts) {
| var eccnth = 0.5 * eccent;
| var con, dphi;
| var phi = HALF_PI - 2 * Math.atan(ts);
| for (var i = 0; i <= 15; i++) {
| con = eccent * Math.sin(phi);
| dphi = HALF_PI - 2 * Math.atan(ts * (Math.pow(((1 - con) / (1 + con)), eccnth))) - phi;
| phi += dphi;
| if (Math.abs(dphi) <= 0.0000000001) {
| return phi;
| }
| }
| //console.log("phi2z has NoConvergence");
| return -9999;
| };
|
|