iteration   0 rho =  0.668733 enthalpy =     -15.797466
iteration   1 rho =  0.188680 enthalpy =     -15.871148
iteration   2 rho =  0.078088 enthalpy =     -15.874713
iteration   3 rho =  0.048497 enthalpy =     -15.875793
iteration   4 rho =  0.030832 enthalpy =     -15.876304
iteration   5 rho =  0.008297 enthalpy =     -15.876404