anchor.diff:    !!      write(6,*) 'WARNING : (claw_out2.f ) Setting q to 0'
     DO j = 1,my
        DO i = 1,mx
  -        DO mq = 1,mfields
  +        DO mq = 1,meqn
              IF (ABS(q(i,j,mq)) .LT. 1d-99) THEN
                 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)
  +        WRITE(fid_com,120) (q(i,j,mq),mq=1,meqn)
        ENDDO
        WRITE(fid_com,*) ' '
     ENDDO
   !! Print the first field assuming it is the MPI rank (i.e. an integer)
  -120 FORMAT (I5,10E26.16)
  +120 FORMAT (10E26.16)
   
     CLOSE(fid_com)
   
  
positive.code:      READ(20,*) ngrid2
       READ(20,*) level2
       READ(20,*) blockno2
       READ(20,*) mx2
       READ(20,*) my2
       READ(20,*) xlow2
  
negative.code:   ENDIF
    t = t1
    ngrids = ngrids1
    meqn = meqn1
  
    WRITE(6,'(A,F24.8)') 'Time of simulation : t = ',t
  
    dir3_fname3 = TRIM(dir3)//'/'//TRIM(fname3)
    WRITE(6,*) 'A summary of the results will be stored in : '
    WRITE(6,*) dir3_fname3
  
  
    ALLOCATE(max_grids(0:ngrids-1,meqn))
  
    CALL write_tfile(iframe,t,meqn,ngrids,dir3)
    CALL new_qfile(iframe,dir3)
  
    ALLOCATE(grid_max(meqn),global_max(meqn))
    grid_max = 0
    global_max = 0
  
  
  
    OPEN(10,FILE=dir1_fname1,ERR=100)
    OPEN(20,FILE=dir2_fname1,ERR=200)
    DO n = 0, ngrids-1
       READ(10,*) ngrid1
       READ(10,*) level1
       READ(10,*) blockno1
       READ(10,*) mpirank1
       READ(10,*) mx1
       READ(10,*) my1
       READ(10,*) xlow1
  