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?