// Creation Date : Wed May 28 17:57:15 1997
// Last Modified : Tue Aug  5 15:43:17 1997
// $Id$

#ifndef _SPEC_H
#define _SPEC_H
//=========================
// User specification file
//=========================

//_APPLICATION: Appear as the title of output
#define _APPLICATION "Test"

//_2DIM: Defined if only use 2 dimension (Only support 2d now!!)
#define _2DIM

//_NP: Number of particles.
//_NSTEP: Number of timesteps (include rescale if any)
//_NDUMP: #of timestesp to write output to `stdout'
#define _NP 2000
#define _NSTEP 12000000
#define _NDUMP 10000

//_NRESCALE: Defined if want to rescale temperature for certain timesteps
//           default behaviour of rescale is bulk rescale to _TR,
//           this behaviour can be changed, see below user part
#define _NRESCALE 2000000

//_NPROC: Number of processes
#define _NPROC 4

//_LCFG: Defined if want to load config file, from _LCFG_READ,
//       if _LCFG_READ is not defined, default is "config.dat"
//       if the file does not exist, program will do normal initialization
//       just as _LCFG is not defined.
// This program always write configuration to output file _LCFG_WRITE
//       if _LCFG_WRITE is not defined, default is "config.dat"
//#define _LCFG 
//#define _LCFG_READ "config.dat"

//_RPM: define if RPM is enabled, value is the probability
//    (should be between 0 and 1, use 1 to maximize flow rate)
//_RPMAX, _RPMIN: the bin boudarys which act as RPM plane(parallel to y axis)
//     _RPMIN must be small than _RPMAX, the boundary is inclusive, i.e.
//     _RPMIN==1, _RPMAX==2 means two RPM planes:xbin0/xbin1 and xbin1/xbin2
#define _RPM 1
#define _RPMAX 1
#define _RPMIN 0

//_NNL: Neighbor list size,
//_NBL: Bin list size, NB:don't modify them unless you know what you are doing
#define _NBL 80
#define _NNL 50

//_DT, _RC: time step in reduced time, potential cutoff in reduced length
#define _DT 0.0008
#define _RC 2.5

//_TR: reduced temperature
//_CUBEX,_CUBEY: geometry, only PBC now (easy to add free, bounce, etc.)
#define _TR 1.8
#define _CUBEX 150.0
#define _CUBEY 50.0

//_PROP_MESH: defined if want mesh output, to file
//   now only include density, Vx, Vy.(_PROP_MESH_VX,_VY not used now)
//_MESHX,_MESHY: number of meshes on geometry
#define _PROP_MESH "mesh.out"
#define _PROP_MESH_VX
#define _PROP_MESH_VY
#define _MESHX (53*3*3)
#define _MESHY (53*3)

//_PROP_POSITION: defined if want position output.
#define _PROP_POSITION "pos.out"

//_PROP_PARTICLE: defined if want particle data output
//_NCORR: no correlation after _NCORR timesteps(depend on _DT)
#define _PROP_PARTICLE "particle.dat"
#define _NCORR 5000

//_CR: defined to enable center wall
//_CR_X0,etc: geometry of wall(inclusive), every virtaul wall particle sits on
//       interger points,
//       e.g.:X0=1,X1=3;Y0=1,Y1=2 means 6 wall particles at
//            [1,1][1,2][2,1][2,2][3,1][3,2](in reduced length)
#define _CR
#define _CR_X0 74
#define _CR_X1 76
#define _CR_Y0 15
#define _CR_Y1 35

//-------------------------
// Available user functions
//-------------------------
// _USER_RESCALE
// _USER_INIT
// _USER_FINISH
// _USER_STEP
// _USER_SINGLE
// _USER_PAIR

#endif // _SPEC_H