Average and output rho and volume in plotSDNoAverageCesaro-fit-ave2.m

The upperhalf rho is now called rho2
Change dataBaseName to dataFilename in plotNoAverageCesaro-fit-ave.m

And make it runnable with only one MD; that is no averaging.
Add volume to the output of average-sdCorr.m

and fix a bug of reading rho incorrectly