- Timestamp:
- May 3, 2011 1:06:10 PM (14 years ago)
- Location:
- trunk
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified trunk/fsource/histogram2d.for ΒΆ
r271 r274 30 30 INTEGER*4 I,J,K 31 31 REAL*8 DX,DY 32 33 32 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) 43 39 END IF 44 40 END DO
Note: See TracChangeset
for help on using the changeset viewer.