source: palm/trunk/UTIL/compare_palm_logs.f90 @ 253

 PROGRAM compare_palm_logs

!------------------------------------------------------------------------------!
! Actual revisions:
! -----------------
! 
!
! Former revisions:
! -----------------
!
! Description:
! ------------
! This routine compares the log files from two different PALM runs.
!
! This routine must be compiled with:
! decalpha:
!    f95 -cpp -fast -r8
! IBM-Regatta:
!    xlf95 -qsuffix=cpp=f90 -qrealsize=8 -q64 -qmaxmem=-1 -Q -O3
! IMUK:
!    ifort compare...f90 -o compare...x
!    -cpp -axW -r8 -nbs -Vaxlib
! NEC-SX6:
!    sxf90 compare...f90 -o compare...x 
!    -C hopt -Wf '-A idbl4'
!------------------------------------------------------------------------------!

    IMPLICIT NONE

!
!-- Local variables
    CHARACTER (LEN=5) ::  id_char
    CHARACTER (LEN=80), DIMENSION(2)  ::  directory, log_message
    CHARACTER (LEN=100), DIMENSION(2) ::  filename

    INTEGER ::  count=0, i, id, i1(2), i2(2), j, j1(2), j2(2), k, k1(2), k2(2), &
                n_err, n_files(2)

    LOGICAL ::  found

    REAL    ::  simtime(2)

    INTEGER, DIMENSION(:,:),   ALLOCATABLE ::  array_2d_i_1, array_2d_i_2

    REAL, DIMENSION(:,:),   ALLOCATABLE ::  array_2d_1, array_2d_2
    REAL, DIMENSION(:,:,:), ALLOCATABLE ::  array_1, array_2

    directory(1) = 'const_log/'
    directory(2) = 'const_log.1/'

!
!-- Check, if file from PE0 exists on directory 1. Stop, if it does not exist.
    n_files(1) = 0

    WRITE (id_char,'(''_'',I4.4)')  n_files(1)
    INQUIRE ( FILE=TRIM( directory(1) )//id_char, EXIST=found )
!
!-- Find out the number of files (equal to the number of PEs which
!-- have been used in PALM) and open them
    DO  WHILE ( found )

       OPEN ( n_files(1)+100, FILE=TRIM( directory(1) )//id_char, &
              FORM='UNFORMATTED' )
       n_files(1) = n_files(1) + 1
       WRITE (id_char,'(''_'',I4.4)')  n_files(1)
       INQUIRE ( FILE=TRIM( directory(1) )//id_char, EXIST=found )

    ENDDO

    IF ( n_files(1) == 0 )  THEN
       PRINT*, '+++ no file _0000 in directory "', TRIM( directory(1) ), '"'
       STOP
    ELSE
       PRINT*, '*** directory "', TRIM( directory(1) ), '": ', n_files(1), &
               ' files found'
    ENDIF

!
!-- Same for the second directory
    n_files(2) = 0

    WRITE (id_char,'(''_'',I4.4)')  n_files(2)
    INQUIRE ( FILE=TRIM( directory(2) )//id_char, EXIST=found )

    DO  WHILE ( found )

       OPEN ( n_files(2)+200, FILE=TRIM( directory(2) )//id_char, &
              FORM='UNFORMATTED' )
       n_files(2) = n_files(2) + 1
       WRITE (id_char,'(''_'',I4.4)')  n_files(2)
       INQUIRE ( FILE=TRIM( directory(2) )//id_char, EXIST=found )

    ENDDO

!
!-- Number of files must be identical
    IF ( n_files(1) /= n_files(2) )  THEN
       PRINT*, '+++ file number mismatch'
       PRINT*, '    ', TRIM( directory(1) ), ': ', n_files(1), ' files'
       PRINT*, '    ', TRIM( directory(2) ), ': ', n_files(2), ' files'
       STOP
    ENDIF

!
!-- Compare the data file by file
    DO  id = 0, n_files(1)-1

       count = 0

       WRITE (filename(1),'(A,''_'',I4.4)')  TRIM( directory(1) ), id
       WRITE (filename(2),'(A,''_'',I4.4)')  TRIM( directory(2) ), id

       PRINT*, '*** comparing files "', TRIM( filename(1) ),'" "', &
               TRIM( filename(2) ), '"'
       DO
          PRINT*,' '
          READ ( id+100, END=100 )  log_message(1)
          PRINT*,'    --- ', TRIM( log_message(1) )
          READ ( id+200, END=900 )  log_message(2)

          IF ( TRIM( log_message(1) ) /= TRIM( log_message(2) ) )  THEN
             PRINT*,'    +++ log message on file 2 does not match:'
             PRINT*,'        ', TRIM( log_message(2) )
          ENDIF

          count = count + 1
          IF ( log_message(1)(1:2) == '3d' )  THEN
             PRINT*,'    *** reading 3d array'
             READ ( id+100, END=901 )  simtime(1), i1(1), i2(1), j1(1), j2(1), &
                                       k1(1), k2(1)
             PRINT*,'        time=', simtime(1)
             PRINT*,'        array size=(',i1(1),':',i2(1), &
                                       ',',j1(1),':',j2(1),',',k1(1),':',k2(1),')'
             READ ( id+200, END=902 )  simtime(2), i1(2), i2(2), j1(2), j2(2), &
                                       k1(2), k2(2)
             IF ( simtime(1) /= simtime(2) .OR. i1(1) /= i1(2) .OR. &
                  i2(1) /= i2(2) .OR. j1(1) /= j1(2) .OR. j2(1) /= j2(2) .OR. &
                  k1(1) /= k1(2) .OR. k2(1) /= k2(2) )  THEN
                PRINT*,'    +++ time/indices on file 2 does not match:'
                PRINT*,'        time=', simtime(2)
                PRINT*,'        array size=(',i1(2),':', &
                                i2(2), ',',j1(2),':',j2(2),',',k1(2),':',k2(2),')'
                STOP
             ENDIF

             ALLOCATE( array_1(i1(1):i2(1),j1(1):j2(1),k1(1):k2(1)), &
                       array_2(i1(2):i2(2),j1(2):j2(2),k1(2):k2(2)) )

             READ ( id+100, END=903 )  array_1
             READ ( id+200, END=904 )  array_2

             n_err = 0
loop:        DO  k = k1(1), k2(1)
loop1:           DO  j = j1(1), j2(1)
                   DO  i = i1(1), i2(1)
                      IF ( array_1(i,j,k) /= array_2(i,j,k) )  THEN
                         PRINT*,'+++ data mismatch on element (',i,',',j,',',k,')'
                         PRINT*,'    array_1: ', array_1(i,j,k)
                         PRINT*,'    array_2: ', array_2(i,j,k)
                         n_err = n_err + 1
                         IF ( n_err > 5 )  EXIT loop
                      ENDIF
                   ENDDO
                ENDDO loop1
             ENDDO loop

             DEALLOCATE( array_1, array_2 )

          ELSEIF ( log_message(1)(1:2) == '2d' )  THEN
             PRINT*,'    *** reading 2d array'
             READ ( id+100, END=901 )  simtime(1), i1(1), i2(1), j1(1), j2(1)
             PRINT*,'        time=', simtime(1)
             PRINT*,'        array size=(',i1(1),':',i2(1), &
                                       ',',j1(1),':',j2(1),')'
             READ ( id+200, END=902 )  simtime(2), i1(2), i2(2), j1(2), j2(2)
             IF ( simtime(1) /= simtime(2) .OR. i1(1) /= i1(2) .OR. &
                  i2(1) /= i2(2) .OR. j1(1) /= j1(2) .OR. j2(1) /= j2(2) )  THEN
                PRINT*,'    +++ time/indices on file 2 does not match:'
                PRINT*,'        time=', simtime(2)
                PRINT*,'        array size=(',i1(2),':', &
                                i2(2), ',',j1(2),':',j2(2),')'
             ENDIF

             ALLOCATE( array_2d_1(i1(1):i2(1),j1(1):j2(1)), &
                       array_2d_2(i1(2):i2(2),j1(2):j2(2)) )

             READ ( id+100, END=903 )  array_2d_1
             READ ( id+200, END=904 )  array_2d_2

             IF ( i1(1) /= i1(2) )  i1(1) = i1(2)
             IF ( i2(1) /= i2(2) )  i2(1) = i2(2)
             IF ( j1(1) /= j1(2) )  j1(1) = j1(2)
             IF ( j2(1) /= j2(2) )  j2(1) = j2(2)

             n_err = 0
loop2:       DO  j = j1(1), j2(1)
                DO  i = i1(1), i2(1)
                   IF ( array_2d_1(i,j) /= array_2d_2(i,j) )  THEN
                      PRINT*,'+++ data mismatch on element (',i,',',j,')'
                      PRINT*,'    array_1: ', array_2d_1(i,j)
                      PRINT*,'    array_2: ', array_2d_2(i,j)
                      n_err = n_err + 1
                      IF ( n_err > 5 )  EXIT loop2
                   ENDIF
                ENDDO
             ENDDO loop2

             DEALLOCATE( array_2d_1, array_2d_2 )

          ELSE
             PRINT*,'    *** reading 2d int array'
             READ ( id+100, END=901 )  simtime(1), i1(1), i2(1), j1(1), j2(1)
             PRINT*,'        time=', simtime(1)
             PRINT*,'        array size=(',i1(1),':',i2(1), &
                                       ',',j1(1),':',j2(1),')'
             READ ( id+200, END=902 )  simtime(2), i1(2), i2(2), j1(2), j2(2)
             IF ( simtime(1) /= simtime(2) .OR. i1(1) /= i1(2) .OR. &
                  i2(1) /= i2(2) .OR. j1(1) /= j1(2) .OR. j2(1) /= j2(2) )  THEN
                PRINT*,'    +++ time/indices on file 2 does not match:'
                PRINT*,'        time=', simtime(2)
                PRINT*,'        array size=(',i1(2),':', &
                                i2(2), ',',j1(2),':',j2(2),')'
             ENDIF

             ALLOCATE( array_2d_i_1(i1(1):i2(1),j1(1):j2(1)), &
                       array_2d_i_2(i1(2):i2(2),j1(2):j2(2)) )

             READ ( id+100, END=903 )  array_2d_i_1
             READ ( id+200, END=904 )  array_2d_i_2

             IF ( i1(1) /= i1(2) )  i1(1) = i1(2)
             IF ( i2(1) /= i2(2) )  i2(1) = i2(2)
             IF ( j1(1) /= j1(2) )  j1(1) = j1(2)
             IF ( j2(1) /= j2(2) )  j2(1) = j2(2)

             n_err = 0
loop3:       DO  j = j1(1), j2(1)
                DO  i = i1(1), i2(1)
                   IF ( array_2d_i_1(i,j) /= array_2d_i_2(i,j) )  THEN
                      PRINT*,'+++ data mismatch on element (',i,',',j,')'
                      PRINT*,'    array_1: ', array_2d_i_1(i,j)
                      PRINT*,'    array_2: ', array_2d_i_2(i,j)
                      n_err = n_err + 1
                      IF ( n_err > 5 )  EXIT loop3
                   ENDIF
                ENDDO
             ENDDO loop3

             DEALLOCATE( array_2d_i_1, array_2d_i_2 )

          ENDIF

!          IF ( count > 8 )  STOP
       ENDDO

100    PRINT*, '*** end of data on file "', TRIM( filename(1) ), '"'
       PRINT*, '*** files seem to be identical'
       PRINT*, ' '
    ENDDO

    STOP

900 PRINT*,'+++ unexpected end on file "', TRIM( filename(2) ), '"'
    STOP
901 PRINT*,'+++ unexpected end on file "', TRIM( filename(1) ), '"'
    PRINT*,'    while reading indices'
    STOP
902 PRINT*,'+++ unexpected end on file "', TRIM( filename(2) ), '"'
    PRINT*,'    while reading indices'
    STOP
903 PRINT*,'+++ unexpected end on file "', TRIM( filename(1) ), '"'
    PRINT*,'    while reading array data'
    STOP
904 PRINT*,'+++ unexpected end on file "', TRIM( filename(2) ), '"'
    PRINT*,'    while reading array data'
    STOP

 END PROGRAM compare_palm_logs