将UTCG转换为MJD（IAU）

%--------------------------------------------------------------------------
% Purpose:
%
%   Modified Julian Date from calendar date and time
%
% Input(s):
%
%   Year      
%   Month
%   Day
%   Hour      
%   Min
%   Sec
% 
% output(s):  
%    
%    djm0  
%    djm     Modified Julian Date
% 
% Reference: O. Montenbruck, E. Gill (2000), "Satellite Orbits"
% 
%--------------------------------------------------------------------------

function [djm0 djm] = iauCal2jd(Year, Month, Day, Hour, Min, Sec)

  if (nargin == 3)
      Hour = 0;
      Min = 0;
      Sec = 0;
  end

  djm0 = 2400000.5;
   
  if ( Month<=2 ) 
      Month = Month+12;
      Year = Year-1;
  end
  
  if ( (10000*Year+100*Month+Day) <= 15