From a9d952f7a4ba07351155bc17b3bdcfa1617e4422 Mon Sep 17 00:00:00 2001
From: lkugler <lukas.kugler@gmail.com>
Date: Tue, 14 Sep 2021 11:34:57 +0200
Subject: [PATCH] bugfix error assimilate
---
scripts/assim_synth_obs.py | 9 +++++++--
1 file changed, 7 insertions(+), 2 deletions(-)
diff --git a/scripts/assim_synth_obs.py b/scripts/assim_synth_obs.py
index df29afc..e28376c 100755
--- a/scripts/assim_synth_obs.py
+++ b/scripts/assim_synth_obs.py
@@ -405,8 +405,13 @@ if __name__ == "__main__":
run_perfect_model_obs() # actually create observations that are used to assimilate
################################################
- print(' 3) assimilate with adapted observation-errors ')
- error_assimilate = np.zeros(n_obs_3d) + obscfg['error_assimilate'] # the obs-error we assume for assimilation
+ print(' 3) assimilate with observation-errors for assimilation')
+
+ error_assimilate = []
+ for i, obscfg in enumerate(exp.observations):
+ err_this_type = np.zeros(n_obs_3d) + obscfg['error_assimilate'] # the obs-error we assume for assimilation
+ error_assimilate.extend(err_this_type)
+
replace_errors_obsseqout(cluster.dartrundir+'/obs_seq.out', error_assimilate)
t = time_module.time()
assimilate()
--
GitLab