function nd = days2(mo,da,leap)

nd = 0;

m = 1;

while m < mo
        
    if m == 1 | m == 3 | m == 5 | ...
       m == 7 | m == 8 | m == 10
       nday = 31;
    elseif m == 4 | m == 6 | m == 9 | m == 11
       nday = 30;
    else % m == 2
       nday = 28+leap;
    end
    
    nd = nd+nday;
    
    m = m+1;
end

nd = nd+da;