Commit 3183a467 authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

Renamed streamXXX() to pstreamXXX().

parent eae92c35
......@@ -152,8 +152,8 @@ void write_data(int streamID, int nvars, double *data[], int *nmiss)
{
for ( int varID = 0; varID < nvars; ++varID )
{
streamDefRecord(streamID, varID, 0);
streamWriteRecord(streamID, data[varID], nmiss[varID]);
pstreamDefRecord(streamID, varID, 0);
pstreamWriteRecord(streamID, data[varID], nmiss[varID]);
}
}
......@@ -190,7 +190,7 @@ void *Importamsr(void *argument)
int vdate = getDate(cdoStreamName(0)->args);
if ( vdate <= 999999 ) vdate = vdate*100 + 1;
int streamID = streamOpenWrite(cdoStreamName(1), cdoFiletype());
int streamID = pstreamOpenWrite(cdoStreamName(1), cdoFiletype());
/*
Longitude is 0.25*xdim-0.125 degrees east
......@@ -219,7 +219,7 @@ void *Importamsr(void *argument)
init_amsr_day(vlistID, gridID, zaxisID, nvars);
streamDefVlist(streamID, vlistID);
pstreamDefVlist(streamID, vlistID);
vtime = 13000; /* 1:30:00 */
for ( tsID = 0; tsID < 2; ++tsID )
......@@ -227,7 +227,7 @@ void *Importamsr(void *argument)
taxisDefVdate(taxisID, vdate);
taxisDefVtime(taxisID, vtime);
vtime += 120000; /* 13:30:00 */
streamDefTimestep(streamID, tsID);
pstreamDefTimestep(streamID, tsID);
processDefTimesteps(streamID);
read_amsr(fp, vlistID, nvars, data, nmiss);
......@@ -245,12 +245,12 @@ void *Importamsr(void *argument)
init_amsr_averaged(vlistID, gridID, zaxisID, nvars);
/* vlistDefNtsteps(vlistID, 0);*/
streamDefVlist(streamID, vlistID);
pstreamDefVlist(streamID, vlistID);
taxisDefVdate(taxisID, vdate);
taxisDefVtime(taxisID, vtime);
tsID = 0;
streamDefTimestep(streamID, tsID);
pstreamDefTimestep(streamID, tsID);
processDefTimesteps(streamID);
read_amsr(fp, vlistID, nvars, data, nmiss);
......@@ -264,7 +264,7 @@ void *Importamsr(void *argument)
processDefVarNum(vlistNvars(vlistID), streamID);
streamClose(streamID);
pstreamClose(streamID);
fclose(fp);
......
......@@ -338,9 +338,9 @@ void *Importbinary(void *argument)
vlistDefTaxis(vlistID, taxisID);
int streamID = streamOpenWrite(cdoStreamName(1), cdoFiletype());
int streamID = pstreamOpenWrite(cdoStreamName(1), cdoFiletype());
streamDefVlist(streamID, vlistID);
pstreamDefVlist(streamID, vlistID);
gridsize = pfi.dnum[0]*pfi.dnum[1];
......@@ -459,7 +459,7 @@ void *Importbinary(void *argument)
taxisDefVdate(taxisID, vdate);
taxisDefVtime(taxisID, vtime);
streamDefTimestep(streamID, tsID);
pstreamDefTimestep(streamID, tsID);
for ( int recID = 0; recID < nrecs; ++recID )
{
......@@ -544,8 +544,8 @@ void *Importbinary(void *argument)
*/
varID = recVarID[recID];
levelID = recLevelID[recID];
streamDefRecord(streamID, varID, levelID);
streamWriteRecord(streamID, array, nmiss);
pstreamDefRecord(streamID, varID, levelID);
pstreamWriteRecord(streamID, array, nmiss);
}
}
......@@ -557,7 +557,7 @@ void *Importbinary(void *argument)
processDefVarNum(vlistNvars(vlistID), streamID);
streamClose(streamID);
pstreamClose(streamID);
vlistDestroy(vlistID);
gridDestroy(gridID);
......
......@@ -1549,9 +1549,9 @@ void *Importcmsaf(void *argument)
vdate = get_vdate(vlistID);
if ( vdate == 0 ) vdate = 10101;
streamID = streamOpenWrite(cdoStreamName(1), cdoFiletype());
streamID = pstreamOpenWrite(cdoStreamName(1), cdoFiletype());
streamDefVlist(streamID, vlistID);
pstreamDefVlist(streamID, vlistID);
for ( tsID = 0; tsID < nt; ++tsID )
{
......@@ -1559,7 +1559,7 @@ void *Importcmsaf(void *argument)
vtime = 0;
if ( vtimes ) vtime = vtimes[tsID];
taxisDefVtime(taxisID, vtime);
streamDefTimestep(streamID, tsID);
pstreamDefTimestep(streamID, tsID);
for ( ivar = 0; ivar < dsets.nsets; ++ivar )
{
......@@ -1600,8 +1600,8 @@ void *Importcmsaf(void *argument)
cdoWarning(" Missval is inside of valid values! Name: %s Range: %g - %g Missval: %g\n",
dsets.obj[ivar].name, minval, maxval, missval);
*/
streamDefRecord(streamID, varID, levelID);
streamWriteRecord(streamID, array, nmiss);
pstreamDefRecord(streamID, varID, levelID);
pstreamWriteRecord(streamID, array, nmiss);
}
}
}
......@@ -1612,7 +1612,7 @@ void *Importcmsaf(void *argument)
processDefVarNum(vlistNvars(vlistID), streamID);
streamClose(streamID);
pstreamClose(streamID);
vlistDestroy(vlistID);
gridDestroy(gridID);
......
......@@ -64,13 +64,13 @@ void write_data(int streamID, int vlistID, int nvars, double *data[])
int gridsize = gridInqSize(vlistInqVarGrid(vlistID, varID));
double missval = vlistInqVarMissval(vlistID, varID);
streamDefRecord(streamID, varID, 0);
pstreamDefRecord(streamID, varID, 0);
int nmiss = 0;
for ( int i = 0; i < gridsize; ++i )
if ( DBL_IS_EQUAL(data[varID][i], missval) ) nmiss++;
streamWriteRecord(streamID, data[varID], nmiss);
pstreamWriteRecord(streamID, data[varID], nmiss);
}
}
......@@ -142,7 +142,7 @@ void *Importobs(void *argument)
int vdate = getDate(cdoStreamName(0)->args);
if ( vdate <= 999999 ) vdate = vdate*100 + 1;
int streamID = streamOpenWrite(cdoStreamName(1), cdoFiletype());
int streamID = pstreamOpenWrite(cdoStreamName(1), cdoFiletype());
int zaxisID = zaxisCreate(ZAXIS_SURFACE, 1);
......@@ -155,7 +155,7 @@ void *Importobs(void *argument)
init_vars(vlistID, gridID, zaxisID, nvars);
streamDefVlist(streamID, vlistID);
pstreamDefVlist(streamID, vlistID);
int vdate0 = 0;
int vtime0 = 0;
......@@ -179,7 +179,7 @@ void *Importobs(void *argument)
*/
taxisDefVdate(taxisID, vdate);
taxisDefVtime(taxisID, vtime);
streamDefTimestep(streamID, tsID);
pstreamDefTimestep(streamID, tsID);
init_data(vlistID, nvars, data);
......@@ -233,7 +233,7 @@ void *Importobs(void *argument)
processDefVarNum(vlistNvars(vlistID), streamID);
streamClose(streamID);
pstreamClose(streamID);
fclose(fp);
......
......@@ -292,9 +292,9 @@ void *Info(void *argument)
for ( int indf = 0; indf < cdoStreamCnt(); indf++ )
{
int streamID = streamOpenRead(cdoStreamName(indf));
int streamID = pstreamOpenRead(cdoStreamName(indf));
int vlistID = streamInqVlist(streamID);
int vlistID = pstreamInqVlist(streamID);
int taxisID = vlistInqTaxis(vlistID);
int nvars = vlistNvars(vlistID);
......@@ -309,7 +309,7 @@ void *Info(void *argument)
int indg = 0;
int tsID = 0;
while ( (nrecs = streamInqTimestep(streamID, tsID)) )
while ( (nrecs = pstreamInqTimestep(streamID, tsID)) )
{
dtlist_taxisInqTimestep(dtlist, taxisID, 0);
int vdate = dtlist_get_vdate(dtlist, 0);
......@@ -337,8 +337,8 @@ void *Info(void *argument)
fprintf(stdout, "\n" );
}
streamInqRecord(streamID, &varID, &levelID);
streamReadRecord(streamID, array, &nmiss);
pstreamInqRecord(streamID, &varID, &levelID);
pstreamReadRecord(streamID, array, &nmiss);
infostat_type *infostatp = &infostat[varID];
indg = (operatorID == XINFON) ? varID+1 : indg + 1;
......@@ -525,7 +525,7 @@ void *Info(void *argument)
tsID++;
}
streamClose(streamID);
pstreamClose(streamID);
if ( array ) Free(array);
if ( infostat ) Free(infostat);
......
......@@ -227,27 +227,27 @@ void *Input(void *argument)
taxisID = taxisCreate(TAXIS_RELATIVE);
vlistDefTaxis(vlistID, taxisID);
streamID = streamOpenWrite(cdoStreamName(0), output_filetype);
streamID = pstreamOpenWrite(cdoStreamName(0), output_filetype);
streamDefVlist(streamID, vlistID);
pstreamDefVlist(streamID, vlistID);
}
int vdate = date;
int vtime = time;
taxisDefVdate(taxisID, vdate);
taxisDefVtime(taxisID, vtime);
streamDefTimestep(streamID, tsID);
pstreamDefTimestep(streamID, tsID);
for ( int levelID = 0; levelID < nlevs; levelID++ )
{
streamDefRecord(streamID, varID, levelID);
pstreamDefRecord(streamID, varID, levelID);
int offset = gridsize*levelID;
int nmiss = 0;
for ( i = 0; i < gridsize; ++i )
if ( DBL_IS_EQUAL(array[offset+i], missval) ) nmiss++;
streamWriteRecord(streamID, array+offset, nmiss);
pstreamWriteRecord(streamID, array+offset, nmiss);
}
nrecs++;
......@@ -256,7 +256,7 @@ void *Input(void *argument)
if ( streamID >= 0 )
{
streamClose(streamID);
pstreamClose(streamID);
vlistDestroy(vlistID);
}
......
......@@ -302,7 +302,7 @@ void *Intgrid(void *argument)
int operatorID = cdoOperatorID();
int streamID1 = streamOpenRead(cdoStreamName(0));
int streamID1 = pstreamOpenRead(cdoStreamName(0));
if ( operatorID == INTGRIDBIL || operatorID == INTERPOLATE )
{
......@@ -329,7 +329,7 @@ void *Intgrid(void *argument)
yinc = parameter2int(operatorArgv()[1]);
}
int vlistID1 = streamInqVlist(streamID1);
int vlistID1 = pstreamInqVlist(streamID1);
int vlistID2 = vlistDuplicate(vlistID1);
int taxisID1 = vlistInqTaxis(vlistID1);
......@@ -369,31 +369,30 @@ void *Intgrid(void *argument)
vlistChangeGridIndex(vlistID2, index, gridID2);
}
int streamID2 = streamOpenWrite(cdoStreamName(1), cdoFiletype());
streamDefVlist(streamID2, vlistID2);
int streamID2 = pstreamOpenWrite(cdoStreamName(1), cdoFiletype());
pstreamDefVlist(streamID2, vlistID2);
int gridsize = vlistGridsizeMax(vlistID1);
double *array1 = (double*) Malloc(gridsize*sizeof(double));
double *array1 = (double*) Malloc(gridsize*sizeof(double));
gridsize = gridInqSize(gridID2);
double *array2 = (double*) Malloc(gridsize*sizeof(double));
double *array2 = (double*) Malloc(gridsize*sizeof(double));
field_type field1, field2;
field_init(&field1);
field_init(&field2);
int tsID = 0;
while ( (nrecs = streamInqTimestep(streamID1, tsID)) )
while ( (nrecs = pstreamInqTimestep(streamID1, tsID)) )
{
taxisCopyTimestep(taxisID2, taxisID1);
streamDefTimestep(streamID2, tsID);
pstreamDefTimestep(streamID2, tsID);
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID1, &varID, &levelID);
streamReadRecord(streamID1, array1, &nmiss);
pstreamInqRecord(streamID1, &varID, &levelID);
pstreamReadRecord(streamID1, array1, &nmiss);
gridID1 = vlistInqVarGrid(vlistID1, varID);
missval = vlistInqVarMissval(vlistID1, varID);
......@@ -417,14 +416,14 @@ void *Intgrid(void *argument)
nmiss = field2.nmiss;
streamDefRecord(streamID2, varID, levelID);
streamWriteRecord(streamID2, array2, nmiss);
pstreamDefRecord(streamID2, varID, levelID);
pstreamWriteRecord(streamID2, array2, nmiss);
}
tsID++;
}
streamClose(streamID2);
streamClose(streamID1);
pstreamClose(streamID2);
pstreamClose(streamID1);
if ( array2 ) Free(array2);
if ( array1 ) Free(array1);
......
......@@ -77,9 +77,9 @@ void *Intgridtraj(void *argument)
juldate_t juldate;
readnextpos(fp, calendar, &juldate, &xpos, &ypos);
int streamID1 = streamOpenRead(cdoStreamName(0));
int streamID1 = pstreamOpenRead(cdoStreamName(0));
int vlistID1 = streamInqVlist(streamID1);
int vlistID1 = pstreamInqVlist(streamID1);
field_type field1, field2;
field_init(&field1);
......@@ -129,33 +129,32 @@ void *Intgridtraj(void *argument)
int taxisID2 = taxisDuplicate(taxisID1);
vlistDefTaxis(vlistID2, taxisID2);
int streamID2 = streamOpenWrite(cdoStreamName(1), cdoFiletype());
streamDefVlist(streamID2, vlistID2);
int streamID2 = pstreamOpenWrite(cdoStreamName(1), cdoFiletype());
pstreamDefVlist(streamID2, vlistID2);
int tsID = 0;
int nrecs = streamInqTimestep(streamID1, tsID++);
int nrecs = pstreamInqTimestep(streamID1, tsID++);
juldate_t juldate1 = juldate_encode(calendar, taxisInqVdate(taxisID1), taxisInqVtime(taxisID1));
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID1, &varID, &levelID);
pstreamInqRecord(streamID1, &varID, &levelID);
gridsize = gridInqSize(vlistInqVarGrid(vlistID1, varID));
offset = gridsize*levelID;
single1 = vardata1[varID] + offset;
streamReadRecord(streamID1, single1, &nmiss);
pstreamReadRecord(streamID1, single1, &nmiss);
if ( nmiss ) cdoAbort("Missing values unsupported for this operator!");
}
int tsIDo = 0;
while ( juldate_to_seconds(juldate1) <= juldate_to_seconds(juldate) )
{
nrecs = streamInqTimestep(streamID1, tsID++);
nrecs = pstreamInqTimestep(streamID1, tsID++);
if ( nrecs == 0 ) break;
juldate_t juldate2 = juldate_encode(calendar, taxisInqVdate(taxisID1), taxisInqVtime(taxisID1));
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID1, &varID, &levelID);
pstreamInqRecord(streamID1, &varID, &levelID);
recVarID[recID] = varID;
recLevelID[recID] = levelID;
......@@ -163,7 +162,7 @@ void *Intgridtraj(void *argument)
gridsize = gridInqSize(vlistInqVarGrid(vlistID1, varID));
offset = gridsize*levelID;
single2 = vardata2[varID] + offset;
streamReadRecord(streamID1, single2, &nmiss);
pstreamReadRecord(streamID1, single2, &nmiss);
if ( nmiss ) cdoAbort("Missing values unsupported for this operator!");
}
......@@ -175,7 +174,7 @@ void *Intgridtraj(void *argument)
juldate_decode(calendar, juldate, &vdate, &vtime);
taxisDefVdate(taxisID2, vdate);
taxisDefVtime(taxisID2, vtime);
streamDefTimestep(streamID2, tsIDo++);
pstreamDefTimestep(streamID2, tsIDo++);
double fac1 = juldate_to_seconds(juldate_sub(juldate2, juldate)) /
juldate_to_seconds(juldate_sub(juldate2, juldate1));
......@@ -210,8 +209,8 @@ void *Intgridtraj(void *argument)
intgridbil(&field1, &field2);
streamDefRecord(streamID2, varID, levelID);
streamWriteRecord(streamID2, &point, nmiss);
pstreamDefRecord(streamID2, varID, levelID);
pstreamWriteRecord(streamID2, &point, nmiss);
}
}
if ( readnextpos(fp, calendar, &juldate, &xpos, &ypos) == EOF ) break;
......@@ -229,8 +228,8 @@ void *Intgridtraj(void *argument)
}
fclose(fp);
streamClose(streamID2);
streamClose(streamID1);
pstreamClose(streamID2);
pstreamClose(streamID1);
for ( varID = 0; varID < nvars; varID++ )
{
......
......@@ -186,9 +186,9 @@ void *Intlevel(void *argument)
if ( cdoVerbose ) for ( i = 0; i < nlev2; ++i ) printf("lev2 %d: %g\n", i, lev2[i]);
int streamID1 = streamOpenRead(cdoStreamName(0));
int streamID1 = pstreamOpenRead(cdoStreamName(0));
int vlistID1 = streamInqVlist(streamID1);
int vlistID1 = pstreamInqVlist(streamID1);
int vlistID2 = vlistDuplicate(vlistID1);
int taxisID1 = vlistInqTaxis(vlistID1);
......@@ -281,9 +281,8 @@ void *Intlevel(void *argument)
if ( zaxisID1 == vlistZaxis(vlistID1, i) )
vlistChangeZaxisIndex(vlistID2, i, zaxisID2);
int streamID2 = streamOpenWrite(cdoStreamName(1), cdoFiletype());
streamDefVlist(streamID2, vlistID2);
int streamID2 = pstreamOpenWrite(cdoStreamName(1), cdoFiletype());
pstreamDefVlist(streamID2, vlistID2);
int nvars = vlistNvars(vlistID1);
......@@ -321,22 +320,22 @@ void *Intlevel(void *argument)
int tsID = 0;
while ( (nrecs = streamInqTimestep(streamID1, tsID)) )
while ( (nrecs = pstreamInqTimestep(streamID1, tsID)) )
{
for ( varID = 0; varID < nvars; ++varID ) vars[varID] = false;
taxisCopyTimestep(taxisID2, taxisID1);
streamDefTimestep(streamID2, tsID);
pstreamDefTimestep(streamID2, tsID);
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID1, &varID, &levelID);
pstreamInqRecord(streamID1, &varID, &levelID);
gridsize = gridInqSize(vlistInqVarGrid(vlistID1, varID));
offset = gridsize*levelID;
single1 = vardata1[varID] + offset;
streamReadRecord(streamID1, single1, &varnmiss[varID][levelID]);
pstreamReadRecord(streamID1, single1, &varnmiss[varID][levelID]);
vars[varID] = true;
}
......@@ -374,8 +373,8 @@ void *Intlevel(void *argument)
gridsize = gridInqSize(vlistInqVarGrid(vlistID2, varID));
offset = gridsize*levelID;
single2 = vardata2[varID] + offset;
streamDefRecord(streamID2, varID, levelID);
streamWriteRecord(streamID2, single2, varnmiss[varID][levelID]);
pstreamDefRecord(streamID2, varID, levelID);
pstreamWriteRecord(streamID2, single2, varnmiss[varID][levelID]);
}
}
}
......@@ -383,8 +382,8 @@ void *Intlevel(void *argument)
tsID++;
}
streamClose(streamID2);
streamClose(streamID1);
pstreamClose(streamID2);
pstreamClose(streamID1);
for ( varID = 0; varID < nvars; varID++ )
{
......
......@@ -69,20 +69,20 @@ void *Intlevel3d(void *argument)
operatorInputArg("icoordinate");
int streamID1 = streamOpenRead(cdoStreamName(0)); /* input data */
int streamID2 = streamOpenRead(cdoStreamName(1)); /* 3d target vertical coordinate */
int streamID3 = streamOpenWrite(cdoStreamName(2),cdoFiletype()); /* output stream */
int streamID1 = pstreamOpenRead(cdoStreamName(0)); /* input data */
int streamID2 = pstreamOpenRead(cdoStreamName(1)); /* 3d target vertical coordinate */
int streamID3 = pstreamOpenWrite(cdoStreamName(2),cdoFiletype()); /* output stream */
/* Read filename from Parameter */
operatorInputArg("filename for vertical source coordinates variable");
operatorCheckArgc(1);
argument_t *fileargument = file_argument_new(operatorArgv()[0]);
int streamID0 = streamOpenRead(fileargument); /* 3d vertical input coordinate */
int streamID0 = pstreamOpenRead(fileargument); /* 3d vertical input coordinate */
file_argument_free(fileargument);
int vlistID0 = streamInqVlist(streamID0);
int vlistID1 = streamInqVlist(streamID1); taxisID1 = vlistInqTaxis(vlistID1);
int vlistID2 = streamInqVlist(streamID2);
int vlistID0 = pstreamInqVlist(streamID0);
int vlistID1 = pstreamInqVlist(streamID1); taxisID1 = vlistInqTaxis(vlistID1);
int vlistID2 = pstreamInqVlist(streamID2);
int vlistID3 = vlistDuplicate(vlistID1); taxisID3 = taxisDuplicate(taxisID1);
vlistDefTaxis(vlistID3, taxisID1);
......@@ -110,16 +110,16 @@ void *Intlevel3d(void *argument)
zlevels_in = (double*) Malloc(gridsize*(nlevel+2)*sizeof(double));
nlevi = nlevel; /* number of input levels for later use */
gridsizei = gridsize; /* horizontal gridsize of input z coordinate */
nrecs = streamInqTimestep(streamID0, 0);
nrecs = pstreamInqTimestep(streamID0, 0);
if (cdoVerbose) cdoPrint("%d records input 3d vertical height",nrecs);
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID0, &varID, &levelID);
pstreamInqRecord(streamID0, &varID, &levelID);
gridsize = gridInqSize(vlistInqVarGrid(vlistID0, varID));
offset = gridsize + gridsize*levelID;
single1 = zlevels_in + offset;
streamReadRecord(streamID0, single1, &zlevels_in_miss);
pstreamReadRecord(streamID0, single1, &zlevels_in_miss);
}
}
......@@ -146,16 +146,16 @@ void *Intlevel3d(void *argument)
zlevels_out = (double*) Malloc(gridsize*nlevel*sizeof(double));
nlevo = nlevel; /* number of output levels for later use */
gridsizeo = gridsize;/* horizontal gridsize of output z coordinate */
nrecs = streamInqTimestep(streamID2, 0);
nrecs = pstreamInqTimestep(streamID2, 0);
if (cdoVerbose) cdoPrint("%d records target 3d vertical height and gridsize %d",nrecs,gridsize);
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID2, &varID, &levelID);
pstreamInqRecord(streamID2, &varID, &levelID);
gridsize = gridInqSize(vlistInqVarGrid(vlistID2, varID));
offset = gridsize*levelID;
single1 = zlevels_out + offset;
streamReadRecord(streamID2, single1, &zlevels_out_miss);
pstreamReadRecord(streamID2, single1, &zlevels_out_miss);
}
}
......@@ -311,7 +311,7 @@ void *Intlevel3d(void *argument)
if ( str[0] ) vlistDefVarUnits(vlistID3,oz3dvarID, str);
}
streamDefVlist(streamID3, vlistID3);
pstreamDefVlist(streamID3, vlistID3);
maxlev = nlevi > nlevo ? nlevi : nlevo;
nvars = vlistNvars(vlistID1);
......@@ -384,25 +384,25 @@ void *Intlevel3d(void *argument)
if ( varID == nvars ) cdoAbort("No processable variable found!");
tsID = 0;
while ( (nrecs = streamInqTimestep(streamID1, tsID)) )
while ( (nrecs = pstreamInqTimestep(streamID1, tsID)) )
{
for ( varID = 0; varID < nvars; ++varID ) vars[varID] = false;
taxisCopyTimestep(taxisID3, taxisID1);
streamDefTimestep(streamID3, tsID);
pstreamDefTimestep(streamID3, tsID);
/*
* Read the whole 3d data field
*/
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID1, &varID, &levelID);
pstreamInqRecord(streamID1, &varID, &levelID);
vlistInqVarName(vlistID1, varID, varname);
gridsize = gridInqSize(vlistInqVarGrid(vlistID1, varID));
nlevel = zaxisInqSize(vlistInqVarZaxis(vlistID1, varID));
offset = gridsize*levelID;
single1 = vardata1[varID] + offset;
streamReadRecord(streamID1, single1, &varnmiss[varID][levelID]);
pstreamReadRecord(streamID1, single1, &varnmiss[varID][levelID]);
vars[varID] = true;
}
......@@ -448,9 +448,8 @@ void *Intlevel3d(void *argument)
offset = gridsize*levelID;