Contact information for heat/gran/conduction

Hello,

I have a LIGGGHTS model with a fix heat/gran/conduction. I want to output the conduction contact area information using the dump command. The following properties are mentioned in the online documentation - https://www.cfdem.com/media/DEM/docu/fix_heat_gran_conduction.html

"You can visualize the heat sources by accessing f_heatSource[0], and the heatFluxes by f_heatFlux[0] . With f_directionalHeatFlux[0], f_directionalHeatFlux[1] and f_directionalHeatFlux[2] you can access the conductive heat fluxes in x,y,z directions. The conductive heat fluxes are calculated per-contact and half the value is stored in each atom participating in the contact. With f_Temp[0] you can access the per-particle temperature. You can also access the total thermal energy of the fix group (useful for the thermo command) with f_id .

In case the store_contact_data option is set to yes, additional output of per-atom properties is available. The property ids are nContactsConduction, contactAreaConduction and wallHeatFlux which represent the number of contacts a particle has, the combined contact area with other particles and the wall heat flux, respectively. These ids can be used as usual to output them in a dump command or to calculate atom properties."

But when I try to output them in the following format:

fix heattransfer all heat/gran/conduction initial_temperature 320 &
contact_area overlap area_correction yes store_contact_data yes

dump atomwise all custom 50000 post/atoms_*.vtk f_Temp[0] f_heatFlux[0] &
f_directionalHeatFlux[0] f_directionalHeatFlux[1] f_directionalHeatFlux[2] &
f_nContactsConduction[0] f_contactAreaConduction[0] f_wallHeatFlux[0]

I get the following error: ERROR: Dump particle fix does not compute per-atom vector (../dump_particle.cpp:596). I am not sure how to output the stored contact area information. Any help is appreciated. Thank you.