Skip to content
GitLab
Commit 6d9469c6 authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

Merge declaration and definition.

parent b3aee032
Hide whitespace changes
Inline Side-by-side
src/stream_cgribex.c
View file @ 6d9469c6
......@@ -753,14 +753,13 @@ int cgribexScanTimestep1(stream_t *streamptr)
size_t buffersize = 0;
int leveltype = 0, level1 = 0, level2 = 0, vdate = 0, vtime = 0;
DateTime datetime, datetime0 = { LONG_MIN, LONG_MIN };
unsigned nrecords, recID;
unsigned recID;
int nrecs_scanned = 0;
bool warn_time = true;
bool warn_numavg = true;
int taxisID = -1;
int rdate = 0, rtime = 0, tunit = 0;
bool fcast = false;
int vlistID;
char paramstr[32];
streamptr->curTsID = 0;
......@@ -892,10 +891,10 @@ int cgribexScanTimestep1(stream_t *streamptr)
taxis->vdate = (int)datetime0.date;
taxis->vtime = (int)datetime0.time;
vlistID = streamptr->vlistID;
int vlistID = streamptr->vlistID;
vlistDefTaxis(vlistID, taxisID);
nrecords = (unsigned)streamptr->tsteps[0].nallrecs;
unsigned nrecords = (unsigned)streamptr->tsteps[0].nallrecs;
if ( nrecords < (unsigned)streamptr->tsteps[0].recordSize )
{
streamptr->tsteps[0].recordSize = (int)nrecords;
......@@ -970,7 +969,7 @@ int cgribexScanTimestep2(stream_t * streamptr)
if ( rindex > nrecords ) break;
size_t recsize = gribGetSize(fileID);
recpos = fileGetPos(fileID);
recpos = fileGetPos(fileID);
if ( recsize == 0 )
{
streamptr->ntsteps = 2;
......
src/stream_ext.c
View file @ 6d9469c6
......@@ -141,30 +141,30 @@ void extCopyRecord(stream_t *streamptr2, stream_t *streamptr1)
void extDefRecord(stream_t *streamptr)
{
Record *record = streamptr->record;
int pdis, pcat, pnum;
cdiDecodeParam(streamptr->record->param, &pnum, &pcat, &pdis);
cdiDecodeParam(record->param, &pnum, &pcat, &pdis);
int header[4];
header[0] = streamptr->record->date;
header[0] = record->date;
header[1] = pnum;
header[2] = streamptr->record->level;
int gridID = streamptr->record->gridID;
header[2] = record->level;
int gridID = record->gridID;
cdi_check_gridsize_int_limit("EXTRA", gridInqSize(gridID));
header[3] = (int) gridInqSize(gridID);
extrec_t *extp = (extrec_t*) streamptr->record->exsep;
extDefDatatype(streamptr->record->prec, &extp->prec, &extp->number);
extrec_t *extp = (extrec_t*) record->exsep;
extDefDatatype(record->prec, &extp->prec, &extp->number);
extDefHeader(extp, header);
}
void extWriteRecord(stream_t *streamptr, const double *data)
{
int fileID = streamptr->fileID;
void *extp = streamptr->record->exsep;
extDefDataDP(extp, data);
extWrite(fileID, extp);
extWrite(streamptr->fileID, extp);
}
static
......@@ -189,16 +189,15 @@ void extAddRecord(stream_t *streamptr, int param, int level, size_t xysize,
struct addIfNewRes gridAdded = cdiVlistAddGridIfNew(vlistID, grid, 0);
int gridID = gridAdded.Id;
if (!gridAdded.isNew) Free(grid);
/*
if ( level == 0 ) leveltype = ZAXIS_SURFACE;
else leveltype = ZAXIS_GENERIC;
*/
int leveltype = ZAXIS_GENERIC;
int datatype = extInqDatatype(prec, number);
int varID, levelID = 0;
varAddRecord(recID, param, gridID, leveltype, 0, level, 0, 0, 0,
extInqDatatype(prec, number), &varID, &levelID, TSTEP_INSTANT, 0, 0, -1,
datatype, &varID, &levelID, TSTEP_INSTANT, 0, 0, -1,
NULL, NULL, NULL, NULL, NULL, NULL, 0);
record->varID = (short)varID;
record->levelID = (short)levelID;
......@@ -214,9 +213,7 @@ void extScanTimestep1(stream_t *streamptr)
{
int header[4];
DateTime datetime0 = { LONG_MIN, LONG_MIN };
int varID;
off_t recpos = 0;
long recsize;
int recID;
extcompvar_t compVar, compVar0;
extrec_t *extp = (extrec_t*) streamptr->record->exsep;
......@@ -224,10 +221,9 @@ void extScanTimestep1(stream_t *streamptr)
streamptr->curTsID = 0;
int tsID = tstepsNewEntry(streamptr);
taxis_t *taxis = &streamptr->tsteps[tsID].taxis;
if ( tsID != 0 )
Error("Internal problem! tstepsNewEntry returns %d", tsID);
taxis_t *taxis = &streamptr->tsteps[tsID].taxis;
int fileID = streamptr->fileID;
......@@ -241,7 +237,7 @@ void extScanTimestep1(stream_t *streamptr)
streamptr->ntsteps = 1;
break;
}
recsize = fileGetPos(fileID) - recpos;
size_t recsize = (size_t)(fileGetPos(fileID) - recpos);
extInqHeader(extp, header);
......@@ -264,7 +260,7 @@ void extScanTimestep1(stream_t *streamptr)
compVar.level = rlevel;
for ( recID = 0; recID < nrecs; recID++ )
{
compVar0.param = streamptr->tsteps[0].records[recID].param;
compVar0.param = streamptr->tsteps[0].records[recID].param;
compVar0.level = streamptr->tsteps[0].records[recID].ilevel;
if ( memcmp(&compVar0, &compVar, sizeof(extcompvar_t)) == 0 ) break;
......@@ -280,7 +276,7 @@ void extScanTimestep1(stream_t *streamptr)
if ( CDI_Debug )
Message("%4d%8d%4d%8d%8d%6d", nrecs, (int)recpos, rcode, rlevel, vdate, vtime);
extAddRecord(streamptr, param, rlevel, rxysize, (size_t)recsize, recpos, extp->prec, extp->number);
extAddRecord(streamptr, param, rlevel, rxysize, recsize, recpos, extp->prec, extp->number);
}
streamptr->rtsteps = 1;
......@@ -309,12 +305,12 @@ void extScanTimestep1(stream_t *streamptr)
streamptr->tsteps[0].recIDs = (int *) Malloc((size_t)nrecords * sizeof (int));
streamptr->tsteps[0].nrecs = nrecords;
for ( recID = 0; recID < nrecords; recID++ )
for ( int recID = 0; recID < nrecords; recID++ )
streamptr->tsteps[0].recIDs[recID] = recID;
if ( streamptr->ntsteps == -1 )
{
tsID = tstepsNewEntry(streamptr);
int tsID = tstepsNewEntry(streamptr);
if ( tsID != streamptr->rtsteps )
Error("Internal error. tsID = %d", tsID);
......@@ -327,7 +323,7 @@ void extScanTimestep1(stream_t *streamptr)
if ( taxis->vdate == 0 && taxis->vtime == 0 )
{
streamptr->ntsteps = 0;
for ( varID = 0; varID < streamptr->nvars; varID++ )
for ( int varID = 0; varID < streamptr->nvars; varID++ )
vlistDefVarTimetype(vlistID, varID, TIME_CONSTANT);
}
}
......@@ -337,7 +333,6 @@ static
int extScanTimestep2(stream_t *streamptr)
{
int header[4];
int varID;
off_t recpos = 0;
extcompvar_t compVar, compVar0;
void *extp = streamptr->record->exsep;
......@@ -348,8 +343,7 @@ int extScanTimestep2(stream_t *streamptr)
int vlistID = streamptr->vlistID;
int tsID = streamptr->rtsteps;
if ( tsID != 1 )
Error("Internal problem! unexpected timestep %d", tsID+1);
if ( tsID != 1 ) Error("Internal problem! unexpected timestep %d", tsID+1);
taxis_t *taxis = &streamptr->tsteps[tsID].taxis;
......@@ -365,7 +359,6 @@ int extScanTimestep2(stream_t *streamptr)
for ( int recID = 0; recID < nrecords; recID++ )
{
varID = streamptr->tsteps[0].records[recID].varID;
streamptr->tsteps[tsID].records[recID].position =
streamptr->tsteps[0].records[recID].position;
streamptr->tsteps[tsID].records[recID].size =
......@@ -456,7 +449,7 @@ int extScanTimestep2(stream_t *streamptr)
{
if ( ! streamptr->tsteps[tsID].records[recID].used )
{
varID = streamptr->tsteps[tsID].records[recID].varID;
int varID = streamptr->tsteps[tsID].records[recID].varID;
vlistDefVarTimetype(vlistID, varID, TIME_CONSTANT);
}
else
......@@ -470,7 +463,7 @@ int extScanTimestep2(stream_t *streamptr)
if ( streamptr->ntsteps == -1 )
{
tsID = tstepsNewEntry(streamptr);
int tsID = tstepsNewEntry(streamptr);
if ( tsID != streamptr->rtsteps )
Error("Internal error. tsID = %d", tsID);
......@@ -488,11 +481,9 @@ int extInqContents(stream_t *streamptr)
extScanTimestep1(streamptr);
int status = 0;
if ( streamptr->ntsteps == -1 ) status = extScanTimestep2(streamptr);
int status = (streamptr->ntsteps == -1) ? extScanTimestep2(streamptr) : 0;
int fileID = streamptr->fileID;
fileSetPos(fileID, 0, SEEK_SET);
fileSetPos(streamptr->fileID, 0, SEEK_SET);
return status;
}
......@@ -506,15 +497,6 @@ long extScanTimestep(stream_t *streamptr)
int nrecs = 0;
extcompvar_t compVar, compVar0;
void *extp = streamptr->record->exsep;
/*
if ( CDI_Debug )
{
Message("streamID = %d", streamptr->self);
Message("cts = %d", streamptr->curTsID);
Message("rts = %d", streamptr->rtsteps);
Message("nts = %d", streamptr->ntsteps);
}
*/
int tsID = streamptr->rtsteps;
taxis_t *taxis = &streamptr->tsteps[tsID].taxis;
......
src/stream_ieg.c
View file @ 6d9469c6
......@@ -57,8 +57,7 @@ void iegReadRecord(stream_t *streamptr, double *data, size_t *nmiss)
void *iegp = streamptr->record->exsep;
int status = iegRead(fileID, iegp);
if ( status != 0 )
Error("Could not read IEG record!");
if ( status != 0 ) Error("Could not read IEG record!");
iegInqDataDP(iegp, data);
......@@ -448,7 +447,6 @@ void iegCopyRecord(stream_t *streamptr2, stream_t *streamptr1)
void iegDefRecord(stream_t *streamptr)
{
Record *record = streamptr->record;
iegrec_t *iegp = (iegrec_t*) record->exsep;
int vlistID = streamptr->vlistID;
int byteorder = streamptr->byteorder;
......@@ -460,6 +458,7 @@ void iegDefRecord(stream_t *streamptr)
int gridID = vlistInqVarGrid(vlistID, varID);
int zaxisID = vlistInqVarZaxis(vlistID, varID);
iegrec_t *iegp = (iegrec_t*) record->exsep;
iegInitMem(iegp);
for ( int i = 0; i < 37; i++ ) iegp->ipdb[i] = -1;
......@@ -616,8 +615,7 @@ void iegAddRecord(stream_t *streamptr, int param, int *pdb, int *gdb, double *vc
int datatype = iegInqDatatype(prec);
int varID;
int levelID = 0;
int varID, levelID = 0;
varAddRecord(recID, param, gridID, leveltype, lbounds, level1, level2, 0, 0,
datatype, &varID, &levelID, TSTEP_INSTANT, 0, 0, -1,
NULL, NULL, NULL, NULL, NULL, NULL, 0);
......@@ -687,10 +685,9 @@ void iegScanTimestep1(stream_t *streamptr)
streamptr->curTsID = 0;
int tsID = tstepsNewEntry(streamptr);
taxis_t *taxis = &streamptr->tsteps[tsID].taxis;
if ( tsID != 0 )
Error("Internal problem! tstepsNewEntry returns %d", tsID);
taxis_t *taxis = &streamptr->tsteps[tsID].taxis;
int fileID = streamptr->fileID;
......@@ -781,7 +778,7 @@ void iegScanTimestep1(stream_t *streamptr)
if ( streamptr->ntsteps == -1 )
{
tsID = tstepsNewEntry(streamptr);
int tsID = tstepsNewEntry(streamptr);
if ( tsID != streamptr->rtsteps )
Error("Internal error. tsID = %d", tsID);
......@@ -937,7 +934,7 @@ int iegScanTimestep2(stream_t *streamptr)
if ( streamptr->ntsteps == -1 )
{
tsID = tstepsNewEntry(streamptr);
int tsID = tstepsNewEntry(streamptr);
if ( tsID != streamptr->rtsteps )
Error("Internal error. tsID = %d", tsID);
......@@ -951,16 +948,13 @@ int iegScanTimestep2(stream_t *streamptr)
int iegInqContents(stream_t *streamptr)
{
int fileID = streamptr->fileID;
streamptr->curTsID = 0;
iegScanTimestep1(streamptr);
int status = 0;
if ( streamptr->ntsteps == -1 ) status = iegScanTimestep2(streamptr);
int status = (streamptr->ntsteps == -1) ? iegScanTimestep2(streamptr) : 0;
fileSetPos(fileID, 0, SEEK_SET);
fileSetPos(streamptr->fileID, 0, SEEK_SET);
return status;
}
......
src/stream_srv.c
View file @ 6d9469c6
......@@ -92,8 +92,7 @@ void srvReadRecord(stream_t *streamptr, double *data, size_t *nmiss)
void *srvp = streamptr->record->exsep;
int status = srvRead(fileID, srvp);
if ( status != 0 )
Error("Failed to read record from SRV file");
if ( status != 0 ) Error("Failed to read record from SRV file");
int header[8];
srvInqHeader(srvp, header);
......@@ -124,7 +123,6 @@ void srvCopyRecord(stream_t *streamptr2, stream_t *streamptr1)
void srvDefRecord(stream_t *streamptr)
{
Record *record = streamptr->record;
srvrec_t *srvp = (srvrec_t*) record->exsep;
int pdis, pcat, pnum;
cdiDecodeParam(record->param, &pnum, &pcat, &pdis);
......@@ -154,8 +152,8 @@ void srvDefRecord(stream_t *streamptr)
header[6] = 0;
header[7] = 0;
int datatype = record->prec;
srvp->dprec = srvDefDatatype(datatype);
srvrec_t *srvp = (srvrec_t*) record->exsep;
srvp->dprec = srvDefDatatype(record->prec);
srvDefHeader(srvp, header);
}
......@@ -163,11 +161,9 @@ void srvDefRecord(stream_t *streamptr)
void srvWriteRecord(stream_t *streamptr, const double *data)
{
int fileID = streamptr->fileID;
void *srvp = streamptr->record->exsep;
srvDefDataDP(srvp, data);
srvWrite(fileID, srvp);
srvWrite(streamptr->fileID, srvp);
}
static
......@@ -192,16 +188,11 @@ void srv_add_record(stream_t *streamptr, int param, int level, size_t xsize, siz
struct addIfNewRes gridAdded = cdiVlistAddGridIfNew(vlistID, grid, 0);
int gridID = gridAdded.Id;
if (!gridAdded.isNew) Free(grid);
/*
if ( level == 0 ) leveltype = ZAXIS_SURFACE;
else leveltype = ZAXIS_GENERIC;
*/
int leveltype = ZAXIS_GENERIC;
int leveltype = ZAXIS_GENERIC;
int datatype = srvInqDatatype(prec);
int levelID = 0;
int varID;
int varID, levelID = 0;
varAddRecord(recID, param, gridID, leveltype, 0, level, 0, 0, 0,
datatype, &varID, &levelID, TSTEP_INSTANT, 0, 0, -1,
NULL, NULL, NULL, NULL, NULL, NULL, 0);
......@@ -222,17 +213,14 @@ void srvScanTimestep1(stream_t *streamptr)
{
DateTime datetime0 = { LONG_MIN, LONG_MIN };
off_t recpos;
taxis_t *taxis;
srvrec_t *srvp = (srvrec_t*) streamptr->record->exsep;
streamptr->curTsID = 0;
{
int tsID = tstepsNewEntry(streamptr);
if ( tsID != 0 )
Error("Internal problem! tstepsNewEntry returns %d", tsID);
taxis = &streamptr->tsteps[tsID].taxis;
}
int tsID = tstepsNewEntry(streamptr);
if ( tsID != 0 )
Error("Internal problem! tstepsNewEntry returns %d", tsID);
taxis_t *taxis = &streamptr->tsteps[tsID].taxis;
int fileID = streamptr->fileID;
......@@ -407,7 +395,7 @@ int srvScanTimestep2(stream_t *streamptr)
int recID;
for ( recID = 0; recID < nrecords; recID++ )
{
compVar0.param = streamptr->tsteps[tsID].records[recID].param;
compVar0.param = streamptr->tsteps[tsID].records[recID].param;
compVar0.level = streamptr->tsteps[tsID].records[recID].ilevel;
if ( memcmp(&compVar0, &compVar, sizeof(srvcompvar_t)) == 0 )
......@@ -437,7 +425,7 @@ int srvScanTimestep2(stream_t *streamptr)
streamptr->tsteps[tsID].records[recID].size = recsize;
compVar0.param = streamptr->tsteps[tsID].records[recID].param;
compVar0.param = streamptr->tsteps[tsID].records[recID].param;
compVar0.level = streamptr->tsteps[tsID].records[recID].ilevel;
if ( memcmp(&compVar0, &compVar, sizeof(srvcompvar_t)) != 0 )
......@@ -471,7 +459,7 @@ int srvScanTimestep2(stream_t *streamptr)
if ( streamptr->ntsteps == -1 )
{
tsID = tstepsNewEntry(streamptr);
int tsID = tstepsNewEntry(streamptr);
if ( tsID != streamptr->rtsteps )
Error("Internal error. tsID = %d", tsID);
......@@ -489,11 +477,9 @@ int srvInqContents(stream_t *streamptr)
srvScanTimestep1(streamptr);
int status = 0;
if ( streamptr->ntsteps == -1 ) status = srvScanTimestep2(streamptr);
int status = (streamptr->ntsteps == -1) ? srvScanTimestep2(streamptr) : 0;
int fileID = streamptr->fileID;
fileSetPos(fileID, 0, SEEK_SET);
fileSetPos(streamptr->fileID, 0, SEEK_SET);
return status;
}
......@@ -502,20 +488,10 @@ static
long srvScanTimestep(stream_t *streamptr)
{
int header[8];
/* int rxsize = 0, rysize = 0; */
off_t recpos = 0;
int recID;
int nrecs = 0;
void *srvp = streamptr->record->exsep;
/*
if ( CDI_Debug )
{
Message("streamID = %d", streamptr->self);
Message("cts = %d", streamptr->curTsID);
Message("rts = %d", streamptr->rtsteps);
Message("nts = %d", streamptr->ntsteps);
}
*/
int tsID = streamptr->rtsteps;
taxis_t *taxis = &streamptr->tsteps[tsID].taxis;
......
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