We model the L1551 IRS5 source as a system containing two protostars, each surrounded by a circumstellar disk, both encircled by a circumbinary disk, and all disks surrounded by an extended infalling flattened envelope With this composite model, we can calculate self-consistently the spectral energy distribution of the source the shape of the ice and silicate features, and the spatial intensity distributions of the envelope and disks. We compare our model results with the observations, determining the physical parameters of the disks and the envelope. We find that flattened envelope collapse models are required in order to explain the observations and that the infall rate of the envelope is much larger than the accretion rates of the binary disks.
