# -------------------------------------------------------------
# Analysis from File
# -------------------------------------------------------------
configuration = nlread('Bi2Se3_slab.nc', object_id='gID001')[0]

# -------------------------------------------------------------
# Density Of States
# -------------------------------------------------------------
density_of_states = DensityOfStates(
    configuration=configuration,
    kpoints=MonkhorstPackGrid(21,21,1),
    energy_zero_parameter=FermiLevel,
    bands_above_fermi_level=None,
    )
nlsave('dos.nc', density_of_states)
nlprint(density_of_states)