Changeset 451 for trunk/fsource/histogram2d.for
 Timestamp:
 Jan 13, 2012 11:48:53 AM (12 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/fsource/histogram2d.for
r450 r451 30 30 INTEGER*4 I,J,K 31 31 REAL*8 DX,DY 32 REAL*4 DDX,DDY 32 33 DO K=0,N1 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) ) THEN39 C if ( j.ge.500 ) 40 C 1 print *,i,j,nst(i,j),hst(i,j),z(k)34 C if ( mod(k,8000) .eq. 0 ) 35 C 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) 41 42 NST(I,J) = NST(I,J)+1.0 42 43 HST(I,J) = HST(I,J)+Z(K) 44 C if ( mod(k,8000) .eq. 0 ) 45 C 1 print *,i,j,nst(i,j),hst(i,j) 43 46 END IF 44 47 END DO
Note: See TracChangeset
for help on using the changeset viewer.