Create test for `interpolate_eddy`

Create main program `test_interpolate_eddy`. Add corresponding test in
`short_tests.json`. Add corresponding input files.

Tests/CMakeLists.txt

@@ -207,3 +207,16 @@ target_include_directories(test_spher_polyline_area PRIVATE
   ${contour_531_INCLUDE_DIR})
 
 target_link_libraries(test_spher_polyline_area ${geometry_LIBRARY})
+
+# test_interpolate_eddy
+
+add_executable(test_interpolate_eddy derived_types.f90 init_shapefiles.f90
+  read_eddy.f90 write_eddy.f90 read_field_indices.f90 interpolate_eddy.f90
+  ${CMAKE_CURRENT_LIST_DIR}/test_interpolate_eddy.f90)
+
+target_include_directories(test_interpolate_eddy PRIVATE
+  ${contour_531_INCLUDE_DIR} ${fortrangis_INCLUDE_DIR})
+
+target_link_libraries(test_interpolate_eddy ${contour_531_LIBRARY}
+  ${shapelib_03_LIBRARY} ${fortrangis_LIBRARY} ${fortranc_LIBRARY}
+  ${shapelib_LIBRARY} ${GPC_F_LIBRARY} ${GPC_LIBRARY})

Tests/Input/Snapshot_region_1/mkdir_run_stdout.txt

++ mkdir Snapshot
++ /home/guez/build/Ifort/Detection_eddies_debug/extraction_eddies
+ &MAIN_NML
+ MIN_AMP = 1.0000000E-03  ,
+ MAX_RADIUS      = 2*20,
+ MIN_RADIUS      = 25.00000       
+ /
+ Enter namelist main_nml.
+ &MAIN_NML
+ MIN_AMP = 1.0000000E-03  ,
+ MAX_RADIUS      = 2*20,
+ MIN_RADIUS      = 25.00000       
+ /
+ periodic =  F
+ Number of extrema: 8
+ Created shapefile "Snapshot/extremum".
+ Created shapefile "Snapshot/outermost_contour".
+ Created shapefile "Snapshot/max_speed_contour".

Tests/short_tests.json

@@ -318,5 +318,12 @@
 	"title" : "Successive_overlap_global",
 	"description": "Global grid, normal 0.25° resolution.",
 	"stdin_filename": "$src_dir/Tests/Input/successive_overlap_global_nml.txt"
+    },
+    {
+	"args": ["$src_dir/mkdir_run.sh", "$build_dir/test_interpolate_eddy"],
+	"required":
+	[["$src_dir/Tests/Input/Snapshot_region_1", "Snapshot_old"]],
+	"title": "Interpolate_eddy",
+	"input" : "&main_nml ISHAPE_2 = 2, k2 = 20457, k = 20456/\n"
     }
 ]

Tests/test_interpolate_eddy.f90

+program test_interpolate_eddy
+
+  use, intrinsic:: ISO_FORTRAN_ENV
+
+  ! Libraries:
+  use shapelib, only: shpfileobject, shpclose
+  use shapelib_03, only: shp_open_03
+
+  use derived_types, only: eddy
+  use init_shapefiles_m, only: init_shapefiles
+  use interpolate_eddy_m, only: interpolate_eddy
+  use read_eddy_m, only: read_eddy
+  use read_field_indices_m, only: read_field_indices
+  use write_eddy_m, only: write_eddy
+
+  implicit none
+
+  type(eddy) e1, e2
+  integer k1, i
+  integer:: k2 = 90001, k = 90000, ishape_1 = 0, ishape_2 = 1
+  TYPE(shpfileobject) hshp_extremum ! shapefile extremum
+  TYPE(shpfileobject) hshp_outermost ! shapefile outermost_contour
+  TYPE(shpfileobject) hshp_max_speed ! shapefile max_speed_contour
+  namelist /main_nml/ ishape_1, ishape_2, k2, k
+
+  !-------------------------------------------------------------------------
+
+  write(unit = error_unit, nml = main_nml)
+  write(unit = error_unit, fmt = *) "Enter namelist main_nml."
+  read(unit = *, nml = main_nml)
+  write(unit = *, nml = main_nml)
+  
+  call shp_open_03(hshp_extremum, "Snapshot_old/extremum", pszaccess = "rb")
+  call shp_open_03(hshp_outermost, "Snapshot_old/outermost_contour", &
+       pszaccess = "rb")
+  call shp_open_03(hshp_max_speed, "Snapshot_old/max_speed_contour", &
+       pszaccess = "rb")
+  call read_field_indices(hshp_extremum, hshp_outermost, hshp_max_speed)
+  call read_eddy(e1, k1, i, hshp_extremum, hshp_outermost, hshp_max_speed, &
+       ishape_1)
+  call read_eddy(e2, k1, i, hshp_extremum, hshp_outermost, hshp_max_speed, &
+       ishape_2)
+  CALL shpclose(hshp_extremum)
+  CALL shpclose(hshp_outermost)
+  CALL shpclose(hshp_max_speed)
+
+  call init_shapefiles(hshp_extremum, hshp_outermost, hshp_max_speed)
+
+  call write_eddy(interpolate_eddy(e1, e2, k1, k2, k), hshp_extremum, &
+       hshp_outermost, hshp_max_speed, k, i = 1)
+  CALL shpclose(hshp_extremum)
+  print *, 'Created shapefile "extremum".'
+  CALL shpclose(hshp_outermost)
+  print *, 'Created shapefile "outermost_contour".'
+  CALL shpclose(hshp_max_speed)
+  print *, 'Created shapefile "max_speed_contour".'
+
+end program test_interpolate_eddy