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

Allow missing `outside_points.csv`

parent 7902a18d
No related branches found
No related tags found
No related merge requests found
...@@ -59,6 +59,7 @@ program test_get_1_outerm ...@@ -59,6 +59,7 @@ program test_get_1_outerm
TYPE(shpfileobject) hshp_cont_list TYPE(shpfileobject) hshp_cont_list
namelist /main_nml/ ind_targ_extr, innermost_level, cyclone, date, & namelist /main_nml/ ind_targ_extr, innermost_level, cyclone, date, &
grid_lon_lat grid_lon_lat
logical exist
!---------------------------------------------------------------- !----------------------------------------------------------------
...@@ -68,8 +69,16 @@ program test_get_1_outerm ...@@ -68,8 +69,16 @@ program test_get_1_outerm
write(unit = *, nml = main_nml) write(unit = *, nml = main_nml)
call config call config
call input_ssh(corner, step, periodic, ssh, u, v) call input_ssh(corner, step, periodic, ssh, u, v)
call csvread("outside_points.csv", outside_points, first_r = 3) inquire(file = "outside_points.csv", exist = exist)
outside_points = transpose(outside_points) * deg_to_rad
if (exist) then
call csvread("outside_points.csv", outside_points, first_r = 3)
outside_points = transpose(outside_points) * deg_to_rad
else
print *, '"outside_points.csv" not found. Assuming no outside points...'
allocate(outside_points(2, 0))
end if
e%extr%coord_proj = ind_targ_extr e%extr%coord_proj = ind_targ_extr
e%extr%coord = corner + (ind_targ_extr - 1) * step e%extr%coord = corner + (ind_targ_extr - 1) * step
call get_1_outerm(e%out_cont, cont_list, n_cont, cyclone, e%extr%coord, & call get_1_outerm(e%out_cont, cont_list, n_cont, cyclone, e%extr%coord, &
......
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