Heat &cntrl imin=0, ntx=1, irest=0, nstlim=1000000, dt=0.002, ntf=2, ntc=2, tempi=0.0, temp0=300.0, ntpr=2000, ntwx=20000, cut=8.0, ntb=1, ntp=0, ntt=3, gamma_ln=2.0, ig=-1, ioutfm=1, ntwv=-1, iwrap=1, /