moved working fn to BMLfn, started report - plotting p, pc red vs velocity