Add point to test the case when the test point coincides with

one of the polygon's corner points.

Fixes ufs-community#1000.

tests/orog/ftst_inside_polygon.F90

@@ -145,6 +145,27 @@ program inside_polygon
 
  if (inside) stop 12 ! Test point should be outside polygon.
 
+! Test the case when the test point coincides with a corner point
+! of the polygon. Result should be 'inside'.
+
+ print*, "Test point 7"
+
+ lon1 = 90.0 * D2R      ! Test point.
+ lat1 = 0.0 * D2R
+
+ lon2(1) = 90.0 * D2R   ! Polygon.
+ lat2(1) =  0.0 * D2R
+ lon2(2) = 90.0 * D2R
+ lat2(2) =  1.0 * D2R
+ lon2(3) = 91.0 * D2R
+ lat2(3) =  1.0 * D2R
+ lon2(4) = 91.0 * D2R
+ lat2(4) =  0.0 * D2R
+
+ inside=inside_a_polygon(lon1, lat1, npts, lon2, lat2)
+
+ if (.not.inside) stop 14  ! Test point should be inside polygon.
+
  print*,"OK"
  print*,"SUCCSSS"