Production MD &cntrl imin=0, irest=1, ntx=5, nstlim=10000,dt=0.001, ntc=1, ntpr=20, ntwx=1000, cut=16, ntb=0, igb=1, ntt=3, gamma_ln=1.0, tempi=300.0, temp0=300.0, /