calc_vis

input: int qt, double *e, double *mu, double vis, double e_0, double e_ref, double *wk

output:

Description: This function calculates viscosity using sutherlands law in the form: mu = vis (e/e_0)(3/2) (e_0+e_ref)/(e + e_ref) if e_ref is not defined the just fill with vis