
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">


<HTML><HEAD>


<META http-equiv=Content-Type content="text/html; charset=gb2312">


<META content="MSHTML 6.00.2900.5969" name=GENERATOR>


<STYLE></STYLE>


</HEAD>


<BODY bgColor=#c5e4cd>


<DIV><FONT size=2><FONT size=3>Dear all,</FONT><BR>&nbsp;&nbsp; </FONT>I am 


trying&nbsp; to analysis the force of a atom.So I change "fr.x[n][XX]" to 


"fr.f[n][XX]" in the template.c .But it turns up "Segmentation fault" .(The 


pointer "fr.f" is NULL) How could I get the force?</DIV>


<DIV>&nbsp;</DIV>


<DIV>My mdp file:</DIV>


<DIV><FONT size=2></FONT>&nbsp;</DIV>


<DIV><FONT 


size=2>title&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


Yo<BR>cpp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


/usr/bin/cpp<BR>constraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


all-bonds<BR>integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


md<BR>dt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 0.002&nbsp;; ps 


!<BR>nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 2500000&nbsp;; total 5ns.<BR>comm_mode&nbsp;&nbsp;&nbsp;&nbsp; = 


NonE<BR>nstcomm&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


<BR>nstxout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


50<BR>nstvout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


0<BR>nstfout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


50<BR>nstxtcout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


0&nbsp;<BR>nstlog&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


5000<BR>nstenergy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


5000<BR>nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


10<BR>ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


grid<BR>coulombtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=PME<BR>rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


1.4<BR>rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


1.4<BR>rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


1.4<BR>pbc&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; xyz<BR>; Berendsen temperature coupling is on in two 


groups<BR>Tcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


berendsen<BR>tc-grps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


 =&nbsp; CNT 


SOL&nbsp;<BR>tau_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 0.1&nbsp;0.1 </FONT></DIV>


<DIV><FONT 


size=2>ref_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 300 300 <BR>; Energy 


monitoring<BR>energygrps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; CNT SOL<BR>; Non-equilibrium MD</FONT></DIV>


<DIV>&nbsp;</DIV>


<DIV><FONT size=2>freezegrps&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp;CNP&nbsp;&nbsp;<BR>freezedim&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; Y Y 


Y&nbsp;&nbsp;&nbsp; <BR>; Isotropic pressure coupling is now 


on<BR>Pcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; no<BR>;Pcoupltype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


= 


isotropic<BR>;tau_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 0.5<BR>;compressibility&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 


4.5e-5<BR>;ref_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 1.0<BR>; Generate velocites is off at 300 


K.<BR>gen_vel&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


no<BR>gen_temp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 


300.0<BR>gen_seed&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 


=&nbsp; 537129</FONT></DIV>


<DIV><FONT size=2></FONT>&nbsp;</DIV>


<DIV><FONT size=2></FONT>&nbsp;</DIV>


<DIV><FONT size=2><FONT size=3>Many thanks in advance,<BR>Chuan</FONT></DIV>


<DIV><BR></DIV></FONT></BODY></HTML>