molecule m;
float 	x[7500], f[7500], v[7500];
float 	dgrad, fret, time;
int 	ier, mytaskid, rank, size; 

mytaskid = mpiinit( argc, argv, rank, size );

m = getpdb("mb.pdb");
readparm(m, "mb.prmtop");
mm_options( "nsnb=10, ntpr=20, ntpr_md=20, gb=1, kappa=0.14701, cut=20.0, diel=C, rgbmax=15., rattle=1, tautp=999999.");
mme_init( m, NULL, "::Z", x, NULL);
getxv( "mb.inpcrd", m.natoms, time, x, v );

md(3*m.natoms, 100, x, f, v, mme );
mme_timer();

mpifinalize();