Skip to content

Commit 1584480

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 44d63f6 commit 1584480

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
@@ -81,8 +81,8 @@ contains
8181
perturb_alpha = q_prim_vf(E_idx + perturb_flow_fluid)%sf(i, j, k)
8282
call random_number(rand_real)
8383
rand_real = rand_real*perturb_flow_mag
84-
q_prim_vf(mom_idx%beg)%sf(i, j, k) = (1._wp + rand_real)*q_prim_vf(mom_idx%beg)%sf(i, j, k)
8584
q_prim_vf(mom_idx%end)%sf(i, j, k) = rand_real*q_prim_vf(mom_idx%beg)%sf(i, j, k)
85+
q_prim_vf(mom_idx%beg)%sf(i, j, k) = (1._wp + rand_real)*q_prim_vf(mom_idx%beg)%sf(i, j, k)
8686
if (bubbles_euler) then
8787
q_prim_vf(alf_idx)%sf(i, j, k) = (1._wp + rand_real)*q_prim_vf(alf_idx)%sf(i, j, k)
8888
end if

0 commit comments

Comments
 (0)