read_3d and coe_sum

Add a new read_3d function, and modified coe_sum to output pdf in the
same folder as those coefficients output.