load curvemat

figure;
c = [0 0 1];
wh = [1 1 1];

cs = 0.5*(wh+c);
plot(curve.brupt,curve.brupN,'Color',c);

hold on;
plot(curve.brupt,curve.brupN+curve.brupNe,'Color',cs);
plot(curve.brupt,curve.brupN-curve.brupNe,'Color',cs);


c = [0 1 1];
wh = [1 1 1];

cs = 0.5*(wh+c);
plot(curve.brupt,curve.brupC,'Color',c);

plot(curve.brupt,curve.brupC+curve.brupCe,'Color',cs);
plot(curve.brupt,curve.brupC-curve.brupCe,'Color',cs);



c = [1 0 0];
wh = [1 1 1];

cs = 0.5*(wh+c);
plot(curve.mgupt,curve.mgupN,'Color',c);

plot(curve.mgupt,curve.mgupN+curve.mgupNe,'Color',cs);
plot(curve.mgupt,curve.mgupN-curve.mgupNe,'Color',cs);


c = [1 0 1];
wh = [1 1 1];

cs = 0.5*(wh+c);
plot(curve.mgupt,curve.mgupC,'Color',c);

plot(curve.mgupt,curve.mgupC+curve.mgupCe,'Color',cs);
plot(curve.mgupt,curve.mgupC-curve.mgupCe,'Color',cs);

xlabel('Time (hours)');
ylabel('Cell width');

print('-dpdf','upshift-growth.pdf');

figure;

c = [0 0 1];
wh = [1 1 1];

cs = 0.5*(wh+c);
plot(curve.brdownt,curve.brdownN,'Color',c);

hold on;
plot(curve.brdownt,curve.brdownN+curve.brdownNe,'Color',cs);
plot(curve.brdownt,curve.brdownN-curve.brdownNe,'Color',cs);


c = [0 1 1];
wh = [1 1 1];

cs = 0.5*(wh+c);
plot(curve.brdownt,curve.brdownC,'Color',c);

plot(curve.brdownt,curve.brdownC+curve.brdownCe,'Color',cs);
plot(curve.brdownt,curve.brdownC-curve.brdownCe,'Color',cs);



c = [1 0 0];
wh = [1 1 1];

cs = 0.5*(wh+c);
plot(curve.mgdownt,curve.mgdownN,'Color',c);

plot(curve.mgdownt,curve.mgdownN+curve.mgdownNe,'Color',cs);
plot(curve.mgdownt,curve.mgdownN-curve.mgdownNe,'Color',cs);


c = [1 0 1];
wh = [1 1 1];

cs = 0.5*(wh+c);
plot(curve.mgdownt,curve.mgdownC,'Color',c);

plot(curve.mgdownt,curve.mgdownC+curve.mgdownCe,'Color',cs);
plot(curve.mgdownt,curve.mgdownC-curve.mgdownCe,'Color',cs);

xlabel('Time (hours)');
ylabel('Cell width');

print('-dpdf','downshift-growth.pdf');
