#include "stdio.h"
#include <string.h>
#include <stdlib.h>
#include <time.h>
#include "dcp3rd.h"
void p3_read(char *bullx, int *lenbx, int *ibxptr, int *iunbfo, int *iret)
/************************************************************************
 * p3_read(bullx, lenbx, ibxptr, iret)                                  *
 *                                                                      *
 * This program reads binary p3 radar data and write them into BUFR     *
 * output.                                                              *
 * Input parameters:                                                    *
 *      bullx           char*           Bulletin                        *
 *      lenbx           int             Length of bullx                 *
 *      ibxptr          int             Pointer within bullx            *
 *                                                                      *
 * Output parameters:                                                   *
 *      iret            int             Return code                     *
 *                                        0 = normal return             *
 *                                       -1 = no more reports in        *
 *                                            bulletin                  *
 *                                                                      *
 **                                                                     *
 * Log:                                                                 *
 * S. Guan/NCEP         06/09  Initial version                          *
 ***********************************************************************/

{
   FILE *in;
   char buf[80], message[120];
   P3_header_t *p3h;
   double *dpt;
   float *fpt;
   short *spt;
   int  *ipt;
   P3_header_t vp3h;
   double time_stamp;
   struct tm  *current;
   time_t time_secs; 
   int lunout, tref, tspeed, twidth;
   int    hrs;
   int    mins;
   int    secs;
   int    year;
   int    month;
   int    day;

/*   array[5] = Radar Latitude, Radar Longitude, Radar Altitude, 
                Earth Relative Azimuth, Earth Relative Incidence  */
   float lat, lon, altitude, azimuth, incidence, array[5];
   char acid[14], acrn[10], stmid[10]; 
   int lacid, lacrn, lstmid;
   short ref[NMSGSIZE], refws[3][NMSGSIZE];
   int i, j, ier, blength, totlenth;
   short  vshort;

   *iret = 0;
   totlenth = *lenbx;
   blength = *ibxptr;
   if ( (blength + 22) > totlenth)
   {
      *iret = 0;
      sprintf(message, "End of this report \n");
      dc_wclg ( 0, "DC", 2, message, &ier );
      return;
   }

   p3h = &vp3h;
   dpt = &time_stamp;
   fpt = &lat;
   spt = &vshort;
   ipt = &i;

   lunout = *iunbfo;

   ipt = (int *) (bullx + blength );
   vp3h.flight_date = *ipt;
   spt = ( short * ) (bullx + blength + 4 );
   vp3h.aircraft_id = *spt;
   spt = ( short * ) (bullx + blength + 6 );
   vp3h.storm_no = *spt;
   spt = ( short * ) (bullx + blength + 8 );
   vp3h.bins  = *spt;
   fpt = (float * ) (bullx + blength + 10 );
   vp3h.range_delay = *fpt;
   fpt = (float * ) (bullx + blength + 14 );
   vp3h.bin_resolution = *fpt;
   spt = ( short * ) (bullx + blength + 18 );
   vp3h.data_type = *spt; 
   spt = ( short * ) (bullx + blength + 20 );
   vp3h.antenna_no  = *spt;

   blength =  blength + 22;

   vp3h.flight_date = reverse_int4(vp3h.flight_date);
   vp3h.aircraft_id =  reverse_int2(vp3h.aircraft_id);
   vp3h.storm_no =  reverse_int2(vp3h.storm_no);
   vp3h.bins  =  reverse_int2(vp3h.bins);
   vp3h.range_delay = reverse_float(vp3h.range_delay);
   vp3h.bin_resolution = reverse_float(vp3h.bin_resolution);
   vp3h.data_type =  reverse_int2(vp3h.data_type);
   vp3h.antenna_no =  reverse_int2(vp3h.antenna_no);
   if ( (vp3h.aircraft_id == 42) || (vp3h.aircraft_id == 43) || (vp3h.aircraft_id == 49))
   {
      if ( (vp3h.antenna_no < 0) || (vp3h.antenna_no > 10))
      {
         sprintf(message, "Wrong Antenna Number. \n");
         dc_wclg ( 0, "DC", 2, message, &ier );
         *iret = -1;
         return;
      }
   }
   else
   {
      sprintf(message, "Data are not from N42RF  or N43RF or N49RF. \n");
      dc_wclg ( 0, "DC", 2, message, &ier );
      *iret = -1;
      return; 
   }

/*  Decide data type from Radial Data Type Flag */
   tref = 0;
   tspeed = 0;
   twidth = 0;
   if (vp3h.data_type == 1 || vp3h.data_type == 3 || vp3h.data_type == 5 || vp3h.data_type == 7)
   {
      tref = 1;
   }
   if (vp3h.data_type == 2 || vp3h.data_type == 3 || vp3h.data_type == 6 || vp3h.data_type == 7)
   {
      tspeed = 1;
   }
   if (vp3h.data_type >= 4 )
   {
      twidth = 1;
   }

   for(i=0; i<3; i++) 
   for (j=0; j<vp3h.bins; j++)
   {
       refws[i][j] = -8888;
   }

   while (  (blength + 34) < totlenth ) 
   {
      dpt = ( double * ) ( bullx + blength );
      time_stamp = *dpt;
      blength = blength + 8;
      time_stamp = reverse_double(time_stamp);
      time_secs = (time_t) time_stamp;
      current = gmtime ( &time_secs );
      year = (*current).tm_year + 1900;
      month = (*current).tm_mon + 1;
      day = (*current).tm_mday;
      hrs = (*current).tm_hour;
      mins =(*current).tm_min;
      secs =(*current).tm_sec;

      for(i=0; i<5; i++)
      {
         fpt = (float * ) (bullx + blength);
         array[i] = *fpt;
         blength = blength + 4;
      }
      swp_float(5, array); 

      if ( tref == 1 )
      {
         for(i=0; i< vp3h.bins; i++)
         {
            spt = ( short * ) (bullx + blength);
            ref[i] = *spt;
            blength = blength + 2;
         }
         swp_int2(vp3h.bins, ref); 
         for (j=0; j<vp3h.bins; j++)
             refws[0][j] = ref[j];
      }  
      if ( tspeed == 1 )
      {
         for(i=0; i< vp3h.bins; i++)
         {
            spt = ( short * ) (bullx + blength);
            blength = blength + 2;
            ref[i] = *spt;
         }
         swp_int2(vp3h.bins, ref);
         for (j=0; j<vp3h.bins; j++)
         {
             refws[1][j] = ref[j];
         }
      }
      if ( twidth == 1 )
      {
         for(i=0; i< vp3h.bins; i++)
         {
            spt = ( short * ) (bullx + blength);
            blength = blength + 2;
            ref[i] = *spt;
         }
         swp_int2(vp3h.bins, ref);
         for (j=0; j<vp3h.bins; j++)
             refws[2][j] = ref[j];
      }

      lat = array[0];
      lon = array[1];
      altitude = array[2];
      azimuth = array[3];
      incidence = array[4];
      sprintf(acid, "%d",  vp3h.flight_date);
      sprintf(acrn, "%hd",  vp3h.aircraft_id);
      sprintf(stmid, "%hd",  vp3h.storm_no);
      lacid = strlen(acid);
      lacrn = strlen(acrn);
      lstmid = strlen(stmid);
      sprintf(message, "Bin = %d;  Range_delay = %f;   Resolution = %f\n", vp3h.bins, vp3h.range_delay, vp3h.bin_resolution);
      dc_wclg ( 3, "DC", 2, message, &ier ); 
      sprintf(message, "Year = %d; Month = %d; Day = %d; Hour = %d; Minute = %d; Second = %d\n", year, month, day, hrs, mins, secs);
      dc_wclg ( 3, "DC", 2, message, &ier );
      sprintf(message, "Latitude = %f; Longitude = %f;  Altitude = %f\n", lat, lon, altitude);
      dc_wclg ( 3, "DC", 2, message, &ier );
      sprintf(message, "Azimuth = %f; Incidence = %f\n", azimuth, incidence);
      dc_wclg ( 3, "DC", 2, message, &ier );
      p3_bout(&lunout, refws, &vp3h.bins, &vp3h.range_delay, &vp3h.bin_resolution, &year,
          &month, &day, &hrs, &mins, &secs, &lat, &lon, &altitude, &azimuth, &incidence, 
         &vp3h.antenna_no, &lacid, &lacrn, &lstmid, acid, acrn, stmid, strlen(acid), strlen(acrn), strlen(stmid));
   }
   sprintf(message, "Finished processing bulletin (%d bytes).\n", blength);
   dc_wclg ( 0, "DC", 2, message, &ier );
   *iret = -1;   

   return;
}