function out = w03reflF1(x,bc,n)
% 
%  initial data functions (djm: 21 sept 2010)
% 

switch n
    case{1}
        slope = (2+sign(x-6));
        out1 = (1/2 - abs(x-6)./slope);
        slope = (2+sign(x+2));
        out2 = (1/2 - abs(x+2)./slope);
        out = out1.*(out1>0) + bc*out2.*(out2>0);
    case{2}
        x = -x;

        slope = (2+sign(x-6));
        out1 = (1/2 - abs(x-6)./slope);
        slope = (2+sign(x+2));
        out2 = (1/2 - abs(x+2)./slope);
        out = bc*out1.*(out1>0) + out2.*(out2>0);
end

