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