def METHOD_NAME(
    create_skymodel_in_phase_center, bda_predict_step_type
):
    check_call(
        [
            tcf.DP3EXE,
            "checkparset=1",
            f"msin={MSIN_REGULAR}",
            "msout=bdapredict.MS",
            "msout.overwrite=true",
            "steps=[bdaaverager, predict]",
            "bdaaverager.timebase=600",
            "bdaaverager.frequencybase=1000",
            f"predict.type={bda_predict_step_type}",
            "predict.sourcedb=test.sourcedb",
            "predict.usebeammodel=F",
        ]
    )
    taql_command = (
        "select from bdapredict.MS "
        "where all(near(abs(DATA[,0]),10,1e-6))"
        "  and all(near(abs(DATA[,1]), 0,1e-6))"
        "  and all(near(abs(DATA[,2]), 0,1e-6))"
        "  and all(near(abs(DATA[,3]),10,1e-6))"
    )
    assert_taql(taql_command, 153)