%% Initialize

load = (vload/N_load)*ones(N_load,1);

unload_1  = - sc.Vmax1  / N_1_un .* (cp(Loc_1_un)./( cp(Loc_1_un) + sc.Km1));
unload_2  = - sc.Vmax2  / N_2_un .* (cp(Loc_2_un)./( cp(Loc_2_un) + sc.Km2));


%% Loading zone
i=2:Loc_load(end);
JS_load = [load];


%Transport equations sink1 
i=Loc_1;
JS_load   =   [JS_load ; zeros(N_1,1)];

i=Loc_1_un;
JS_load   =   [JS_load ; unload_1 .* ones(N_1_un,1)];

%Transport equations sink2 
i=Loc_2;
JS_load   =   [JS_load ; zeros(N_2,1)];

i=Loc_2_un;
JS_load   =   [JS_load ; unload_2 .* ones(N_2_un,1)];
