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