Commit 9366bbf4 authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

Rename streamXXX() to pstreamXXX().

parent b9409888
......@@ -35,18 +35,17 @@ void *Deltat(void *argument)
cdoInitialize(argument);
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);
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);
field_type **vars = field_malloc(vlistID1, FIELD_PTR);
......@@ -55,25 +54,25 @@ void *Deltat(void *argument)
double *array2 = (double*) Malloc(gridsizemax*sizeof(double));
int tsID = 0;
int nrecs = streamInqTimestep(streamID1, tsID);
int nrecs = pstreamInqTimestep(streamID1, tsID);
for ( int recID = 0; recID < nrecs; ++recID )
{
streamInqRecord(streamID1, &varID, &levelID);
streamReadRecord(streamID1, vars[varID][levelID].ptr, &nmiss);
pstreamInqRecord(streamID1, &varID, &levelID);
pstreamReadRecord(streamID1, vars[varID][levelID].ptr, &nmiss);
vars[varID][levelID].nmiss = nmiss;
}
tsID++;
int tsID2 = 0;
while ( (nrecs = streamInqTimestep(streamID1, tsID)) )
while ( (nrecs = pstreamInqTimestep(streamID1, tsID)) )
{
taxisCopyTimestep(taxisID2, taxisID1);
streamDefTimestep(streamID2, tsID2);
pstreamDefTimestep(streamID2, tsID2);
for ( int recID = 0; recID < nrecs; ++recID )
{
streamInqRecord(streamID1, &varID, &levelID);
streamReadRecord(streamID1, array1, &nmiss);
pstreamInqRecord(streamID1, &varID, &levelID);
pstreamReadRecord(streamID1, array1, &nmiss);
double missval = vars[varID][levelID].missval;
double *array0 = vars[varID][levelID].ptr;
......@@ -100,8 +99,8 @@ void *Deltat(void *argument)
for ( int i = 0; i < gridsize; ++i ) array0[i] = array1[i];
streamDefRecord(streamID2, varID, levelID);
streamWriteRecord(streamID2, array2, nmiss);
pstreamDefRecord(streamID2, varID, levelID);
pstreamWriteRecord(streamID2, array2, nmiss);
}
tsID++;
......@@ -112,8 +111,8 @@ void *Deltat(void *argument)
Free(array2);
field_free(vars, vlistID1);
streamClose(streamID2);
streamClose(streamID1);
pstreamClose(streamID2);
pstreamClose(streamID1);
cdoFinish();
......
......@@ -78,9 +78,9 @@ void *Deltime(void *argument)
if ( cdoVerbose ) cdoPrint("delete day %d%s", dday, cmons[dmon]);
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);
......@@ -88,9 +88,8 @@ void *Deltime(void *argument)
taxisDefCalendar(taxisID2, CALENDAR_365DAYS);
vlistDefTaxis(vlistID2, taxisID2);
int streamID2 = streamOpenWrite(cdoStreamName(1), cdoFiletype());
streamDefVlist(streamID2, vlistID2);
int streamID2 = pstreamOpenWrite(cdoStreamName(1), cdoFiletype());
pstreamDefVlist(streamID2, vlistID2);
if ( ! lcopy )
{
......@@ -101,7 +100,7 @@ void *Deltime(void *argument)
int nfound = 0;
int tsID = 0;
int tsID2 = 0;
while ( (nrecs = streamInqTimestep(streamID1, tsID)) )
while ( (nrecs = pstreamInqTimestep(streamID1, tsID)) )
{
vdate = taxisInqVdate(taxisID1);
// vtime = taxisInqVtime(taxisID1);
......@@ -122,20 +121,20 @@ void *Deltime(void *argument)
{
taxisCopyTimestep(taxisID2, taxisID1);
streamDefTimestep(streamID2, tsID2++);
pstreamDefTimestep(streamID2, tsID2++);
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID1, &varID, &levelID);
streamDefRecord(streamID2, varID, levelID);
pstreamInqRecord(streamID1, &varID, &levelID);
pstreamDefRecord(streamID2, varID, levelID);
if ( lcopy )
{
streamCopyRecord(streamID2, streamID1);
pstreamCopyRecord(streamID2, streamID1);
}
else
{
streamReadRecord(streamID1, array, &nmiss);
streamWriteRecord(streamID2, array, nmiss);
pstreamReadRecord(streamID1, array, &nmiss);
pstreamWriteRecord(streamID2, array, nmiss);
}
}
}
......@@ -143,8 +142,8 @@ void *Deltime(void *argument)
tsID++;
}
streamClose(streamID2);
streamClose(streamID1);
pstreamClose(streamID2);
pstreamClose(streamID1);
if ( nfound == 0 )
cdoWarning("Day %d%s not found!", dday, cmons[dmon]);
......
......@@ -63,9 +63,9 @@ void *Derivepar(void *argument)
int operatorID = cdoOperatorID();
int streamID1 = streamOpenRead(cdoStreamName(0));
int streamID1 = pstreamOpenRead(cdoStreamName(0));
int vlistID1 = streamInqVlist(streamID1);
int vlistID1 = pstreamInqVlist(streamID1);
int gridID = vlistGrid(vlistID1, 0);
if ( gridInqType(gridID) == GRID_SPECTRAL )
......@@ -288,24 +288,23 @@ void *Derivepar(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;
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);
pstreamInqRecord(streamID1, &varID, &levelID);
zaxisID = vlistInqVarZaxis(vlistID1, varID);
nlevel = zaxisInqSize(zaxisID);
offset = gridsize*levelID;
streamReadRecord(streamID1, array, &nmiss);
pstreamReadRecord(streamID1, array, &nmiss);
if ( zaxisIDh != -1 )
{
......@@ -393,8 +392,8 @@ void *Derivepar(void *argument)
nlevel = nhlevf;
for ( levelID = 0; levelID < nlevel; levelID++ )
{
streamDefRecord(streamID2, varID, levelID);
streamWriteRecord(streamID2, gheight+levelID*gridsize, nmissout);
pstreamDefRecord(streamID2, varID, levelID);
pstreamWriteRecord(streamID2, gheight+levelID*gridsize, nmissout);
}
}
else if ( operatorID == SEALEVELPRESSURE )
......@@ -403,8 +402,8 @@ void *Derivepar(void *argument)
extra_P(sealevelpressure, half_press+gridsize*(nhlevf), full_press+gridsize*(nhlevf-1), sgeopot, temp+gridsize*(nhlevf-1), gridsize);
streamDefRecord(streamID2, 0, 0);
streamWriteRecord(streamID2, sealevelpressure, 0);
pstreamDefRecord(streamID2, 0, 0);
pstreamWriteRecord(streamID2, sealevelpressure, 0);
}
else
cdoAbort("Internal error");
......@@ -412,8 +411,8 @@ void *Derivepar(void *argument)
tsID++;
}
streamClose(streamID2);
streamClose(streamID1);
pstreamClose(streamID2);
pstreamClose(streamID1);
vlistDestroy(vlistID2);
......
......@@ -82,23 +82,22 @@ void *Detrend(void *argument)
cdoInitialize(argument);
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);
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 nvars = vlistNvars(vlistID1);
int tsID = 0;
while ( (nrecs = streamInqTimestep(streamID1, tsID)) )
while ( (nrecs = pstreamInqTimestep(streamID1, tsID)) )
{
if ( tsID >= nalloc )
{
......@@ -112,11 +111,11 @@ void *Detrend(void *argument)
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID1, &varID, &levelID);
pstreamInqRecord(streamID1, &varID, &levelID);
gridID = vlistInqVarGrid(vlistID1, varID);
gridsize = gridInqSize(gridID);
vars[tsID][varID][levelID].ptr = (double*) Malloc(gridsize*sizeof(double));
streamReadRecord(streamID1, vars[tsID][varID][levelID].ptr, &nmiss);
pstreamReadRecord(streamID1, vars[tsID][varID][levelID].ptr, &nmiss);
vars[tsID][varID][levelID].nmiss = nmiss;
}
......@@ -168,7 +167,7 @@ void *Detrend(void *argument)
for ( int tsID = 0; tsID < nts; tsID++ )
{
dtlist_taxisDefTimestep(dtlist, taxisID2, tsID);
streamDefTimestep(streamID2, tsID);
pstreamDefTimestep(streamID2, tsID);
for ( varID = 0; varID < nvars; varID++ )
{
......@@ -178,8 +177,8 @@ void *Detrend(void *argument)
if ( vars[tsID][varID][levelID].ptr )
{
nmiss = vars[tsID][varID][levelID].nmiss;
streamDefRecord(streamID2, varID, levelID);
streamWriteRecord(streamID2, vars[tsID][varID][levelID].ptr, nmiss);
pstreamDefRecord(streamID2, varID, levelID);
pstreamWriteRecord(streamID2, vars[tsID][varID][levelID].ptr, nmiss);
Free(vars[tsID][varID][levelID].ptr);
vars[tsID][varID][levelID].ptr = NULL;
}
......@@ -193,8 +192,8 @@ void *Detrend(void *argument)
dtlist_delete(dtlist);
streamClose(streamID2);
streamClose(streamID1);
pstreamClose(streamID2);
pstreamClose(streamID1);
cdoFinish();
......
......@@ -55,11 +55,11 @@ void *Diff(void *argument)
if ( operatorArgc() == 2 ) rellim = parameter2double(operatorArgv()[1]);
if ( rellim < -1.e33 || rellim > 1.e+33 ) cdoAbort("Rel. limit out of range!");
int streamID1 = streamOpenRead(cdoStreamName(0));
int streamID2 = streamOpenRead(cdoStreamName(1));
int streamID1 = pstreamOpenRead(cdoStreamName(0));
int streamID2 = pstreamOpenRead(cdoStreamName(1));
int vlistID1 = streamInqVlist(streamID1);
int vlistID2 = streamInqVlist(streamID2);
int vlistID1 = pstreamInqVlist(streamID1);
int vlistID2 = pstreamInqVlist(streamID2);
vlistCompare(vlistID1, vlistID2, CMP_ALL);
......@@ -73,21 +73,21 @@ void *Diff(void *argument)
int taxisID = vlistInqTaxis(vlistID1);
while ( TRUE )
{
nrecs = streamInqTimestep(streamID1, tsID);
nrecs = pstreamInqTimestep(streamID1, tsID);
if ( nrecs > 0 )
{
date2str(taxisInqVdate(taxisID), vdatestr, sizeof(vdatestr));
time2str(taxisInqVtime(taxisID), vtimestr, sizeof(vtimestr));
}
nrecs2 = streamInqTimestep(streamID2, tsID);
nrecs2 = pstreamInqTimestep(streamID2, tsID);
if ( nrecs == 0 || nrecs2 == 0 ) break;
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID1, &varID1, &levelID);
streamInqRecord(streamID2, &varID2, &levelID);
pstreamInqRecord(streamID1, &varID1, &levelID);
pstreamInqRecord(streamID2, &varID2, &levelID);
indg += 1;
......@@ -104,8 +104,8 @@ void *Diff(void *argument)
cdiParamToString(param, paramstr, sizeof(paramstr));
streamReadRecord(streamID1, array1, &nmiss1);
streamReadRecord(streamID2, array2, &nmiss2);
pstreamReadRecord(streamID1, array1, &nmiss1);
pstreamReadRecord(streamID2, array2, &nmiss2);
int ndiff = 0;
bool dsgn = false;
......@@ -234,8 +234,8 @@ void *Diff(void *argument)
if ( nrecs > 0 && nrecs2 == 0 )
cdoWarning("stream1 has more time steps than stream2!");
streamClose(streamID1);
streamClose(streamID2);
pstreamClose(streamID1);
pstreamClose(streamID2);
if ( array1 ) Free(array1);
if ( array2 ) Free(array2);
......
......@@ -225,9 +225,9 @@ void *Distgrid(void *argument)
if ( nxblocks <= 0 ) cdoAbort("nxblocks has to be greater than 0!");
if ( nyblocks <= 0 ) cdoAbort("nyblocks has to be greater than 0!");
int streamID1 = streamOpenRead(cdoStreamName(0));
int streamID1 = pstreamOpenRead(cdoStreamName(0));
int vlistID1 = streamInqVlist(streamID1);
int vlistID1 = pstreamInqVlist(streamID1);
int ngrids = vlistNgrids(vlistID1);
......@@ -319,7 +319,7 @@ void *Distgrid(void *argument)
const char *refname = cdoStreamName(0)->argv[cdoStreamName(0)->argc-1];
filesuffix[0] = 0;
cdoGenFileSuffix(filesuffix, sizeof(filesuffix), streamInqFiletype(streamID1), vlistID1, refname);
cdoGenFileSuffix(filesuffix, sizeof(filesuffix), pstreamInqFiletype(streamID1), vlistID1, refname);
for ( int index = 0; index < nsplit; index++ )
{
......@@ -328,23 +328,23 @@ void *Distgrid(void *argument)
sprintf(filename+nchars+5, "%s", filesuffix);
argument_t *fileargument = file_argument_new(filename);
streamIDs[index] = streamOpenWrite(fileargument, cdoFiletype());
streamIDs[index] = pstreamOpenWrite(fileargument, cdoFiletype());
file_argument_free(fileargument);
streamDefVlist(streamIDs[index], vlistIDs[index]);
pstreamDefVlist(streamIDs[index], vlistIDs[index]);
}
if ( ngrids > 1 ) cdoPrint("Baustelle: number of different grids > 1!");
int tsID = 0;
while ( (nrecs = streamInqTimestep(streamID1, tsID)) )
while ( (nrecs = pstreamInqTimestep(streamID1, tsID)) )
{
for ( int index = 0; index < nsplit; index++ )
streamDefTimestep(streamIDs[index], tsID);
pstreamDefTimestep(streamIDs[index], tsID);
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID1, &varID, &levelID);
streamReadRecord(streamID1, array1, &nmiss);
pstreamInqRecord(streamID1, &varID, &levelID);
pstreamReadRecord(streamID1, array1, &nmiss);
double missval = vlistInqVarMissval(vlistID1, varID);
......@@ -352,25 +352,25 @@ void *Distgrid(void *argument)
{
i = 0;
window_cell(array1, array2, grids[i].gridsize[index], grids[i].gridindex[index]);
streamDefRecord(streamIDs[index], varID, levelID);
pstreamDefRecord(streamIDs[index], varID, levelID);
if ( nmiss > 0 )
{
nmiss = 0;
for ( int k = 0; k < grids[i].gridsize[index]; ++k )
if ( DBL_IS_EQUAL(array2[k], missval) ) nmiss++;
}
streamWriteRecord(streamIDs[index], array2, nmiss);
pstreamWriteRecord(streamIDs[index], array2, nmiss);
}
}
tsID++;
}
streamClose(streamID1);
pstreamClose(streamID1);
for ( int index = 0; index < nsplit; index++ )
{
streamClose(streamIDs[index]);
pstreamClose(streamIDs[index]);
vlistDestroy(vlistIDs[index]);
}
......
......@@ -41,9 +41,9 @@ void *Duplicate(void *argument)
if ( cdoVerbose ) cdoPrint("ndup = %d", ndup);
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);
......@@ -67,14 +67,13 @@ void *Duplicate(void *argument)
vlistDefVarTsteptype(vlistID2, varID, TSTEP_INSTANT);
}
int streamID2 = streamOpenWrite(cdoStreamName(1), cdoFiletype());
streamDefVlist(streamID2, vlistID2);
int streamID2 = pstreamOpenWrite(cdoStreamName(1), cdoFiletype());
pstreamDefVlist(streamID2, vlistID2);
nvars = vlistNvars(vlistID1);
int tsID = 0;
while ( (nrecs = streamInqTimestep(streamID1, tsID)) )
while ( (nrecs = pstreamInqTimestep(streamID1, tsID)) )
{
if ( tsID >= nalloc )
{
......@@ -91,11 +90,11 @@ void *Duplicate(void *argument)
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID1, &varID, &levelID);
pstreamInqRecord(streamID1, &varID, &levelID);
int gridID = vlistInqVarGrid(vlistID1, varID);
int gridsize = gridInqSize(gridID);
vars[tsID][varID][levelID].ptr = (double*) Malloc(gridsize*sizeof(double));
streamReadRecord(streamID1, vars[tsID][varID][levelID].ptr, &nmiss);
pstreamReadRecord(streamID1, vars[tsID][varID][levelID].ptr, &nmiss);
vars[tsID][varID][levelID].nmiss = nmiss;
}
......@@ -110,7 +109,7 @@ void *Duplicate(void *argument)
{
taxisDefVdate(taxisID2, vdate[tsID]);
taxisDefVtime(taxisID2, vtime[tsID]);
streamDefTimestep(streamID2, idup*nts+tsID);
pstreamDefTimestep(streamID2, idup*nts+tsID);
for ( varID = 0; varID < nvars; varID++ )
{
......@@ -120,8 +119,8 @@ void *Duplicate(void *argument)
if ( vars[tsID][varID][levelID].ptr )
{
nmiss = vars[tsID][varID][levelID].nmiss;
streamDefRecord(streamID2, varID, levelID);
streamWriteRecord(streamID2, vars[tsID][varID][levelID].ptr, nmiss);
pstreamDefRecord(streamID2, varID, levelID);
pstreamWriteRecord(streamID2, vars[tsID][varID][levelID].ptr, nmiss);
}
}
}
......@@ -134,8 +133,8 @@ void *Duplicate(void *argument)
if ( vdate ) Free(vdate);
if ( vtime ) Free(vtime);
streamClose(streamID2);
streamClose(streamID1);
pstreamClose(streamID2);
pstreamClose(streamID1);
cdoFinish();
......
......@@ -211,8 +211,8 @@ void *EOFs(void * argument)
enum T_EIGEN_MODE eigen_mode = get_eigenmode();
enum T_WEIGHT_MODE weight_mode = get_weightmode();
int streamID1 = streamOpenRead(cdoStreamName(0));
int vlistID1 = streamInqVlist(streamID1);
int streamID1 = pstreamOpenRead(cdoStreamName(0));
int vlistID1 = pstreamInqVlist(streamID1);
int taxisID1 = vlistInqTaxis(vlistID1);
int gridID1 = vlistInqVarGrid(vlistID1, 0);
int gridsize = vlistGridsizeMax(vlistID1);
......@@ -248,17 +248,17 @@ void *EOFs(void * argument)
if ( nts == -1 )
{
nts = 0;
while ( streamInqTimestep(streamID1, nts) ) nts++;
while ( pstreamInqTimestep(streamID1, nts) ) nts++;
if ( cdoVerbose ) cdoPrint("Counted %i timeSteps", nts);
}
else
if ( cdoVerbose ) cdoPrint("Found %i timeSteps", nts);
streamClose(streamID1);
pstreamClose(streamID1);
streamID1 = streamOpenRead(cdoStreamName(0));
vlistID1 = streamInqVlist(streamID1);
streamID1 = pstreamOpenRead(cdoStreamName(0));
vlistID1 = pstreamInqVlist(streamID1);
taxisID1 = vlistInqTaxis(vlistID1);
if ( nts < gridsize || operfunc == EOF_TIME )
......@@ -350,13 +350,13 @@ void *EOFs(void * argument)
/* read the data and create covariance matrices for each var & level */
while ( TRUE )
{
nrecs = streamInqTimestep(streamID1, tsID);
nrecs = pstreamInqTimestep(streamID1, tsID);
if ( nrecs == 0 ) break;
for ( int recID = 0; recID < nrecs; recID++ )
{
streamInqRecord(streamID1, &varID, &levelID);
streamReadRecord(streamID1, in, &nmiss);
pstreamInqRecord(streamID1, &varID, &levelID);
pstreamReadRecord(streamID1, in, &nmiss);
gridsize = gridInqSize(vlistInqVarGrid(vlistID1, varID));
missval = vlistInqVarMissval(vlistID1, varID);
......@@ -440,7 +440,7 @@ void *EOFs(void * argument)
/* write files with eigenvalues (ID3) and eigenvectors (ID2) */
/* eigenvalues */
int streamID2 = streamOpenWrite(cdoStreamName(1), cdoFiletype());
int streamID2 = pstreamOpenWrite(cdoStreamName(1), cdoFiletype());
int vlistID2 = vlistDuplicate(vlistID1);
int taxisID2 = taxisDuplicate(taxisID1);
......@@ -458,7 +458,7 @@ void *EOFs(void * argument)
vlistChangeGridIndex(vlistID2, i, gridID2);
/* eigenvectors */
int streamID3 = streamOpenWrite(cdoStreamName(2), cdoFiletype());
int streamID3 = pstreamOpenWrite(cdoStreamName(2), cdoFiletype());
int vlistID3 = vlistDuplicate(vlistID1);
int taxisID3 = taxisDuplicate(taxisID1);
......@@ -466,8 +466,8 @@ void *EOFs(void * argument)
taxisDefRtime(taxisID3, 0);
vlistDefTaxis(vlistID3, taxisID3);
streamDefVlist(streamID2, vlistID2);
streamDefVlist(streamID3, vlistID3);
pstreamDefVlist(streamID2, vlistID2);
pstreamDefVlist(streamID3, vlistID3);
int vdate = 10101;
int vtime = 0;
......@@ -486,13 +486,13 @@ void *EOFs(void * argument)
taxisDefVdate(taxisID2, vdate);
taxisDefVtime(taxisID2, vtime);
streamDefTimestep(streamID2, tsID);
pstreamDefTimestep(streamID2, tsID);
if ( tsID < n_eig )
{
taxisDefVdate(taxisID3, vdate);
taxisDefVtime(taxisID3, vtime);
streamDefTimestep(streamID3, tsID);
pstreamDefTimestep(streamID3, tsID);
}
for ( varID = 0; varID < nvars; varID++ )
......@@ -615,14 +615,14 @@ void *EOFs(void * argument)
nmiss = 0;
for ( int i = 0; i < gridsize; i++ ) if ( DBL_IS_EQUAL(out[i], missval) ) nmiss++;
streamDefRecord(streamID3, varID, levelID);
streamWriteRecord(streamID3, out, nmiss);
pstreamDefRecord(streamID3, varID, levelID);
pstreamWriteRecord(streamID3, out, nmiss);
} // loop n_eig
nmiss = 0;
if ( DBL_IS_EQUAL(eig_val[tsID], missval) ) nmiss = 1;
streamDefRecord(streamID2, varID, levelID);
streamWriteRecord(streamID2, &eig_val[tsID], nmiss);
pstreamDefRecord(streamID2, varID, levelID);
pstreamWriteRecord(streamID2, &eig_val[tsID], nmiss);
} // loop nlevs
} // loop nvars
}
......@@ -652,9 +652,9 @@ void *EOFs(void * argument)
Free(pack);
Free(weight);