Changeset 451 for trunk/fsource/histogram2d.for

Timestamp:

Jan 13, 2012 11:48:53 AM (12 years ago)

Author:

vondreele

Message:

 

File:

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

trunk/fsource/histogram2d.for

@@ -30 +30 @@
       INTEGER*4   I,J,K
       REAL*8      DX,DY
+      REAL*4      DDX,DDY
       DO K=0,N-1
-        I = INT((X(K)-XLIM(0))/DX)
-        J = INT((Y(K)-YLIM(0))/DY)
-C        if ( j.ge.500 )
-C     1    print *,k,x(k),y(k),i,j,xlim(0),dx,ylim(0),dy,nxbins,nybins
-        IF ( (I .GE. 0 .AND. I .LT. NXBINS) .AND.
-     1     (J .GE. 0 .AND. J .LT. NYBINS) ) THEN
-C          if ( j.ge.500 )
-C     1      print *,i,j,nst(i,j),hst(i,j),z(k)
+C        if ( mod(k,8000) .eq. 0 )
+C     1    print *,k,x(k),xlim,y(k),ylim
+        IF ( ( X(K).GE.XLIM(0) .AND. X(K).LT.XLIM(1) ) .AND.
+     1    ( Y(K).GE.YLIM(0) .AND. Y(K).LT.YLIM(1) )) THEN
+          DDX = (X(K)-XLIM(0))/DX
+          I = INT(DDX)
+          DDY = (Y(K)-YLIM(0))/DY
+          J = INT(DDY)
           NST(I,J) = NST(I,J)+1.0
           HST(I,J) = HST(I,J)+Z(K)
+C          if ( mod(k,8000) .eq. 0 )
+C     1      print *,i,j,nst(i,j),hst(i,j)
         END IF
       END DO