Paul Richard Thomas wrote: + real :: fnres(numclusters) ! change to REAL and it works!! That's a good start to 2008! The comment is wrong and should read: + real :: fnres(numclusters) ! change to INTEGER and it works!! Paul