1 lat1,lon1=70,10 
   2 lat2,lon2=75,15
   3 xpc1,ypc1=mapll(lat1,lon1,1) 
   4 xpc2,ypc2=mapll(lat2,lon2,1)
   5 
   6 shiftx=(xpc1-xpc2) #in [km]
   7 shifty=(ypc1-ypc2) #in [km]
   8 
   9 A=sqrt(shiftx**2+shifty**2)