Polish

Inst_eddies/inst_eddies.f90

@@ -89,6 +89,7 @@ program inst_eddies
      call shpc_open(hshpc_anti, shpc_dir = "SHPC", cyclone = .false., &
           slice = slice, pszaccess = "rb+")
   else
+     print *, "inst_eddies: Did not find slice, will create it..."
      call shpc_create(hshpc_cyclo, shpc_dir = "SHPC", cyclone = .true., &
           slice = slice)
      call shpc_create(hshpc_anti, shpc_dir = "SHPC", cyclone = .false., &

Overlap/eddy_graph.f90

@@ -99,6 +99,17 @@ program eddy_graph
      call assert(max_delta >= 1, "eddy_graph max_delta")
      call assert(n_dates == huge(0) .or. n_slices == 1, "If we have several " &
           // "slices then we read all the dates from them")
+
+     ! As we are requiring the grid spacing to be uniform, the value of
+     ! "periodic" may be deduced from the values of step_deg(1) and nlon:
+     periodic = nint(360. / step_deg(1)) == nlon
+     print *, "periodic = ", periodic
+     if (periodic) call assert(2 * dist_lim * step_deg(1) < 180., &
+          "eddy_graph dist_lim")
+     copy = merge(dist_lim, 0, periodic)
+
+     corner = corner_deg * deg_to_rad
+     step = step_deg * deg_to_rad
   end if
 
   call ezmpi_bcast(n_slices, root = 0)
@@ -114,14 +125,6 @@ program eddy_graph
   end do
 
   if (rank == 0) then
-     ! As we are requiring the grid spacing to be uniform, the value of
-     ! "periodic" may be deduced from the values of step_deg(1) and nlon:
-     periodic = nint(360. / step_deg(1)) == nlon
-     print *, "periodic = ", periodic
-     if (periodic) call assert(2 * dist_lim * step_deg(1) < 180., &
-          "eddy_graph dist_lim")
-     copy = merge(dist_lim, 0, periodic)
-
      if (n_dates == huge(0)) &
           n_dates = sum([(size(ssm(i)%ishape_last), i = 1, n_slices)])
      print *, "n_dates = ", n_dates
@@ -134,8 +137,6 @@ program eddy_graph
      open(unit, file = "e_overestim.txt", status = "replace", action = "write")
      write(unit, fmt = *) e_overestim
      close(unit)
-     corner = corner_deg * deg_to_rad
-     step = step_deg * deg_to_rad
   end if
 
   call ezmpi_bcast(max_delta, root = 0)