Skip to content
GitLab
Commit 8b3df29b authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

merge declaration and definition

parent 9d2fb4f1
Hide whitespace changes
Inline Side-by-side
src/stream_cgribex.c
View file @ 8b3df29b
......@@ -578,18 +578,15 @@ cgribexScanTsConstAdjust(stream_t *streamptr, taxis_t *taxis)
int cgribexScanTimestep1(stream_t * streamptr)
{
int *isec0, *isec1, *isec2, *isec3, *isec4;
double fsec2[512], fsec3[2], *fsec4 = NULL;
int lmv = 0, iret = 0;
off_t recpos = 0;
void *gribbuffer = NULL;
size_t buffersize = 0;
int rstatus;
int fileID;
int param = 0;
int level1 = 0, level2 = 0, vdate = 0, vtime = 0;
DateTime datetime, datetime0 = { LONG_MIN, LONG_MIN };
int tsID;
size_t readsize;
unsigned nrecords, recID;
int nrecs_scanned = 0;
......@@ -599,7 +596,6 @@ int cgribexScanTimestep1(stream_t * streamptr)
int warn_numavg = TRUE;
int taxisID = -1;
int rdate = 0, rtime = 0, tunit = 0, fcast = 0;
taxis_t *taxis;
int vlistID;
int comptype;
long unzipsize;
......@@ -608,19 +604,19 @@ int cgribexScanTimestep1(stream_t * streamptr)
streamptr->curTsID = 0;
isec0 = streamptr->record->sec0;
isec1 = streamptr->record->sec1;
isec2 = streamptr->record->sec2;
isec3 = streamptr->record->sec3;
isec4 = streamptr->record->sec4;
int *isec0 = streamptr->record->sec0;
int *isec1 = streamptr->record->sec1;
int *isec2 = streamptr->record->sec2;
int *isec3 = streamptr->record->sec3;
int *isec4 = streamptr->record->sec4;
tsID = tstepsNewEntry(streamptr);
taxis = &streamptr->tsteps[tsID].taxis;
int tsID = tstepsNewEntry(streamptr);
taxis_t *taxis = &streamptr->tsteps[tsID].taxis;
if ( tsID != 0 )
Error("Internal problem! tstepsNewEntry returns %d", tsID);
fileID = streamptr->fileID;
int fileID = streamptr->fileID;
while ( nskip-- > 0 )
{
......@@ -700,9 +696,7 @@ int cgribexScanTimestep1(stream_t * streamptr)
{
datetime.date = vdate;
datetime.time = vtime;
compvar_t compVar = cgribexVarSet(param, level1, level2, ISEC1_LevelType, ISEC1_TimeRange);
for ( recID = 0; recID < nrecs; recID++ )
{
if ( cgribexVarCompare(compVar, streamptr->tsteps[0].records[recID], 0) == 0 ) break;
......@@ -805,56 +799,47 @@ int cgribexScanTimestep1(stream_t * streamptr)
int cgribexScanTimestep2(stream_t * streamptr)
{
int rstatus = 0;
int *isec0, *isec1, *isec2, *isec3, *isec4;
double fsec2[512], fsec3[2], *fsec4 = NULL;
int lmv = 0, iret = 0;
off_t recpos = 0;
void *gribbuffer = NULL;
size_t buffersize = 0;
int fileID;
int param = 0;
int level1 = 0, level2 = 0, vdate = 0, vtime = 0;
DateTime datetime, datetime0 = { LONG_MIN, LONG_MIN };
int tsID;
int varID, gridID;
size_t readsize;
int nrecords, nrecs, recID, rindex;
int nrecs_scanned = 0;
int nrecs, recID;
long recsize = 0;
int warn_numavg = TRUE;
int tsteptype;
int taxisID = -1;
taxis_t *taxis;
int vlistID;
long unzipsize;
char paramstr[32];
streamptr->curTsID = 1;
isec0 = streamptr->record->sec0;
isec1 = streamptr->record->sec1;
isec2 = streamptr->record->sec2;
isec3 = streamptr->record->sec3;
isec4 = streamptr->record->sec4;
int *isec0 = streamptr->record->sec0;
int *isec1 = streamptr->record->sec1;
int *isec2 = streamptr->record->sec2;
int *isec3 = streamptr->record->sec3;
int *isec4 = streamptr->record->sec4;
fileID = streamptr->fileID;
vlistID = streamptr->vlistID;
taxisID = vlistInqTaxis(vlistID);
int fileID = streamptr->fileID;
int vlistID = streamptr->vlistID;
int taxisID = vlistInqTaxis(vlistID);
gribbuffer = streamptr->record->buffer;
buffersize = streamptr->record->buffersize;
void *gribbuffer = streamptr->record->buffer;
size_t buffersize = streamptr->record->buffersize;
tsID = streamptr->rtsteps;
int tsID = streamptr->rtsteps;
if ( tsID != 1 )
Error("Internal problem! unexpected timestep %d", tsID+1);
taxis = &streamptr->tsteps[tsID].taxis;
taxis_t *taxis = &streamptr->tsteps[tsID].taxis;
fileSetPos(fileID, streamptr->tsteps[tsID].position, SEEK_SET);
cdi_create_records(streamptr, tsID);
nrecords = streamptr->tsteps[tsID].nallrecs;
int nrecords = streamptr->tsteps[tsID].nallrecs;
if ( nrecords ) streamptr->tsteps[1].recIDs = (int *) Malloc((size_t)nrecords * sizeof(int));
streamptr->tsteps[1].nrecs = 0;
for ( recID = 0; recID < nrecords; recID++ )
......@@ -867,8 +852,8 @@ int cgribexScanTimestep2(stream_t * streamptr)
streamptr->tsteps[tsID].records[recID].size = streamptr->tsteps[0].records[recID].size;
}
nrecs_scanned = nrecords;
rindex = 0;
int nrecs_scanned = nrecords;
int rindex = 0;
while ( TRUE )
{
if ( rindex > nrecords ) break;
......
src/stream_gribapi.c
View file @ 8b3df29b
......@@ -427,29 +427,23 @@ void gribapiAddRecord(stream_t * streamptr, int param, grib_handle *gh,
int leveltype1, int leveltype2, int lbounds, int level1, int level2, int level_sf, int level_unit,
const var_tile_t *tiles, int lread_additional_keys)
{
int zaxistype;
int gridID = CDI_UNDEFID, varID;
int varID;
int levelID = 0;
int tsID, recID;
int numavg;
int tsteptype;
record_t *record;
grid_t grid;
int vlistID;
long lpar;
int status;
char stdname[CDI_MAX_NAME], longname[CDI_MAX_NAME], units[CDI_MAX_NAME];
size_t vlen;
long ens_index = 0, ens_count = 0, ens_forecast_type = 0;
vlistID = streamptr->vlistID;
tsID = streamptr->curTsID;
recID = recordNewEntry(streamptr, tsID);
record = &streamptr->tsteps[tsID].records[recID];
int vlistID = streamptr->vlistID;
int tsID = streamptr->curTsID;
int recID = recordNewEntry(streamptr, tsID);
record_t *record = &streamptr->tsteps[tsID].records[recID];
tsteptype = gribapiGetTsteptype(gh);
int tsteptype = gribapiGetTsteptype(gh);
// numavg = ISEC1_AvgNum;
numavg = 0;
int numavg = 0;
// fprintf(stderr, "param %d %d %d %d\n", param, level1, level2, leveltype1);
......@@ -477,9 +471,9 @@ void gribapiAddRecord(stream_t * streamptr, int param, grib_handle *gh,
gribapiGetGrid(gh, &grid);
gridID = varDefGrid(vlistID, &grid, 0);
int gridID = varDefGrid(vlistID, &grid, 0);
zaxistype = gribapiGetZaxisType(gribEditionNumber(gh), leveltype1);
int zaxistype = gribapiGetZaxisType(gribEditionNumber(gh), leveltype1);
switch (zaxistype)
{
......@@ -1006,9 +1000,6 @@ int gribapiScanTimestep2(stream_t * streamptr)
{
int rstatus = 0;
off_t recpos = 0;
unsigned char *gribbuffer = NULL;
size_t buffersize = 0;
int fileID;
DateTime datetime0 = { LONG_MIN, LONG_MIN };
// int gridID;
int recID;
......@@ -1017,12 +1008,12 @@ int gribapiScanTimestep2(stream_t * streamptr)
streamptr->curTsID = 1;
fileID = streamptr->fileID;
int fileID = streamptr->fileID;
int vlistID = streamptr->vlistID;
int taxisID = vlistInqTaxis(vlistID);
gribbuffer = (unsigned char *) streamptr->record->buffer;
buffersize = streamptr->record->buffersize;
unsigned char *gribbuffer = (unsigned char *) streamptr->record->buffer;
size_t buffersize = streamptr->record->buffersize;
int tsID = streamptr->rtsteps;
if ( tsID != 1 )
......@@ -1471,8 +1462,7 @@ int gribapiDecode(unsigned char *gribbuffer, int gribsize, double *data, int gri
int status = 0;
long lpar;
long numberOfPoints;
size_t datasize, dummy, recsize;
grib_handle *gh = NULL;
size_t datasize;
UNUSED(vlistID);
UNUSED(varID);
......@@ -1488,8 +1478,8 @@ int gribapiDecode(unsigned char *gribbuffer, int gribsize, double *data, int gri
}
}
recsize = (size_t)gribsize;
gh = grib_handle_new_from_message(NULL, (void *) gribbuffer, recsize);
size_t recsize = (size_t)gribsize;
grib_handle *gh = grib_handle_new_from_message(NULL, (void *) gribbuffer, recsize);
GRIB_CHECK(my_grib_set_double(gh, "missingValue", missval), 0);
/* get the size of the values array*/
......@@ -1500,12 +1490,11 @@ int gribapiDecode(unsigned char *gribbuffer, int gribsize, double *data, int gri
if ( gridsize != (int) datasize )
Error("Internal problem: gridsize(%d) != datasize(%d)!", gridsize, datasize);
dummy = datasize;
size_t dummy = datasize;
GRIB_CHECK(grib_get_double_array(gh, "values", data, &dummy), 0);
int gridtype;
GRIB_CHECK(grib_get_long(gh, "gridDefinitionTemplateNumber", &lpar), 0);
gridtype = (int) lpar;
int gridtype = (int) lpar;
*nmiss = 0;
if ( gridtype < 50 || gridtype > 53 )
......@@ -1735,7 +1724,6 @@ int gribapiDefDateTimeRel(int editionNumber, grib_handle *gh, int rdate, int rti
int status = -1;
int year, month, day, hour, minute, second;
int julday1, secofday1, julday2, secofday2, days, secs;
int factor;
long startStep = 0, endStep;
cdiDecodeDate(rdate, &year, &month, &day);
......@@ -1750,7 +1738,7 @@ int gribapiDefDateTimeRel(int editionNumber, grib_handle *gh, int rdate, int rti
(void) julday_sub(julday1, secofday1, julday2, secofday2, &days, &secs);
factor = getTimeunitFactor(timeunit);
int factor = getTimeunitFactor(timeunit);
if ( !(int) fmod(days*86400.0 + secs, factor) )
{
......@@ -1828,13 +1816,12 @@ void gribapiDefTime(int editionNumber, int productDefinitionTemplate, int typeOf
static
void gribapiDefGrid(int editionNumber, grib_handle *gh, int gridID, int comptype, int lieee, int datatype, int nmiss, int gcinit)
{
int gridtype;
int status;
static short lwarn = TRUE;
UNUSED(nmiss);
gridtype = gridInqType(gridID);
int gridtype = gridInqType(gridID);
if ( editionNumber <= 1 )
if ( gridtype == GRID_GME || gridtype == GRID_UNSTRUCTURED )
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment

Imprint | Privacy Policy