harmonic1d(100);
system('lmp_daily -in thermal.lammps');
