anchor.diff:  
         !$OMP CRITICAL
         call map_update(map_c, key, value, idx,1.d-15) 
  -!      call map_append(map_c, key, value, idx)
         !$OMP END CRITICAL
   
     !WRITE OUTPUT
  
positive.code:   integer*8 :: tempspace
  
    tempspace = (new_size * 16_8) / (1024_8 * 1024_8)
    npass = min(l_end-l_start,1 + tempspace / 2048)   ! 2 GiB of scratch space
    l_block = (l_end-l_start+1)/npass
  
    ipass = 0
  
negative.code:     deallocate(key,value,V,T)
      !$OMP END PARALLEL
      call map_merge(map_c)
  
      deallocate(l_pointer)
    enddo
  