Changeset 111
 Timestamp:
 Jul 16, 2010 10:10:40 AM (11 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/fsource/histogram2d.for
r106 r111 1 SUBROUTINE HISTOGRAM2D(N,X,Y,Z,NXBINS,NYBINS,XLIM,YLIM, 2 1 NST,HST ,HSTX,HSTY)1 SUBROUTINE HISTOGRAM2D(N,X,Y,Z,NXBINS,NYBINS,XLIM,YLIM,DX,DY, 2 1 NST,HST) 3 3 4 4 Cf2py intent(in) N … … 13 13 Cf2py intent(in) XLIM 14 14 Cf2py intent(in) YLIM 15 Cf2py intent(in) DX 16 Cf2py intent(in) DY 15 17 Cf2py intent(in,out) NST 16 18 Cf2py depend(NXBINS,NYBINS) NST 17 19 Cf2py intent(in,out) HST 18 20 Cf2py depend(NXBINS,NYBINS) HST 19 Cf2py intent(in,out) HSTX20 Cf2py depend(NXBINS) HSTX21 Cf2py intent(in,out) HSTY22 Cf2py depend(NYBINS) HSTY23 21 24 22 IMPLICIT NONE … … 27 25 INTEGER*4 NXBINS,NYBINS 28 26 REAL*8 XLIM(0:1),YLIM(0:1) 27 REAL*8 DX,DY 29 28 INTEGER*8 NST(0:NXBINS1,0:NYBINS1) 30 29 REAL*4 HST(0:NXBINS1,0:NYBINS1) 31 REAL*8 HSTX(0:NXBINS),HSTY(0:NYBINS)32 30 33 31 INTEGER*4 I,J,K 34 REAL*8 DX,DY35 32 36 DX = (XLIM(1)XLIM(0))/FLOAT(NXBINS)37 DY = (YLIM(1)YLIM(0))/FLOAT(NYBINS)38 39 DO I=0,NXBINS40 HSTX(I) = XLIM(0)+FLOAT(I)*DX41 END DO42 DO J=0,NYBINS43 HSTY(J) = YLIM(0)+FLOAT(J)*DY44 END DO45 46 33 DO K=0,N1 47 34 IF ( ( X(K) .GE. XLIM(0) .AND. X(K) .LE. XLIM(1)) .AND.
Note: See TracChangeset
for help on using the changeset viewer.