Spatial temperature profile

Hi LAMMPS developers,

I would like to know whether it is possible in LAMMPS to bin the domain, subtract out the average velocity in each bin and compute the temperature in each bin and output a file just like fix ave/spatial. Compute temp/profile is very close to what I want but unfortunately does not output the temperature for each bin. Thank you for your help.

Marcus

You can use fix ave/spatial with an atom-style
variable. If that variable calculates the KE of
each atom, then summing and averaging
it on a per-bin basis as fix ave/spatial does
is effectively temperature. I think you could
also write the variable to subtract out a streaming
velocity from a compute, such as compute temp/profile.

Steve