add new fprintf for debuging viscosity