351     double lati, latj, longi, longj;  in geographic_edgelen()  local
361     latj = GH_PI * (deg + 5.0 * min / 3.0) / 180.0;  in geographic_edgelen()
371     q2 = cos (lati - latj);  in geographic_edgelen()
372     q3 = cos (lati + latj);  in geographic_edgelen()
380     double lati, latj, longi, longj;  in geom_edgelen()  local
384     latj = M_PI * dat->x[j] / 180.0;  in geom_edgelen()
389     q1 = cos (latj) * sin(longi - longj);  in geom_edgelen()
392     q2 = sin(lati + latj) * q3 * q3 - sin(lati - latj) * q4 * q4;  in geom_edgelen()
393     q5 = cos(lati - latj) * q4 * q4 - cos(lati + latj) * q3 * q3;  in geom_edgelen()
400     double lati, latj, longi, longj;
405     latj = M_PI * (dat->x[j] / 180.0);
411     q2 = cos (lati - latj);
412     q3 = cos (lati + latj);