anchor.diff:  c.... read the temporary iBC array
   c
         ione=1
  -      write (temp1,"('(''bc codes array@'',i',i1,',A1)')") itmp
  -      write (fname2,temp1) (myrank+1),'?'
  -      call readheader(igeom,fname2 // char(0) ,numpbc,ione,
  +      call readheader(igeom,'bc codes array' // char(0) ,numpbc,ione,
        & 'integer' // char(0),iotype)
   
   c      ione = 1
  
positive.code:       if(numpe > 1) then
  
  cc MR CHANGE
           write (temp1,"('(''size of ilwork array@'',i',i1,',A1)')") itmp
           write (fname2,temp1) (myrank+1),'?'
           call readheader(igeom,fname2 // char(0),nlwork,ione,
       &   'integer' // char(0) ,iotype)
  
           write (temp1,"('(''ilwork@'',i',i1,',A1)')") itmp
           write (fname2,temp1) (myrank+1),'?'
           call readheader(igeom,fname2 //char(0) ,nlwork,ione,
       &   'integer' // char(0) , iotype)
  
           allocate( point2ilwork(nlwork) )
           allocate( ilworkread(nlwork) )
           call readdatablock(igeom,fname2 // char(0),ilworkread,
       &                      nlwork,'integer' // char(0) , iotype)
  
  c      call closefile( igeom, "read" )
  
negative.code:          allocate( acread(nshg,ndof2) )
           acread=zero
           iacsiz=nshg*ndof2
           call readdatablock(descriptor,
       &    'time derivative of solution' // char(0),acread,
       &    iacsiz, 'double' // char(0),iotype)
           acold(:,1:ndof)=acread(:,1:ndof)
           deallocate(acread)
  