<P>
&nbsp; Dear Xavier and Jochen,<BR>
<BR>
Thank you for the suggestions: it works well now.<BR>
<BR>
Jes.<BR>
<BR>
<BR>
On Mon, 27 Oct 2008 Xavier Periole wrote :<BR>
&gt;On 27 Oct 2008 06:36:09 -0000<BR>
&gt;&nbsp; &quot;JMandumpal&quot; &lt;jesbman@rediffmail.com&gt; wrote:<BR>
&gt;&gt;&nbsp; Dear GROMACS users,<BR>
&gt;&gt;<BR>
&gt;&gt;I have got a problem while running a simulated anealing simulations, using GROMACS 3.3.3 version, of pure water:<BR>
&gt;&gt;<BR>
&gt;&gt;the task:<BR>
&gt;&gt;<BR>
&gt;&gt;I want to cool the system (water- 500 molecules) from 360 K to 0K in 100ps time ( fast cooling). For that I used the following script:<BR>
&gt;&gt;<BR>
&gt;&gt;title&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = Water annealing simulation<BR>
&gt;&gt;cpp&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = /usr/bin/cpp<BR>
&gt;&gt;; RUN CONTROL PARAMETERS<BR>
&gt;&gt;integrator&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = md<BR>
&gt;&gt;; Start time and timestep in ps<BR>
&gt;&gt;dt&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 0.001; 1fs<BR>
&gt;&gt;nsteps&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 100000&nbsp; ; 1ps<BR>
&gt;&gt;; For exact run continuation or redoing part of a run<BR>
&gt;&gt;init_step&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0<BR>
&gt;&gt;; number of steps for center of mass motion removal<BR>
&gt;&gt;nstcomm&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1<BR>
&gt;&gt;; group(s) for center of mass motion removal<BR>
&gt;&gt;comm-grps&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = SOL<BR>
&gt;&gt;; OUTPUT CONTROL OPTIONS<BR>
&gt;&gt;; Output frequency for coords (x), velocities (v) and forces (f)<BR>
&gt;&gt;nstxout&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 10000 ; 10ps<BR>
&gt;&gt;nstvout&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 10000 ; 10ps<BR>
&gt;&gt;nstfout&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0<BR>
&gt;&gt;; Checkpointing helps you continue after crashes<BR>
&gt;&gt;nstcheckpoint&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 100<BR>
&gt;&gt;; Output frequency for energies to log file and energy file<BR>
&gt;&gt;nstlog&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 0<BR>
&gt;&gt;nstenergy&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1000 ; 1ps<BR>
&gt;&gt;; Output frequency and precision for xtc file<BR>
&gt;&gt;;nstxtcout&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0<BR>
&gt;&gt;;xtc-precision&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1000<BR>
&gt;&gt;; select multiple groups. By default all atoms will be written.<BR>
&gt;&gt;; NEIGHBORSEARCHING PARAMETERS<BR>
&gt;&gt;nstlist&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 10<BR>
&gt;&gt;; ns algorithm (simple or grid)<BR>
&gt;&gt;ns_type&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = grid<BR>
&gt;&gt;; Periodic boundary conditions:<BR>
&gt;&gt;pbc&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = xyz<BR>
&gt;&gt;; nblist cut-off<BR>
&gt;&gt;rlist&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1.0<BR>
&gt;&gt;; OPTIONS FOR ELECTROSTATICS AND VDW<BR>
&gt;&gt;; Method for doing electrostatics<BR>
&gt;&gt;coulombtype&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = pme<BR>
&gt;&gt;rcoulomb&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 1.0<BR>
&gt;&gt;ewald_rtol&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 1.0e-5&nbsp; ; since erfc(sigma*rcutoff)=ewald_rtol<BR>
&gt;&gt;optimize_fft&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = yes<BR>
&gt;&gt;; Relative dielectric constant for the medium<BR>
&gt;&gt;epsilon_r&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 1&nbsp; ; for water<BR>
&gt;&gt;; Method for doing Van der Waals<BR>
&gt;&gt;vdw-type&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = cut-off<BR>
&gt;&gt;rvdw&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 1.0<BR>
&gt;&gt;; Spacing for the PME/PPPM FFT grid<BR>
&gt;&gt;fourierspacing&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 0.12 ; check for the range 0.1 to 0.15<BR>
&gt;&gt;; EWALD/PME/PPPM parameters<BR>
&gt;&gt;pme_order&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 4<BR>
&gt;&gt;ewald_geometry&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 3d<BR>
&gt;&gt;epsilon_surface&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0<BR>
&gt;&gt;; OPTIONS FOR WEAK COUPLING ALGORITHMS<BR>
&gt;&gt;; OPTIONS FOR WEAK COUPLING ALGORITHMS<BR>
&gt;in what follows you have commented out the temperature coupling details.<BR>
&gt;I have never used SA but this for sure will not help the program controling<BR>
&gt;the temperature! Put these back on should help! :))<BR>
&gt;&gt;; Temperature coupling<BR>
&gt;&gt;;Tcoupl&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = berendsen<BR>
&gt;&gt;; Groups to couple separately<BR>
&gt;&gt;;tc-grps&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = system<BR>
&gt;&gt;; Time constant (ps) and reference temperature (K)<BR>
&gt;&gt;;tau-t&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 0.1<BR>
&gt;&gt;;ref-t&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 360<BR>
&gt;&gt;; simulated anealing<BR>
&gt;&gt;annealing&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = single<BR>
&gt;&gt;annealing_npoints&nbsp; &nbsp; &nbsp; &nbsp; = 2<BR>
&gt;&gt;annealing_time&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 0&nbsp;  0<BR>
&gt;&gt;annealing_temp&nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 360 0<BR>
&gt;&gt;; GENERATE VELOCITIES FOR STARTUP RUN<BR>
&gt;&gt;gen_vel&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = yes<BR>
&gt;&gt;gen-temp&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 360<BR>
&gt;&gt;gen-seed&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  = 173529<BR>
&gt;&gt;; OPTIONS FOR BONDS<BR>
&gt;&gt;constraints&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = all-bonds<BR>
&gt;&gt;constraint_algorithm&nbsp; &nbsp;  = shake<BR>
&gt;&gt;shake_tol&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = 0.0001<BR>
&gt;&gt;morse&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; = no<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  -----------------------------------------------------------------------<BR>
&gt;&gt;<BR>
&gt;&gt;Problem:<BR>
&gt;&gt;<BR>
&gt;&gt;Having done the simulations, the output shows ( ener.edr , Temperature)&nbsp; a high temperature value around 850K. Clearly I did something wrong in my parameter file. <BR>
&gt;&gt;-------------------------------------------------------<BR>
&gt;&gt;the temperature extracted after the simulations<BR>
&gt;&gt;-------------------------------------------------------<BR>
&gt;&gt;s0 legend &quot;Temperature&quot;<BR>
&gt;&gt;&nbsp; &nbsp; 0.000000&nbsp; 376.058044<BR>
&gt;&gt;&nbsp; &nbsp; 1.000000&nbsp; 864.033020<BR>
&gt;&gt;&nbsp; &nbsp; 2.000000&nbsp; 873.607605<BR>
&gt;&gt;&nbsp; &nbsp; 3.000000&nbsp; 851.308838<BR>
&gt;&gt;&nbsp; &nbsp; 4.000000&nbsp; 847.223877<BR>
&gt;&gt;&nbsp; &nbsp; 5.000000&nbsp; 860.511841<BR>
&gt;&gt;&nbsp; &nbsp; 6.000000&nbsp; 813.547241<BR>
&gt;&gt;&nbsp; &nbsp; 7.000000&nbsp; 845.261658<BR>
&gt;&gt;&nbsp; &nbsp; 8.000000&nbsp; 872.090027<BR>
&gt;&gt;&nbsp; &nbsp; 9.000000&nbsp; 852.079468<BR>
&gt;&gt;&nbsp;  10.000000&nbsp; 849.152466<BR>
&gt;&gt;&nbsp;  11.000001&nbsp; 850.055237<BR>
&gt;&gt;&nbsp;  12.000001&nbsp; 834.468323<BR>
&gt;&gt;&nbsp;  13.000001&nbsp; 861.107300<BR>
&gt;&gt;&nbsp;  14.000001&nbsp; 828.380493<BR>
&gt;&gt;&nbsp;  15.000001&nbsp; 850.166809<BR>
&gt;&gt;&nbsp;  16.000000&nbsp; 850.382019<BR>
&gt;&gt;&nbsp;  17.000000&nbsp; 847.214478<BR>
&gt;&gt;&nbsp;  18.000000&nbsp; 851.201538<BR>
&gt;&gt;&nbsp;  19.000000&nbsp; 857.533691<BR>
&gt;&gt;&nbsp;  20.000000&nbsp; 849.760071<BR>
&gt;&gt;&nbsp;  21.000002&nbsp; 864.907715<BR>
&gt;&gt;&nbsp;  22.000002&nbsp; 837.307495<BR>
&gt;&gt;&nbsp;  23.000002&nbsp; 827.407593<BR>
&gt;&gt;&nbsp;  24.000002&nbsp; 841.464111<BR>
&gt;&gt;&nbsp;  25.000002&nbsp; 863.330688<BR>
&gt;&gt;&nbsp;  26.000002&nbsp; 841.096069<BR>
&gt;&gt;&nbsp;  27.000002&nbsp; 856.469910<BR>
&gt;&gt;&nbsp;  28.000002&nbsp; 854.449768<BR>
&gt;&gt;&nbsp;  29.000002&nbsp; 845.534973<BR>
&gt;&gt;&nbsp;  30.000002&nbsp; 853.421448<BR>
&gt;&gt;&nbsp;  31.000002&nbsp; 825.198059<BR>
&gt;&gt;&nbsp;  32.000000&nbsp; 857.317566<BR>
&gt;&gt;&nbsp;  33.000000&nbsp; 823.510925<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  34.000000 858.875122<BR>
&gt;&gt;&nbsp;  35.000000&nbsp; 880.710266<BR>
&gt;&gt;&nbsp;  36.000000&nbsp; 853.692749<BR>
&gt;&gt;&nbsp;  37.000000&nbsp; 847.079468<BR>
&gt;&gt;&nbsp;  38.000000&nbsp; 851.004822<BR>
&gt;&gt;&nbsp;  39.000000&nbsp; 858.750000<BR>
&gt;&gt;&nbsp;  40.000000&nbsp; 872.787048<BR>
&gt;&gt;&nbsp;  41.000004&nbsp; 843.033813<BR>
&gt;&gt;&nbsp;  42.000004&nbsp; 874.984009<BR>
&gt;&gt;&nbsp;  43.000004&nbsp; 845.272583<BR>
&gt;&gt;&nbsp;  44.000004&nbsp; 822.993713<BR>
&gt;&gt;&nbsp;  45.000004&nbsp; 863.116455<BR>
&gt;&gt;&nbsp;  46.000004&nbsp; 865.640198<BR>
&gt;&gt;&nbsp;  47.000004&nbsp; 873.489563<BR>
&gt;&gt;&nbsp;  48.000004&nbsp; 849.170776<BR>
&gt;&gt;&nbsp;  49.000004&nbsp; 861.631958<BR>
&gt;&gt;&nbsp;  50.000004&nbsp; 860.999939<BR>
&gt;&gt;&nbsp;  51.000004&nbsp; 834.341003<BR>
&gt;&gt;&nbsp;  52.000004&nbsp; 864.601196<BR>
&gt;&gt;&nbsp;  53.000004&nbsp; 838.268372<BR>
&gt;&gt;&nbsp;  54.000004&nbsp; 842.860046<BR>
&gt;&gt;&nbsp;  55.000004&nbsp; 837.590332<BR>
&gt;&gt;&nbsp;  56.000004&nbsp; 855.985046<BR>
&gt;&gt;&nbsp;  57.000004&nbsp; 836.168030<BR>
&gt;&gt;&nbsp;  58.000004&nbsp; 862.147644<BR>
&gt;&gt;&nbsp;  59.000004&nbsp; 869.183044<BR>
&gt;&gt;&nbsp;  60.000004&nbsp; 873.057556<BR>
&gt;&gt;&nbsp;  61.000004&nbsp; 871.051636<BR>
&gt;&gt;&nbsp;  62.000004&nbsp; 857.647217<BR>
&gt;&gt;&nbsp;  63.000004&nbsp; 857.834106<BR>
&gt;&gt;&nbsp;  64.000000&nbsp; 871.380981<BR>
&gt;&gt;&nbsp;  65.000000&nbsp; 850.045654<BR>
&gt;&gt;&nbsp;  66.000000&nbsp; 880.142456<BR>
&gt;&gt;&nbsp;  67.000000&nbsp; 844.183350<BR>
&gt;&gt;&nbsp;  68.000000&nbsp; 850.980042<BR>
&gt;&gt;&nbsp;  69.000000&nbsp; 854.431091<BR>
&gt;&gt;&nbsp;  70.000000&nbsp; 847.898743<BR>
&gt;&gt;&nbsp;  71.000000&nbsp; 859.750488<BR>
&gt;&gt;&nbsp;  72.000000&nbsp; 852.219910<BR>
&gt;&gt;&nbsp;  73.000000&nbsp; 854.833984<BR>
&gt;&gt;&nbsp;  74.000000&nbsp; 851.676514<BR>
&gt;&gt;&nbsp;  75.000000&nbsp; 866.519287<BR>
&gt;&gt;&nbsp;  76.000000&nbsp; 861.120605<BR>
&gt;&gt;&nbsp;  77.000000&nbsp; 857.873230<BR>
&gt;&gt;&nbsp;  78.000000&nbsp; 845.287720<BR>
&gt;&gt;&nbsp;  79.000000&nbsp; 874.603333<BR>
&gt;&gt;&nbsp;  80.000000&nbsp; 854.244751<BR>
&gt;&gt;&nbsp;  81.000008&nbsp; 864.049561<BR>
&gt;&gt;&nbsp;  82.000008&nbsp; 892.515991<BR>
&gt;&gt;&nbsp;  83.000008&nbsp; 855.164734<BR>
&gt;&gt;&nbsp;  84.000008&nbsp; 867.054077<BR>
&gt;&gt;&nbsp;  85.000008&nbsp; 862.560059<BR>
&gt;&gt;&nbsp;  86.000008&nbsp; 829.368469<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; 87.000008 864.729553<BR>
&gt;&gt;&nbsp;  88.000008&nbsp; 882.954102<BR>
&gt;&gt;&nbsp;  89.000008&nbsp; 865.048401<BR>
&gt;&gt;&nbsp;  90.000008&nbsp; 864.581848<BR>
&gt;&gt;&nbsp;  91.000008&nbsp; 870.596619<BR>
&gt;&gt;&nbsp;  92.000008&nbsp; 856.474976<BR>
&gt;&gt;&nbsp;  93.000008&nbsp; 861.473450<BR>
&gt;&gt;&nbsp;  94.000008&nbsp; 838.902344<BR>
&gt;&gt;&nbsp;  95.000008&nbsp; 896.907715<BR>
&gt;&gt;&nbsp;  96.000008&nbsp; 865.042480<BR>
&gt;&gt;&nbsp;  97.000008&nbsp; 877.761292<BR>
&gt;&gt;&nbsp;  98.000008&nbsp; 877.356750<BR>
&gt;&gt;&nbsp;  99.000008&nbsp; 855.421753<BR>
&gt;&gt;&nbsp; 100.000008&nbsp; 863.720398<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  ------------------------------------------------------<BR>
&gt;&gt;<BR>
&gt;&gt;Could anyone tell me what changes I should make in the parameter file?<BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;Many thanks,<BR>
&gt;&gt;Jes.<BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;<BR>
&gt;&gt;&nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp; &nbsp;  <BR>
&gt;<BR>
&gt;-----------------------------------------------------<BR>
&gt;XAvier Periole - PhD<BR>
&gt;<BR>
&gt;- Molecular Dynamics Group -<BR>
&gt;NMR and Computation<BR>
&gt;University of Groningen<BR>
&gt;The Netherlands<BR>
&gt;-----------------------------------------------------<BR>

</P>
<br><br>
<Table border=0 Width=644 Height=57 cellspacing=0 cellpadding=0 style='font-family:Verdana;font-size:11px;line-height:15px;'><TR><td><a href='http://adworks.rediff.com/cgi-bin/AdWorks/click.cgi/www.rediff.com/signature-home.htm/1050715198@Middle5/2606998_2599290/2602379/1?PARTNER=3&OAS_QUERY=null' target=new ><img src ='http://imadworks.rediff.com/cgi-bin/AdWorks/adimage.cgi/2606998_2599290/creative_2602379.gif'  alt='Black578x38_banner2.gif'  border=0></a></td></TR></Table>