timestep             2.0
numsteps            250000000 # 500 000 ps
firsttimestep       200000000

useConstantRatio      yes
useflexiblecell      yes
set temperatur       310

structure            system.psf
parameters           lib/par_all27_prot_lipid_c36.prm
parameters           lib/chol_new.str 
paraTypeCharmm       on

# INPUT
#binvelocities        prod2.vel
velocities           prod4.vel
coordinates          prod4.coor
#bincoordinates       prod2.coor.old
extendedsystem       prod4.xsc
#extendedsystem       equil3.xsc


# OUTPUT

outputenergies       500
outputtiming         500
outputpressure       500
binaryoutput         no

outputname           prod5
restartname          prod5
restartfreq          5000
binaryrestart        yes

XSTFreq              5000

dcdfile              prod5.dcd
dcdfreq              5000

# CUT-OF
splitpatch           hydrogen
hgroupcutoff         2.8
stepspercycle        8
switching            on
vdwForceSwitching on
switchdist           11.0
cutoff               12.0
pairlistdist         14.0
#LJcorrection         no

# CONSTANT-T
langevin             on
langevintemp         $temperatur
langevindamping      1.0

# CONSTANT-P
langevinpiston       on
langevinpistontarget 1
langevinpistonperiod 200
langevinpistondecay  50
langevinpistontemp   $temperatur

strainrate           0.  0.  0.
usegrouppressure     yes

# PME
PME                  yes
PMETolerance         10e-6
PMEInterpOrder       4

PMEGridSpacing 1.0

# MULTIPLE TIME-STEP
fullelectfrequency   2
nonbondedfreq        1


# SHAKE/RATTLE
rigidbonds           all
rigidtolerance       0.00001
rigiditerations      400


# 1-4's
exclude              scaled1-4
1-4scaling           1.0