Changeset 3522 for trunk/fsource/histogram2d.for

Timestamp:

Aug 1, 2018 12:36:58 PM (3 years ago)

Author:

vondreele

Message:

fix array overflow problem in histogram2d.for
clean up some stuff in ImageIntegrate?

File:

• trunk/fsource/histogram2d.for (modified) (1 diff)

trunk/fsource/histogram2d.for

@@ -37 +37 @@
      1    ( Y(K).GE.YLIM(0) .AND. Y(K).LT.YLIM(1) )) THEN
           DDX = (X(K)-XLIM(0))/DX
-          I = INT(DDX)
+          I = MIN(MAX(0,INT(DDX)),NXBINS-1)
           DDY = (Y(K)-YLIM(0))/DY
-          J = INT(DDY)
+          J = MIN(MAX(0,INT(DDY)),NYBINS-1)
           NST(I,J) = NST(I,J)+1.0
           HST(I,J) = HST(I,J)+Z(K)