Skip to content

Commit 4166fc2

Browse files
wilfonbawilfonba
committed
precheck and bug fixes
1 parent 6fc430a commit 4166fc2

5 files changed

Lines changed: 56 additions & 56 deletions

File tree

‎src/simulation/m_compute_levelset.fpp‎

Lines changed: 23 additions & 23 deletions
Original file line numberDiff line numberDiff line change
@@ -102,8 +102,8 @@ contains
102102

103103
radius = patch_ib(ib_patch_id)%radius
104104

105-
dist_vec(1) = x_cc(i) - patch_ib(ib_patch_id)%x_centroid - real(gp%x_periodicity, wp)*(x_domain%end - x_domain%beg)
106-
dist_vec(2) = y_cc(j) - patch_ib(ib_patch_id)%y_centroid - real(gp%y_periodicity, wp)*(y_domain%end - y_domain%beg)
105+
dist_vec(1) = x_cc(i) - patch_ib(ib_patch_id)%x_centroid - real(gp%x_periodicity, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
106+
dist_vec(2) = y_cc(j) - patch_ib(ib_patch_id)%y_centroid - real(gp%y_periodicity, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
107107
dist_vec(3) = 0._wp
108108
dist = sqrt(sum(dist_vec**2))
109109

@@ -137,8 +137,8 @@ contains
137137
i = gp%loc(1)
138138
j = gp%loc(2)
139139

140-
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(x_domain%end - x_domain%beg)
141-
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(y_domain%end - y_domain%beg)
140+
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
141+
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
142142
inverse_rotation(:, :) = patch_ib(ib_patch_id)%rotation_matrix_inverse(:, :)
143143
rotation(:, :) = patch_ib(ib_patch_id)%rotation_matrix(:, :)
144144
offset(:) = patch_ib(ib_patch_id)%centroid_offset(:)
@@ -223,9 +223,9 @@ contains
223223
j = gp%loc(2)
224224
l = gp%loc(3)
225225
226-
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(x_domain%end - x_domain%beg)
227-
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(y_domain%end - y_domain%beg)
228-
center(3) = patch_ib(ib_patch_id)%z_centroid + real(gp%z_periodicity, wp)*(z_domain%end - z_domain%beg)
226+
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
227+
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
228+
center(3) = patch_ib(ib_patch_id)%z_centroid + real(gp%z_periodicity, wp)*(glb_bounds(3)%end - glb_bounds(3)%beg)
229229
lz = patch_ib(ib_patch_id)%length_z
230230
inverse_rotation(:, :) = patch_ib(ib_patch_id)%rotation_matrix_inverse(:, :)
231231
rotation(:, :) = patch_ib(ib_patch_id)%rotation_matrix(:, :)
@@ -328,8 +328,8 @@ contains
328328

329329
length_x = patch_ib(ib_patch_id)%length_x
330330
length_y = patch_ib(ib_patch_id)%length_y
331-
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(x_domain%end - x_domain%beg)
332-
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(y_domain%end - y_domain%beg)
331+
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
332+
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
333333
inverse_rotation(:, :) = patch_ib(ib_patch_id)%rotation_matrix_inverse(:, :)
334334
rotation(:, :) = patch_ib(ib_patch_id)%rotation_matrix(:, :)
335335

@@ -399,8 +399,8 @@ contains
399399

400400
length_x = patch_ib(ib_patch_id)%length_x
401401
length_y = patch_ib(ib_patch_id)%length_y
402-
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(x_domain%end - x_domain%beg)
403-
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(y_domain%end - y_domain%beg)
402+
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
403+
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
404404
inverse_rotation(:, :) = patch_ib(ib_patch_id)%rotation_matrix_inverse(:, :)
405405
rotation(:, :) = patch_ib(ib_patch_id)%rotation_matrix(:, :)
406406

@@ -453,9 +453,9 @@ contains
453453
length_y = patch_ib(ib_patch_id)%length_y
454454
length_z = patch_ib(ib_patch_id)%length_z
455455

456-
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(x_domain%end - x_domain%beg)
457-
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(y_domain%end - y_domain%beg)
458-
center(3) = patch_ib(ib_patch_id)%z_centroid + real(gp%z_periodicity, wp)*(z_domain%end - z_domain%beg)
456+
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
457+
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
458+
center(3) = patch_ib(ib_patch_id)%z_centroid + real(gp%z_periodicity, wp)*(glb_bounds(3)%end - glb_bounds(3)%beg)
459459

460460
inverse_rotation(:, :) = patch_ib(ib_patch_id)%rotation_matrix_inverse(:, :)
461461
rotation(:, :) = patch_ib(ib_patch_id)%rotation_matrix(:, :)
@@ -535,9 +535,9 @@ contains
535535
k = gp%loc(3)
536536
537537
radius = patch_ib(ib_patch_id)%radius
538-
periodicity(1) = real(gp%x_periodicity, wp)*(x_domain%end - x_domain%beg)
539-
periodicity(2) = real(gp%y_periodicity, wp)*(y_domain%end - y_domain%beg)
540-
periodicity(3) = real(gp%z_periodicity, wp)*(z_domain%end - z_domain%beg)
538+
periodicity(1) = real(gp%x_periodicity, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
539+
periodicity(2) = real(gp%y_periodicity, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
540+
periodicity(3) = real(gp%z_periodicity, wp)*(glb_bounds(3)%end - glb_bounds(3)%beg)
541541
center(1) = patch_ib(ib_patch_id)%x_centroid
542542
center(2) = patch_ib(ib_patch_id)%y_centroid
543543
center(3) = patch_ib(ib_patch_id)%z_centroid
@@ -579,9 +579,9 @@ contains
579579
k = gp%loc(3)
580580
581581
radius = patch_ib(ib_patch_id)%radius
582-
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(x_domain%end - x_domain%beg)
583-
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(y_domain%end - y_domain%beg)
584-
center(3) = patch_ib(ib_patch_id)%z_centroid + real(gp%z_periodicity, wp)*(z_domain%end - z_domain%beg)
582+
center(1) = patch_ib(ib_patch_id)%x_centroid + real(gp%x_periodicity, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
583+
center(2) = patch_ib(ib_patch_id)%y_centroid + real(gp%y_periodicity, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
584+
center(3) = patch_ib(ib_patch_id)%z_centroid + real(gp%z_periodicity, wp)*(glb_bounds(3)%end - glb_bounds(3)%beg)
585585
length(1) = patch_ib(ib_patch_id)%length_x
586586
length(2) = patch_ib(ib_patch_id)%length_y
587587
length(3) = patch_ib(ib_patch_id)%length_z
@@ -658,10 +658,10 @@ contains
658658
total_vertices = gpu_total_vertices(patch_id)
659659

660660
center = 0._wp
661-
if (.not. f_is_default(patch_ib(patch_id)%x_centroid)) center(1) = patch_ib(patch_id)%x_centroid + real(gp%x_periodicity, wp)*(x_domain%end - x_domain%beg)
662-
if (.not. f_is_default(patch_ib(patch_id)%y_centroid)) center(2) = patch_ib(patch_id)%y_centroid + real(gp%y_periodicity, wp)*(y_domain%end - y_domain%beg)
661+
if (.not. f_is_default(patch_ib(patch_id)%x_centroid)) center(1) = patch_ib(patch_id)%x_centroid + real(gp%x_periodicity, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
662+
if (.not. f_is_default(patch_ib(patch_id)%y_centroid)) center(2) = patch_ib(patch_id)%y_centroid + real(gp%y_periodicity, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
663663
if (p > 0) then
664-
if (.not. f_is_default(patch_ib(patch_id)%z_centroid)) center(3) = patch_ib(patch_id)%z_centroid + real(gp%z_periodicity, wp)*(z_domain%end - z_domain%beg)
664+
if (.not. f_is_default(patch_ib(patch_id)%z_centroid)) center(3) = patch_ib(patch_id)%z_centroid + real(gp%z_periodicity, wp)*(glb_bounds(3)%end - glb_bounds(3)%beg)
665665
end if
666666

667667
inverse_rotation(:, :) = patch_ib(patch_id)%rotation_matrix_inverse(:, :)

‎src/simulation/m_global_parameters.fpp‎

Lines changed: 6 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -713,8 +713,8 @@ contains
713713
bub_pp%gam_g = dflt_real; gam_g = dflt_real
714714
bub_pp%M_v = dflt_real; M_v = dflt_real
715715
bub_pp%M_g = dflt_real; M_g = dflt_real
716-
bub_pp%k_v = dflt_real;
717-
bub_pp%k_g = dflt_real;
716+
bub_pp%k_v = dflt_real;
717+
bub_pp%k_g = dflt_real;
718718
bub_pp%cp_v = dflt_real; cp_v = dflt_real
719719
bub_pp%cp_g = dflt_real; cp_g = dflt_real
720720
bub_pp%R_v = dflt_real; R_v = dflt_real
@@ -1377,15 +1377,15 @@ contains
13771377
@:PREFER_GPU(x_cc)
13781378
@:PREFER_GPU(dx)
13791379
1380-
if (n == 0) return;
1380+
if (n == 0) return;
13811381
@:ALLOCATE(y_cb(-1 - buff_size:n + buff_size))
13821382
@:ALLOCATE(y_cc(-buff_size:n + buff_size))
13831383
@:ALLOCATE(dy(-buff_size:n + buff_size))
13841384
@:PREFER_GPU(y_cb)
13851385
@:PREFER_GPU(y_cc)
13861386
@:PREFER_GPU(dy)
13871387
1388-
if (p == 0) return;
1388+
if (p == 0) return;
13891389
@:ALLOCATE(z_cb(-1 - buff_size:p + buff_size))
13901390
@:ALLOCATE(z_cc(-buff_size:p + buff_size))
13911391
@:ALLOCATE(dz(-buff_size:p + buff_size))
@@ -1476,10 +1476,10 @@ contains
14761476
! Deallocating grid variables for the x-, y- and z-directions
14771477
@:DEALLOCATE(x_cb, x_cc, dx)
14781478
1479-
if (n == 0) return;
1479+
if (n == 0) return;
14801480
@:DEALLOCATE(y_cb, y_cc, dy)
14811481
1482-
if (p == 0) return;
1482+
if (p == 0) return;
14831483
@:DEALLOCATE(z_cb, z_cc, dz)
14841484
14851485
end subroutine s_finalize_global_parameters_module

‎src/simulation/m_ib_patches.fpp‎

Lines changed: 25 additions & 25 deletions
Original file line numberDiff line numberDiff line change
@@ -143,8 +143,8 @@ contains
143143

144144
! Transferring the circular patch's radius, centroid, smearing patch
145145
! identity and smearing coefficient information
146-
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(x_domain%end - x_domain%beg)
147-
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(y_domain%end - y_domain%beg)
146+
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
147+
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
148148
radius = patch_ib(patch_id)%radius
149149
150150
! encode the periodicity information into the patch_id
@@ -197,8 +197,8 @@ contains
197197
real(wp), dimension(1:2) :: center !< x and y coordinates in local IB frame
198198
real(wp), dimension(1:3, 1:3) :: inverse_rotation
199199
200-
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(x_domain%end - x_domain%beg)
201-
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(y_domain%end - y_domain%beg)
200+
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
201+
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
202202
ca_in = patch_ib(patch_id)%c
203203
pa = patch_ib(patch_id)%p
204204
ma = patch_ib(patch_id)%m
@@ -360,9 +360,9 @@ contains
360360
real(wp), dimension(1:3) :: xyz_local, center, offset !< x, y, z coordinates in local IB frame
361361
real(wp), dimension(1:3, 1:3) :: inverse_rotation
362362

363-
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(x_domain%end - x_domain%beg)
364-
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(y_domain%end - y_domain%beg)
365-
center(3) = patch_ib(patch_id)%z_centroid + real(zp, wp)*(z_domain%end - z_domain%beg)
363+
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
364+
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
365+
center(3) = patch_ib(patch_id)%z_centroid + real(zp, wp)*(glb_bounds(3)%end - glb_bounds(3)%beg)
366366
lz = patch_ib(patch_id)%length_z
367367
ca_in = patch_ib(patch_id)%c
368368
pa = patch_ib(patch_id)%p
@@ -533,8 +533,8 @@ contains
533533
real(wp), dimension(1:3, 1:3) :: inverse_rotation
534534
535535
! Transferring the rectangle's centroid and length information
536-
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(x_domain%end - x_domain%beg)
537-
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(y_domain%end - y_domain%beg)
536+
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
537+
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
538538
length(1) = patch_ib(patch_id)%length_x
539539
length(2) = patch_ib(patch_id)%length_y
540540
inverse_rotation(:, :) = patch_ib(patch_id)%rotation_matrix_inverse(:, :)
@@ -604,9 +604,9 @@ contains
604604

605605
! Transferring spherical patch's radius, centroid, smoothing patch
606606
! identity and smoothing coefficient information
607-
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(x_domain%end - x_domain%beg)
608-
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(y_domain%end - y_domain%beg)
609-
center(3) = patch_ib(patch_id)%z_centroid + real(zp, wp)*(z_domain%end - z_domain%beg)
607+
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
608+
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
609+
center(3) = patch_ib(patch_id)%z_centroid + real(zp, wp)*(glb_bounds(3)%end - glb_bounds(3)%beg)
610610
radius = patch_ib(patch_id)%radius
611611
612612
! encode the periodicity information into the patch_id
@@ -675,9 +675,9 @@ contains
675675
real(wp) :: corner_distance
676676
677677
! Transferring the cuboid's centroid and length information
678-
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(x_domain%end - x_domain%beg)
679-
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(y_domain%end - y_domain%beg)
680-
center(3) = patch_ib(patch_id)%z_centroid + real(zp, wp)*(z_domain%end - z_domain%beg)
678+
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
679+
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
680+
center(3) = patch_ib(patch_id)%z_centroid + real(zp, wp)*(glb_bounds(3)%end - glb_bounds(3)%beg)
681681
length(1) = patch_ib(patch_id)%length_x
682682
length(2) = patch_ib(patch_id)%length_y
683683
length(3) = patch_ib(patch_id)%length_z
@@ -760,9 +760,9 @@ contains
760760
real(wp) :: corner_distance
761761

762762
! Transferring the cylindrical patch's centroid, length, radius,
763-
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(x_domain%end - x_domain%beg)
764-
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(y_domain%end - y_domain%beg)
765-
center(3) = patch_ib(patch_id)%z_centroid + real(zp, wp)*(z_domain%end - z_domain%beg)
763+
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
764+
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
765+
center(3) = patch_ib(patch_id)%z_centroid + real(zp, wp)*(glb_bounds(3)%end - glb_bounds(3)%beg)
766766
length(1) = patch_ib(patch_id)%length_x
767767
length(2) = patch_ib(patch_id)%length_y
768768
length(3) = patch_ib(patch_id)%length_z
@@ -845,8 +845,8 @@ contains
845845
real(wp), dimension(1:3, 1:3) :: inverse_rotation
846846

847847
! Transferring the ellipse's centroid and length information
848-
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(x_domain%end - x_domain%beg)
849-
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(y_domain%end - y_domain%beg)
848+
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
849+
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
850850
ellipse_coeffs(1) = 0.5_wp*patch_ib(patch_id)%length_x
851851
ellipse_coeffs(2) = 0.5_wp*patch_ib(patch_id)%length_y
852852
inverse_rotation(:, :) = patch_ib(patch_id)%rotation_matrix_inverse(:, :)
@@ -905,8 +905,8 @@ contains
905905
real(wp), dimension(1:3, 1:3) :: inverse_rotation, rotation
906906
907907
center = 0._wp
908-
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(x_domain%end - x_domain%beg)
909-
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(y_domain%end - y_domain%beg)
908+
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
909+
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
910910
inverse_rotation(:, :) = patch_ib(patch_id)%rotation_matrix_inverse(:, :)
911911
rotation(:, :) = patch_ib(patch_id)%rotation_matrix(:, :)
912912
offset(:) = patch_ib(patch_id)%centroid_offset(:)
@@ -987,9 +987,9 @@ contains
987987
real(wp), dimension(1:3) :: bbox_min, bbox_max, local_corner, world_corner
988988
989989
center = 0._wp
990-
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(x_domain%end - x_domain%beg)
991-
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(y_domain%end - y_domain%beg)
992-
center(3) = patch_ib(patch_id)%z_centroid + real(zp, wp)*(z_domain%end - z_domain%beg)
990+
center(1) = patch_ib(patch_id)%x_centroid + real(xp, wp)*(glb_bounds(1)%end - glb_bounds(1)%beg)
991+
center(2) = patch_ib(patch_id)%y_centroid + real(yp, wp)*(glb_bounds(2)%end - glb_bounds(2)%beg)
992+
center(3) = patch_ib(patch_id)%z_centroid + real(zp, wp)*(glb_bounds(3)%end - glb_bounds(3)%beg)
993993
inverse_rotation(:, :) = patch_ib(patch_id)%rotation_matrix_inverse(:, :)
994994
offset(:) = patch_ib(patch_id)%centroid_offset(:)
995995
spc = patch_ib(patch_id)%model_spc

‎src/simulation/m_ibm.fpp‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -706,7 +706,7 @@ contains
706706
if (p /= 0) then
707707
k = gp%ip_grid(3)
708708
else
709-
k = 0;
709+
k = 0;
710710
end if
711711
712712
! get the distance to a cell in each direction

‎src/simulation/m_sim_helpers.fpp‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -222,7 +222,7 @@ contains
222222
Rc = min(dx(j)*(abs(vel(1)) + c), &
223223
dy(k)*(abs(vel(2)) + c)) &
224224
/maxval(1._wp/Re_l)
225-
else !1D
225+
else !1D
226226
vcfl = maxval(dt/Re_l/rho)/dx(j)**2._wp
227227
Rc = dx(j)*(abs(vel(1)) + c)/maxval(1._wp/Re_l)
228228
end if

0 commit comments

Comments
 (0)