Last change
on this file since 111 was
111,
checked in by vondreel, 13 years ago
|
modified for refactoring of integration
|
File size:
1.1 KB
|
Line | |
---|
1 | SUBROUTINE HISTOGRAM2D(N,X,Y,Z,NXBINS,NYBINS,XLIM,YLIM,DX,DY, |
---|
2 | 1 NST,HST) |
---|
3 | |
---|
4 | Cf2py intent(in) N |
---|
5 | Cf2py intent(in) X |
---|
6 | Cf2py depend(N) X |
---|
7 | Cf2py intent(in) Y |
---|
8 | Cf2py depend(N) Y |
---|
9 | Cf2py intent(in) Z |
---|
10 | Cf2py depend(N) Z |
---|
11 | Cf2py intent(in) NXBINS |
---|
12 | Cf2py intent(in) NYBINS |
---|
13 | Cf2py intent(in) XLIM |
---|
14 | Cf2py intent(in) YLIM |
---|
15 | Cf2py intent(in) DX |
---|
16 | Cf2py intent(in) DY |
---|
17 | Cf2py intent(in,out) NST |
---|
18 | Cf2py depend(NXBINS,NYBINS) NST |
---|
19 | Cf2py intent(in,out) HST |
---|
20 | Cf2py depend(NXBINS,NYBINS) HST |
---|
21 | |
---|
22 | IMPLICIT NONE |
---|
23 | INTEGER*4 N |
---|
24 | REAL*4 X(0:N-1),Y(0:N-1),Z(0:N-1) |
---|
25 | INTEGER*4 NXBINS,NYBINS |
---|
26 | REAL*8 XLIM(0:1),YLIM(0:1) |
---|
27 | REAL*8 DX,DY |
---|
28 | INTEGER*8 NST(0:NXBINS-1,0:NYBINS-1) |
---|
29 | REAL*4 HST(0:NXBINS-1,0:NYBINS-1) |
---|
30 | |
---|
31 | INTEGER*4 I,J,K |
---|
32 | |
---|
33 | 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+0.5)-1 |
---|
37 | J = NINT((Y(K)-YLIM(0))/DY+0.5)-1 |
---|
38 | NST(I,J) = NST(I,J)+1 |
---|
39 | HST(I,J) = HST(I,J)+Z(K) |
---|
40 | END IF |
---|
41 | END DO |
---|
42 | RETURN |
---|
43 | END |
---|
Note: See
TracBrowser
for help on using the repository browser.