rewired collection of obs into optimizer arrays

da/tools/pipeline.py

@@ -222,7 +222,7 @@ def SampleOneCycle(DaCycle,Samples,StateVector, ObservationOperator,lag):
     # Read forecast model samples that were written to NetCDF files by each member, also add the obs to the statevector
     # Note that obs will only be added to the statevector is either this is the first step (restart=False), or lag==nlag
 
-    if lag == nlag-1 or DaCycle['time.restart']==True :
+    if lag == DaCycle['nlag']-1 or DaCycle['time.restart']==True :
 
         silent = False
         members = StateVector.EnsembleMembers[lag]