If all you want is the time average of the global pressure components
then you don’t need per-atom stress and compute reduce. Just:
fix 2 all ave/time 1000 100 200000 c_thermo_pressure file PT.data
Thermo_pressure is defined for you by the thermo command and is
global pressure.
Steve