open stream before calling cdoDefineGrid().

src/Enlarge.cc

@@ -36,6 +36,7 @@ void *Enlarge(void *argument)
 
   operatorCheckArgc(1);
 
+  // open stream before calling cdoDefineGrid!!!
   int streamID1 = pstreamOpenRead(cdoStreamName(0));
 
   int gridID2 = cdoDefineGrid(operatorArgv()[0]);

src/Enlargegrid.cc

@@ -166,6 +166,7 @@ void *Enlargegrid(void *argument)
   if ( operatorArgc() < 1 ) cdoAbort("Too few arguments!");
   if ( operatorArgc() > 2 ) cdoAbort("Too many arguments!");
 
+  // open stream before calling cdoDefineGrid!!!
   int streamID1 = pstreamOpenRead(cdoStreamName(0));
 
   int gridID2 = cdoDefineGrid(operatorArgv()[0]);

src/Intgrid.cc

@@ -304,6 +304,7 @@ void *Intgrid(void *argument)
 
   int operatorID = cdoOperatorID();
 
+  // open stream before calling cdoDefineGrid!!!
   int streamID1 = pstreamOpenRead(cdoStreamName(0));
 
   if ( operatorID == INTGRIDBIL || operatorID == INTERPOLATE )

src/MapReduce.cc

@@ -79,6 +79,7 @@ void *MapReduce(void *argument)
 
   cdoInitialize(argument);
 
+  // open stream before calling cdoDefineGrid!!!
   int streamID1 = pstreamOpenRead(cdoStreamName(0));
 
   /* check input grid type and size - this will be used for selecting relevant

src/Remap.cc

@@ -840,6 +840,7 @@ void *Remap(void *argument)
 	cdoPrint("Extrapolation disabled!");
     }
 
+  // open stream before calling cdoDefineGrid!!!
   int streamID1 = pstreamOpenRead(cdoStreamName(0));
 
   if ( lremapxxx )

src/Setgrid.cc

@@ -67,6 +67,9 @@ void *Setgrid(void *argument)
 
   int operatorID = cdoOperatorID();
 
+  // open stream before calling cdoDefineGrid!!!
+  int streamID1 = pstreamOpenRead(cdoStreamName(0));
+
   if ( operatorID != UNSETGRIDMASK )
     operatorInputArg(cdoOperatorEnter(operatorID));  
 
@@ -177,8 +180,6 @@ void *Setgrid(void *argument)
       griduri = operatorArgv()[0];
     }
 
-  int streamID1 = pstreamOpenRead(cdoStreamName(0));
-
   int vlistID1 = pstreamInqVlist(streamID1);
   int vlistID2 = vlistDuplicate(vlistID1);
 

src/YAR.cc

@@ -714,6 +714,7 @@ void *YAR(void *argument)
   field_init(&field1);
   field_init(&field2);
 
+  // open stream before calling cdoDefineGrid!!!
   int streamID1 = pstreamOpenRead(cdoStreamName(0));
 
   int gridID2 = cdoDefineGrid(operatorArgv()[0]);