Please continue to post on this thread to the mail list, so others
Parallel tempering is a bit confusing with respect to what you are asking.
If you weren't doing tempering, but just a bunch of independent
runs at different Temps, then what you have below would make sense,
and you would get 1 dump file for a run at each Temp.
However, when you do tempering, the Temps are swapped between
worlds, and so for the simulation running on one world of processors
(and its dump file), the temperature will change continuously during
the simulation, as will its fix nvt params (these are changed behind
the scenes. So you will confuse yourself if you name it's dump file
snap.Temp.lammpstrj. You will see this by looking at the log
file produced by that proc - the temperature of the run will vary
all over the place as tempering proceeds.
If you want stats on all snapshots at a given Temp, then they
are scattered through all the dump files produced by all the sets
of procs. At any point in time, exactly one snapshot in one of
the files has that Temp, but you will need to look at the hi-level
tempering log file to figure out which one it is, and extract it
in a post-processing step. Like the Pizza.py suggestion I made before.