def analyze_eos_flow(flow, **kwargs): work = flow[0] etotals = work.read_etotals(unit="eV") eos_fit = EOS.Birch_Murnaghan().fit(flow.volumes, etotals) return eos_fit.plot(**kwargs)