Skip to content

Commit 5f75aa6

Browse files
sbryngelsonclaude
andcommitted
Fix y-velocity perturbation using already-modified x-velocity
Co-Authored-By: Claude Opus 4.6 <noreply@anthropic.com>
1 parent 1fca6b6 commit 5f75aa6

1 file changed

Lines changed: 1 addition & 1 deletion

File tree

src/pre_process/m_perturbation.fpp

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -83,8 +83,8 @@ contains
8383
perturb_alpha = q_prim_vf(E_idx + perturb_flow_fluid)%sf(i, j, k)
8484
call random_number(rand_real)
8585
rand_real = rand_real*perturb_flow_mag
86-
q_prim_vf(mom_idx%beg)%sf(i, j, k) = (1._wp + rand_real)*q_prim_vf(mom_idx%beg)%sf(i, j, k)
8786
q_prim_vf(mom_idx%end)%sf(i, j, k) = rand_real*q_prim_vf(mom_idx%beg)%sf(i, j, k)
87+
q_prim_vf(mom_idx%beg)%sf(i, j, k) = (1._wp + rand_real)*q_prim_vf(mom_idx%beg)%sf(i, j, k)
8888
if (bubbles_euler) then
8989
q_prim_vf(alf_idx)%sf(i, j, k) = (1._wp + rand_real)*q_prim_vf(alf_idx)%sf(i, j, k)
9090
end if

0 commit comments

Comments
 (0)