Error in `populating_lvls_example.py`

In the code section corresponding to plotting the populating fraction against temperature in the file `populating_lvls_example.py`, there appears to be a misunderstanding regarding the definition of `level 1`. If I correctly understand, the index should be, `be.data['processed']['pop_lvl'][0,1,0,:,0]/np.sum(be.data['processed']['pop_lvl'][ 0,0,0,:,0], axis=0))` instead of `be.data['processed']['pop_lvl][10,10,0,:,0]/np.sum(be.data['processed']['pop_lvl'][10,:,0,:,0], axis=0))`. I believe the latter refers to `level number 10` instead of `level number 1`. I might be wrong in understanding. A clarification on this is needed.