legacy fix for haversine formula

Author: jjokella

src/PDAFomi_obs_l.F90

@@ -3103,10 +3103,7 @@ SUBROUTINE PDAFomi_check_dist2(thisobs, thisobs_l, coordsA, coordsB, distance2,
                 distflag = .FALSE.
              ELSE
                 ! full squared distance
-                distance2 = 0.0
-                DO k = 1, thisobs%ncoord
-                   distance2 = distance2 + dists(k)*dists(k)
-                END DO
+                distance2 = dists(2)*dists(2)
              END IF
           END IF
        END IF
@@ -3496,10 +3493,7 @@ SUBROUTINE PDAFomi_check_dist2_noniso(thisobs, thisobs_l, coordsA, coordsB, dist
                 distflag = .FALSE.
              ELSE
                 ! full squared distance
-                distance2 = 0.0
-                DO k = 1, thisobs%ncoord
-                   distance2 = distance2 + dists(k)*dists(k)
-                END DO
+                distance2 = dists(2)*dists(2)
              END IF
           END IF
        END IF