alban's code now in cic_velpower_dmnu.f90