diff --git a/Inst_eddies/set_max_speed.f90 b/Inst_eddies/set_max_speed.f90 index 10e3b003d9253b34e5825363cb0bd045970b849f..eea012289e5e94a4af3a0025ff6f828e052ceb8d 100644 --- a/Inst_eddies/set_max_speed.f90 +++ b/Inst_eddies/set_max_speed.f90 @@ -6,20 +6,15 @@ contains subroutine set_max_speed(e, ssh, u, v, corner, step, cont_list) - ! This procedure defines the components speed_cont, max_speed and - ! radius4 of argument e. On return, e%speed_cont may be a null ssh - ! contour but e%max_speed is never set to missing_speed. However, - ! e%max_speed might be NaN if speed_outerm is NaN (with any value - ! of e%radius4). The value of e%radius4 set here is always >= - ! 1. On return, if e%speed_cont is not a null ssh contour then - ! e%max_speed is the speed on e%speed_cont. + ! This procedure defines the components speed_cont and max_speed + ! of argument e. On return, e%speed_cont may be a null ssh contour + ! and e%max_speed may be set to missing_speed. However, + ! e%max_speed is never NaN. On return, if e%speed_cont is not a + ! null ssh contour then e%max_speed is the speed on e%speed_cont. ! The length of the interval of longitudes of the domain, step(1) ! * (size(ssh, 1) - 1), should be < 180 degrees, so that the - ! geometrical processing done here is non-ambiguous. The - ! longitudes of outside points should be in the interval - ! [corner(1), corner(1) + step(1) * (size(ssh, 1) - 1)] as there - ! will be no automatic shifting by a multiple of 360 degrees. + ! geometrical processing done here is non-ambiguous. use, intrinsic:: IEEE_ARITHMETIC, only: IEEE_IS_NAN