bcf_call_glfgen(n_plp[i], plp[i], ref16, bca, bcr + i);
bcf_call_combine(n, bcr, ref16, &bc);
bcf_call2bcf(tid, pos, &bc, b);
bcf_call_glfgen(n_plp[i], plp[i], ref16, bca, bcr + i);
bcf_call_combine(n, bcr, ref16, &bc);
bcf_call2bcf(tid, pos, &bc, b);