It is not clear what you expect as result or why you are using the temp array, the constant 8888, please clarify.   ... View more

Hello @scad,   I agree with Reeza that you should use the GEODIST function if your coordinates are latitudes and longitudes (as the variable names suggest). It looks like you interchanged latitudes and longitudes, though: Latitudes must be between -90 (South Pole) and 90 (North Pole). After this correction all your locations would be in South Korea. Does this sound plausible? If so, you can use the code provided by Mike Zdeb in Using Proc Distance to Calculate Geodesic Distance (the title refers to the OP's question in that old thread) to compute the desired distances. Omit the last argument 'm' in the GEODIST function call if you want distances in kilometres rather than miles.   The DISTANCE procedure does not support geodesic distance, so the results you would obtain (with PROC DISTANCE options method=euclid nostd and VAR statement var interval(longti lati)), i.e., the "distances" 1.50073, 1.68964, etc., would be questionable. ... View more

SAS has a function to calculate Euclidian distance   distance = EUCLID(dist1, dist2);   https://documentation.sas.com/?docsetId=lefunctionsref&docsetTarget=p01uvzvdpw514in18j8xy1hwsd8d.htm&docsetVersion=9.4&locale=en  ... View more

Look at function geodist   https://documentation.sas.com/?docsetId=lefunctionsref&docsetTarget=n1korpfg2e18lon1nwpow9qijdxe.htm&docsetVersion=9.4&locale=en  ... View more