/[MITgcm]/MITgcm/eesupp/src/different_multiple.F
ViewVC logotype

Diff of /MITgcm/eesupp/src/different_multiple.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by cnh, Tue Nov 3 15:28:03 1998 UTC revision 1.3 by adcroft, Mon May 24 14:57:55 1999 UTC
# Line 22  CEndofinterface Line 22  CEndofinterface
22    
23  C     o Do easy cases first.  C     o Do easy cases first.
24        DIFFERENT_MULTIPLE = .FALSE.        DIFFERENT_MULTIPLE = .FALSE.
25        IF ( freq .EQ. 0. ) RETURN  
26        IF ( ABS(val1-val2) .GT. freq ) THEN        IF ( freq .NE. 0. ) THEN
27         DIFFERENT_MULTIPLE = .TRUE.          IF ( ABS(val1-val2) .GT. freq ) THEN
28         RETURN           DIFFERENT_MULTIPLE = .TRUE.
29        ENDIF          ELSE
30    
31  C     o This case is more complex because of round-off error  C         o This case is more complex because of round-off error
32        f = freq            f = freq
33        v1 = val1            v1 = val1
34        v2 = val2            v2 = val2
35        step = v1-v2            step = v1-v2
36    
37  C     Test v1 to see if its a "closest multiple"  C         Test v1 to see if its a "closest multiple"
38        v3 = v1 + step            v3 = v1 + step
39        v4 = NINT(v1/f)*f            v4 = NINT(v1/f)*f
40        d1 = v1-v4            d1 = v1-v4
41        d2 = v2-v4            d2 = v2-v4
42        d3 = v3-v4            d3 = v3-v4
43        IF ( ABS(d1) .LE. ABS(d2) .AND.            IF ( ABS(d1) .LE. ABS(d2) .AND. ABS(d1) .LE. ABS(d3) )
44       &     ABS(d1) .LE. ABS(d3) ) DIFFERENT_MULTIPLE = .TRUE.       &        DIFFERENT_MULTIPLE = .TRUE.
45    
46            ENDIF ! |val1-val2| > freq
47          ENDIF ! freq != 0
48    
49        END        END
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22