--- MITgcm/eesupp/src/different_multiple.F	2001/09/21 03:54:34	1.5
+++ MITgcm/eesupp/src/different_multiple.F	2005/05/15 02:58:40	1.8
@@ -1,18 +1,20 @@
-C     $Header: /home/ubuntu/mnt/e9_copy/MITgcm/eesupp/src/different_multiple.F,v 1.5 2001/09/21 03:54:34 cnh Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/eesupp/src/different_multiple.F,v 1.8 2005/05/15 02:58:40 jmc Exp $
 C $Name:  $
+
 #include "CPP_EEOPTIONS.h"
 
 CBOP
 C     !ROUTINE: DIFFERENT_MULTIPLE
 
 C     !INTERFACE:
-      LOGICAL FUNCTION DIFFERENT_MULTIPLE( freq, val1, val2 )
+      LOGICAL FUNCTION DIFFERENT_MULTIPLE( freq, val1, step )
       IMPLICIT NONE
 
 C     !DESCRIPTION:
 C     *==========================================================*
-C     | LOGICAL FUNCTION DIFFERENT_MULTIPLE                       
-C     | o Checks two numbers multiple of a third number.          
+C     | LOGICAL FUNCTION DIFFERENT\_MULTIPLE                       
+C     | o Checks if a multiple of freq exist
+C     |   around val1 +/- step/2
 C     *==========================================================*
 C     | This routine is used for diagnostic and other periodic    
 C     | operations. It is very sensitive to arithmetic precision. 
@@ -24,43 +26,43 @@
 
 C     !INPUT PARAMETERS:
 C     == Routine arguments ==
-C     val1, val2 :: Two times that are checked 
-C     freq       :: Frequency by which times are divided.
-      _RL  freq, val1, val2
+C     freq       :: Frequency by which time is divided.
+C     val1       :: time that is checked 
+C     step       :: length of time interval (around val1) that is checked 
+      _RL  freq, val1, step
+
+C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
 
 C     !LOCAL VARIABLES:
 C     == Local variables ==
-C     f              :: Temp. for holding freq
 C     v1, v2, v3, v4 :: Temp. for holding time
 C     d1, d2, d3     :: Temp. for hold difference
-C     step           :: Temp. for hold difference used as increment
-      _RL  f, v1, v2, v3, v4, d1, d2, d3, step
+      _RL  v1, v2, v3, v4, d1, d2, d3
 CEOP
 
 C     o Do easy cases first.
       DIFFERENT_MULTIPLE = .FALSE.
 
       IF ( freq .NE. 0. ) THEN
-        IF ( ABS(val1-val2) .GT. freq ) THEN
+        IF ( ABS(step) .GT. freq ) THEN
          DIFFERENT_MULTIPLE = .TRUE.
         ELSE
 
 C         o This case is more complex because of round-off error
-          f = freq
           v1 = val1
-          v2 = val2
-          step = v1-v2
+          v2 = val1 - step
+          v3 = val1 + step
 
 C         Test v1 to see if its a "closest multiple"
-          v3 = v1 + step
-          v4 = NINT(v1/f)*f
+          v4 = NINT(v1/freq)*freq
           d1 = v1-v4
           d2 = v2-v4
           d3 = v3-v4
           IF ( ABS(d1) .LE. ABS(d2) .AND. ABS(d1) .LE. ABS(d3) )
      &        DIFFERENT_MULTIPLE = .TRUE.
 
-        ENDIF ! |val1-val2| > freq
-      ENDIF ! freq != 0
+        ENDIF
+      ENDIF
 
+      RETURN
       END