@@ -63,7 +63,7 @@ contains
6363 patch_ib(i)%geometry == 12 ) then
6464 call s_check_model_ib_patch_geometry(i)
6565 else if (patch_ib(i)%geometry == 6 ) then
66- print * , " Ellipse Patch "
66+ call s_check_ellipse_ib_patch_geometry(i)
6767 else
6868 call s_prohibit_abort(" Invalid IB patch" , &
6969 " patch_ib(" // trim (iStr)// " )%geometry must be " // &
@@ -94,6 +94,25 @@ contains
9494 .or. f_is_default(patch_ib(patch_id)%y_centroid), &
9595 ' in circle IB patch ' // trim (iStr))
9696
97+ end subroutine s_check_circle_ib_patch_geometry
98+
99+ !> This subroutine verifies that the geometric parameters of
100+ !! the circle patch have consistently been inputted by the
101+ !! user.
102+ !! @param patch_id Patch identifier
103+ impure subroutine s_check_ellipse_ib_patch_geometry (patch_id )
104+
105+ integer , intent (in ) :: patch_id
106+
107+ call s_int_to_str(patch_id, iStr)
108+
109+ @:PROHIBIT(n == 0 .or. p > 0 &
110+ .or. patch_ib(patch_id)%length_x <= 0._wp &
111+ .or. patch_ib(patch_id)%length_y <= 0._wp &
112+ .or. f_is_default(patch_ib(patch_id)%x_centroid) &
113+ .or. f_is_default(patch_ib(patch_id)%y_centroid), &
114+ ' in ellipse IB patch ' // trim (iStr))
115+
97116 end subroutine s_check_circle_ib_patch_geometry
98117
99118 !> This subroutine verifies that the geometric parameters of
0 commit comments