loop like a native
renaming
fix fortran source (was printing for each species)
add viscosity regression test