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

Read `grid_lon_lat` and cylone from namelist

parent c95f8197
No related branches found
No related tags found
No related merge requests found
......@@ -17,18 +17,19 @@ program test_read_eddy
integer:: ishape = 0
TYPE(shpc_slice_handler) hshp_in, hshp_out
character(len = :), allocatable:: shpc_dir
namelist /main_nml/ ishape
logical:: cyclone = .false., grid_lon_lat = .true.
namelist /main_nml/ ishape, cyclone, grid_lon_lat
!-------------------------------------------------------------------------
call get_command_arg_dyn(1, shpc_dir, "Required argument: SHPC-directory")
print *, "Enter namelist main_nml."
read(unit = *, nml = main_nml)
call shpc_open(hshp_in, shpc_dir, cyclone = .false., slice = 0, &
call shpc_open(hshp_in, shpc_dir, cyclone = cyclone, slice = 0, &
pszaccess = "rb")
call read_eddy(e, k, eddy_i, hshp_in, ishape)
call shpc_create(hshp_out, shpc_dir = "SHPC", cyclone = hshp_in%cyclone, &
slice = 0, grid_lon_lat = .true.)
slice = 0, grid_lon_lat = grid_lon_lat)
CALL shpc_close(hshp_in)
call write_eddy(e, hshp_out, k, eddy_i)
write(hshp_out%unit, fmt = *) 0
......
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