 Input file for pmd 
&cntrl
  imin = 0, irest = 1, igb = 1, dielc = 1.0,
  ntx  = 5, ntb   = 0, ntt = 3, gamma_ln = 1.0, ntp = 0,
  temp0 = 300.0, tautp = 2.0,
  tol = 0.00001, ntc = 2, ntf = 2,
  cut = 32.0, ibelly = 0, dt = 0.001,
  ntpr = 100,
  ntwr = 100,
  ntwx = 100,
  nstlim = 100000,
  nscm = 100,
  ntr = 0, ipol = 0
&end

ncsu_pmd

  output_file = 'pmd.XXX.txt'
  output_freq = 100
  
  variable
    type = TORSION
    i=( 5, 7, 9, 15 )
    anchor_position = YYY
    anchor_strength = 100.0
  end variable

  variable
    type = TORSION
    i=( 7, 9, 15, 17 )
    anchor_position = ZZZ
    anchor_strength = 100.0
  end variable

end ncsu_pmd
