Polish.

successive_overlap.f

@@ -25,6 +25,7 @@ contains
     integer, intent(in):: j, k
 
     ! Local:
+    
     integer i1, i1_lon, i1_lat, m, l, i2
     type(polyline) polyline_1, polyline_2
     type(polygon) res_pol

weight.f

@@ -29,13 +29,13 @@ contains
     ! a quiet NaN, and then the result is a quiet NaN.
 
     use derived_types, only: eddy
-    use nr_util, only: pi
+    use nr_util, only: pi, twopi
 
     type(eddy), intent(in):: e
 
     ! Local:
 
-    real, parameter:: omega = 2 * pi / 86164.
+    real, parameter:: omega = twopi / 86164.
     ! angular speed of rotation of the Earth, in rad s-1
 
     real radius ! of the disk of same area