add function from dominick for zbi calculation
add check -- if signif is too big, just return 8
add pujetid option