anchor.diff:  subroutine dirichlet(n,counts,nc,alpha,estimate)
     ! posterior mean entropy (averaged over Dirichlet distribution) given alpha
     use iso_fortran_env
  -  use dirichlet_mod, only: initialize_dirichlet, compute_multiplicities, dirichlet_finalize
  +  use dirichlet_mod, only: initialize_dirichlet, compute_multiplicities, finalize
     use dirichlet_mod, only: h_dir
     implicit none
   
  
positive.code: 
  subroutine phony_4(n,counts,nc,estimate,err_estimate)
    use iso_fortran_env
    use dirichlet_mod, only: initialize_dirichlet, compute_multiplicities, dirichlet_finalize
    use nsb_mod, only: hnsb
    use nsb_mod, only: compute_integration_range
    implicit none
  
negative.code: 
    ! call hnsb(estimate,err_estimate, err)
  
    call finalize()
  
    call cpu_time(finish)
  
  