coi_um.c/h: function to estimate intensity