anchor.diff:  
     tempspace = (new_size * 16_8) / (1024_8 * 1024_8)
     npass = min(l_end-l_start,1 + tempspace / 2048)   ! 2 GiB of scratch space
  -  npass = 1
  -  l_block = (l_end-l_start)/npass
  +  l_block = (l_end-l_start+1)/npass
   
     ipass = 0
     do l_start_block = l_start, l_end, l_block
  
positive.code:       do b=b_start,d
          do c=c_start,c_end
            do a=a_start,min(b,c)
              if (dabs(U(a,c,b)) < 1.d-15) then
                cycle
              endif
  
negative.code:     deallocate(key,value,V,T)
      !$OMP END PARALLEL
      call map_merge(map_c)
  !    call map_sort(map_c)
  
      deallocate(l_pointer)
    enddo
  