anchor.diff:    ENDDO
   
     dir1_fname1 = TRIM(dir1)//'/'//trim(fname1)
  -  dir2_fname1 = TRIM(dir2)//'/'//trim(fname1)
  +   dir2_fname1 = TRIM(dir2)//'/'//trim(fname1)
   
     dir1_fname2 = TRIM(dir1)//'/'//trim(fname2)
     dir2_fname2 = TRIM(dir2)//'/'//trim(fname2)
  
positive.code:   character(100) :: dir3_fname3
  
    DOUBLE PRECISION, DIMENSION(:,:), ALLOCATABLE :: max_grids
    DOUBLE PRECISION :: global_max
    INTEGER :: grid_max, field_max
  
    !! CHARACTER :: c
    INTEGER :: nstp, ipos, idigit, l1, l2
  
negative.code:               q(i,j,mq) = 0.d0
             ENDIF
          END DO
          WRITE(fid_com,120) INT(q(i,j,1)), (q(i,j,mq),mq=2,mfields)
       ENDDO
       WRITE(fid_com,*) ' '
    ENDDO
  