getTransportRxns

 getTransportRxns
   Retrieves the transport reactions from a model

   model           a model structure

   transportRxns   logical array with true if the corresponding 
                   reaction is a transport reaction

   Transport reactions are defined as reactions involving (at least) one
   metabolite name in more than one compartment.

   Usage: transportRxns=getTransportRxns(model)

0001 function transportRxns=getTransportRxns(model)
0002 % getTransportRxns
0003 %   Retrieves the transport reactions from a model
0004 %
0005 %   model           a model structure
0006 %
0007 %   transportRxns   logical array with true if the corresponding
0008 %                   reaction is a transport reaction
0009 %
0010 %   Transport reactions are defined as reactions involving (at least) one
0011 %   metabolite name in more than one compartment.
0012 %
0013 %   Usage: transportRxns=getTransportRxns(model)
0014 
0015 transportRxns=false(numel(model.rxns),1);
0016 
0017 for i=1:numel(model.rxns)
0018     %Get the involved metabolites in each reaction
0019     mets=model.metNames(model.S(:,i)~=0);
0020     transportRxns(i)=numel(mets)~=numel(unique(mets));
0021 end
0022 end