anchor.diff:                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
  
positive.code:            iiBCBtpsiz=neltp*ndiBCB
             call phio_readdatablock(fhandle, fname2 // char(0),
       &      c_loc(iBCBtp),iiBCBtpsiz,dataInt,iotype)
  
  !MR CHANGE
  !           print *, "ndiBCB is ",ndiBCB
  !           write (temp1,"('nbcCodesDatablock_',i1,'_',i1)")
  !     &                                              iblk,myrank
  !           temp1=trim(temp1)
  !
  !           open(unit=13,file=temp1,status='unknown')
  !           do i=1,neltp
  !              do j=1,ndiBCB
  !                write(13,*) iBCBtp(i,j)
  !              enddo
  !           enddo
  !           write(13,*) iiBCBtpsiz
  !           close(13)
  !!           iBCBtp(:,:) = 0 ! JUST FOR TEST
  !MR CHANGE END
  
  c     
  c.... read the boundary condition data
  c     
             
  ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
             call MPI_BARRIER(MPI_COMM_WORLD, ierr)
  ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
              if(nfiles.gt.0)then
                 write (fname2,"('nbc values',i1)") iblk
              else
  
negative.code:            iBCBtpsiz=neltp*ndBCB
             call phio_readdatablock(fhandle, fname2 // char(0),
       &      c_loc(BCBtp),iBCBtpsiz,dataDbl,iotype)
  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
             if(writeLock==0) then
                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
  