anchor.diff:             iBCBtpsiz=neltp*ndBCB
              call phio_readdatablock(fhandle, fname2 // char(0),
        &      c_loc(BCBtp),iBCBtpsiz,dataDbl,iotype)
  -
  -!MR CHANGE
  -!           write (temp1,"('nbcValuesDatablock_',i1,'_',i1)")
  -!     &                                              iblk,myrank
  -!           temp1=trim(temp1)
  -!           open(unit=12,file=temp1,status='unknown')
  -!           do i=1,neltp
  -!              do j=1,ndBCB
  -!                write(12,*) BCBtp(i,j)
  -!              enddo
  -!           enddo
  -!           write(12,*) iBCBtpsiz
  -!           close(12)
  -!MR CHANGE END
  -
   c
   c This is a temporary fix until NSpre properly zeros this array where it
   c is not set.  DEC has indigestion with these arrays though the
   c result is never used (never effects solution).
   c
  -
  -
  -!MR CHANGE
              if(writeLock==0) then
  -!MR CHANGE
  -!              print *,"In ASSIGN ASSIGN",myrank
                 where(.not.btest(iBCBtp(:,1),0)) BCBtp(:,1)=zero
                 where(.not.btest(iBCBtp(:,1),1)) BCBtp(:,2)=zero
                 where(.not.btest(iBCBtp(:,1),3)) BCBtp(:,6)=zero
  
positive.code:               writeLock=1;
             endif
  
  !           print *, "neltp is ", neltp
  
             call phio_readdatablock(fhandle, fname2 // char(0),
       &      c_loc(ientp),iientpsiz,dataInt,iotype)
  
  !MR CHANGE
  !           write (temp1,"('connectivityBoundaryDatablock_',i1,'_',i1)")
  !     &                                              iblk,myrank
  !           temp1=trim(temp1)
  
  !           open(unit=14,file=temp1,status='unknown')
  !           do i=1,neltp
  !              do j=1,nshl
  !                write(14,*) ientp(i,j)
  !              enddo
  !           enddo
  !           write(14,*) iientpsiz
  !           close(14)
  !MR CHANGE END
  
  
  c     
  c.... Read the boundary flux codes
  c
       
  !           print *,"connectivity [] is ", trim(fname2),ientp(0,0)
  
  
  ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
             call MPI_BARRIER(MPI_COMM_WORLD, ierr)
             
  ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
              if(nfiles.gt.0)then
                 write (fname2,"('nbc codes',i1)") iblk
              else
  
negative.code:               do n=1,neltp,ibksz 
                   nelblb=nelblb+1
                   npro= min(IBKSZ, neltp - n + 1)
                   lcblkb(1,nelblb)  = iel
                   lcblkb(3,nelblb)  = lcsyst
                   lcblkb(4,nelblb)  = ipordl
                   lcblkb(5,nelblb)  = nenl
  