Chris, Thanks - we've added the fix to the main code in checkpoint33. The new impldiff.F (different to your version) should be easily insertible into your code; it uses 3d arrays (a,b,c) rather then the 2d arrays in your version but other wise is identical. A.