function outvar = Source_Term(x,y,t)
if t<1
    if sqrt((x-.5)*(x-.5)+(y-.75)*(y-.75)) <.1
        outvar = 5;
    else
        if sqrt((x+.5)*(x+.5)+(y+.75)*(y+.75)) <.2
            outvar=-2;
        else
            outvar = 0;
        end
    end
else
    if sqrt((x-.5)*(x-.5)+(y-.75)*(y-.75)) <.1
        outvar = -3;
    else
        if sqrt((x+.5)*(x+.5)+(y+.75)*(y+.75)) <.2
            outvar=1;
        else
            outvar = 0;
        end
    end
end
