anchor.diff:                writeLock=1;
              endif
   
  -           call readdatablock(igeom,fname2 // char(0),ientp,iientpsiz,
  +           call phio_readdatablock(igeom,fname2 // char(0),ientp,iientpsiz,
        &                     "integer" // char(0), iotype)
   
   !            call closefile( igeom, "read" // char(0) )
  
positive.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
          write(temp1,
       &   "('(''total number of interior tpblocks@'',i',i1,',A1)')") itmp
  
          write (fname2,temp1) (myrank+1),'?'
          call phio_readheader(igeom,fname2 // char(0) ,itpblktot,ione,
       &  'integer' // char(0),iotype) 
  
  !        write (*,*) 'Rank: ',myrank,' interior itpblktot intermediate:',
  !     &               itpblktot
  
negative.code: 
              intfromfile(:)=-1
              iblk = iblk+1
              write (fname2,"('connectivity interior',i1)") iblk
  
              !write(*,*) 'rank, fname2',myrank, trim(adjustl(fname2))
              call phio_readheader(igeom,fname2 // char(0),intfromfile,
  