%Set initial conditions for landscape

function z = setelev(x,y) 


    for i=1:x
        for j=1:y
            z(i,j) = 0;
        end
    end

    
    for i=1:x
        for j=1:y
            f = 1.0 - abs(j-y/2)*2/y;
             if (f>0) 
                   z(i,j) =  100 * f;
             end
        end
   end

    
    % If I want randomness, but it seems to bias flow
    %A = 1.2;
    %for i=1:x
    %    for j=1:y
    %        z(i,j) = z(i,j) + A*rand();
    %    end
    %end