Dear Colleague,

Two bugs in the Dissipative Particle Dynamics code in DL_MESO 2.3 have been 
discovered:

1. The search_list routines in bond_module.f90 and bond_module_ser.f90 do not 
   detect duplicate entries in the global/local particle number array lblclst 
   for every single case. This is due to a bug in the logic to discover if there 
   are any duplicates before or after the entry found by the routine. This bug 
   may have caused DL_MESO_DPD to not find the permutation of local particle 
   numbers giving the minimum distance and thus give overly long bond lengths.

2. The serial utility exportimage_ser.f90 as distributed with DL_MESO produces 
   an "End of file" error message, due to the READ statement in line 107 
   attempting to read data from the wrong input channel.

The first bug may be fixed by modifying lines 250-254 in bond_module.f90 and 
lines 118-122 in bond_module_ser.f90 to the following:

      IF (ind>1) THEN
        IF (lblclst (ind-1, 1)==aim) ldupl = .true.
      END IF
      IF (ind<nlist) THEN
        IF (lblclst (ind+1, 1)==aim) ldupl = .true.
      END IF

Note that this bug is fixed by separating the conditions for searching both 
before and after the array index discovered by the binary search.

The second bug can be fixed by changing line 107 in exportimage_ser.f90 to:

        READ (nsave+j-1) beads, bonds, angles, diheds, told, tsold, vold

These bug fixes have been applied to the current DL_MESO release and a corrected 
version may be downloaded by registered users of version 2.3 without re-
registering and decrypted using the same password.

Michael Seaton
2010-08-25