#running MD in 20ps with restraint on the model
 &cntrl
  imin=0,
  irest=0,
  ntx=1,
  ntb=1,
  cut=10.0,
  ntr=1,
  ntc=2,
  ntf=2,
  tempi=0.0,
  temp0=300.0,
  ntt=3,
  gamma_ln=1.0,
  nstlim=10000, dt=0.002,
  ntpr=100, ntwx=100,ntwr=1000
 /
Keep the model fixed with weak restraints
10.0
RES 1 200
END
END