How to couple density distribution function and total energy function in DDF in palabos

Hello everyone, I write two dynamics and D2Q13 descripter, then I need to write a data process to couple density distribution function and total energy function. Could anyone give me some advice?

Thanks, best wishes.