#include #include #include #include #include #include #include "nrutil.c" #include "UVStruct.c" #include "UVFile.c" #include "UVCommon.c" #include "SERC.c" void main(void) { int i, j; int n_col, n_file, n_row, ok; int hh, mm, ss, dd, yy; int day, year; double fraction, integer; float *wl, dark, time; char st[128], c; char sch_file[20], in_file[20], out_file[20], wl_file[20]; struct tm date; struct scan_header SERC; struct SERC_data S; FILE *fp_sch, *fp_in, *fp_out, *fp_data; wl = vector(1, 18); /*printf("Scan Schedule File: "); // schedule of scans scanf("%s", &sch_file); */ strcpy(sch_file, "96S_Sch.SER"); /*printf("Wavelength File: "); //wavelength file scanf("%s", &wl_file);*/ strcpy(wl_file, "96WL01.SER"); fp_in = fopen(wl_file, "r"); do{ fscanf(fp_in, "%s", &st); } while(strcmp(st, "Number_of_rows") != 0); fscanf(fp_in, "%d", &n_row); do{ fscanf(fp_in, "%s", &st); } while(strcmp(st, "Data:") != 0); for(i = 1; i <= n_row; i++){ fscanf(fp_in, "%s%f", &st, wl + i); printf("%d\t%f\n", i, wl[i]); } fclose(fp_in); fp_sch = fopen(sch_file, "r"); fp_data = fopen("96FE.SER", "w"); // compilation of header data fprintf(fp_data, "Header:\n"); fprintf(fp_data, "Number_of_columns\t22\n"); fprintf(fp_data, "Labels:\n"); fprintf(fp_data, "File\tInstrument\tSource\tLocation\tEnvironment\t"); fprintf(fp_data, "Year\tDay\tDate\tOperation\tMeasurement_type\t"); fprintf(fp_data, "Measurement_purpose\tNumber_of_scans\tNumber_of_columns\t"); fprintf(fp_data, "Number_of_rows\tNumber_of_channels\tStart_time\tStop_time\tIntegration_time\t"); fprintf(fp_data, "Comments\tData_file\tWavelength_file\tAnalysis_file\n"); fprintf(fp_data, "\t\t\t\t\t\t\t\t\t\t\t\t\t\t\t"); fprintf(fp_data, "[UTC]\t[UTC]\t[s]\n"); fprintf(fp_data, "Data:\n"); fscanf(fp_sch, "%s %s %d", &st, &st, &n_col); // read schedule file header do{ fscanf(fp_sch, "%s", &st); } while(strcmp(st, "Data:") != 0); while(fscanf(fp_sch, "%s", &in_file) != EOF){ // read schedule file info. strcpy(SERC.data_file, in_file); fscanf(fp_sch, "%s", SERC.inst); fscanf(fp_sch, "%s", SERC.source); fscanf(fp_sch, "%s", SERC.locat); fscanf(fp_sch, "%s", SERC.envir); fscanf(fp_sch, "%s", SERC.operat); fscanf(fp_sch, "%s", SERC.meas_type); fscanf(fp_sch, "%s", SERC.meas_purp); fscanf(fp_sch, "%s", SERC.comments); fp_in = fopen(in_file, "r"); printf("%s\t", in_file); fp_out = fopen("temp.dat", "w"); strcpy(S.inst, "UI"); n_row = 1; while(read_SERC_line(fp_in, &S) != 0){ // read SERC line sscanf(S.time, "%2d%2d%2d", &hh, &mm, &ss); //calculate time time = (float)hh + (float)mm/60.0 + (float)ss/3600.0; fprintf(fp_out, "%.5f\t", time); dark = (S.data[10]+S.data[20])/2.0; //calculate dark for(i = 1; i <= 9; i++){ //subtract dark fprintf(fp_out, "%.6f\t", S.data[i] - dark); } for(i = 11; i <= 19; i++){ fprintf(fp_out, "%.6f\t", S.data[i] - dark); } fprintf(fp_out,"%.6f\t%.6f\t", S.data[10], S.data[20]); //dark signals fprintf(fp_out, "%.2f\t", S.temp); //temperature fprintf(fp_out, "%d\n", S.fsc); //full-scale-counts if(n_row == 1) SERC.start_time = time; SERC.stop_time = time; n_row++; } n_row--; fclose(fp_in); fclose(fp_out); sscanf(in_file, "I%2d%2d%2d%s", &yy, &mm, &dd, &st); SERC.year = 1900 + yy; date.tm_sec = 0; date.tm_min = 0; date.tm_hour = 0; date.tm_mday = dd; date.tm_mon = mm - 1; date.tm_year = yy; mktime(&date); SERC.day = date.tm_yday + 1; day_to_date(SERC.year, SERC.day, SERC.date); n_file = 1; do{ // check if file exists sprintf(out_file, "96%03dF%02d.SER", SERC.day, n_file); ok = file_exist(out_file); n_file++; } while(ok == 1); fp_in = fopen("temp.dat", "r"); fp_out = fopen(out_file, "w"); // store header info. printf("%s\n", out_file); fprintf(fp_data, "%s\t", out_file); fprintf(fp_out, "Header:\n"); fprintf(fp_out, "Instrument\t%s\n", SERC.inst); fprintf(fp_data, "%s\t", SERC.inst); fprintf(fp_out, "Source\t%s\n", SERC.source); fprintf(fp_data, "%s\t", SERC.source); fprintf(fp_out, "Location\t%s\n", SERC.locat); fprintf(fp_data, "%s\t", SERC.locat); fprintf(fp_out, "Environment\t%s\n", SERC.envir); fprintf(fp_data, "%s\t", SERC.envir); fprintf(fp_out, "Year\t%d\n", SERC.year); fprintf(fp_data, "%d\t", SERC.year); fprintf(fp_out, "Day\t%d\n", SERC.day); fprintf(fp_data, "%d\t", SERC.day); fprintf(fp_out, "Date\t%s\n", SERC.date); fprintf(fp_data, "%s\t", SERC.date); fprintf(fp_out, "Operation\t%s\n", SERC.operat); fprintf(fp_data, "%s\t", SERC.operat); fprintf(fp_out, "Measurement_type\t%s\n", SERC.meas_type); fprintf(fp_data, "%s\t", SERC.meas_type); fprintf(fp_out, "Measurement_purpose\t%s\n", SERC.meas_purp); fprintf(fp_data, "%s\t", SERC.meas_purp); fprintf(fp_out, "Number_of_scans\t1\n"); fprintf(fp_data, "1\t"); fprintf(fp_out, "Number_of_columns\t23\n"); fprintf(fp_data, "23\t"); fprintf(fp_out, "Number_of_rows\t%d\n", n_row); fprintf(fp_data, "%d\t", n_row); fprintf(fp_out, "Number_of_channels\t%d\n", 18); fprintf(fp_data, "%d\t", 18); fprintf(fp_out, "Start_time\t%.5f\t[UTC]\n", SERC.start_time); fprintf(fp_data, "%.5f\t", SERC.start_time); fprintf(fp_out, "Stop_time\t%.5f\t[UTC]\n", SERC.stop_time); fprintf(fp_data, "%.5f\t", SERC.stop_time); fprintf(fp_out, "Integration_time\t60\t[s]\n"); fprintf(fp_data, "60\t"); fprintf(fp_out, "Comments\t%s\n", SERC.comments); fprintf(fp_data, "%s\t", SERC.comments); fprintf(fp_out, "Data_file\t%s\n", SERC.data_file); fprintf(fp_data, "%s\t", SERC.data_file); fprintf(fp_out, "Wavelength_file\t%s\n", wl_file); fprintf(fp_data, "%s\t", wl_file); fprintf(fp_out, "Analysis_program\tFmt_96SERC_S.c\n"); fprintf(fp_data, "Fmt_96SERC_S.c\n"); fprintf(fp_out, "Labels:\n"); fprintf(fp_out, "Time\t"); for(i = 1; i <= 18; i++){ fprintf(fp_out, "Signal_%3.2f[nm]\t", wl[i]); } fprintf(fp_out, "Dark_J\t"); fprintf(fp_out, "Dark_T\t"); fprintf(fp_out, "Filter_temperature\tFull_scale_counts\n"); fprintf(fp_out, "[UTC]\t"); for(i = 1; i <= 20; i++){ fprintf(fp_out, "[V]\t"); } fprintf(fp_out, "[C]\n"); fprintf(fp_out, "Data:\n"); while((c = fgetc(fp_in)) != EOF){ fputc(c, fp_out); } fclose(fp_in); fclose(fp_out); } fclose(fp_sch); fclose(fp_data); }