title               =  ezq_drg position restraining
warnings            =  10
cpp                 =  /lib/cpp
define              =  -DPOSRES
constraints         =  all-bonds
integrator          =  md
dt                  =  0.002	; ps !
nsteps              =  10000	; total 20.0 ps.
nstcomm             =  1
nstxout             =  250      ; collect data every 0.5 ps
nstvout             =  1000
nstfout             =  0
nstlog              =  10
nstenergy           =  10
nstlist             =  10
ns_type             =  grid
rlist               =  0.9
coulombtype         =  PME
rcoulomb            =  0.9
rvdw                =  0.9
fourierspacing      =  0.12
fourier_nx          =  0
fourier_ny          =  0
fourier_nz          =  0
pme_order           =  4
ewald_rtol          =  le-5
optimize_fft        =  yes
; Berendsen temperature coupling is on in five groups
Tcoupl              =  berendsen
tau_t               =  0.1                0.1         0.1   0.1    0.1
tc-grps		    =  Protein            Ca          sol   RPR    Cl
ref_t               =  300                300         300   300    300
; Pressure coupling is on
Pcoupl              =  berendsen
pcoupltype          =  isotropic
tau_p               =  1.0
compressibility     =  4.5e-5
ref_p               =  1.0
; Generate velocites is on at 300 K.
gen_vel             =  yes
gen_temp            =  300.0
gen_seed            =  173529