nmxp_base.c 19.2 KB
Newer Older
1 2
/*! \file
 *
Matteo Quintiliani's avatar
Matteo Quintiliani committed
3
 * \brief Base for Nanometrics Protocol Library
4 5 6 7 8 9
 *
 * Author:
 * 	Matteo Quintiliani
 * 	Istituto Nazionale di Geofisica e Vulcanologia - Italy
 *	quintiliani@ingv.it
 *
10
 * $Id: nmxp_base.c,v 1.71 2008-03-31 11:32:36 mtheo Exp $
Matteo Quintiliani's avatar
Matteo Quintiliani committed
11
 *
12 13
 */

14
#include "config.h"
15
#include "nmxp_base.h"
16
#include "nmxp_memory.h"
Matteo Quintiliani's avatar
Matteo Quintiliani committed
17
#ifdef HAVE_WINDOWS_H
18
#include "nmxp_win.h"
Matteo Quintiliani's avatar
Matteo Quintiliani committed
19
#endif
20 21 22 23 24 25

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <stdarg.h>

26 27 28 29
#ifdef HAVE_WINDOWS_H
#include "winsock2.h"
#warning You are compiling on Windows MinGW...
#else
30 31 32 33 34 35 36
#include <sys/socket.h>
#include <netinet/in.h>
#include <arpa/inet.h>
#include <netdb.h>
#include <errno.h>
#include <time.h>
#include <unistd.h>
37
#endif
Matteo Quintiliani's avatar
Matteo Quintiliani committed
38

Matteo Quintiliani's avatar
Matteo Quintiliani committed
39 40 41
#define MAX_OUTDATA 4096


42 43 44 45 46 47 48
int nmxp_openSocket(char *hostname, int portNum)
{
  static int sleepTime = 1;
  int isock = 0;
  struct hostent *hostinfo = NULL;
  struct sockaddr_in psServAddr;
  struct in_addr hostaddr;
49

50
#ifdef HAVE_WINDOWS_H
51
  nmxp_initWinsock();
52
#endif
53

54 55
  if (!hostname)
  {
56
    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Empty host name?\n");
57 58 59 60
    return -1;
  }

  if ( (hostinfo = gethostbyname(hostname)) == NULL) {
61 62
    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Cannot lookup host: %s\n",
	    NMXP_LOG_STR(hostname));
63 64 65 66 67 68
    return -1;
  }

  while(1)
  {
    isock = socket (AF_INET, SOCK_STREAM, 0);
69 70 71 72 73 74 75 76

#ifdef HAVE_WINDOWS_H
    if (isock == INVALID_SOCKET) {
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW,"Error at socket()\n");
	    WSACleanup();
	    exit(1);
    }
#else
77 78
    if (isock < 0)
    {
79
      nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Can't open stream socket\n");
80 81
      exit(1);
    }
82
#endif
83 84 85 86 87 88

    /* Fill in the structure "psServAddr" with the address of server
       that we want to connect with */
    memset (&psServAddr, 0, sizeof(psServAddr));
    psServAddr.sin_family = AF_INET;
    psServAddr.sin_port = htons((unsigned short) portNum);
89 90 91 92 93
#ifdef HAVE_WINDOWS_H
    unsigned long address;
    memcpy(&address, hostinfo->h_addr, (size_t) hostinfo->h_length);
    psServAddr.sin_addr.s_addr = address;
#else
94
    psServAddr.sin_addr = *(struct in_addr *)hostinfo->h_addr_list[0];
95
#endif
96 97 98 99

    /* Report action and resolved address */
    memcpy(&hostaddr.s_addr, *hostinfo->h_addr_list,
	   sizeof (hostaddr.s_addr));
100
    nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW, "Attempting to connect to %s port %d\n",
101
	    NMXP_LOG_STR(inet_ntoa(hostaddr)), portNum);
102 103 104 105

    if (connect(isock, (struct sockaddr *)&psServAddr, sizeof(psServAddr)) >= 0)
    {
      sleepTime = 1;
106
      nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW, "Connection established: socket=%i,IP=%s,port=%d\n",
107
	      isock, NMXP_LOG_STR(inet_ntoa(hostaddr)), portNum);
108 109 110 111
      return isock;
    }
    else
    {
112
      nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Connecting to %s port %d. Trying again after %d seconds...\n",
113
	      NMXP_LOG_STR(inet_ntoa(hostaddr)), portNum, sleepTime);
114
      nmxp_closeSocket(isock);
115
      nmxp_sleep (sleepTime);
116 117 118 119 120 121 122 123
      sleepTime *= 2;
      if (sleepTime > NMXP_SLEEPMAX)
        sleepTime = NMXP_SLEEPMAX;
    }
  }
}


124 125 126 127 128 129 130 131 132 133
int nmxp_closeSocket(int isock) {
	int ret;
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_CONNFLOW, "Closed connection.\n");
#ifdef HAVE_WINDOWS_H
	ret = closesocket(isock);
	WSACleanup();
#else
	ret = close(isock);
#endif
	return ret;
134 135 136 137 138 139 140 141 142 143 144 145 146
}


int nmxp_send_ctrl(int isock, void* buffer, int length)
{
  int sendCount = send(isock, (char*) buffer, length, 0);
  
  if (sendCount != length)
    return NMXP_SOCKET_ERROR;

  return NMXP_SOCKET_OK;
}

147 148 149 150 151 152 153 154
#ifdef HAVE_BROKEN_SO_RCVTIMEO
#warning Managing non-blocking I/O using select()
int nmxp_recv_select_timeout(int s, char *buf, int len, int timeout)
{
    fd_set fds;
    int n;
    struct timeval tv;

155
    /* set up the file descriptor set*/
156 157 158
    FD_ZERO(&fds);
    FD_SET(s, &fds);

159
    /* set up the struct timeval for the timeout*/
160 161 162
    tv.tv_sec = timeout;
    tv.tv_usec = 0;

163
    /* wait until timeout or data received*/
Matteo Quintiliani's avatar
Matteo Quintiliani committed
164
    errno = 0;
165
    n = select(s+1, &fds, NULL, NULL, &tv);
166
    if (n == 0) return -2; /* timeout!*/
167
    if(errno == EINTR) return -2; /* timeout! "Interrupted system call" */
168
    if (n == -1) return -1; /* error*/
169
 
170
    /* data must be here, so do a normal recv()*/
171 172 173
    return recv(s, buf, len, 0);
}
#endif
174

175 176
int nmxp_setsockopt_RCVTIMEO(int isock, int timeoutsec) {
    int ret = 0;
177
#ifdef HAVE_WINDOWS_H
178
    int timeos;
179
#else
180 181 182

#ifndef HAVE_BROKEN_SO_RCVTIMEO
    struct timeval timeo;
183
#endif
184

185
#endif
186

187
    if(timeoutsec > 0) {
188
#ifdef HAVE_WINDOWS_H
189 190 191 192 193 194
	timeos  = timeoutsec * 1000;
	ret = setsockopt(isock, SOL_SOCKET, SO_RCVTIMEO, (char *) &timeos, sizeof(timeos));
	if (ret < 0)
	{
	    perror("setsockopt SO_RCVTIMEO");
	}
195
#else
196 197

#ifndef HAVE_BROKEN_SO_RCVTIMEO
198 199 200 201 202 203 204 205
	timeo.tv_sec  = timeoutsec;
	timeo.tv_usec = 0;
	ret = setsockopt(isock, SOL_SOCKET, SO_RCVTIMEO, &timeo, sizeof(timeo));
	if (ret < 0) {
	    perror("setsockopt SO_RCVTIMEO");
	}
#else
#warning nmxp_setsockopt_RCVTIMEO() do nothing for your system.
206
#endif
207

208
#endif
209
    }
210

211 212
    return ret;
}
213

Matteo Quintiliani's avatar
Matteo Quintiliani committed
214

215 216 217 218
#define MAXLEN_RECV_ERRNO_STR 200

char *nmxp_strerror(int errno_value) {
    static char ret_recv_errno_str[MAXLEN_RECV_ERRNO_STR];
219
#ifdef HAVE_WINDOWS_H
220
    char *recv_errno_str;
221
#else
222

223
#ifdef HAVE_STRERROR_R
224 225 226 227 228 229 230 231 232 233 234 235 236 237 238
    char recv_errno_str[MAXLEN_RECV_ERRNO_STR];
#else
    char *recv_errno_str;
#endif

#endif

    ret_recv_errno_str[0] = 0;

#ifdef HAVE_WINDOWS_H
    recv_errno_str = WSAGetLastErrorMessage(errno_value);
#else

#ifdef HAVE_STRERROR_R
    strerror_r(errno_value, recv_errno_str, MAXLEN_RECV_ERRNO_STR);
239
#else
240
    recv_errno_str = strerror(errno_value);
241 242
#endif

243
#endif
244

245 246 247 248 249 250 251 252 253 254 255 256 257 258 259
    if(recv_errno_str) {
	strncpy(ret_recv_errno_str, recv_errno_str, MAXLEN_RECV_ERRNO_STR);
    }

    return ret_recv_errno_str;
}


int nmxp_recv_ctrl(int isock, void *buffer, int length, int timeoutsec, int *recv_errno )
{
  int recvCount;
  int cc;
  char *buffer_char = buffer;
  char *recv_errno_str = NULL;

260
  nmxp_setsockopt_RCVTIMEO(isock, timeoutsec);
261
  
262
  cc = 1;
263 264
  *recv_errno  = 0;
  recvCount = 0;
265
  while(cc > 0 && *recv_errno == 0  && recvCount < length) {
266

Matteo Quintiliani's avatar
Matteo Quintiliani committed
267
      /* TODO some operating system could not reset errno */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
268
      errno = 0;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
269

270 271 272 273 274 275 276 277 278 279 280 281
#ifdef HAVE_BROKEN_SO_RCVTIMEO
      if(timeoutsec == 0) {
#endif

	  cc = recv(isock, buffer_char + recvCount, length - recvCount, 0);

#ifdef HAVE_BROKEN_SO_RCVTIMEO
      } else {
	  cc = nmxp_recv_select_timeout(isock, buffer_char + recvCount, length - recvCount, timeoutsec);
      }
#endif

282 283 284
#ifdef HAVE_WINDOWS_H
      *recv_errno  = WSAGetLastError();
#else
285 286 287 288 289
      if(cc == -2) {
	  *recv_errno  = EWOULDBLOCK;
      } else {
	  *recv_errno  = errno;
      }
290
#endif
291
      if(cc <= 0) {
292
	  /*
293 294
	  nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "nmxp_recv_ctrl(): (cc=%d <= 0) errno=%d  recvCount=%d  length=%d\n",
	  cc, *recv_errno, recvCount, length);
295
	  */
296 297 298
      } else {
	  recvCount += cc;
      }
299
  }
300

301
  nmxp_setsockopt_RCVTIMEO(isock, 0);
302

303
  if (recvCount != length  ||  *recv_errno != 0  ||  cc <= 0) {
304

305
      recv_errno_str = nmxp_strerror(*recv_errno);
306

307
#ifdef HAVE_WINDOWS_H
308
      if(*recv_errno != WSAEWOULDBLOCK  &&  *recv_errno != WSAETIMEDOUT)
309 310 311 312
#else
      if(*recv_errno != EWOULDBLOCK)
#endif
      {
313 314
	  nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "nmxp_recv_ctrl(): %s (errno=%d recvCount=%d length=%d cc=%d)\n",
		  NMXP_LOG_STR(recv_errno_str), *recv_errno, recvCount, length, cc);
315 316 317 318 319 320 321 322
      }

      /* TO IMPROVE 
       * Fixed bug receiving zero byte from recv() 'TCP FIN or EOF received'
       * */
      if(cc == 0  &&  *recv_errno == 0) {
	  *recv_errno = -100;
      }
323

324
#ifdef HAVE_WINDOWS_H
325
      if(recvCount != length || (*recv_errno != WSAEWOULDBLOCK  &&  *recv_errno != WSAETIMEDOUT))
326
#else
327
      if(recvCount != length || *recv_errno != EWOULDBLOCK)
328
#endif
329 330 331
      {
	  return NMXP_SOCKET_ERROR;
      }
332 333 334 335 336 337
  }
  
  return NMXP_SOCKET_OK;
}


338
int nmxp_sendHeader(int isock, NMXP_MSG_CLIENT type, int32_t length)
339
{  
Matteo Quintiliani's avatar
Matteo Quintiliani committed
340
    NMXP_MESSAGE_HEADER msg;
341 342 343 344 345

    msg.signature = htonl(NMX_SIGNATURE);
    msg.type      = htonl(type);
    msg.length    = htonl(length);

Matteo Quintiliani's avatar
Matteo Quintiliani committed
346
    return nmxp_send_ctrl(isock, &msg, sizeof(NMXP_MESSAGE_HEADER));
347 348 349
}


350
int nmxp_receiveHeader(int isock, NMXP_MSG_SERVER *type, int32_t *length, int timeoutsec, int *recv_errno )
351 352
{  
    int ret ;
Matteo Quintiliani's avatar
Matteo Quintiliani committed
353
    NMXP_MESSAGE_HEADER msg;
354

355
    ret = nmxp_recv_ctrl(isock, &msg, sizeof(NMXP_MESSAGE_HEADER), timeoutsec, recv_errno);
356 357 358 359 360 361 362 363 364 365 366 367

    *type = 0;
    *length = 0;

    if(ret == NMXP_SOCKET_OK) {
	msg.signature = ntohl(msg.signature);
	msg.type      = ntohl(msg.type);
	msg.length    = ntohl(msg.length);

	if (msg.signature != NMX_SIGNATURE)
	{
	    ret = NMXP_SOCKET_ERROR;
368 369
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW,
		    "nmxp_receiveHeader(): signature mismatches. signature = %d, type = %d, length = %d\n",
370 371 372 373 374 375 376 377 378 379 380
		    msg.signature, msg.type, msg.length);
	} else {
	    *type = msg.type;
	    *length = msg.length;
	}
    }

    return ret;
}


381
int nmxp_sendMessage(int isock, NMXP_MSG_CLIENT type, void *buffer, int32_t length) {
382 383 384 385 386 387 388 389 390 391 392
    int ret;
    ret = nmxp_sendHeader(isock, type, length);
    if( ret == NMXP_SOCKET_OK) {
	if(buffer && length > 0) {
	    ret = nmxp_send_ctrl(isock, buffer, length);
	}
    }
    return ret;
}


393 394 395 396 397 398 399 400 401 402 403 404 405 406
int32_t nmxp_display_error_from_server(char *buffer, int32_t length) {
    /* NMXP_MSG_TERMINATESUBSCRIPTION */
    char *str_msg = NULL;
    int32_t reason;
    memcpy(&reason, buffer, sizeof(reason));
    reason = ntohl(reason);
    str_msg = buffer + sizeof(reason);
    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_ANY, "%d %s shutdown: %s\n",
	    reason,
	    (reason == 0)? "Normal" : (reason == 1)? "Error" : (reason == 2)? "Timeout" : "Unknown",
	    str_msg);
    return reason;
}

407
int nmxp_receiveMessage(int isock, NMXP_MSG_SERVER *type, void **buffer, int32_t *length, int timeoutsec, int *recv_errno ) {
408 409 410 411
    int ret;
    *buffer = NULL;
    *length = 0;

412
    ret = nmxp_receiveHeader(isock, type, length, timeoutsec, recv_errno);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
413

Matteo Quintiliani's avatar
Matteo Quintiliani committed
414
    if( ret == NMXP_SOCKET_OK  ) {
Matteo Quintiliani's avatar
Matteo Quintiliani committed
415
	if (*length > 0) {
416
	    *buffer = NMXP_MEM_MALLOC(*length);
417
	    ret = nmxp_recv_ctrl(isock, *buffer, *length, 0, recv_errno);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
418

419 420 421 422
	    if(*type == NMXP_MSG_TERMINATESUBSCRIPTION) {
		nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Received TerminateSubscritption.\n");
		nmxp_display_error_from_server(*buffer, *length);
	    } else if(*type == NMXP_MSG_ERROR) {
423
		nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Received ErrorMessage: %s\n", NMXP_LOG_STR(*buffer));
Matteo Quintiliani's avatar
Matteo Quintiliani committed
424
	    } else {
425
		nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_PACKETMAN, "Received message type: %d  length=%d\n", *type, *length);
Matteo Quintiliani's avatar
Matteo Quintiliani committed
426 427 428
	    }

	}
Matteo Quintiliani's avatar
Matteo Quintiliani committed
429 430 431
    }

    if(*recv_errno != 0) {
432
#ifdef HAVE_WINDOWS_H
433
	if(*recv_errno == WSAEWOULDBLOCK  ||  *recv_errno == WSAETIMEDOUT) {
434
#else
435
	if(*recv_errno == EWOULDBLOCK) {
436
#endif
437
	    nmxp_log(NMXP_LOG_WARN, NMXP_LOG_D_DOD, "Timeout receiving in nmxp_receiveMessage()\n");
Matteo Quintiliani's avatar
Matteo Quintiliani committed
438
	} else {
439 440
	    /* Log message is not necessary because managed by nmxp_recv_ctrl() */
	    /* nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_CONNFLOW, "Error in nmxp_receiveMessage()\n"); */
441
	}
442
    }
Matteo Quintiliani's avatar
Matteo Quintiliani committed
443

444 445 446
    return ret;
}

447

448
NMXP_DATA_PROCESS *nmxp_processDecompressedData(char* buffer_data, int length_data, NMXP_CHAN_LIST_NET *channelList, const char *network_code_default)
449 450 451 452 453 454
{
  int32_t   netInt    = 0;
  int32_t   pKey      = 0;
  double    pTime     = 0.0;
  int32_t   pNSamp    = 0;
  int32_t   pSampRate = 0;
455
  int32_t  *pDataPtr  = NULL;
456 457
  int       swap      = 0;
  int       idx;
458
  static int32_t outdata[MAX_OUTDATA];
459

460 461
  char station_code[20];
  char channel_code[20];
462
  char network_code[20];
463

464
  char *nmxp_channel_name = NULL;
465 466 467 468 469 470 471
  static NMXP_DATA_PROCESS pd;

  /* copy the header contents into local fields and swap */
  memcpy(&netInt, &buffer_data[0], 4);
  pKey = ntohl(netInt);
  if ( pKey != netInt ) { swap = 1; }

472 473 474 475 476 477
  nmxp_data_init(&pd);

  nmxp_channel_name = nmxp_chan_lookupName(pKey, channelList);

  if(nmxp_channel_name) {

478
  memcpy(&pTime, &buffer_data[4], 8);
479
  if ( swap ) { nmxp_data_swap_8b(&pTime); }
480 481 482 483 484 485 486

  memcpy(&netInt, &buffer_data[12], 4);
  pNSamp = ntohl(netInt);
  memcpy(&netInt, &buffer_data[16], 4);
  pSampRate = ntohl(netInt);

  /* There should be (length_data - 20) bytes of data as 32-bit ints here */
Matteo Quintiliani's avatar
Matteo Quintiliani committed
487 488
  memcpy(outdata , (int32_t *) &buffer_data[20], length_data - 20);
  pDataPtr = outdata;
489 490 491 492 493 494 495

  /* Swap the data samples to host order */
  for ( idx=0; idx < pNSamp; idx++ ) {
      netInt = ntohl(pDataPtr[idx]);
      pDataPtr[idx] = netInt;
  }

496
  if(!nmxp_chan_cpy_sta_chan(nmxp_channel_name, station_code, channel_code, network_code)) {
497 498
    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not in STA.CHAN format: %s\n",
	    NMXP_LOG_STR(nmxp_channel_name));
499 500 501
  }
  
  pd.key = pKey;
502
  if(network_code[0] != 0) {
503
      strcpy(pd.network, network_code);
504 505
  } else {
      strcpy(pd.network, network_code_default);
506
  }
507 508
  if(station_code[0] != 0) {
      strncpy(pd.station, station_code, STATION_LENGTH);
509
  }
510 511
  if(channel_code[0] != 0) {
      strncpy(pd.channel, channel_code, CHANNEL_LENGTH);
512
  }
513
  pd.packet_type = NMXP_MSG_DECOMPRESSED;
514 515 516
  pd.x0 = -1;
  pd.xn = -1;
  pd.x0n_significant = 0;
517 518 519
  /* TODO*/
  /* pd.oldest_seq_no = ;*/
  /* pd.seq_no = ;*/
520 521 522 523 524 525 526
  pd.time = pTime;
  pd.buffer = buffer_data;
  pd.length = length_data;
  pd.nSamp = pNSamp;
  pd.pDataPtr = pDataPtr;
  pd.sampRate = pSampRate;

527 528
  } else {
      nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not found for key %d\n", pKey);
529 530
  }

531 532 533 534
  return &pd;
}


535
NMXP_DATA_PROCESS *nmxp_processCompressedData(char* buffer_data, int length_data, NMXP_CHAN_LIST_NET *channelList, const char *network_code_default)
536 537 538 539 540
{
    int32_t   pKey      = 0;
    double    pTime     = 0.0;
    int32_t   pNSamp    = 0;
    int32_t   pSampRate = 0;
541
    int32_t  *pDataPtr  = NULL;
542

543 544
    char station_code[20];
    char channel_code[20];
545
    char network_code[20];
546 547 548

    static NMXP_DATA_PROCESS pd;

549
    int32_t nmx_rate_code_to_sample_rate[32] = {
550 551 552 553 554
	0,1,2,5,10,20,40,50,
	80,100,125,200,250,500,1000,25,
	120,0,0,0,0,0,0,0,
	0,0,0,0,0,0,0,0};

555
	int32_t nmx_oldest_sequence_number;
556 557
	char nmx_hdr[25];
	unsigned char nmx_ptype;
558
	int32_t nmx_seconds;
559
	double nmx_seconds_double;
560 561
	int16_t nmx_ticks, nmx_instr_id;
	int32_t nmx_seqno;
562
	unsigned char nmx_sample_rate;
563 564
	int32_t nmx_x0;
	int32_t rate_code, chan_code, this_sample_rate;
565

566
	int32_t comp_bytecount;
567
	unsigned char *indata;
568
	static int32_t outdata[MAX_OUTDATA];
569 570
	int32_t nout, i, k;
	int32_t prev_xn;
571 572
	const uint32_t high_scale = 4096 * 2048;
	const uint32_t high_scale_p = 4096 * 4096;
573

574 575
	char *nmxp_channel_name = NULL;

576
	/* TOREMOVE int my_order = get_my_wordorder();*/
Matteo Quintiliani's avatar
Matteo Quintiliani committed
577
	int my_host_is_bigendian = nmxp_data_bigendianhost();
578
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "my_host_is_bigendian %d\n", my_host_is_bigendian);
579 580

	memcpy(&nmx_oldest_sequence_number, buffer_data, 4);
581 582 583
	if (my_host_is_bigendian) {
	    nmxp_data_swap_4b (&nmx_oldest_sequence_number);
	}
584
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "Oldest sequence number = %d\n", nmx_oldest_sequence_number);
585 586 587 588 589 590

	memcpy(nmx_hdr, buffer_data+4, 17);
	/* Decode the Nanometrics packet header bundle. */
	memcpy (&nmx_ptype, nmx_hdr+0, 1);
	if ( (nmx_ptype & 0xf) == 9) {
	    /* Filler packet.  Discard entire packet.   */
591
	    nmxp_log (NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Filler packet - discarding\n");
592
	    /*m continue;*/
593 594 595 596 597 598 599 600 601 602 603
	    exit(0);
	}

	nmx_x0 = 0;
	memcpy (&nmx_seconds, nmx_hdr+1, 4);
	memcpy (&nmx_ticks, nmx_hdr+5, 2);
	memcpy (&nmx_instr_id, nmx_hdr+7, 2);
	memcpy (&nmx_seqno, nmx_hdr+9, 4);
	memcpy (&nmx_sample_rate, nmx_hdr+13, 1);
	memcpy (&nmx_x0, nmx_hdr+14, 3);

604
	if (my_host_is_bigendian) {
605
	    nmxp_data_swap_4b ((int32_t *)&nmx_seconds);
606 607 608 609
	    nmxp_data_swap_2b (&nmx_ticks);
	    nmxp_data_swap_2b (&nmx_instr_id);
	    nmxp_data_swap_4b (&nmx_seqno);
	    nmxp_data_swap_4b (&nmx_x0);
610
	}
611 612 613 614 615 616 617

	/* check if nmx_x0 is negative like as signed 3-byte int */
	if( (nmx_x0 & high_scale) ==  high_scale) {
	    /* nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "WARNING: changed nmx_x0, old value = %d\n",  nmx_x0);*/
	    nmx_x0 -= high_scale_p;
	}

618 619 620 621 622
	nmx_seconds_double = (double) nmx_seconds + ( (double) nmx_ticks / 10000.0 );
	rate_code = nmx_sample_rate>>3;
	chan_code = nmx_sample_rate&7;
	this_sample_rate = nmx_rate_code_to_sample_rate[rate_code];

623 624 625
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_ptype          = %d\n", nmx_ptype);
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_seconds        = %d\n", nmx_seconds);
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_ticks          = %d\n", nmx_ticks);
626

627 628
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_seconds_double = %f\n", nmx_seconds_double);
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_x0             = %d\n", nmx_x0);
629

630 631 632
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_instr_id       = %d\n", nmx_instr_id);
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_seqno          = %d\n", nmx_seqno);
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "nmx_sample_rate    = %d\n", nmx_sample_rate);
633
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "this_sample_rate   = %d\n", this_sample_rate);
634

635 636 637 638 639 640
	pKey = (nmx_instr_id << 16) | ( 1 << 8) | ( chan_code);

	pTime = nmx_seconds_double;

	pSampRate = this_sample_rate;

641 642 643 644 645 646 647
	nmxp_data_init(&pd);

	nmxp_channel_name = nmxp_chan_lookupName(pKey, channelList);

	if(nmxp_channel_name) {

	if(!nmxp_chan_cpy_sta_chan(nmxp_channel_name, station_code, channel_code, network_code)) {
648 649
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not in STA.CHAN format: %s\n",
		    NMXP_LOG_STR(nmxp_channel_name));
650 651
	}
  
652 653
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "Channel key %d for %s.%s\n",
		pKey, NMXP_LOG_STR(station_code), NMXP_LOG_STR(channel_code));
654

655 656 657
	comp_bytecount = length_data-21;
	indata = (unsigned char *) buffer_data + 21;

658 659
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "comp_bytecount     = %d  (N = %.2f)\n", comp_bytecount, (double) comp_bytecount / 17.0);

660 661 662 663 664 665
	/* Unpack the data bundles, each 17 bytes long. */
	prev_xn = nmx_x0;
	outdata[0] = nmx_x0;
	nout = 1;
	for (i=0; i<comp_bytecount; i+=17) {
	    if (i+17>comp_bytecount) {
666
		nmxp_log (NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "comp_bytecount = %d, i+17 = %d\n",
667 668 669 670
			comp_bytecount, i+17);
		exit(1);
	    }
	    if (nout+16 > MAX_OUTDATA)  {
671
		nmxp_log (NMXP_LOG_ERR,  NMXP_LOG_D_PACKETMAN, "Output buffer size too small\n");
672 673 674
		exit(1);
	    }
	    k = nmxp_data_unpack_bundle (outdata+nout,indata+i,&prev_xn);
675 676 677 678 679
	    if (k < 0) nmxp_log (NMXP_LOG_WARN, NMXP_LOG_D_PACKETMAN, "Null bundle: %s.%s.%s (k=%d) %s %d\n",
		    NMXP_LOG_STR(network_code),
		    NMXP_LOG_STR(station_code),
		    NMXP_LOG_STR(channel_code), k,
		    __FILE__,  __LINE__);
680 681 682 683 684 685
	    if (k < 0) break;
	    nout += k;
	    /* prev_xn = outdata[nout-1]; */
	}
	nout--;

686
	nmxp_log(NMXP_LOG_NORM, NMXP_LOG_D_PACKETMAN, "Unpacked %d samples.\n", nout);
687 688 689 690 691 692

	pDataPtr = outdata;

	pNSamp = nout;

	pd.key = pKey;
693
	if(network_code[0] != 0) {
694
	    strcpy(pd.network, network_code);
695 696
	} else {
	    strcpy(pd.network, network_code_default);
697
	}
698 699
	if(station_code[0] != 0) {
	    strncpy(pd.station, station_code, STATION_LENGTH);
700
	}
701 702
	if(channel_code[0] != 0) {
	    strncpy(pd.channel, channel_code, CHANNEL_LENGTH);
703
	}
704 705
	pd.packet_type = nmx_ptype;
	pd.x0 = nmx_x0;
706
	pd.xn = pDataPtr[nout];
707
	pd.x0n_significant = 1;
708
	pd.oldest_seq_no = nmx_oldest_sequence_number;
709 710 711 712 713 714 715 716
	pd.seq_no = nmx_seqno;
	pd.time = pTime;
	pd.buffer = buffer_data;
	pd.length = length_data;
	pd.nSamp = pNSamp;
	pd.pDataPtr = pDataPtr;
	pd.sampRate = pSampRate;

717 718
	} else {
	    nmxp_log(NMXP_LOG_ERR, NMXP_LOG_D_PACKETMAN, "Channel name not found for key %d\n", pKey);
719 720
	}

721 722 723
	return &pd;
}

724 725 726

unsigned int nmxp_sleep(unsigned int sleep_time) {
#ifdef HAVE_WINDOWS_H
727
    Sleep(sleep_time * 1000);
728
    return 0;
729 730 731 732 733
#else
    return sleep(sleep_time);
#endif
}