|
|
/* $Id: splitCbfFile.c 3304 2011-01-17 15:25:59Z brideout $ */
|
|
|
/*
|
|
|
modification history
|
|
|
--------------------
|
|
|
00a,02apr02 original
|
|
|
*/
|
|
|
|
|
|
/*
|
|
|
USAGE: splitCbfFile cbfFileName
|
|
|
e.g.: splitCbfFile sts660125a.cbf
|
|
|
|
|
|
splitCbfFile splits the files in a multi-file CBF file into individual
|
|
|
files. This needs to be done prior to entering CBF files from NCAR
|
|
|
into the Madrigal Database. This is a rare example of a production program
|
|
|
which uses the low-level cedarIO routines. The reason is that above
|
|
|
this level the cbf distiction between end-of-file and end-of-data is
|
|
|
lost.
|
|
|
|
|
|
*/
|
|
|
|
|
|
#include <stdlib.h>
|
|
|
#include <string.h>
|
|
|
#include <sys/types.h>
|
|
|
#include <sys/stat.h>
|
|
|
#include <fcntl.h>
|
|
|
#include <stdio.h>
|
|
|
#include <cedarIO.h>
|
|
|
#include <cedar.h>
|
|
|
|
|
|
int
|
|
|
main (argc, argv)
|
|
|
int argc;
|
|
|
char *argv[];
|
|
|
{
|
|
|
|
|
|
int ind1, ind2, nrec;
|
|
|
char infile[100];
|
|
|
FILE *foutp = NULL;
|
|
|
FILE *finp, *flp;
|
|
|
Int16 *cedarRecordp=(Int16 *)NULL;
|
|
|
int fileType;
|
|
|
|
|
|
/* madrigal state variables */
|
|
|
int madrigalBlockSize=6720;
|
|
|
|
|
|
/* Cbf state parameters */
|
|
|
int cbfBlockSize=4096;
|
|
|
int forceCosRead=0, initPos8, initPos=0, initLCosBlock=0,
|
|
|
lCosBlock=0, gpos=0, initFwi=0, fwi=0;
|
|
|
Int16 *cosRecordp=(Int16 *) NULL;
|
|
|
|
|
|
int n, ibyr, ibdt, ibymd, newfile, kinst, month, day,
|
|
|
fileCategory, action;
|
|
|
|
|
|
char outfile[128], tag[6], staName[4], expDir[128], expName[128],
|
|
|
y[3], m[3], d[3], fileList[100];
|
|
|
char months[12][4];
|
|
|
kinst = 0;
|
|
|
sprintf(months[0], "%s", "jan");
|
|
|
sprintf(months[1], "%s", "feb");
|
|
|
sprintf(months[2], "%s", "mar");
|
|
|
sprintf(months[3], "%s", "apr");
|
|
|
sprintf(months[4], "%s", "may");
|
|
|
sprintf(months[5], "%s", "jun");
|
|
|
sprintf(months[6], "%s", "jul");
|
|
|
sprintf(months[7], "%s", "aug");
|
|
|
sprintf(months[8], "%s", "sep");
|
|
|
sprintf(months[9], "%s", "oct");
|
|
|
sprintf(months[10], "%s", "nov");
|
|
|
sprintf(months[11], "%s", "dec");
|
|
|
|
|
|
if (argc == 3) {
|
|
|
strcpy(infile, argv[1]);
|
|
|
strcpy(fileList, argv[2]);
|
|
|
} else {
|
|
|
(void) printf("Usage: splitCbfFile inputFile outputFile\n");
|
|
|
return(1);
|
|
|
}
|
|
|
flp = fopen(fileList, "w");
|
|
|
|
|
|
/* Check file type of start file */
|
|
|
(void) printf("Input file = %s\n", infile);
|
|
|
fileType = cedarFileType(infile, madrigalBlockSize, cbfBlockSize);
|
|
|
(void) printf("%s - fileType = %d\n", infile,fileType);
|
|
|
|
|
|
|
|
|
/*** Cbf to Ascii ***/
|
|
|
finp = fopen(infile, "r");
|
|
|
n = 0;
|
|
|
lCosBlock = 0;
|
|
|
gpos = 0;
|
|
|
fwi = 0;
|
|
|
cosRecordp = (Int16 *) NULL;
|
|
|
nrec = 0;
|
|
|
newfile = 1;
|
|
|
while((ind1 = getNextCedarCbfRecord(finp, &cedarRecordp,
|
|
|
forceCosRead,
|
|
|
cbfBlockSize,
|
|
|
&initPos8,
|
|
|
&initFwi,
|
|
|
&initPos,
|
|
|
&initLCosBlock,
|
|
|
&lCosBlock,
|
|
|
&gpos,
|
|
|
&fwi,
|
|
|
&cosRecordp)) != -2) {
|
|
|
if (ind1 == 0) {
|
|
|
kinst = cedarGetKinst(cedarRecordp);
|
|
|
ibyr = cedarGetIbyr(cedarRecordp);
|
|
|
ibdt = cedarGetIbdt(cedarRecordp);
|
|
|
month = ibdt/100;
|
|
|
day = ibdt - 100*month;
|
|
|
sprintf(m, "%2.2d", month);
|
|
|
sprintf(d, "%2.2d", day);
|
|
|
sprintf(y, "%2.2d", ibyr - 100*(ibyr/100));
|
|
|
ibymd = 10000*(ibyr - 100*(ibyr/100)) + ibdt;
|
|
|
if (newfile == 1) {
|
|
|
n++;
|
|
|
sprintf(tag,"%6.6d", ibymd);
|
|
|
staName[0] = infile[0];
|
|
|
staName[1] = infile[1];
|
|
|
staName[2] = infile[2];
|
|
|
staName[3] = '\0';
|
|
|
strcpy(outfile, staName);
|
|
|
strcat(outfile, tag);
|
|
|
strcat(outfile, "g.001");
|
|
|
foutp = fopen(outfile, "w");
|
|
|
newfile = 0;
|
|
|
nrec = 0;
|
|
|
|
|
|
sprintf(expDir, "%4.4d", ibyr);
|
|
|
strcat(expDir, "/");
|
|
|
strcat(expDir, staName);
|
|
|
strcat(expDir, "/");
|
|
|
strcat(expDir, d);
|
|
|
strcat(expDir, months[month-1]);
|
|
|
strcat(expDir, y);
|
|
|
|
|
|
/* sts660125g.001 1966/sts/25jan66 41 "St Santin I.S. Radar data" 1 1 */
|
|
|
|
|
|
}
|
|
|
ind2 = putNextCedarAsciiRecord(foutp, &cedarRecordp);
|
|
|
nrec++;
|
|
|
} else if (ind1 == -1) {
|
|
|
/* printf("Wrote %d records to %s\n", nrec,outfile); */
|
|
|
strcpy(expName, "\"St Santin Data Archive\"");
|
|
|
fileCategory = 1;
|
|
|
action = 1;
|
|
|
printf("%s %s %d %s %d %d\n",
|
|
|
outfile, expDir, kinst, expName, fileCategory, action);
|
|
|
fprintf(flp, "%s %s %d %s %d %d\n",
|
|
|
outfile, expDir, kinst, expName, fileCategory, action);
|
|
|
ibyr = 0;
|
|
|
ibdt = 0;
|
|
|
ibymd = 0;
|
|
|
(void) fclose(foutp);
|
|
|
newfile = 1;
|
|
|
} else {
|
|
|
printf("Fatal error: ind1 = %d\n", ind1);
|
|
|
return(-1);
|
|
|
}
|
|
|
}
|
|
|
|
|
|
(void) fclose(finp);
|
|
|
(void) fclose(foutp);
|
|
|
|
|
|
return(0);
|
|
|
}
|
|
|
|
