diff --git a/Common/shpc_open.f90 b/Common/shpc_open.f90
index b943c2d7a616d82489380a8eb5bcbfd5bb5d87de..01da636c27e857bd6e7291aad11d8bb74330c8f1 100644
--- a/Common/shpc_open.f90
+++ b/Common/shpc_open.f90
@@ -7,8 +7,8 @@ contains
subroutine shpc_open(hshp, shpc_dir, pszaccess)
! Libraries:
- USE jumble, only: new_unit, read_column
- use shapelib_03, only: shp_open_03, dbf_read_attribute_03
+ USE jumble, only: new_unit
+ use shapelib_03, only: shp_open_03
use derived_types, only: shpc
use read_field_indices_m, only: read_field_indices
@@ -36,10 +36,6 @@ contains
close(unit)
hshp%cyclone = trim(adjustl(orientation)) == "cyclones"
hshp%dir = shpc_dir
- call dbf_read_attribute_03(hshp%d0, hshp%extremum, hshp%extr_date, &
- ishape = 0)
- call read_column(hshp%ishape_last, file = shpc_dir // "/ishape_last.txt", &
- my_lbound = hshp%d0)
end subroutine shpc_open
diff --git a/Inst_eddies/Tests/test_nearby_extr.f90 b/Inst_eddies/Tests/test_nearby_extr.f90
index 705b8b019487616682049d78ebbf59a1a7e02b64..ab3754b6c1d03e1ac7ee5f06f185aaaba0d3aa24 100644
--- a/Inst_eddies/Tests/test_nearby_extr.f90
+++ b/Inst_eddies/Tests/test_nearby_extr.f90
@@ -1,7 +1,7 @@
program test_nearby_extr
! Libraries:
- use jumble, only: get_command_arg_dyn, new_unit
+ use jumble, only: get_command_arg_dyn, new_unit, read_column
use jumble, only: deg_to_rad, rad_to_deg
use shapelib_03, only: dbf_read_attribute_03
@@ -44,6 +44,9 @@ program test_nearby_extr
close(unit)
call shpc_open(hshp, trim(shpc_dir), pszaccess = "rb")
+ call dbf_read_attribute_03(hshp%d0, hshp%extremum, hshp%extr_date, ishape = 0)
+ call read_column(hshp%ishape_last, file = shpc_dir // "/ishape_last.txt", &
+ my_lbound = hshp%d0)
call read_snapshot(s, [hshp], nlon, nlat, k = hshp%d0, &
corner = corner_deg * deg_to_rad, step = step_deg * deg_to_rad, &
copy = 0)
diff --git a/Overlap/Tests/test_get_dispatch_snap.f90 b/Overlap/Tests/test_get_dispatch_snap.f90
index 6d5edb68571b7e215e259fd681d231759b95dd54..38377a9cb6f5bfbaf423a7e3f3005f124e5a4b1a 100644
--- a/Overlap/Tests/test_get_dispatch_snap.f90
+++ b/Overlap/Tests/test_get_dispatch_snap.f90
@@ -73,6 +73,9 @@ program test_get_dispatch_snap
end if
call shpc_open(hshp, trim(shpc_dir), pszaccess = "rb")
+ call dbf_read_attribute_03(hshp%d0, hshp%extremum, hshp%extr_date, ishape = 0)
+ call read_column(hshp%ishape_last, file = shpc_dir // "/ishape_last.txt", &
+ my_lbound = hshp%d0)
if (rank == 0) then
n_dates = size(hshp%ishape_last)
diff --git a/Overlap/Tests/test_overlap.f90 b/Overlap/Tests/test_overlap.f90
index c95e675bce4b3b71a55c3b7b774cca47f9abfb25..35d919a18082edce753718273f601d2052eb52a0 100644
--- a/Overlap/Tests/test_overlap.f90
+++ b/Overlap/Tests/test_overlap.f90
@@ -65,6 +65,10 @@ program test_overlap
step = step_deg * deg_to_rad
allocate(flow(max_delta + 1))
call shpc_open(hshp, trim(shpc_dir), pszaccess = "rb")
+ call dbf_read_attribute_03(hshp%d0, hshp%extremum, hshp%extr_date, &
+ ishape = 0)
+ call read_column(hshp%ishape_last, file = shpc_dir // "/ishape_last.txt", &
+ my_lbound = hshp%d0)
n_dates = size(hshp%ishape_last)
call assert(hshp%d0 <= [k_test_1, k_test_2] .and. [k_test_1, k_test_2] &
< hshp%d0 + n_dates, "test_overlap k_test_1, k_test_2")
diff --git a/Overlap/Tests/test_read_eddy.f90 b/Overlap/Tests/test_read_eddy.f90
index 7f63edced9d20e2117d3df664d6ef5ac1e5fa552..eae08f0d2b47ccb9d510ecb484a711654afd2db1 100644
--- a/Overlap/Tests/test_read_eddy.f90
+++ b/Overlap/Tests/test_read_eddy.f90
@@ -1,7 +1,8 @@
program test_read_eddy
! Library:
- use jumble, only: get_command_arg_dyn
+ use jumble, only: get_command_arg_dyn, read_column
+ use shapelib_03, only: dbf_read_attribute_03
use derived_types, only: eddy, shpc
use read_eddy_m, only: read_eddy
@@ -25,6 +26,9 @@ program test_read_eddy
print *, "Enter namelist main_nml."
read(unit = *, nml = main_nml)
call shpc_open(hshp, shpc_dir, pszaccess = "rb")
+ call dbf_read_attribute_03(hshp%d0, hshp%extremum, hshp%extr_date, ishape = 0)
+ call read_column(hshp%ishape_last, file = shpc_dir // "/ishape_last.txt", &
+ my_lbound = hshp%d0)
call read_eddy(e, k, eddy_i, hshp, ishape)
CALL shpc_close(hshp)
diff --git a/Overlap/Tests/test_read_snapshot.f90 b/Overlap/Tests/test_read_snapshot.f90
index ef7107c721e0229c9193602a2fd8fab3303127b3..f540c8e52802c82fc1f297553ab02c4da0f7d77f 100644
--- a/Overlap/Tests/test_read_snapshot.f90
+++ b/Overlap/Tests/test_read_snapshot.f90
@@ -3,7 +3,8 @@ program test_read_snapshot
use, intrinsic:: ISO_FORTRAN_ENV
! Libraries:
- use jumble, only: get_command_arg_dyn, new_unit, deg_to_rad, assert
+ use jumble, only: get_command_arg_dyn, new_unit, deg_to_rad, assert, &
+ read_column
use shapelib_03, only: dbf_read_attribute_03
use derived_types, only: snapshot, shpc
@@ -48,6 +49,10 @@ program test_read_snapshot
do i = 1, n_shpc
call get_command_arg_dyn(i, shpc_dir)
call shpc_open(hshp(i), trim(shpc_dir), pszaccess = "rb")
+ call dbf_read_attribute_03(hshp(i)%d0, hshp(i)%extremum, &
+ hshp(i)%extr_date, ishape = 0)
+ call read_column(hshp(i)%ishape_last, &
+ file = shpc_dir // "/ishape_last.txt", my_lbound = hshp(i)%d0)
end do
! Assuming grid_nml.txt is the same in all input SHPC:
diff --git a/Overlap/Tests/test_send_recv.f90 b/Overlap/Tests/test_send_recv.f90
index e097f374b8bc875ff74e7fad64f6b3f194f43007..75397cbd5693004d685a46cad4184072298b3818 100644
--- a/Overlap/Tests/test_send_recv.f90
+++ b/Overlap/Tests/test_send_recv.f90
@@ -83,6 +83,9 @@ program test_send_recv
end if
call shpc_open(hshp, trim(shpc_dir), pszaccess = "rb")
+ call dbf_read_attribute_03(hshp%d0, hshp%extremum, hshp%extr_date, ishape = 0)
+ call read_column(hshp%ishape_last, file = shpc_dir // "/ishape_last.txt", &
+ my_lbound = hshp%d0)
if (rank == 0) then
n_dates = size(hshp%ishape_last)
diff --git a/Overlap/eddy_graph.f90 b/Overlap/eddy_graph.f90
index 7e3f9bf7ce230be86c96fd34070f42ce9d50019b..76e2222c8084f976232c70400fdb8403326aefe8 100644
--- a/Overlap/eddy_graph.f90
+++ b/Overlap/eddy_graph.f90
@@ -4,9 +4,11 @@ program eddy_graph
! Libraries:
use ezmpi, only: ezmpi_bcast
- use jumble, only: get_command_arg_dyn, new_unit, ediff1d, assert, deg_to_rad
+ use jumble, only: get_command_arg_dyn, new_unit, ediff1d, assert, &
+ deg_to_rad, read_column
use mpi_f08, only: mpi_init, mpi_finalize, MPI_Comm_rank, MPI_Comm_world, &
MPI_Comm_size
+ use shapelib_03, only: dbf_read_attribute_03
use derived_types, only: snapshot, shpc
use dispatch_snapshot_m, only: dispatch_snapshot
@@ -77,6 +79,10 @@ program eddy_graph
do i = 1, n_shpc
call get_command_arg_dyn(i, shpc_dir)
call shpc_open(hshpc(i), trim(shpc_dir), pszaccess = "rb")
+ call dbf_read_attribute_03(hshpc(i)%d0, hshpc(i)%extremum, &
+ hshpc(i)%extr_date, ishape = 0)
+ call read_column(hshpc(i)%ishape_last, &
+ file = shpc_dir // "/ishape_last.txt", my_lbound = hshpc(i)%d0)
end do
if (rank == 0) then