added uni initial velocity field

p-costa · Nov 29, 2021 · 00a1cae · 00a1cae
1 parent 2cdfe59
commit 00a1cae
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 0 deletions.
diff --git a/src/INFO_INPUT.md b/src/INFO_INPUT.md
@@ -251,6 +251,7 @@ This line sets the initial velocity field for the block in question.
 `initvel` **chooses the initial velocity field**. The following options are available:
 
 * `zer`: zero velocity field
+* `uni`: uniform velocity field equal to `uref` ; streamwise direction in `x`
 * `cou`: plane Couette flow profile with symmetric wall velocities equal to `uref/2`; streamwise direction in `x`
 * `poi`: plane Poiseuille flow profile with mean velocity `uref` ; streamwise direction in `x`
 * `log`: logarithmic profile with mean velocity `uref` ; streamwise direction in `x`

diff --git a/src/initflow.f90 b/src/initflow.f90
@@ -41,6 +41,8 @@ subroutine initflow(inivel,is_wallturb,lo,hi,lo_g,hi_g,l,uref,lref,visc,bforce,
  is_mean=.true.
  case('zer')
  u1d(:) = 0._rp
+ case('uni')
+ u1d(:) = uref
  case('log')
  call log_profile(lo(3),hi(3),zc,l(3),uref,lref,visc,u1d)
  is_noise = .true.