@@ -46,6 +46,17 @@ contains
4646 if (i <= num_patches) then
4747 ! call s_check_patch_geometry(i)
4848 call s_int_to_str(i, iStr)
49+ @:PROHIBIT(patch_icpp(i)%geometry == 6 , " Invalid patch geometry number. " // &
50+ " patch_icpp(" // trim (iStr)// " )%geometry is deprecated." )
51+ @:PROHIBIT(patch_icpp(i)%geometry == 7 , " Invalid patch geometry number. " // &
52+ " patch_icpp(" // trim (iStr)// " )%geometry is deprecated." )
53+ @:PROHIBIT(patch_icpp(i)%geometry == 13 , " Invalid patch geometry number. " // &
54+ " patch_icpp(" // trim (iStr)// " )%geometry is deprecated." )
55+ @:PROHIBIT(patch_icpp(i)%geometry == 15 , " Invalid patch geometry number. " // &
56+ " patch_icpp(" // trim (iStr)// " )%geometry is deprecated." )
57+ @:PROHIBIT(patch_icpp(i)%geometry == dflt_int, " Invalid patch geometry number. " // &
58+ " patch_icpp(" // trim (iStr)// " )%geometry must be set." )
59+
4960 ! Constraints on the geometric initial condition patch parameters
5061 if (patch_icpp(i)%geometry == 1 ) then
5162 call s_check_line_segment_patch_geometry(i)
@@ -57,12 +68,6 @@ contains
5768 call s_check_line_sweep_patch_geometry(i)
5869 elseif (patch_icpp(i)%geometry == 5 ) then
5970 call s_check_ellipse_patch_geometry(i)
60- elseif (patch_icpp(i)%geometry == 6 ) then
61- call s_prohibit_abort(" Invalid patch geometry number" , " patch_icpp(" // trim (iStr)// " )%geometry " // &
62- " is deprecated" )
63- elseif (patch_icpp(i)%geometry == 7 ) then
64- call s_prohibit_abort(" Invalid patch geometry number" , " patch_icpp(" // trim (iStr)// " )%geometry " // &
65- " is deprecated" )
6671 elseif (patch_icpp(i)%geometry == 8 ) then
6772 call s_check_sphere_patch_geometry(i)
6873 elseif (patch_icpp(i)%geometry == 9 ) then
@@ -73,31 +78,21 @@ contains
7378 call s_check_plane_sweep_patch_geometry(i)
7479 elseif (patch_icpp(i)%geometry == 12 ) then
7580 call s_check_ellipsoid_patch_geometry(i)
76- elseif (patch_icpp(i)%geometry == 13 ) then
77- call s_prohibit_abort(" Invalid patch geometry number" , " patch_icpp(" // trim (iStr)// " )%geometry " // &
78- " is deprecated" )
7981 elseif (patch_icpp(i)%geometry == 14 ) then
8082 call s_check_spherical_harmonic_patch_geometry(i)
81- elseif (patch_icpp(i)%geometry == 15 ) then
82- call s_prohibit_abort(" Invalid patch geometry number" , " patch_icpp(" // trim (iStr)// " )%geometry " // &
83- " is deprecated" )
8483 elseif (patch_icpp(i)%geometry == 20 ) then
8584 call s_check_2D_TaylorGreen_vortex_patch_geometry(i)
8685 elseif (patch_icpp(i)%geometry == 21 ) then
8786 call s_check_model_geometry(i)
88- elseif (patch_icpp(i)%geometry == dflt_int) then
89- call s_prohibit_abort(" Active patch undefined" , " patch_icpp(" // trim (iStr)// " )%geometry must be set" )
9087 else
9188 call s_prohibit_abort(" Invalid patch geometry number" , " patch_icpp(" // trim (iStr)// " )%geometry " // &
9289 " must be between 1 and 21" )
9390 end if
9491 else
95- if (patch_icpp(i)%geometry == dflt_int) then
96- call s_check_inactive_patch_geometry(i)
97- else
98- call s_prohibit_abort(" Inactive patch defined" , " patch_icpp(" // trim (iStr)// " )%geometry must not be set. " // &
99- " Patch " // trim (iStr)// " is inactive as the number of patches is " // trim (num_patches_str))
100- end if
92+ @:PROHIBIT(patch_icpp(i)%geometry /= dflt_int, " Inactive patch defined. " // &
93+ " patch_icpp(" // trim (iStr)// " )%geometry not be set for inactive patches. " // &
94+ " Patch " // trim (iStr)// " is inactive as the number of patches is " // trim (num_patches_str))
95+ call s_check_inactive_patch_geometry(i)
10196 end if
10297 end do
10398
0 commit comments