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

<html>

  <head>

    <meta content="text/html; charset=ISO-8859-1"

      http-equiv="Content-Type">

  </head>

  <body text="#000000" bgcolor="#ffffff">

    On 21/05/2011 10:45 AM, Elisabeth wrote:

    <blockquote

      cite="mid:BANLkTinO7fTwpeB1JiEcASj75pih25R=ng@mail.gmail.com"

      type="cite">Dear&nbsp; experts,<br>

      <br>

      I have a box of pure hexane with density of around 50 SI (size

      7nm). With the settings below I get error: <br>

      <br>

      vol 0.65&nbsp; imb F&nbsp; 3% step 42100, will finish Fri May 20 18:54:05

      2011<br>

      Warning: 1-4 interaction between 1246 and 1249 at distance 3.580

      which is larger than the 1-4 table size 2.250 nm<br>

      These are ignored for the rest of the simulation<br>

      This usually means your system is exploding,<br>

      if not, you should increase table-extension in your mdp file<br>

      or with user tables increase the table size<br>

      Warning: 1-4 interaction between 1249 and 1252 at distance 3.693

      which is larger than the 1-4 table size 2.250 nm<br>

      These are ignored for the rest of the simulation<br>

      This usually means your system is exploding,<br>

      if not, you should increase table-extension in your mdp file<br>

      or with user tables increase the table size<br>

      <br>

      A list of missing interactions:<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Bond of&nbsp;&nbsp; 2375 missing&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Angle of&nbsp;&nbsp; 4500 missing&nbsp;&nbsp;&nbsp;&nbsp; 11<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ryckaert-Bell. of&nbsp;&nbsp; 5625 missing&nbsp;&nbsp;&nbsp;&nbsp; 21<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ-14 of&nbsp;&nbsp; 5625 missing&nbsp;&nbsp;&nbsp;&nbsp; 15<br>

      <br>

      Molecule type 'Hexane'<br>

      the first 10 missing interactions, except for exclusions:<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ryckaert-Bell. atoms&nbsp;&nbsp;&nbsp; 5&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp; 11&nbsp;&nbsp; 14 global&nbsp; 1245&nbsp; 1248&nbsp;

      1251&nbsp; 1254<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ-14 atoms&nbsp;&nbsp;&nbsp; 5&nbsp;&nbsp; 14&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; global&nbsp; 1245&nbsp; 1254<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Angle atoms&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp; 11&nbsp;&nbsp; 14&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; global&nbsp; 1248&nbsp; 1251&nbsp;

      1254<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ryckaert-Bell. atoms&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp; 11&nbsp;&nbsp; 14&nbsp;&nbsp; 15 global&nbsp; 1248&nbsp; 1251&nbsp;

      1254&nbsp; 1255<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ryckaert-Bell. atoms&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp; 11&nbsp;&nbsp; 14&nbsp;&nbsp; 16 global&nbsp; 1248&nbsp; 1251&nbsp;

      1254&nbsp; 1256<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ryckaert-Bell. atoms&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp; 11&nbsp;&nbsp; 14&nbsp;&nbsp; 17 global&nbsp; 1248&nbsp; 1251&nbsp;

      1254&nbsp; 1257<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ-14 atoms&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp; 15&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; global&nbsp; 1248&nbsp; 1255<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ-14 atoms&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp; 16&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; global&nbsp; 1248&nbsp; 1256<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; Ryckaert-Bell. atoms&nbsp;&nbsp;&nbsp; 9&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp; 11&nbsp;&nbsp; 14 global&nbsp; 1249&nbsp; 1248&nbsp;

      1251&nbsp; 1254<br>

      &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; LJ-14 atoms&nbsp;&nbsp;&nbsp; 9&nbsp;&nbsp; 14&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; global&nbsp; 1249&nbsp; 1254<br>

      <br>

      -------------------------------------------------------<br>

      Program mdrun_mpi, VERSION 4.5.4<br>

      Source code file: domdec_top.c, line: 356<br>

      <br>

      Fatal error:<br>

      49 of the 18125 bonded interactions could not be calculated

      because some atoms involved moved further apart than the

      multi-body cut-off distance (1.25 nm) or the two-body cut-off

      distance (1.25 nm), see option -rdd, for pairs and tabulated bonds

      also see option -ddcheck<br>

      For more information and tips for troubleshooting, please check

      the GROMACS<br>

      website at <a moz-do-not-send="true"

        href="http://www.gromacs.org/Documentation/Errors">http://www.gromacs.org/Documentation/Errors</a><br>

    </blockquote>

    <br>

    OK, your system blew up.<br>

    <br>

    <blockquote

      cite="mid:BANLkTinO7fTwpeB1JiEcASj75pih25R=ng@mail.gmail.com"

      type="cite"><br>

      mdp file:<br>

      <br>

      ;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; Bonds<br>

      constraints&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp; = none&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>

      constraint-algorithm = lincs&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; <br>

      <br>

      ;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; Run control&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br>

      integrator&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; md&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>

      dt&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.001 &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp; <br>

      nsteps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 500000 ;5000 &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>

      nstcomm&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 100&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp; <br>

      <br>

      ;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; Output control<br>

      nstenergy&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 100 &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;; frequency to write

      energies to energy file. i.e., energies and other statistical data

      are stored every 10 steps<br>

      nstxout&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 100;1&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; <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; 0<br>

      nstlog&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1000&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;; frequency to write

      energies to log file<br>

      nstxtcout&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; =&nbsp; 1000&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br>

      <br>

      ;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; Neighbor searching<br>

      nstlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 10&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br>

      ns_type&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; grid&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br>

      <br>

      ;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; Electrostatics/VdW <br>

      coulombtype&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; Shift&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp; <br>

      vdw-type&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; Shift&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp; <br>

      rcoulomb-switch&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br>

      rvdw-switch&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.9 ;0&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>

      <br>

      ;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; Cut-offs<br>

      rlist&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1.25&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp; <br>

      rcoulomb&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1.0&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; <br>

      rvdw&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1.0&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br>

      <br>

      ;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; Temperature coupling&nbsp;&nbsp;&nbsp; <br>

      Tcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; v-rescale&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br>

      tc-grps&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; System <br>

      tau_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 0.1&nbsp;&nbsp; <br>

      ref_t&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 300&nbsp;&nbsp;&nbsp; <br>

      <br>

      ;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; Pressure coupling<br>

      Pcoupl&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; Parrinello-Rahman&nbsp;&nbsp;&nbsp; &nbsp;&nbsp; &nbsp;&nbsp; <br>

      Pcoupltype&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp; =&nbsp; isotropic&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>

      tau_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; <br>

      compressibility&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 4.5e-5 4.5e-5&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br>

      ref_p&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 10&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br>

      <br>

      ;&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; Velocity generation&nbsp;&nbsp;&nbsp; &nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br>

      gen_vel&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; yes&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp;&nbsp; <br>

      gen_temp&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 300.0&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; <br>

      gen_seed&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; =&nbsp; 173529&nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp;&nbsp; &nbsp;&nbsp; <br>

    </blockquote>

    <br>

    Generating velocities and compressing from an unequilibrated initial

    conformation in the same run is asking for trouble, and sometimes

    you're observing it. See <a

href="http://www.gromacs.org/Documentation/How-tos/Steps_to_Perform_a_Simulation">http://www.gromacs.org/Documentation/How-tos/Steps_to_Perform_a_Simulation</a>.

    Equilibrate first. Then slowly change the conditions so that you're

    never all that far from equilibrium, and reach your destination.<br>

    <br>

    <blockquote

      cite="mid:BANLkTinO7fTwpeB1JiEcASj75pih25R=ng@mail.gmail.com"

      type="cite">When I use constraints for all-bonds and dt of 0.002

      compressing the box works well and I get density of above 600 SI

      which is the actual density of alkane. Can any one help why

      setting no constraints and dt=0.001 leads to the error message

      shown above. </blockquote>

    <br>

    You probably got lucky the first time and not the second. Or, your

    atom-atom bond interactions are no good.<br>

    <br>

    <blockquote

      cite="mid:BANLkTinO7fTwpeB1JiEcASj75pih25R=ng@mail.gmail.com"

      type="cite">When no constraints is used PR doesnt work for

      compressing</blockquote>

    <br>

    So you want to restrain the positions of heavy atoms, but compress a

    pure system? This is contradictory. How do you want to decrease the

    volume if the atoms aren't allowed to move much?<br>

    <br>

    <blockquote

      cite="mid:BANLkTinO7fTwpeB1JiEcASj75pih25R=ng@mail.gmail.com"

      type="cite">and berendsen gives desired densty only for high

      pressures above 50 bar. I need to study my system at different

      pressures from 10 bar... please comment on this.. Thanks.<br>

      <br>

      <p class="MsoNormal"><span style="">10 bar</span></p>

      <p class="MsoNormal" style=""><span style="font-family:

          &quot;Courier New&quot;;">;<span style="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>Pressure

          coupling</span></p>

      <p class="MsoNormal" style=""><span style="font-family:

          &quot;Courier New&quot;;">Pcoupl<span style="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>=<span

            style="">&nbsp;

          </span>berendsen</span></p>

      <p class="MsoNormal" style=""><span style="font-family:

          &quot;Courier New&quot;;">Pcoupltype<span style="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span><span

            style="">&nbsp; </span>=<span style="">&nbsp;

          </span>isotropic<span style="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span><span

            style="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</span></span></p>

      <p class="MsoNormal" style=""><span style="font-family:

          &quot;Courier New&quot;;">tau_p<span style="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>=

          <b style=""><span style="">&nbsp;</span>0.1 </b><span style="">&nbsp;&nbsp;

          </span>0.5<span style="">&nbsp;&nbsp;&nbsp; </span><span style="">&nbsp;&nbsp;&nbsp;&nbsp; </span><span

            style="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;</span></span></p>

      <p class="MsoNormal" style=""><span style="font-family:

          &quot;Courier New&quot;;">compressibility<span style="">&nbsp;&nbsp;&nbsp;&nbsp; </span>=<span

            style="">&nbsp;

          </span>4.5e-5 4.5e-5<span style="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span> </span></p>

      <p class="MsoNormal"><span style="font-family: &quot;Courier

          New&quot;;">ref_p<span style="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span>=<span

            style="">&nbsp; </span>10<span style="">&nbsp;&nbsp;&nbsp;&nbsp; </span>10<span

            style="">&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; </span></span><span style=""></span></p>

      <p class="MsoNormal"><span style="">&nbsp;</span></p>

      <p class="MsoNormal"><span style="">does not compress!

          With and without constr.</span></p>

    </blockquote>

    <br>

    We can't tell.<br>

    <br>

    Mark<br>

  </body>

</html>