Skip to content
Snippets Groups Projects
Commit 62d0f068 authored by Lionel GUEZ's avatar Lionel GUEZ
Browse files

Work with date index in tests

Follow-up to commit 3b68dd5.
parent 8f858962
No related branches found
No related tags found
No related merge requests found
......@@ -52,7 +52,7 @@ program test_nearby_extr
call shpc_open(hshp, trim(shpc_dir), rank = 0)
call dbf_read_attribute_03(d_init, hshp%extremum, hshp%extr_date, ishape = 0)
call read_snapshot(s, hshp, nlon, nlat, d_init, k = d_init, &
call read_snapshot(s, hshp, nlon, nlat, d_init, k = 0, &
corner = corner_deg * deg_to_rad, step = step_deg * deg_to_rad, &
copy = 0, ishape_last = [ishape_last])
CALL shpc_close(hshp)
......
......@@ -81,10 +81,10 @@ program test_overlap
call read_column(ishape_last, file = trim(shpc_dir) // "/ishape_last.txt")
call shpc_open(hshp, trim(shpc_dir), rank = 0)
call dbf_read_attribute_03(d_init, hshp%extremum, hshp%extr_date, ishape = 0)
call read_snapshot(flow(1), hshp, nlon, nlat, d_init, k_test_1, corner, &
step, copy, ishape_last)
call read_snapshot(flow(max_delta + 1), hshp, nlon, nlat, d_init, k_test_2, &
call read_snapshot(flow(1), hshp, nlon, nlat, d_init, k_test_1 - d_init, &
corner, step, copy, ishape_last)
call read_snapshot(flow(max_delta + 1), hshp, nlon, nlat, d_init, &
k_test_2 - d_init, corner, step, copy, ishape_last)
CALL shpc_close(hshp)
print *, "Enter flow(1)%list_vis%delta_out (array with ", &
flow(1)%number_vis_extr, " values):"
......
......@@ -70,7 +70,7 @@ program test_read_snapshot
call shpc_open(hshp, trim(shpc_dir), rank = 0)
call dbf_read_attribute_03(d_init, hshp%extremum, hshp%extr_date, ishape = 0)
call read_snapshot(s, hshp, nlon, nlat, d_init, k = d_init, &
call read_snapshot(s, hshp, nlon, nlat, d_init, k = 0, &
corner = corner_deg * deg_to_rad, step = step_deg * deg_to_rad, &
copy = copy, ishape_last = [ishape_last])
CALL shpc_close(hshp)
......
......@@ -102,19 +102,19 @@ program test_send_recv
if (rank == 0) call dbf_read_attribute_03(d_init, hshp%extremum, &
hshp%extr_date, ishape = 0)
call ezmpi_bcast(d_init, root = 0)
if (rank == 1) call read_snapshot(s, hshp, nlon, nlat, d_init, k = d_init, &
if (rank == 1) call read_snapshot(s, hshp, nlon, nlat, d_init, k = 0, &
corner = corner_deg * deg_to_rad, step = step_deg * deg_to_rad, &
copy = copy, ishape_last = ishape_last)
CALL shpc_close(hshp)
if (rank == 1) then
call send_snapshot(s, dest = 0, tag = d_init)
call send_snapshot(s, dest = 0, tag = 0)
else
! rank == 0
call MPI_Comm_get_attr(MPI_Comm_world, MPI_TAG_UB, attribute_val, flag)
call assert(flag, "test_send_recv MPI_Comm_get_attr MPI_TAG_UB")
print *, "MPI_TAG_UB = ", attribute_val
call recv_snapshot(s, nlon, nlat, copy, source = 1, tag = d_init)
call recv_snapshot(s, nlon, nlat, copy, source = 1, tag = 0)
call write_snapshot(s, corner_deg, step_deg, nlon, nlat, copy, d_init)
end if
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment