anchor.diff:    character(100) :: dir3_fname3
   
     DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE :: max_grids
  -  DOUBLE PRECISION :: global_max
  -  INTEGER :: grid_max, field_max
  +  DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: global_max
  +  INTEGER, DIMENSION(:), ALLOCATABLE :: grid_max
   
     !! CHARACTER :: c
     INTEGER :: nstp, ipos, idigit, l1, l2
  
positive.code:   CLOSE(10)
    CLOSE(20)
  
    OPEN(30,file=dir3_fname3)
    WRITE(30,130) global_max, 'Global maximum'
    WRITE(30,130) global_max, 'Global maximum'
    WRITE(30,140) grid_max, 'Global maximum occurs on this grid'
    WRITE(30,140) field_max, 'Global maximum occurs in this field'
    WRITE(30,*) ' '
    DO n = 0,ngrids-1
       WRITE(30,'(I5,10E16.4)') n, (max_grids(ngrid,m),m=1,mfields)
    ENDDO
    CLOSE(30)
  130 FORMAT(E30.16,'     ',A)
  140 FORMAT(I30,   '     ',A)
  
  END PROGRAM compare_files
  
  
negative.code:   ENDDO
  
    dir1_fname1 = TRIM(dir1)//'/'//trim(fname1)
     dir2_fname1 = TRIM(dir2)//'/'//trim(fname1)
  
    dir1_fname2 = TRIM(dir1)//'/'//trim(fname2)
    dir2_fname2 = TRIM(dir2)//'/'//trim(fname2)
  