def print_midplane_tau(self, wavelength):
        for i, disk in enumerate(self.disks):
            if disk.mass > 0.:
                tau = (disk.midplane_cumulative_density(np.array([disk.rmax]))
                       * disk.dust.interp_chi_wav(wavelength))
                print("Disk %i: %.5e" % (i + 1, tau))
dataset/ETHPy150Open hyperion-rt/hyperion/hyperion/model/analytical_yso_model.py/AnalyticalYSOModel.print_midplane_tau/original