anchor.diff:  c     
   c.... read the boundary condition data
   c     
  -           call MPI_BARRIER(MPI_COMM_WORLD, ierr)
  +C           call MPI_BARRIER(MPI_COMM_WORLD, ierr)
               if(input_mode.ge.1)then
                  write (fname2,"('nbc values',i1)") iblk
               else
  
positive.code:            else
               write( fname2,"('material type boundary')")
             endif
             call MPI_BARRIER(MPI_COMM_WORLD, ierr)
             call phio_readheader(fhandle, fname2 // char(0),
       &      c_loc(intfromfile), ione, dataInt, iotype)
             call phio_readdatablock(fhandle, fname2 // char(0),
  
negative.code: c     
  c.... Read the boundary flux codes
  c
  C           call MPI_BARRIER(MPI_COMM_WORLD, ierr)
              if(input_mode.ge.1)then
                 write (fname2,"('nbc codes',i1)") iblk
              else
  