function retval = ahvs(T,ro)
retval=(0.2+0.45*T*ro)*exp(-0.18*T*ro);