Changeset 111


Ignore:
Timestamp:
Jul 16, 2010 10:10:40 AM (11 years ago)
Author:
vondreel
Message:

modified for refactoring of integration

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)
    33
    44Cf2py intent(in) N
     
    1313Cf2py intent(in) XLIM
    1414Cf2py intent(in) YLIM
     15Cf2py intent(in) DX
     16Cf2py intent(in) DY
    1517Cf2py intent(in,out) NST
    1618Cf2py depend(NXBINS,NYBINS) NST
    1719Cf2py intent(in,out) HST
    1820Cf2py depend(NXBINS,NYBINS) HST
    19 Cf2py intent(in,out) HSTX
    20 Cf2py depend(NXBINS) HSTX
    21 Cf2py intent(in,out) HSTY
    22 Cf2py depend(NYBINS) HSTY
    2321
    2422      IMPLICIT NONE
     
    2725      INTEGER*4   NXBINS,NYBINS
    2826      REAL*8      XLIM(0:1),YLIM(0:1)
     27      REAL*8      DX,DY
    2928      INTEGER*8   NST(0:NXBINS-1,0:NYBINS-1)
    3029      REAL*4      HST(0:NXBINS-1,0:NYBINS-1)
    31       REAL*8      HSTX(0:NXBINS),HSTY(0:NYBINS)
    3230
    3331      INTEGER*4   I,J,K
    34       REAL*8      DX,DY
    3532
    36       DX = (XLIM(1)-XLIM(0))/FLOAT(NXBINS)
    37       DY = (YLIM(1)-YLIM(0))/FLOAT(NYBINS)
    38 
    39       DO I=0,NXBINS
    40         HSTX(I) = XLIM(0)+FLOAT(I)*DX
    41       END DO
    42       DO J=0,NYBINS
    43         HSTY(J) = YLIM(0)+FLOAT(J)*DY
    44       END DO
    45      
    4633      DO K=0,N-1
    4734        IF ( ( X(K) .GE. XLIM(0) .AND. X(K) .LE. XLIM(1)) .AND.
Note: See TracChangeset for help on using the changeset viewer.