Finding decay and brancing ratios

Hi!

So I am trying to find the decays and the branching ratios for different benchmark points of the Inert Doublet Model. I already have the UFO.
So, my first idea was to find the widths using e.g. for the h3 particle: (I should add that, due to the Z2 symmetry of the IDM, the h3 has to decay to an h2, which is why I have done > all all all as it can't decay to an onshell Z).

import model InertDoublet_UFO
generate h3 > all all all
output widthh3
launch widthh3
0
set param_card frblock 1 0.007297353
set param_card frblock 2 1.166379e-05
set param_card frblock 4 80.377
set param_card frblock 12 125.25
set param_card frblock 13 65
set param_card frblock 14 71.525
set param_card frblock 15 112.85
set param_card mass 24 80.377
set param_card mass 25 125.25
set param_card mass 35 65
set param_card mass 36 71.525
set param_card mass 37 112.85
set param_card frblock 11 0.3895575
set param_card frblock 10 0.0002
set param_card DECAY 6 auto
set param_card DECAY 23 auto
set param_card DECAY 24 auto
set param_card DECAY 25 auto

Whilst this does get a cross-section output, it does not give me the specific branching ratios.
So my question is:
How do I find the decays and the branching ratios for each benchmark point like above? I've tried using compute_widths but i'm unsure as to how to do this with different benchmark points. The output that I am looking for is similar to the output:

# PDG Width
DECAY 23 2.353339e+00
# BR NDA ID1 ID2 ...
   1.523161e-01 2 1 -1 # 0.358451255228
   1.523153e-01 2 3 -3 # 0.358449523568
   1.510009e-01 2 5 -5 # 0.355356187407
   1.187190e-01 2 2 -2 # 0.279385890968
   1.186130e-01 2 4 -4 # 0.279136567859
   6.798413e-02 2 16 -16 # 0.159989682115
   6.798413e-02 2 14 -14 # 0.159989682115
   6.798413e-02 2 12 -12 # 0.159989682115
   3.438703e-02 2 11 -11 # 0.0809243178553
   3.438675e-02 2 13 -13 # 0.0809236729573
   3.430960e-02 2 15 -15 # 0.0807421157318

Thanks!