DistAutoCorrelation/dist2.f at main · wangdi2016/DistAutoCorrelation · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
      PROGRAM DistAutoCorrelation
!     Declare variables

      implicit none

!     Declare variables
!     PARAMETER (N=2666)
!     PARAMETER (M=2666)
!     PARAMETER (L=1000)
      REAL, DIMENSION(10000) :: W1, W2
      INTEGER, DIMENSION(10000) :: P1, P2
      REAL, DIMENSION(2000) :: OCCURENCE
      INTEGER :: I, J, IMAX, N, M, L, DIST, STATUS, NC, MC
      CHARACTER(len=50) :: FILENAME1, FILENAME2, FILENAME3
      INTEGER :: FILE_UNIT1, FILE_UNIT2, FILE_UNIT3, ierr

      FILE_UNIT1=10
      FILE_UNIT2=20
      FILE_UNIT3=30

      N=10000
      M=10000
      L=1200

!     set NC and MC counters for counting lines in file1 and file2
      NC=0
      MC=0

!     set IMAX = L; L is the number of lines for output
      IMAX=L

!     Prompt the user to enter the input filename1
      WRITE(*,*) 'Enter the filename1: '
      READ(*,*) FILENAME1

!     Prompt the user to enter the input filename2
      WRITE(*,*) 'Enter the filename2: '
      READ(*,*) FILENAME2

!     Prompt the user to enter the output filename
      WRITE(*,*) 'Enter the output filename: '
      READ(*,*) FILENAME3

!CCCCCCCCCCCCCCCCCCCCCCCCCCC
!     Open the file1
      OPEN(FILE_UNIT1, FILE=FILENAME1, STATUS='OLD', IOSTAT=STATUS)
      IF (STATUS.NE.0) THEN
        WRITE(*,*) 'Error: could not open file'
        STOP
      END IF

!     Skip header line
      READ(FILE_UNIT1,*)
!     Read the data from the file into the array
      DO I=1,N
!     set NC to count lines read
        NC=NC+1
!       READ(FILE_UNIT1,*,IOSTAT=STATUS) X(I,1), X(I,2)
        READ(FILE_UNIT1,*,IOSTAT=STATUS) P1(I), W1(I)
        IF (STATUS.NE.0) EXIT
      END DO
!     Close the file
      CLOSE(FILE_UNIT1)
!CCCCCCCCCCCCCCCCCCCCCCCCCCC
!     Open the file2
      OPEN(FILE_UNIT2, FILE=FILENAME2, STATUS='OLD', IOSTAT=STATUS)
      IF (STATUS.NE.0) THEN
        WRITE(*,*) 'Error: could not open file'
        STOP
      END IF

!     Skip header line
      READ(FILE_UNIT2,*)
!     Read the data from the file into the array
      DO I=1,M
!     set MC to count lines read
        MC=MC+1
!       READ(FILE_UNIT2,*,IOSTAT=STATUS) Y(I,1), Y(I,2)
        READ(FILE_UNIT2,*,IOSTAT=STATUS) P2(I), W2(I)
        IF (STATUS.NE.0) EXIT
      END DO
!     Close the file
      CLOSE(FILE_UNIT2)
!CCCCCCCCCCCCCCCCCCCCCCCCCCC

!     Print the contents of the array
      WRITE(*,*) 'The contents of the array are:'
      DO I=1,MC-1
!       WRITE(*,*) X(I,1), X(I,2), Y(I,1), Y(I,2)
        WRITE(*,*) P1(I), W1(I), P2(I), W2(I)
      END DO

      WRITE(*,*) NC-1, MC-1

!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
!C   Calculate Distance and Occurence              CCC
!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC

! Initialize Occurence = 0
      DO I=1,IMAX

         OCCURENCE(I)=0.0

      END DO

!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
!C         Calculate Distance                      CCC
!CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
      DO I=1,NC-1
         DO J=1,MC-1

!           DIST = Y(J,1) - X(I,1) + 1
            DIST = P2(J) - P1(I) + 1

!           WRITE(*,*) DIST

            IF ((DIST .GE. 0) .AND. (DIST .LE. IMAX)) THEN
!                OCCURENCE(DIST) = OCCURENCE(DIST) + X(I,2) * Y(J,2)
                 OCCURENCE(DIST) = OCCURENCE(DIST) + W1(I) * W2(J)
            END IF

         END DO
      END DO

!C
!C Write output array
!C
      OPEN(FILE_UNIT3, FILE=FILENAME3, status='replace', action='write')
      WRITE(FILE_UNIT3,*) "Dist      DAC"
      DO I=1,IMAX
        WRITE(FILE_UNIT3,'(I5, F9.3)') I, OCCURENCE(I)
!       WRITE(FILE_UNIT3,*) I, OCCURENCE(I)
      END DO

      END PROGRAM DistAutoCorrelation