anchor.diff:  c
   
         ione=1
  -      write (temp1,"('(''periodic masters array@'',i',i1,',A1)')") itmp
  -      write (fname2,temp1) (myrank+1),'?'
  -
   c      fname1='periodic masters array?'
   c      call readheader(igeomBAK,fname1,nshg,
   c     &     ione, 'integer', iotype)
  -      call readheader(igeom,fname2 // char(0) ,nshg,
  -     &     ione, 'integer' // char(0), iotype)
  +      call readheader(igeom,'periodic masters array' // char(0) ,nshg,
  +     & ione,'integer' // char(0), iotype)
         allocate( point2iper(nshg) )
         allocate( iperread(nshg) )
   c      call readdatablock(igeomBAK,fname1,iperread,nshg,
   c     &                      'integer',iotype)
  -      call readdatablock(igeom,fname2 // char(0),iperread,nshg,
  -     &                      'integer' // char(0),iotype)
  +      call readdatablock(igeom,'periodic masters array' // char(0),
  +     & iperread,nshg,'integer' // char(0),iotype)
         point2iper=iperread
   
   
  
positive.code: cccccccccccccccccccccccccccccccccccccccccc
  
        itwo=2
        write (temp1,"('(''co-ordinates@'',i',i1,',A1)')") itmp
        write (fname2,temp1) (myrank+1),'?'
  
  c      print *, "fname2 is", fname2
  
  
negative.code: !             fname1='number of father-nodes?'
  !             call readheader(igeomBAK,fname1,nfath,ione,'integer', iotype)
  
              call readheader(igeom,'number of father-nodes' // char(0),
       &       nfath,ione,'integer' // char(0), iotype)
  
  c
  c     fname1='keyword nsons?'
  !             fname1='number of son-nodes for each father?'
  !             call readheader(igeomBAK,fname1,nfath,ione,'integer', iotype)
  
              call readheader(igeom,
       &       'number of son-nodes for each father' // char(0),
       &       nfath,ione,'integer' // char(0), iotype)
  
              allocate (point2nsons(nfath))
  
  !             call readdatablock(igeomBAK,fname1,point2nsons,nfath,
  !      &                      'integer',iotype)
              call readdatablock(igeom,
       &       'number of son-nodes for each father' // char(0),
       &       point2nsons,nfath,'integer' // char(0), iotype)
  
  c
  !             fname1='keyword ifath?'
  !             call readheader(igeomBAK,fname1,nshg,ione,'integer', iotype)
  
              call readheader(igeom,'keyword ifath' // char(0),nshg,ione,
       &      'integer' // char(0), iotype)
  
              allocate (point2ifath(nshg))
  
  !             call readdatablock(igeomBAK,fname1,point2ifath,nshg,
  !      &                      'integer',iotype)
              call readdatablock(igeom,'keyword ifath' // char(0),point2ifath,
       &                      nshg,'integer' // char(0) , iotype)
       
  c     
  