function [mug]=mug(T) 
  %gas dynamic viscosity [Ns/m^2] at T [celsius] 
 
 Tsat=[0:10:370,374.15];
  %saturation temperature in celsius

 muvap=1e-6*[8.105,8.504,8.903,9.305,9.701,10.10,10.50,10.89,11.29,11.67,12.06,12.45,12.83,13.20,13.57,13.94,14.30,14.66,15.02,15.37,15.72,16.07,16.42,16.78,17.14,17.51,17.90,18.31,18.74,19.21,19.73,20.30,20.95,21.70,22.70,24.15,26.45,30.6,41.4];

mug=interp1(Tsat,muvap,T);
