Write `outside_points` to CSV file

Instead of standard output.

Inst_eddies/Tests/examine_eddy.f90

@@ -3,7 +3,7 @@ program examine_eddy
   ! Libraries:
   use contour_531, only: convert_to_ind
   use jumble, only: get_command_arg_dyn, read_column, deg_to_rad, pi, &
-       rad_to_deg, twopi
+       rad_to_deg, twopi, new_unit
   use shapelib, only: shpfileobject, shpclose
   use shapelib_03, only: dbf_read_attribute_03, dbf_get_field_index_03, &
        shp_read_point, shp_open_03
@@ -22,7 +22,7 @@ program examine_eddy
   implicit none
 
   integer:: date = huge(0), eddy_i_target = huge(0)
-  integer ishape_first, d0, ishape, i, l, n1, n2, number_extr
+  integer ishape_first, d0, ishape, i, l, n1, n2, number_extr, unit
   logical:: cyclone = .true. ! orientation of target eddy
   namelist /main_nml/ cyclone, date, eddy_i_target
   TYPE(shpfileobject) hshp
@@ -179,14 +179,18 @@ program examine_eddy
   print *, "corner_window (in degrees):", corner_window * rad_to_deg
   print *, "n1 = ", n1
   print *, "n2 = ", n2
-  print *, "outside_points (in degrees): "
+
+  call new_unit(unit)
+  open(unit, file = "outside_points.csv", status = "replace", action = "write")
 
   do i = 1, n1 + n2
-     print *, outside_points(:, i) * rad_to_deg
+     write(unit, fmt = *) outside_points(:, i) * rad_to_deg
   end do
 
+  close(unit)
   print *, "e%innermost_level = ", e%innermost_level
   print *, "e%extr%coord_proj = ", e%extr%coord_proj
-  print *, "Created shapefile cont_list and shapefile collection SHPC."
+  print *, "Created outside_points.csv, shapefile cont_list and shapefile ", &
+       "collection SHPC."
 
 end program examine_eddy