<div>Hi everyone,</div>
<div>&nbsp;</div>
<div>I use umbrella sample and g_wham to get potential of mean force. However, I also want to compute the mean force. How can I deal with this? Can I get it through calculating the average of data in pullf.xvg file?</div>

<div>&nbsp;</div>
<div>Thanks in advance.</div>
<div>&nbsp;</div>
<div>Wende</div>
<div>&nbsp;</div>
<div>------------------</div>
<div>;spring point is equal to ref_group0+pull_init+time*pull_rate*pull_vec£¬</div>
<div>pull&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Umbrella&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;umbrella,constraint <br>pull_geometry&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Distance&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; ;distance ,direction,cylinder,position<br>pull_dim&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = Y Y Y<br>pull_group0&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = DEN1<br>
pull_group1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = DEN2<br>pull_nstxout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 500<br>pull_nstfout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 50<br>pull_init1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 5.7 <br>pull_k1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; = 1000</div>
<div>-------------------<br clear="all"></div>
<div># This file was created Thu Feb&nbsp; 4 00:51:24 2010<br># by the following command:<br># mdrun_mpi -rdd 1.4 -s npt.tpr -c npt.gro <br>#<br># mdrun_mpi is part of G R O M A C S:<br>#<br># Groningen Machine for Chemical Simulation<br>
#<br>@&nbsp;&nbsp;&nbsp; title &quot;Pull force&quot;<br>@&nbsp;&nbsp;&nbsp; xaxis&nbsp; label &quot;Time (ps)&quot;<br>@&nbsp;&nbsp;&nbsp; yaxis&nbsp; label &quot;Force (kJ/mol/nm)&quot;<br>@TYPE xy<br>0.0000&nbsp;-28.8555<br>1.6000&nbsp;-19.2897<br>3.2000&nbsp;-18.59<br>4.8000&nbsp;-20.7866<br>
6.4000&nbsp;-26.1605<br>8.0000&nbsp;-42.8572<br>9.6000&nbsp;-50.7524<br>11.2000&nbsp;-73.9647<br>12.8000&nbsp;-74.2521<br>14.4000&nbsp;-102.284<br>16.0000&nbsp;-96.5184</div>
<div>&nbsp;</div>
<div><br>-- <br>wende<br></div>