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

Read `grid_lon_lat` from `main_nml`

As `grid_lon_lat` may be false, we define `e%extr%coord_proj`.
parent 06ce5715
No related branches found
No related tags found
No related merge requests found
......@@ -51,7 +51,10 @@ program test_get_1_outerm
! latitude, in rad, of all the significant extrema, except the
! target extremum
namelist /main_nml/ ind_targ_extr, innermost_level, cyclone, date
logical:: grid_lon_lat = .true.
namelist /main_nml/ ind_targ_extr, innermost_level, cyclone, date, &
grid_lon_lat
!----------------------------------------------------------------
......@@ -64,6 +67,7 @@ program test_get_1_outerm
call input_ssh(corner, step, nlon, nlat, periodic, ssh, u, v)
call csvread("outside_points.csv", outside_points, first_r = 3)
outside_points = transpose(outside_points) * deg_to_rad
e%extr%coord_proj = ind_targ_extr
e%extr%coord = corner + (ind_targ_extr - 1) * step
e%out_cont = get_1_outerm(cyclone, e%extr%coord, innermost_level, &
outside_points, ssh, corner, step)
......@@ -76,7 +80,7 @@ program test_get_1_outerm
e%radius4 = - 1
e%valid = .true.
call shpc_create(hshp, shpc_dir = "SHPC", cyclone = cyclone, slice = 0, &
grid_lon_lat = .true.)
grid_lon_lat = grid_lon_lat)
call write_eddy(e, hshp, date, i = 1)
write(hshp%unit, fmt = *) 0
CALL shpc_close(hshp)
......
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