incremental_updates.20220601.tar.gz

/****************************************************************************/ /* program name : RNX2CRX */ /* */ /* RINEX file compression program for UNIX/MS-DOS system. */ /* convert the RINEX format to Compact RINEX format */ /* Created by Yuki HATANAKA / Geospatial Information Authority of Japan */ /* */ /* ver. */ /* 4.0.0 2007-02-05 test version Y. Hatanaka */ /* - CRINEX 1/3 for RINEX 2.x/3.x */ /* 4.0.1 2007-05-08 */ /* - elimination of supports for VMS and SUN OS 4.1.x */ /* - output not to the current directory but the same */ /* directory as the input file. */ /* - the same code for DOS and UNIX */ /* 4.0.2 2007-06-07 */ /* - fixing incompatibility of argument and format */ /* string of printf. */ /* 4.0.3 2007-06-21 */ /* - nothing was changed from 4.0.2 except for the */ /* version string to be the same as crx2rnx.c */ /* 4.0.4 2009-06-31 Y. Hatanaka */ /* - rename function getline to ggetline */ /* - correct typos in error messages */ /* 4.0.5 2012-07-1x Y. Hatanaka */ /* - check length of input file name */ /* - minor changes to suppress warning messages */ /* at compilation. */ /* - initialize "dummy" in main routine */ /* 4.0.6 2014-03-24 Y. Hatanaka */ /* - "include" statements are moved to the top. */ /* - Manipulation of file names in the new file naming */ /* convention (*.rnx/crx) is added. */ /* 4.0.7 2016-04-14 Y. Hatanaka */ /* - increase the following constants */ /* MAXSAT 90 -> 100 */ /* MAXTYPE 50 -> 100 */ /* MAXCLM 1024 -> 2048 */ /* MAX_BUFF_SIZE 131072 -> 204800 */ /* 4.0.8 2019-07-12 Y. Hatanaka */ /* - Correction of a bug which eliminates clock offset */ /* data periodically when used with the option "-e". */ /* - New option "-d" is added to delete input file after */ /* successful conversion. */ /* 4.1.0 2021-12-16 Y. Hatanaka */ /* - RINEX ver.4.xx files are accepted and converted to */ /* Compact RINEX ver. 3 format */ /* - Following bugs are fixed: */ /* + Error in case the number of special records exceeds */ /* 99 in RINEX ver.2 files */ /* + Error in case clock offset is padded with spaces */ /* in RINEX ver. 3 or 4 files. */ /* + Error in case a bad GNSS type is detected even if */ /* option -s is specified. */ /* */ /* Copyright (c) 2007 Geospatial Information Authority of Japan */ /* */ /****************************************************************************/ #include <limits.h> #include <stdio.h> #include <stdlib.h> #include <string.h> #include <ctype.h> #include <time.h> #define VERSION "ver.4.1.0" /**** Exit codes are defined here. ****/ #define EXIT_WARNING 2 #ifndef EXIT_FAILURE #define EXIT_FAILURE 1 #endif #ifndef EXIT_SUCCESS #define EXIT_SUCCESS 0 #endif /*** define macro to flush or clear the output buffer ***/ #define FLUSH_BUFF printf("%s",top_buff), *(p_buff = top_buff)='\0' #define CLEAR_BUFF *(p_buff = top_buff) = '\0' #define CRX_VERSION1 "1.0" /* CRINEX version for RINEX 2.x */ #define CRX_VERSION2 "3.0" /* CRINEX version for RINEX 3.x */ #define PROGNAME "RNX2CRX" #define MAXSAT 100 /* Maximum number of satellites observed at one epoch */ #define MAXTYPE 100 /* Maximum number of data types for a GNSS system */ #define MAXCLM 2048 /* Maximum columns in one line (>MAXTYPE*19+3) */ #define MAX_BUFF_SIZE 204800 /* Maximum size of output buffer (>MAXSAT*(MAXTYPE*19+4)+60 */ #define ARC_ORDER 3 /* order of difference to take */ /* define data structure for fields of clock offset and observation records */ /* Those data will be handled as integers after eliminating decimal points. */ /* Since their size may exceeds range between LONG_MIN and LONG_MAX, */ /* they will be read with being divided properly into upper and lower digits */ typedef struct clock_format{ long u[ARC_ORDER+1]; /* upper X digits */ long l[ARC_ORDER+1]; /* lower 8 digits (can be 9-10 digits for deltas)*/ } clock_format; typedef struct data_format{ long u[ARC_ORDER+1]; /* upper X digits */ long l[ARC_ORDER+1]; /* lower 5 digits (can be 6-7 digits for deltas) */ int order; } data_format; /* define global variables */ long ep_count = 0; long ep_reset = 0; long nl_count = 0; int rinex_version; /* =2, 3 or 4 */ int nsat,ntype,ntype_gnss[UCHAR_MAX],ntype_record[MAXSAT],clk_order = -1; int exit_status = EXIT_SUCCESS; int skip_strange_epoch = 0; /* default : stop with error */ int delete_if_no_error = 0; /* default : not delete */ int n_infile = 0; /* number of input file (must be 0 or 1) */ /* clock_format clk1,clk0 = {0,0,0,0,0,0,0,0}; */ clock_format clk1,clk0 = {{0,0,0,0},{0,0,0,0}}; data_format dy0[MAXSAT][MAXTYPE], dy1[MAXSAT][MAXTYPE]; char flag0[MAXSAT][MAXTYPE*2], flag[MAXSAT][MAXTYPE*2]; char out_buff[MAX_BUFF_SIZE] = {'x','\0'}; /**** a character is put as a stopper to avoid memory overflow ****/ char *top_buff=&out_buff[1],*p_buff; /**** therefore, actual buffer start from the second character ****/ char infile[MAXCLM]; /**** name of input file ****/ size_t C1 = sizeof(""); /* size of one character */ size_t C2 = sizeof(" "); /* size of 2-character string */ size_t C3 = sizeof(" "); /* size of 3-character string */ size_t C14 = sizeof(" "); /* size of 14-character string */ char oldline[MAXCLM] = {'&','\0'}; int nsat_old = 0; /* declaration of functions */ void parse_args(int argc, char *argv[]); void header(void); int get_next_epoch(char *p_line); void skip_to_next(char *p_line); void initialize_all(char *oldline,int *nsat_old, int count); void put_event_data(char *p_line); void read_clock(char *line,int shift_cl); void process_clock(void); int set_sat_table(char *p_new, char *p_old, int nsat_old,int *sa