function jd = julday(y,m,d,h)
% JULDAY Conversion of date as given by
% y ... year (four digits)
% m ... month
% d ... day
% h ... hour and fraction hereof
% The conversion is only valid in the time span
% from March, 1, 1900 to February, 28, 2100
% For further information see
% Meeus, Jean (1991) Astronomical Algorithms,
% Willmann-Bell, Richmond, Virginia, p. 59--62
% Written by Kai Borre
% February 14,2001
if m <= 2, y = y-1; m = m+12; end
jd = floor(365.25*(y+4716))+floor(30.6001*(m+1))+d+h/24-1537.5;
% mjd = jd-2400000.5;
%%%%%%% end julday.m %%%%%%%%%%%%