function zo = cd2zo( cd, z ) % CD2Z0 - Convert Drag coefficient for height z to zo % zo = cd2zo( cd, z ) vk = .408; zo = z ./ exp( vk ./ sqrt(cd ) );