{
const lambda0 = -78.73677825927734;
const phi0 = 42.948787689208984;
const lon = -78.74;
const lat = 42.95;
let phi = lat*(Math.PI/180);
let lambda = lon*(Math.PI/180);
let centerLat = phi0 * (Math.PI/180);
let centerLon = lambda0 * (Math.PI/180);
let rhsf = Math.fround(Math.sin(Math.fround(centerLat))*Math.sin(Math.fround(phi))) + Math.fround(Math.cos(Math.fround(centerLat))*Math.cos(Math.fround(phi))*Math.cos(Math.fround(lambda-centerLon)));
console.log(rhsf);
let c = Math.acos(rhsf);
console.log("c",c)
let kprime = c
console.log("kprim",kprime);
let xout = kprime*Math.cos(phi)*Math.sin(lambda-centerLon);
let yout = kprime*(Math.cos(centerLat)*Math.sin(phi)-Math.sin(centerLat)*Math.cos(phi)*Math.cos(lambda-centerLon));
console.log(xout, yout);
console.log(Math.atan2(yout, xout))
}