Problem with fillGaps function

[JOELROJAS2022]

I was performing a gapfilling procedure on a context-specific model built upon Human1 with ftINIT. The resulting model was closed with the "closeModel" function. I use fillGaps as followed:

[newConnected, cannotConnect, addedRxns, newModel, exitFlag] = fillGaps(model, ihuman, true, false, false, params);

And the error I get is:

WARNING: This algorithm is meant to function on a model with exchange reactions for uptake and excretion of metabolites. The current model still has the "unconstrained" field

Matrix index is out of range for deletion.

Error in removeReactions (line 87)
reducedModel.rxnScores(indexesToDelete,:)=[];
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in simplifyModel (line 98)
reducedModel=removeReactions(reducedModel,rxnsToDelete);
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Error in fillGaps (line 120)
models{i}=simplifyModel(models{i},false,false,true);
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^

Do you have any suggestions?

[edkerk]

Could you share the model so that we could replicate this behaviour?