1:Issue description  2:Solution **Molecular_virial = spec_lj_virial +_coul_virial + spec_coul_virial + ewald_long_virial + bond_virial + angle_virial + dihedral_virial** **if(shake) virial =Molecular_virial+shake_virial** (a ) virial =special ij +coul  (b)virial =Ewald long  (c) virial =special coul  (d)virial =bond  (e)virial =angle  (f)virial =dihedral   3 : result out   **NOSE_HOOVER**   