MtkUtcJdToTaiJd.c

Go to the documentation of this file.

/*===========================================================================
=                                                                           =
=                              MtkUtcJdToTaiJd                              =
=                                                                           =
=============================================================================

                         Jet Propulsion Laboratory
                                   MISR
                               MISR Toolkit

            Copyright 2006, California Institute of Technology.
                           ALL RIGHTS RESERVED.
                 U.S. Government Sponsorship acknowledged.

============================================================================*/

#include "MisrUtil.h"
#include "MisrError.h"

MTKt_status MtkUtcJdToTaiJd(
  double jdUTC[2], 
  double jdTAI[2]  )
{
  MTKt_status status_code; /* Return code of this function */
  double leap_seconds[][2] = LEAP_SECONDS;
  double leapSecs = 0.0;
  int i;
  
  /* Find leap seconds */
  for (i = sizeof(leap_seconds) / sizeof(*leap_seconds) - 1; i >= 0; --i)
    if (leap_seconds[i][0] <= jdUTC[0])
    {
      leapSecs = leap_seconds[i][1];
      break;
    } 

  /* Make sure time isn't before leap second range. */
  if (i == -1)
    MTK_ERR_CODE_JUMP(MTK_BAD_ARGUMENT) 

  jdTAI[0] = jdUTC[0];
  jdTAI[1] = jdUTC[1] + ((leapSecs) / SECONDSperDAY);
    
  if (jdTAI[1] >= 1.0)
  {
    jdTAI[0] +=  1.0;
    jdTAI[1] -=  1.0;
  }   
    
  return MTK_SUCCESS;

ERROR_HANDLE:
  return status_code;
}