title='Simulation at constant T, p'

# ------ initial parameter values ------

ntotal=6.02214E23         # total number of particles (1 mol)
n=50000                   # number of simulated particles
mass=4.64951E-26          # mass of particles [kg] (N2 molecule)
box=1.0,1.0,0.25          # initial box size [m]
tempinit=300.0            # initial temperature [K]
seed=3771                 # random number generator seed
dt=2.0E-5                 # time step [s]
wallmass=0                # boxer: mass of movable wall [kg] (0 for rigid wall)
rheat=1000.0              # heater: heating events per particle and unit time [1/s]
rmix=100.0                # mixer: mixing events per particle and unit time [1/s]
dtprint=0.001             # time interval for reporting [s]
dtave=0.001               # time interval for averaging [s]
report='Time,Step,Temperature,Volume,Pressure,ExtPressure,Energy,Work,Heat,Ideality'
                          # quantities to report
#plotfile='start*.grf'     # plot file (* will be replaced by an incrementing counter)

# ------ stage-specific parameter values ------

stage='Equilibrium'
  duration=0.2            # duration [s]
  tempheater=300.0        # heat bath temperature [K] (0 for adiabatic)
  wallmass=0.2            # mass of movable wall [kg] (0 for rigid wall)
  pressext0=9977.4        # external pressure at start of stage [Pa]
  pressext=9977.4         # external pressure at end of stage [Pa]

