Read index of target eddy

Inst_eddies/Tests/short_tests.json

@@ -431,6 +431,7 @@
             "$build_dir/Inst_eddies/Tests/test_nearby_extr",
 	    "$PWD/Inst_eddies_loop/SHPC"
 	],
+	"input": "&main_nml /\n",
 	"description": "We are using Inst_eddies_loop, which is a previous test."
     },
     {

Inst_eddies/Tests/test_nearby_extr.f90

@@ -22,9 +22,14 @@ program test_nearby_extr
   real, allocatable:: nearby(:, :) ! (2, :) longitude and
   ! latitude, in rad, of extrema near the target extremum   
 
+  integer:: i = 4
+  namelist /main_nml/ i
+
   !-------------------------------------------------------------------------
 
   call get_command_arg_dyn(1, shpc_dir, "Required argument: SHPC-directory")
+  print *, "Enter namelist main_nml."
+  read(unit = *, nml = main_nml)
   call read_grid(rank = 0, shpc_dir = shpc_dir)
   call shpc_open(hshp, trim(shpc_dir), cyclone = .true., slice = 0, &
        pszaccess = "rb")
@@ -33,7 +38,7 @@ program test_nearby_extr
   call read_snapshot(s, [hshp], [ssm], k = ssm%d0, copy = 0)
   CALL shpc_close(hshp)
 
-  nearby = nearby_extr(s%extr_map, s%list, i = 4)
+  nearby = nearby_extr(s%extr_map, s%list, i)
   print *, "nearby:"
 
   do l = 1, size(nearby, 2)