  We revisit the vertical structure of neutrino dominated accretion flows
(NDAFs) in spherical coordinates with a new boundary condition based on the
mechanical equilibrium. The solutions show that NDAF is significantly thick.
The Bernoulli parameter and neutrino trapping are determined by the mass
accretion rate and the viscosity parameter. According to the distribution of
the Bernoulli parameter, the possible outflow may appear in the outer region of
the disk. The neutrino trapping can essentially affect the neutrino radiation
luminosity. The vertical structure of NDAF is like a "sandwich", and the
multilayer accretion may account for the flares in gamma-ray bursts.
