Add fns to generate covariance from flow matrices

[mpadge]

Also include a zero_is_na = FALSE arg, and dump a message if changed about the danger of including zeros in covariance calculations,