anchor.diff:  
               intfromfile(:)=-1
               iblk = iblk+1
  -            write (temp1,"('connectivity boundary',i1)") iblk
  -            temp1 = trim(temp1)
  -            write (temp3,"('(''@'',i',i1,',A1)')") itmp
  -            write (fname2, temp3) (myrank+1), '?'
  -            fname2 = trim(temp1)//trim(fname2)
  -            !write(*,*) 'rank, fname2',myrank, trim(adjustl(fname2))
  -            call readheader(igeom,fname2 // char(0),intfromfile,
  -     &      ieight,'integer' // char(0),iotype)
  +            call readheader(igeom,'connectivity boundary1' // char(0),
  +     &       intfromfile,ieight,'integer' // char(0),iotype)
               neltp = intfromfile(1) ! -1 if fname2 was not found, >=0 otherwise
             end do
             itpblktot = iblk-1   
  
positive.code: !           print *,"connectivity [] is ", trim(fname2),ientp(0,0)
  
  
             fname1='nbc codes?'
  
  ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  
             write (temp1,"('nbc codes',i1)") iblk
             temp1=trim(temp1)
             write (temp3,"('(''@'',i',i1,',A1)')") itmp
             write (fname2, temp3) (myrank+1), '?'
             fname2 = trim(temp1)//trim(fname2)
  !           write(*,*) 'rank, fname2',myrank, trim(adjustl(fname2))
             call MPI_BARRIER(MPI_COMM_WORLD, ierr)
             
  ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  
             call readheader(igeom,fname2 // char(0) ,intfromfile,
       &      ieight,'integer' // char(0),iotype)
             iiBCBtpsiz=neltp*ndiBCB
             call readdatablock(igeom,fname2 // char(0) ,iBCBtp,
       &      iiBCBtpsiz,'integer' // char(0),iotype)
  
  !MR CHANGE
  
negative.code:         ! Get the total number of different interior topologies in the whole domain. 
          ! Try to read from a field. If the field does not exist, scan the geombc file.
          itpblktot=-1
          call readheader(igeom,'total number of boundary tpblocks' // char(0),
       &   itpblktot,ione,'integer' // char(0),iotype)
  
  !        write (*,*) 'Rank: ',myrank,' boundary itpblktot intermediate:',
  !     &               itpblktot
  