anchor.diff:  c
   
         ione=1
  -      write (temp1,"('(''boundary condition array@'',i',i1,',A1)')")
  -     &     itmp
  -      write (fname2,temp1) (myrank+1),'?'
   
   c      ione=1
   c      fname1='boundary condition array?'
   c      call readheader(igeomBAK,fname1,intfromfile,
   c     &     ione, 'double', iotype)
  -      call readheader(igeom,fname2 // char(0),intfromfile,
  -     &     ione, 'double' // char(0), iotype)
  +      call readheader(igeom,'boundary condition array' // char(0),
  +     & intfromfile,ione, 'double' // char(0), iotype)
   c here intfromfile(1) contains (ndof+7)*numpbc
   !MR CHANGE
   !       if ( numpbc > 0 ) then
  
positive.code: cc MR CHANGE
  
        ione=1
        write (temp1,"('(''bc mapping array@'',i',i1,',A1)')") itmp
        write (fname2,temp1) (myrank+1),'?'
        call readheader(igeom,fname2 // char(0),nshg,ione,
       & 'integer' // char(0),iotype)
  
  c      fname1='bc mapping array?'
  
negative.code: 
  c      itmp=1
  c      if (myrank .gt. 0) itmp = int(log10(float(myrank)))+1
        intfromfile=0
        call readheader(descriptor,'time derivative of solution' // char(0),
       & intfromfile,ithree,'integer' // char(0),iotype)
        allocate( acold(nshg,ndof) )
        if(intfromfile(1).ne.0) then
           nshg2=intfromfile(1)
  