In the current workflows I've created these manually to do the rerouting of flows like so:
yakima_topo_scaled = yakima_topo.copy(deep=True)
yakima_topo_scaled['Contrib_Area'] = yakima_topo_scaled['Length'].copy()
yakima_topo_scaled['Contrib_Area'].values[:] = SEC_PER_YEAR * M_PER_KM
yakima_topo_scaled.to_netcdf('../topologies/yakima_topo_scaled_area.nc')
Ideally this should be done behind the scenes since it's always this format (up to a constant for different time steps).