Use `deg_to_rad` from `NR_util`

Inst_eddies/Tests/test_mean_speed.f90

@@ -7,7 +7,7 @@ program test_mean_speed
   use netcdf, only: nf90_nowrite
   use netcdf95, only: nf95_open, nf95_close, nf95_inq_varid, nf95_get_var, &
        find_coord, nf95_inquire_dimension
-  use nr_util, only: pi
+  use nr_util, only: deg_to_rad
   use shapelib, only: shpclose, shpfileobject, shpobject, shpdestroyobject
   use shapelib_03, only: shp_open_03, shp_read_object_03
   
@@ -22,7 +22,6 @@ program test_mean_speed
   real, allocatable:: u(:, :), v(:, :) ! (nlon, nlat) wind, in m s-1
   TYPE(shpfileobject) hshp
   TYPE(shpobject) psobject
-  real, parameter:: deg_over_rad = pi / 180.
   type(polygon) p
 
   namelist /main_nml/ center
@@ -59,9 +58,9 @@ program test_mean_speed
   CALL shpclose(hshp)
   p = shpobj2pol(psobject)
   call shpdestroyobject(psobject)
-  p%part(1)%points = p%part(1)%points * deg_over_rad
-  m = mean_speed(u, v, p%part(1), center * deg_over_rad, &
-       corner * deg_over_rad, step = [0.25, 0.25]  * deg_over_rad)
+  p%part(1)%points = p%part(1)%points * deg_to_rad
+  m = mean_speed(u, v, p%part(1), center * deg_to_rad, &
+       corner * deg_to_rad, step = [0.25, 0.25]  * deg_to_rad)
   print *, "mean azimuthal speed = ", m, "m s-1"
 
 end program test_mean_speed