anchor.diff:             ijunk  =intfromfile(6)
              lcsyst =intfromfile(7)
              allocate (ientp(neltp,nshl))
  -c           read(igeomBAK) ientp
              iientpsiz=neltp*nshl
   
              if (neltp==0) then
  
positive.code:         if (myrank == 0) then
            write(*,*) 'Number of interior topologies: ',itpblktot
          endif
  !        write (*,*) 'Rank: ',myrank,' interior itpblktot final:',
  !     &               itpblktot
  
  !MR CHANGE END
  
          nelblk=0
          mattyp = 0
          ndofl = ndof
          nsymdl = nsymdf
  
  !        call initphmpiio( nfields, nppf, nfiles, igeom )
  !        call openfile( fnamer, 'read', igeom )
  
  !         do iblk = 1, itpblk
          do iblk = 1, itpblktot
             writeLock=0;
  !MR CHANGE END
  c
  c           read(igeomBAK) neltp,nenl,ipordl,nshl, ijunk, ijunk, lcsyst
  c           call creadlist(igeomBAK,iseven,
  c     &          neltp,nenl,ipordl,nshl, ijunk, ijunk, lcsyst)
  
  ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  
              if(nfiles.gt.0) then
                write (fname2,"('connectivity interior',i1)") iblk
              else
                write (fname2,"('connectivity interior linear tetrahedron')") 
              endif
  
  ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
  
  c           fname1='connectivity interior?'
  
             ! Synchronization for performance monitoring, as some parts do not include some topologies
             call MPI_Barrier(MPI_COMM_WORLD,ierr) 
  
negative.code:         call phio_readheader(fhandle,
       &   c_char_'total number of interior tpblocks' // char(0),
       &   c_loc(itpblktot), ione, dataInt, iotype) 
  
          if (itpblktot == -1) then 
            ! The field 'total number of different interior tpblocks' was not found in the geombc file.
  