 equil_md
  &cntrl
  nmropt = 1, ntx = 5, irest = 1, ntrx = 1,
  ntxo = 1, ntr = 1,
  ntpr = 1000, ntwx =1000,
  ntf = 2, ntb = 2, dielc = 1.0, ntc=2,
  cut = 12, nsnb = 10,
  imin = 0,
  nstlim = 10000, dt = 0.002,
  temp0 = 300, tempi = 300,
  ntt = 1,
  ntp = 1,
  &end
  &wt type='END'
  &end
  LISTOUT= rst_st1.lis
  DISANG= rst_st1.dat
 PHOSPHOROUS RESTRAIN
 20.0
 FIND
 P1 * * *
 SEARCH
 RES 190
 END
 END
