coefficients = RANS3PF.Coefficients(
    epsFact=epsFact_viscosity,
    sigma=0.0,
    rho_0=rho_0,
    nu_0=nu_0,
    rho_1=rho_1,
    nu_1=nu_1,
    g=g,
    nd=nd,
    ME_model=V_model,
    PRESSURE_model=PRESSURE_model,
    SED_model=SED_model,
    VOS_model=VOS_model,
    VOF_model=VOF_model,
    LS_model=LS_model,
    Closure_0_model=Closure_0_model,
    Closure_1_model=Closure_1_model,
    epsFact_density=epsFact_density,
    stokes=False,
    useVF=useVF,
    useRBLES=useRBLES,
    useMetrics=useMetrics,
    eb_adjoint_sigma=1.0,
    eb_penalty_constant=weak_bc_penalty_constant,
    forceStrongDirichlet=ns_forceStrongDirichlet,
    turbulenceClosureModel=ns_closure,
    movingDomain=movingDomain,
    dragAlpha=dragAlpha,
    PSTAB=1.0,
    cE=cE,
    cMax=cMax,
    CORRECT_VELOCITY=CORRECT_VELOCITY)
coefficients = RANS3PF.Coefficients(
    epsFact=epsFact_viscosity,
    sigma=0.0,
    rho_0=rho_0,
    nu_0=nu_0,
    rho_1=rho_1,
    nu_1=nu_1,
    g=g,
    nd=nd,
    ME_model=V_model,
    PRESSURE_model=PRESSURE_model,
    SED_model=SED_model,
    VOF_model=VOF_model,
    VOS_model=VOS_model,
    LS_model=LS_model,
    Closure_0_model=Closure_0_model,
    Closure_1_model=Closure_1_model,
    epsFact_density=epsFact_density,
    stokes=False,
    useVF=useVF,
    useRBLES=useRBLES,
    useMetrics=useMetrics,
    eb_adjoint_sigma=1.0,
    eb_penalty_constant=weak_bc_penalty_constant,
    forceStrongDirichlet=ns_forceStrongDirichlet,
    turbulenceClosureModel=ns_closure,
    movingDomain=movingDomain,
    dragAlpha=dragAlpha,
    PSTAB=ct.opts.PSTAB,
    aDarcy=sedClosure.aDarcy,
    betaForch=sedClosure.betaForch,
    grain=sedClosure.grain,
    packFraction=sedClosure.packFraction,
    maxFraction=sedClosure.maxFraction,
    frFraction=sedClosure.frFraction,
    sigmaC=sedClosure.sigmaC,
    C3e=sedClosure.C3e,
    C4e=sedClosure.C4e,
    eR=sedClosure.eR,
    fContact=sedClosure.fContact,
    mContact=sedClosure.mContact,
    nContact=sedClosure.nContact,
    angFriction=sedClosure.angFriction,
    vos_function=ct.vos_function,
    vos_limiter=ct.sedClosure.vos_limiter,
    mu_fr_limiter=ct.sedClosure.mu_fr_limiter,
)
coefficients = RANS3PF.Coefficients(
    epsFact=epsFact_viscosity,
    sigma=0.0,
    rho_0=rho_0,
    nu_0=nu_0,
    rho_1=rho_1,
    nu_1=nu_1,
    g=g,
    nd=nd,
    ME_model=V_model,
    PRESSURE_model=PRESSURE_model,
    SED_model=SED_model,
    VOF_model=VOF_model,
    VOS_model=VOS_model,
    LS_model=LS_model,
    Closure_0_model=Closure_0_model,
    Closure_1_model=Closure_1_model,
    epsFact_density=epsFact_density,
    stokes=False,
    useVF=useVF,
    useRBLES=useRBLES,
    useMetrics=useMetrics,
    eb_adjoint_sigma=1.0,
    eb_penalty_constant=weak_bc_penalty_constant,
    forceStrongDirichlet=ns_forceStrongDirichlet,
    turbulenceClosureModel=ns_closure,
    movingDomain=movingDomain,
    dragAlpha=dragAlpha,
    PSTAB=0.0,
    nParticles=1,
    particle_epsFact=1.5,
    particle_alpha=1e6,
    particle_beta=1e6,
    particle_penalty_constant=100.0,
    particle_sdfList=[particle_sdf],
    particle_velocityList=[particle_vel],
    use_ball_as_particle=1,
    ball_center=numpy.array([[0.2, 0.2, 0.0]]),
    ball_radius=numpy.array([0.05]),
    ball_velocity=numpy.array([[0.0, 0.0, 0.0]]),
    ball_angular_velocity=numpy.array([[0.0, 0.0, 0.0]]),
    use_sbm=USE_SBM)