From d21ba765fc8c52be7b98df1b40acccb0ef6e5947 Mon Sep 17 00:00:00 2001 From: ivar <amvdw95@gmail.com> Date: Fri, 4 Sep 2015 15:32:54 +0000 Subject: [PATCH] bug fix in serial function. KG values were only set zero for the last lag instead of both lags --- da/stilt/optimizer.py | 7 +++++-- 1 file changed, 5 insertions(+), 2 deletions(-) diff --git a/da/stilt/optimizer.py b/da/stilt/optimizer.py index c8f49ca..0ed58f0 100755 --- a/da/stilt/optimizer.py +++ b/da/stilt/optimizer.py @@ -119,9 +119,12 @@ class CO2Optimizer(Optimizer): if 'surface' in self.sitecode[n]: self.KG[-1]=0. + self.KG[3078]=0 + logging.debug("Surface observations do not update BC") if 'aircraft' in self.sitecode[n]: - self.KG[:-1]=0. - + self.KG[0:3078]=0. + self.KG[3079:-1]=0. + logging.debug("Aircraft observations do not update Fluxes") if self.may_localize[n]: logging.debug('Trying to localize observation %s, %i' % (self.sitecode[n], self.obs_ids[n])) -- GitLab