57 real(kind=dp) :: clat1,clat2,slat1,slat2,cdlon,crd
58 real(kind=dp),
parameter :: rerth=6.3712e6_dp
59 real(kind=dp),
parameter :: pi=3.14159265358979_dp
61 if ((abs(rlat1-rlat2).lt.0.0003).and. &
62 (abs(rlon1-rlon2).lt.0.0003))
then
67 clat1=cos(
real(rlat1,kind=dp))
68 slat1=sin(
real(rlat1,kind=dp))
69 clat2=cos(
real(rlat2,kind=dp))
70 slat2=sin(
real(rlat2,kind=dp))
71 cdlon=cos(
real(rlon1-rlon2,kind=dp))
72 crd=slat1*slat2+clat1*clat2*cdlon
real function distance2(rlat1, rlon1, rlat2, rlon2)