bug in calculation of weights for SCRIP format (routine compute_weights_poly)