jd-2*180*asin(sqrt((pow(sin(r/(6378.138*2)),2)-pow((wd*pi()/180-wd*pi()/180)/2,2))/cos(wd*pi()/180)/cos(wd*pi()/180)))/pi()=jd1

def jd(r: Double, jd: Double): Double = {
  val pi = 3.141592653589793
  val radLat1 = r/180*pi
  val radLat2 = jd/180*pi
  val a = radLat1 - radLat2
  val b = (r/180-jd/180)/2
  val s = 2*Asin(Sqrt(Pow(Sin(a/2), 2) + Cos(radLat1)*Cos(radLat2)*Pow(Sin(b), 2)))
  s*6378.138
}