netgraph_external.m

％The Metabolic Networks Toolbox contains functions to create, ％modify, display, and simulate bioche

function [network,p] = netgraph_external(network,draw_only,external)

% [network,p] = netgraph_external(network,draw_only,external)
%
% arguments draw_only,external are optional
% left/centre button: add/ remove external
% right button quit

fprintf('Click on metabolite to switch between internal (yellow) and external (orange) \nRight button: Quit\n');

if exist('external','var'), network.external = external; end

m = network.graphics_par.m;
x = network.graphics_par.x;

p.metvalues=0.2*ones(length(network.metabolites),1);
p.arrowvalues=zeros(size(network.graphics_par.arrowvalues));
p.arrowstyle='fluxes';
p.actvalues=zeros(length(network.actions),1);
p.metvalues(find(network.external))=0.5;
p.metvalues_std=zeros(length(network.metabolites),1);
p.metvaluesmax=1;
p.actvaluesmax=1;
p.metstyle='box';
p.actstyle='none';
p.actprintnames=0;
netgraph_draw(network,p);


if ~exist('draw_only','var'), draw_only = 0; end
if draw_only,  return; end

n_met=length(network.metabolites);
cont = 1;

while cont==1,
     [x_old,y_old,button] = ginput(1);
     dist= sum( (repmat([x_old;y_old],1,size(x,2))-x).^2);
     [dum,i]=min(dist);
     index=i(1);
 switch button
 case 3,
     cont = 0; 
 otherwise
   if index <= n_met, 
     if network.external(index), network.external(index)=0;
    else network.external(index)=1; 
      end
   end
 end
     p.metvalues=0.2*ones(length(network.metabolites),1);
     p.metvalues(find(network.external))=0.5;
     p.metvaluesmax=1;
     p.actvaluesmax=1;
    netgraph_draw(network,p);

end

fprintf('Quitted external metabolites\n');