new fprintf added for debuging viscosity