@@ -46,6 +46,17 @@ contains
46
46
if (i <= num_patches) then
47
47
! call s_check_patch_geometry(i)
48
48
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
+
49
60
! Constraints on the geometric initial condition patch parameters
50
61
if (patch_icpp(i)%geometry == 1 ) then
51
62
call s_check_line_segment_patch_geometry(i)
@@ -57,12 +68,6 @@ contains
57
68
call s_check_line_sweep_patch_geometry(i)
58
69
elseif (patch_icpp(i)%geometry == 5 ) then
59
70
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" )
66
71
elseif (patch_icpp(i)%geometry == 8 ) then
67
72
call s_check_sphere_patch_geometry(i)
68
73
elseif (patch_icpp(i)%geometry == 9 ) then
@@ -73,31 +78,21 @@ contains
73
78
call s_check_plane_sweep_patch_geometry(i)
74
79
elseif (patch_icpp(i)%geometry == 12 ) then
75
80
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" )
79
81
elseif (patch_icpp(i)%geometry == 14 ) then
80
82
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" )
84
83
elseif (patch_icpp(i)%geometry == 20 ) then
85
84
call s_check_2D_TaylorGreen_vortex_patch_geometry(i)
86
85
elseif (patch_icpp(i)%geometry == 21 ) then
87
86
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" )
90
87
else
91
88
call s_prohibit_abort(" Invalid patch geometry number" , " patch_icpp(" // trim (iStr)// " )%geometry " // &
92
89
" must be between 1 and 21" )
93
90
end if
94
91
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)
101
96
end if
102
97
end do
103
98
0 commit comments