A 3D model of doping from the atmosphere for FZ silicon crystal growth is presented. Dopant (phosphorus) transport in gas and melt is modeled. Phosphorus leaves the melt through the free melt surface. A good agreement of the RRV profiles with experiment is obtained. Chemical rate temperature dependence at the free surface is taken into account.