Chun's version of hydro: include evolution of bulk pressure and output
hydro information.
Skip the first line of sd_event_#_block.dat when read that from old
superMC generated profile. In effect when InitialURead==0.
InputFun_1.29RC6: now can read in event_phi2, visbulk and InitialURead.
Hydro: Add more output to hydro.
VISHNew: small bug fix, now can write InitialURead to file correctly.
Test on bulk viscosity.
VBulk resume to temperature dependent one.
Test: output umu, ed, and pl when epsilon_p is large than 100.
VISHNew: ouput total momentum anisotropy for cells inside freeze-out surface.
Add hybrid invertfunction method. Implement it on velocity search.
invertFunctionD: constrain the value of XX2, make it always in the boundary.
Resume to use invertfunctionD() to find velocity. invertfunctionH()
found both f(x_l) and f(x_h) <0.
Stop check large epsilon_p, which causes program to exit.
Limit the value of bulkpi in FYSCOR_velocity(), and use invertfunctionH.
Update the numerical zero in invertfunctionH().
Constrain the boundary of bulkpi in FYSCOR_velocity().
Add initial guess to solving U0_local.
Update invertfunctionH(). now need additional factor varY, the number to
be inverted.
VISHNew: use constant zeta/s.
bug fix: another comparision between Visbulk and 10E-5.
VISHNew bug fix: reset nan output by TEpsP_inside.
FYSCOR_velocity(), modified bulkpi using the minimum value of PL/Ed.