level =  5  
refinementLevel = 1
label = 'channel'
idLabel = 'channel'
refOmega = 1.9
refLevel = 6
u_in = 0.02
-- geometry definition
height =  0.41
length = 5*height
radius = 0.05
usePeriodic = false 
useObstacle = false
testIntersection = false
--
m = 2^(refLevel - level)
omega = 1./(1/m*(1/refOmega -0.5)+0.5)
maxLevel = level+refinementLevel

dx     = length/(2^level)
dxMin  = length/(2^maxLevel)
dxDash = 0.5*dxMin
viscLB = 1/3*(1/omega-1/2)
-- physical reference values
dxPhys = dx
lPhys = length  -- m
csPhys = 300    -- m/s
csLB = 1/math.sqrt(3)
--acoustic scaling with speeds of sound as reference
-- csPhys = csLB * dxPhys / dtPhys
dtPhys = csLB/csPhys*dxPhys
rho0 = 1  -- kg/m^3
viscPhys = viscLB*dxPhys^2/dtPhys


nElemsMax = 2^maxLevel
if useObstacle == true then
  nIters = nElemsMax*50
else
  nIters = nElemsMax*40
end
nIters = 1 --nElemsMax


