#!/bin/bash

###############################################################################

function bytedump2wideenc()
{
 awk -v ptr=0 -v nval=0 \
  '{	for ( i=0;i<NF;i++ )
	 {	d=$(i+1);
		if ( ptr==0 )
		 {	if ( and(d,0xC0)==0 )
			 {	arr[nval+1]=d;
				len[nval+1]=1;
				nval++;
				ptr=0;
			 }
			else if ( and(d,0xE0)==0x40 )
			 {	t=and(d,0x1F);
				l=2;
				ptr=1;
			 }
			else if ( and(d,0xF0)==0x60 )
			 {	t=and(d,0x0F);
				l=3;
				ptr=1;
			 }
			else if ( and(d,0xF8)==0x70 )
			 {	t=and(d,0x07);
				l=4;
				ptr=1;
			 }
			else if ( and(d,0xFC)==0x78 )
			 {	t=and(d,0x07);
				l=5;
				ptr=1;
			 }
			else
				ptr=0;
		 }
		else
		 {	t=or(lshift(t,7),and(d,0x7F));
			ptr++;
			if ( ptr==l )
			 {	arr[nval+1]=t;
				len[nval+1]=l;
				nval++;
				ptr=0;
			 }
		 }
	 }
  } \
 END \
  {	for ( i=0; i<nval; i++ )
	 {	print len[i+1],arr[i+1];
	 }
  }'
} 

###############################################################################

is_array=0

hexdump -v -e '16/1 "%d " "\n"' | \
bytedump2wideenc | \
awk -v offset=0 -v type= -v start_tv_sec=0 -v is_array=0 \
'BEGIN \
  {	if ( is_array )
		printf("[\n");
	is_first_tv=1;
  } \
 function cpredict(time,arr) \
  {	cmd=sprintf("/usr/local/bin/cpredict -f /var/local/tle/47959.tle -t %d --terse",time);
	cmd | getline result;
	close(cmd);
	split(result,arr);
  }
 function flush_data() \
  {	if ( type=="meta" )
	 {	printf(" { \"type\": \"meta\", \"cycle_count\": %d, \"exptime\": %.6f, \"count\": %d, \"cutoff\": %d, \"temperatures\": [ %.2f, %.2f, %.2f ] }%s\n",spectrum_count,exptime/1000000.0,count,cutoff,tmp1,tmp2, tmp3,is_array?",":"");
	 }
	else if ( type=="timing" )
	 {	tstr=strftime("%Y-%m-%dZ%H:%M:%S",tv_sec,1);
		printf(" { \"type\": \"timing\", \"time\": %d.%.6d, \"utc\": \"%s.%.6d\" }%s\n",tv_sec,tv_usec,tstr,tv_usec,is_array?",":"");
		last_tv_sec=tv_sec;
		last_tv_usec=tv_usec;
		had_timing=1;
	 }
	else if ( type=="timing_relative" )
	 {	tv_sec=last_tv_sec;
		tv_usec=last_tv_usec+tv_usec;
		while ( 1000000 <= tv_usec )
		 {	tv_usec -= 1000000;
			tv_sec += 1;
		 }
		tstr=strftime("%Y-%m-%dZ%H:%M:%S",tv_sec,1);
		printf(" { \"type\": \"timing\", \"time\": %d.%.6d, \"utc\": \"%s.%.6d\" }%s\n",tv_sec,tv_usec,tstr,tv_usec,is_array?",":"");
		last_tv_sec=tv_sec;
		last_tv_usec=tv_usec;
		had_timing=1;
	 }
	else if ( type=="spectrum" )
	 {	
		mid_tv_sec = tv_sec;
		mid_tv_usec = tv_usec + int(exptime/2);
		mid_tv_sec += int(mid_tv_usec/1000000);
		mid_tv_usec = mid_tv_usec%1000000;
		tstr=strftime("%Y-%m-%dZ%H:%M:%S",mid_tv_sec,1);

		cpredict(mid_tv_sec,arr);
		lon=arr[5];
		lat=arr[6];
		alt=arr[7];

		if ( 86400 <= tv_sec )
		 {	printf(" { \"type\": \"spectrum\", \"bin_mode\": %d, ",bin_mode);
			printf(" \"mid\": { \"time\": %d.%.6d, \"utc\": \"%s.%.6d\", \"lon\": %g, \"lat\": %g, \"alt\": %g }, ",mid_tv_sec,mid_tv_usec,tstr,mid_tv_usec,lon,lat,alt);
			printf("\"data\": [ ");
			for ( i=0; i<offset; i++ )
			 {	printf("%d",data_arr[i+1]);
				if ( i<offset-1 )
					printf(",");
			 }
			printf(" ] }%s\n",is_array?",":"");
		 }

		tv_usec += int(exptime);
		tv_sec += int(tv_usec/1000000);
		tv_usec = tv_usec%1000000;

	 }
  } \
  {	if ( $1==2 && 32<=$2 && $2<64 )
	 {	flush_data();
	 }

	if ( $1==2 && $2==42 )
	 {	type="meta";
		offset=0;
	 }
	else if ( $1==2 && $2==40 )
	 {	type="timing";
		offset=0;
	 }
	else if ( $1==2 && $2==41 )
	 {	type="timing_relative";
		offset=0;
	 }
	else if ( $1==2 && 32<=$2 && $2<40 )
	 {	type="spectrum";
		offset=0;
		bin_mode=$2-32;
	 }
	else if ( $1==2 && $2<32 && type=="spectrum" )
	 {	nzero=$2;
		for ( i=0; i<nzero ; i++ )
		 {	data_arr[offset+1]=0;
			offset++;
		 }
	 }
	else if ( $1==3 && $2<256 && type=="spectrum" )
	 {	nzero=$2;
		for ( i=0; i<nzero ; i++ )
		 {	data_arr[offset+1]=0;
			offset++;
		 }
	 }
	else if ( ( $1==1 && $2<64 ) || ( $1==2 && 64<=$2 && $2<4096 ) || ( 3<=$1 && 4096<=$2 ) )
	 {	if ( type=="meta" )
		 {	if ( offset==0 )
				spectrum_count=$2;
			else if ( offset==1 )
				exptime=$2;
			else if ( offset==2 )
				count=$2;
			else if ( offset==3 )
				cutoff=$2;
			else if ( offset==4 )
				tmp1=($2<2048?$2:$2-4096)/32.0;
			else if ( offset==5 )
				tmp2=($2<2048?$2:$2-4096)/32.0;
			else if ( offset==6 )
				tmp3=($2<2048?$2:$2-4096)/32.0;
		 }
		else if ( type=="timing" )
		 {	if ( offset==0 )
				tv_sec_2p24=$2;
			else if ( offset==1 )
			 {	tv_sec=tv_sec_2p24*16777216+$2;
				if ( is_first_tv && 0<start_tv_sec )
				 {	diff_tv=start_tv_sec-tv_sec;
					is_first_tv=0;
				 }
				tv_sec += diff_tv;
			 }
			else if ( offset==2 )
				tv_usec=$2;
		 }
		else if ( type=="timing_relative" )
		 {	if ( offset==0 )
				tv_usec=$2;
		 }
		else if ( type=="spectrum" )
			data_arr[offset+1]=$2;

		offset++;
	 }
	else
	 {	printf(" \"unexpected_data\",\n");
	 }
	
  } \
 END \
  {	flush_data();
	if ( is_array )
	 {	printf(" null\n");
		printf("]\n");
	 }
  }'