#include "dbhost.h"     /* PROTOTYPE_ALLOWED */
#include "dbext.h"      /* DMSH_POSITION_TYPE, BAD*, LONG, etc. */
/*****************************************************************************
*                                                                            *
*      COMMON OCEANOGRAPHIC DATA ACCESS SYSTEM (CODAS)                       *
*                                                                            *
*      WRITTEN BY:  RAMON CABRERA, ERIC FIRING, and JULIE RANADA             *
*                   JOINT INSTITUTE FOR MARINE AND ATMOSPHERIC RESEARCH      *
*                   1000 POPE ROAD  MSB 404                                  *
*                   HONOLULU, HI 96822                                       *
*                                                                            *
*      VERSION:     3.00                                                     *
*                                                                            *
*      DATE:        APRIL 1989                                               *
*                                                                            *
*****************************************************************************/
/*-----------------------------------------------------------------------------

       FILE:  posn.c

              POSITION CONVERSION & PRINTING FUNCTIONS

*/
void HUNPOS(DMSH_POSITION_TYPE *p, LONG *h);
LONG POSHUN(DMSH_POSITION_TYPE *p);
/*-----------------------------------------------------------------------------

   FUNCTION:  invalid_position

      It checks a DMSH_POSITION_TYPE structure.

   PARAMETERS:
   
      p = pointer to DMHS_POSITION_TYPE structure

   RETURNS:  0 if position is okay
             INVALID_POSITION (> 0) otherwise

   CALLED FROM:  C

*/
int invalid_position(DMSH_POSITION_TYPE *p, int is_lat)
{
   int test_mask = 0;   /* jr+: must initialize to 0 in case of 0 d, m, s, or h */

   if (p->degree == BADSHORT) return(0);  /* jr+: ok value for no position */
   if ( ( is_lat && (abs_val(p->degree) > 90 ) ) ||
        (!is_lat && (abs_val(p->degree) > 359) ) ||
        (abs_val(p->minute) > 59) ||
        (abs_val(p->second) > 59) ||
        (abs_val(p->hundredth) > 99 )
      )
      return(INVALID_POSITION);
   if (p->degree > 0) 
      test_mask = 1;
   else if (p->degree < 0) 
      test_mask = 2;
   if (p->minute > 0) 
      test_mask = test_mask | 1;
   else if (p->minute < 0) 
      test_mask = test_mask | 2;
   if (p->second > 0) 
      test_mask = test_mask | 1;
   else if (p->second < 0) 
      test_mask = test_mask | 2;
   if (p->hundredth > 0) 
      test_mask = test_mask | 1;
   else if (p->hundredth < 0) 
      test_mask = test_mask | 2;
   if (test_mask == 3)  
      return(INVALID_POSITION);
   return(0);
}

/*-----------------------------------------------------------------------------

   FUNCTION:  POSHUN

      Given a DMSH_POSITION_TYPE structure, it returns the equivalent position
      in hundredths of a second.

   PARAMETER:

      p = pointer to DMSH_POSITION_TYPE structure

   RETURNS:  position in hundredths of a second

   CALLED FROM:  C, FORTRAN

*/
LONG POSHUN(DMSH_POSITION_TYPE *p)
{
   LONG l_360000 = 360000L, l_6000 = 6000, l_100 = 100;

   if (p->degree == BADSHORT) return(BADLONG); /* jr+: bad value */
   return ((p->degree * l_360000)
      + (p->minute * l_6000)
      + (p->second * l_100)
      + (p->hundredth));
}

/*-----------------------------------------------------------------------------

   FUNCTION:  HUNPOS

      Given a position in hundredths of a second, it transforms it into the
      equivalent DMSH_POSITION_TYPE structure.

   PARAMETERS:
   
      p = pointer to DMSH_POSITION_TYPE structure

      h = pointer to position in hundredths of a second

   RETURNS:  VOID

   CALLED FROM:  C, FORTRAN

*/
void HUNPOS(DMSH_POSITION_TYPE *p, LONG *h)
{
   LONG temp, l_360000 = 360000L, l_6000 = 6000, l_100 = 100;

   if (*h != BADLONG)
   {
      p->degree = div_nr(*h,l_360000);
      temp = mod_nr(*h,l_360000);
      p->minute = div_nr(temp,l_6000);
      temp = mod_nr(*h,l_6000);
      p->second = div_nr(temp,l_100);
      p->hundredth = mod_nr(*h,l_100);
   }
   else
      p->degree = p->minute = p->second = p->hundredth = BADSHORT;
}

/*-----------------------------------------------------------------------------

   FUNCTION:  POSDIF

      Given two DMSH_POSITION_TYPE structures, it returns the difference 
      in hundredths of a second.

   PARAMETERS:
   
      p1 = pointer to first DMSH_POSITION_TYPE structure

      p2 = pointer to second DMSH_POSITION_TYPE structure

   RETURNS:  (p2 - p1) in hundredths of a second

   CALLED FROM:  C, FORTRAN

*/
LONG POSDIF(DMSH_POSITION_TYPE *p1, DMSH_POSITION_TYPE *p2)
{
   return(POSHUN(p2) - POSHUN(p1));
}

/*-----------------------------------------------------------------------------

   FUNCTION:  DIFPOS

      Given a reference DMSH_POSITION_TYPE structure and a displacement 
      in hundredths of a second, it calculates a new DMSH_POSITION_TYPE
      structure as the sum.

   PARAMETERS:
   
      p1 = pointer to reference DMSH_POSITION_TYPE structure

      p2 = pointer to resulting DMSH_POSITION_TYPE structure (p1 + s)

      h = pointer to displacement in hundredths of a second

   RETURNS:  VOID

   CALLED FROM:  C, FORTRAN

*/
void DIFPOS(DMSH_POSITION_TYPE *p1, DMSH_POSITION_TYPE *p2, LONG *h)
{
   LONG hh;

   hh = POSHUN(p1) + (*h);
   HUNPOS(p2, &hh);
}

/*-----------------------------------------------------------------------------

   FUNCTION:  POSCMP

      It compares two DMSH_POSITION_TYPE structures.

   PARAMETERS:
   
      p1 = pointer to first DMSH_POSITION_TYPE structure
          
      p2 = pointer to second DMSH_POSITION_TYPE structure

   RETURNS:  0 if p1 = p2
             1 if p1 < p2
             2 if p2 < p1

   CALLED FROM:  C, FORTRAN

*/
LONG POSCMP(DMSH_POSITION_TYPE *p1, DMSH_POSITION_TYPE *p2)
{
   if (p1->degree    < p2->degree   ) return(1);
   if (p1->degree    > p2->degree   ) return(2);
   if (p1->minute    < p2->minute   ) return(1);
   if (p1->minute    > p2->minute   ) return(2);
   if (p1->second    < p2->second   ) return(1);
   if (p1->second    > p2->second   ) return(2);
   if (p1->hundredth < p2->hundredth) return(1);
   if (p1->hundredth > p2->hundredth) return(2);
   return(0);
}