The below lines related to dumping the heat flux components is this correct? and also if this is correct what are the units of this heat flux? please help me out

compute myKE all ke/atom

compute myPE all pe/atom

compute myStress all stress/atom NULL

compute myTemp all temp

compute flux all heat/flux myKE myPE myStress

variable Jx equal c_flux[1]/9000 #9000 is the cross sectional area(yz)

variable Jy equal c_flux[2]/9000 #9000 is the cross sectional area(xz)

variable Jz equal c_flux[3]/3600 #9000 is the cross sectional area(xy)

fix JJ all ave/correlate 1 1000 1000 c_flux[1] c_flux[2] c_flux[3] type auto file J0Jt.dat ave running

If you have the correct syntax, you will find out when you run this.

If this does what you want it to do, only you can tell.

What the units of compute heat/flux are is explained in its documentation.

How many times did I have to remind you to check the documentation so far?

This is becoming a bit of a pattern.

Donâ€™t you think this is a sign that you may need to change your approach?