viscair

 VISCAIR: computes viscosity of air 
 vis=VISCAIR(Ta) computes the kinematic viscosity of dry air as a 
 function of air temperature following Andreas (1989), CRREL Report 
 89-11.

 INPUT:   Ta  -  air temperature  [C]

 OUTPUT:  vis  -  air viscosity  [m^2/s]

0001 function  vis=viscair(Ta)
0002 % VISCAIR: computes viscosity of air
0003 % vis=VISCAIR(Ta) computes the kinematic viscosity of dry air as a
0004 % function of air temperature following Andreas (1989), CRREL Report
0005 % 89-11.
0006 %
0007 % INPUT:   Ta  -  air temperature  [C]
0008 %
0009 % OUTPUT:  vis  -  air viscosity  [m^2/s]
0010 
0011 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0012 % 3/8/97: version 1.0
0013 % 8/5/99: version 2.0
0014 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
0015 
0016 vis = 1.326e-5*(1 + 6.542e-3*Ta + 8.301e-6*Ta.^2 - 4.84e-9*Ta.^3);
0017