boolean mica2_gps_parse_gga(char* str, gps_gga_t* gga)
{
//printf("%s\n", str);
memset(gga, 0, sizeof(gps_gga_t));
// need to pass this as the 'unit' for the horiz. dilution.
char dummy = '\0';
uint16_t value_count = 0;
char* values[20];
split_commas(str, values, &value_count);
// GGA message should have 14(?) fields
if (value_count != 14)
{
//printf("not enough fields in GGA data\n");
return false;
}
else
{
gga->utc = get_timestamp(values[0]);
gga->latitude = get_latlong(values[1], values[2], true);
gga->longitude = get_latlong(values[3], values[4], false);
gga->fix_quality = atoi(values[5]);
gga->satellite_count = atoi(values[6]);
gga->horiz_dilution = get_measurement(values[7], &dummy);
gga->altitude = get_measurement(values[8], values[9]);
gga->geoid_height = get_measurement(values[10], values[11]);
gga->since_last_update = atoi(values[12]);
gga->dgps_id = atoi(values[13]);
return true;
}
}
void main(int argc, char **argv)
{
FILE *fpmet, *fpwnd, *fpll, *fpout;
char met_dir[BUFSIZ+1], wnd_dir[BUFSIZ+1], out_dir[BUFSIZ+1];
char metname[BUFSIZ+1], wndname[BUFSIZ+1], outname[BUFSIZ+1];
char met[BUFSIZ+1], wnd[BUFSIZ+1], out[BUFSIZ+1], ll[BUFSIZ+1];
char fname[BUFSIZ+1];
char str[210];
int i,j,app_flag,tsteps,n=0,maxline=210;
short int prec, tmax, tmin, wind, dummy;
unsigned short int usprec;
float lat,lon;
if (argc!=6) { /* Must be exactly 5 arguments after filename */
printf("Incorrect number of commandline arguments \n");
printf("usage: combine_wind \"met_file_path\" \"wind_file_path\" \"latlong_file\" \"out_dir\" \"APPEND_FLAG\"\n");
exit(EXIT_FAILURE); }
/* read command line agruments */
strcpy(met_dir,argv[1]); printf("prcp,temp file dir: %s \n",met_dir);
strcpy(wnd_dir,argv[2]); printf("wind file dir: %s \n",wnd_dir);
strcpy(ll,argv[3]); printf("file_list: %s \n",ll);
strcpy(out_dir,argv[4]); printf("out_dir: %s \n",out_dir);
app_flag=atoi(argv[5]); printf("append flag = %d \n",app_flag);
/* open lat long list file */
if((fpll = fopen(ll,"r"))==NULL) {
printf("Cannot open file %s \n",ll);exit(0);}
/* count number of cells from lines in latlong file */
while (fgets(str,maxline,fpll) != '\0') n++;
printf("number of files = %d\n",n);
rewind(fpll);
/* Combine the data for each cell file */
for(i=0;i<n;i++){
/* retrieve lat and long from filename of current cell */
if(fgets(fname,maxline,fpll) == '\0')
printf("Latlon file error at %s",fname);
else printf("Cell %d of %d %s",i+1,n,fname);
get_latlong(fname, &lat, &lon );
/* make filenames */
sprintf(metname,"data_%.4f_%.4f",lat,lon);
sprintf(wndname,"wind_%.4f_%.4f",lat,lon);
sprintf(outname,"data_%.4f_%.4f",lat,lon);
strcpy(met,met_dir);
strcat(met,metname);
strcpy(wnd,wnd_dir);
strcat(wnd,wndname);
strcpy(out,out_dir);
strcat(out,outname);
/* open met, wind, output files */
if((fpmet = fopen(met,"rb"))==NULL) {
printf("Cannot open file %s \n",met);exit(0);}
if((fpwnd = fopen(wnd,"rb"))==NULL) {
printf("Cannot open file %s \n",wnd);exit(0);}
if(app_flag!=1){
if((fpout = fopen(out,"wb"))==NULL) {
printf("Cannot open file %s \n",out);exit(0);}
}
else {
if((fpout = fopen(out,"ab"))==NULL) {
printf("Cannot open file %s \n",out);exit(0);}
}
/* count number of timesteps in first wind file */
if(i==0) {
while(fread(&dummy,2,1,fpwnd) != 0 ) tsteps++;
rewind(fpwnd);
printf("Number of timeteps = %d\n",tsteps);}
/* read data from two files and combine into one output file */
for(j=0;j<tsteps;j++){
if (fread(&prec,2,1,fpmet) !=1) printf("Error reading prcp\n");
if (fread(&tmax,2,1,fpmet) !=1) printf("Error reading tmax\n");
if (fread(&tmin,2,1,fpmet) !=1) printf("Error reading tmin\n");
if (fread(&wind,2,1,fpwnd) !=1) printf("Error reading wind\n");
usprec=(unsigned short int) prec;
fwrite(&usprec,sizeof(unsigned short int),1,fpout);
fwrite(&tmax,sizeof(short int),1,fpout);
fwrite(&tmin,sizeof(short int),1,fpout);
fwrite(&wind,sizeof(short int),1,fpout);
}
fclose(fpmet); fclose(fpwnd); fclose(fpout);
} /* end loop for each file */
fclose(fpll);
} /* end main */
