& sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2,
_RL sideDragFactor
& )*drF(k)*sideDragFactor
IF ( sideDragFactor.LE.0. ) THEN
& -sideDragFactor
sideMaskFac = sideDragFactor
CALL WRITE_0D_RL( sideDragFactor, INDEX_NONE, 'sideDragFactor =',
& no_slip_sides, sideDragFactor,
sideDragFactor = 2. _d 0