Adding bond types later

Hello
I am currently using 'LIGGGHTS-flexible fibers'.
My simulation involves several steps, and I create restart files at each step.
I initially defined n_bondtype=2 for my simulations. Later I realized that I have to include one more bond type. To do the same I did the following:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
# Specify what we are simulating
atom_style hybrid granular bond/gran n_bondtypes 3 bonds_per_atom 50 disableNormalContact 1 %%% IMPORTATNT
atom_modify map array
hard_particles yes # <- Needed because of the high Youngs modulus

## newton off is needed due to the tangential history contact model
newton off
communicate single vel yes
#units si

# Use hertz-mindlin contact model
pair_style gran model hertz tangential history
# Use the stiff granular bond model
bond_style gran

processors * * *

# Set neighbor bin sizes and update after each run step if needed
neighbor 0.001 bin
neigh_modify delay 0 one 10000

## Read in the restart file
read_restart restarts/restart4.liggghts
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Here I am getting the error: ERROR: Numeric index is out of bounds (../force.cpp:786)
Because of the extra Bond coefficient defined later

I tried to define atom_style after I read restart, in that case: ERROR: Atom_style command after simulation box is defined (../input.cpp:1226)

Is there any way to include the extra bond type in my simulations ??

Thanks