set outName	        [file rootname [file tail [info script]]]


# MD SETUP

timestep                2.0
numSteps                250000


# INPUT

structure               solvate.psf 
parameters              ../common/par_all22_prot.inp
paraTypeCharmm          on

coordinates             minsolv_0.coor
velocities              minsolv_0.vel 


# OUTPUT

outputenergies          500
outputtiming            500
restartfreq             500

binaryoutput            no
binaryrestart           yes

outputname              [format "%so" $outName]
restartname             $outName


# DCD

DCDfile                 $outName.dcd
DCDfreq                 100
DCDUnitCell             yes


# PBC

cellbasisvector1        24.00  0.00   0.00 
cellbasisvector2         0.00 24.00   0.00 
cellbasisvector3         0.00  0.00  24.00


# PME

PME                     yes
PMETolerance            10e-6
PMEInterpOrder          4

PMEGridSpacing          1.0


# WRAP ATOMS FOR OUTPUT

wrapAll                 on


# CONSTANT-T

langevin                on
langevinTemp            300.0
langevinDamping         1.0


# CONSTANT-P
langevinPiston          on
langevinPistonTarget    1.01325 
langevinPistonPeriod    100.
langevinPistonDecay     50.
langevinPistonTemp      300.0

useGroupPressure        yes 
useFlexibleCell         no


# SPACE PARTITIONING

splitpatch              hydrogen
hgroupcutoff            2.8
stepspercycle           20
margin                  1.0


# CUT-OFFS

switching               on
switchdist               9.0
cutoff                  11.0
pairlistdist            13.0


# RESPA 

fullElectFrequency      2
nonbondedFreq           1


# 1-4 NON-BONDED

exclude                 scaled1-4
1-4scaling              1.0


# COM

commotion               no
zeroMomentum            yes


# SHAKE

rigidbonds              all