Bugfix ra calculation, put limit on minimum wind speed (also in surface layer), put limit on Obukhovlength iteration