Ignore:
Timestamp:
May 3, 2011 1:06:10 PM (11 years ago)
Author:
vondreele
Message:

fix (again!) histogram2d

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/fsource/histogram2d.for

    r271 r274  
    3030      INTEGER*4   I,J,K
    3131      REAL*8      DX,DY
    32 
    3332      DO K=0,N-1
    34         IF ( ( X(K) .GE. XLIM(0) .AND. X(K) .LE. XLIM(1)) .AND.
    35      1    (Y(K) .GE. YLIM(0) .AND. Y(K). LE. YLIM(1)) ) THEN
    36           I = NINT((X(K)-XLIM(0))/DX)
    37           J = NINT((Y(K)-YLIM(0))/DY)
    38           IF ( (I .GE. 0 .AND. I .LT. NXBINS) .AND.
    39      1       (J .GE. 0 .AND. J .LT. NYBINS) ) THEN
    40             NST(I,J) = NST(I,J)+1.0
    41             HST(I,J) = HST(I,J)+Z(K)
    42           END IF
     33        I = INT((X(K)-XLIM(0))/DX)
     34        J = INT((Y(K)-YLIM(0))/DY)
     35        IF ( (I .GE. 0 .AND. I .LT. NXBINS) .AND.
     36     1     (J .GE. 0 .AND. J .LT. NYBINS) ) THEN
     37          NST(I,J) = NST(I,J)+1.0
     38          HST(I,J) = HST(I,J)+Z(K)
    4339        END IF
    4440      END DO
Note: See TracChangeset for help on using the changeset viewer.