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

Rename variable `shp_tr_dir` to `shpc_dir`

parent 70e88286
No related branches found
No related tags found
No related merge requests found
......@@ -25,7 +25,7 @@ program eddy_graph
integer rank, n_proc, copy, delta, j
integer:: n_dates = huge(0) ! number of dates to read
integer k1, k_begin, k_end, k, k_end_main_loop
character(len = :), allocatable:: shp_tr_dir
character(len = :), allocatable:: shpc_dir
integer:: max_delta = 1
real:: corner_deg(2) = [0.125, - 59.875]
......@@ -69,12 +69,11 @@ program eddy_graph
call mpi_init
call MPI_Comm_rank(MPI_Comm_world, rank)
call MPI_Comm_size(MPI_Comm_world, n_proc)
call get_command_arg_dyn(1, shp_tr_dir, &
"Required argument: SHP-triplet-directory")
call get_command_arg_dyn(1, shpc_dir, "Required argument: SHPC-directory")
if (rank == 0) then
call new_unit(unit)
open(unit, file = shp_tr_dir // "/grid_nml.txt", status = "old", &
open(unit, file = shpc_dir // "/grid_nml.txt", status = "old", &
action = "read", position = "rewind")
read(unit, nml = grid_nml)
close(unit)
......@@ -92,7 +91,7 @@ program eddy_graph
"eddy_graph dist_lim")
copy = merge(dist_lim, 0, periodic)
call read_column(ishape_last, &
file = trim(shp_tr_dir) // "/ishape_last.txt", last = n_dates)
file = trim(shpc_dir) // "/ishape_last.txt", last = n_dates)
n_dates = size(ishape_last)
call assert(n_dates >= max_delta + 1, &
"eddy_graph: n_dates should be >= max_delta + 1")
......@@ -114,7 +113,7 @@ program eddy_graph
if (rank /= 0) allocate(ishape_last(n_dates))
call ezmpi_bcast(ishape_last, root = 0)
call shp_tr_open(hshp, trim(shp_tr_dir), rank)
call shp_tr_open(hshp, trim(shpc_dir), rank)
if (rank == 0) call dbf_read_attribute_03(k1, hshp%extremum, hshp%extr_date, &
ishape = 0)
call ezmpi_bcast(k1, root = 0)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment