molecule m;
float 	m_xyz[dynamic], f_xyz[dynamic], v_xyz[dynamic];
float 	dgrad, fret;
int 	ier, i; 

m = getpdb("Na+.pdb");
readparm(m, "Na+.parm7");
allocate m_xyz[3*m.natoms];
allocate f_xyz[3*m.natoms];
allocate v_xyz[3*m.natoms];
mm_options( "e_debug=0, ntpr=1, ntpr_md=1, cut=999.0, dt=0.001, rattle=0, tempi=0, zerov=0");
mm_options( "rism=1, ntpr_rism=1, apply_rism_force=0");
mm_options( "buffer=12, solvcut=999");
mm_options( "progress=1, xvvfile=spc-nacl.xvv");
mme_init( m, NULL, "::Z", m_xyz, NULL);

setxyz_from_mol( m, NULL, m_xyz );

ier = md(3*m.natoms,0,m_xyz,f_xyz,v_xyz, mme );

if(mytaskid==0)
  printf( "md returns %d\n", ier );

mme_timer();