Ignore:
Timestamp:
Jan 13, 2012 11:48:53 AM (10 years ago)
Author:
vondreele
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/fsource/histogram2d.for

    r450 r451  
    3030      INTEGER*4   I,J,K
    3131      REAL*8      DX,DY
     32      REAL*4      DDX,DDY
    3233      DO K=0,N-1
    33         I = INT((X(K)-XLIM(0))/DX)
    34         J = INT((Y(K)-YLIM(0))/DY)
    35 C        if ( j.ge.500 )
    36 C     1    print *,k,x(k),y(k),i,j,xlim(0),dx,ylim(0),dy,nxbins,nybins
    37         IF ( (I .GE. 0 .AND. I .LT. NXBINS) .AND.
    38      1     (J .GE. 0 .AND. J .LT. NYBINS) ) THEN
    39 C          if ( j.ge.500 )
    40 C     1      print *,i,j,nst(i,j),hst(i,j),z(k)
     34C        if ( mod(k,8000) .eq. 0 )
     35C     1    print *,k,x(k),xlim,y(k),ylim
     36        IF ( ( X(K).GE.XLIM(0) .AND. X(K).LT.XLIM(1) ) .AND.
     37     1    ( Y(K).GE.YLIM(0) .AND. Y(K).LT.YLIM(1) )) THEN
     38          DDX = (X(K)-XLIM(0))/DX
     39          I = INT(DDX)
     40          DDY = (Y(K)-YLIM(0))/DY
     41          J = INT(DDY)
    4142          NST(I,J) = NST(I,J)+1.0
    4243          HST(I,J) = HST(I,J)+Z(K)
     44C          if ( mod(k,8000) .eq. 0 )
     45C     1      print *,i,j,nst(i,j),hst(i,j)
    4346        END IF
    4447      END DO
Note: See TracChangeset for help on using the changeset viewer.