From 6c0917a8cc95ff25a54dc9dbad0851ec3048146c Mon Sep 17 00:00:00 2001
From: lkugler <lukas.kugler@gmail.com>
Date: Fri, 9 Jun 2023 11:40:06 +0200
Subject: [PATCH] quickfix error_assimilate
---
dartwrf/assim_synth_obs.py | 17 +++--------------
1 file changed, 3 insertions(+), 14 deletions(-)
diff --git a/dartwrf/assim_synth_obs.py b/dartwrf/assim_synth_obs.py
index b756608..e2132f1 100755
--- a/dartwrf/assim_synth_obs.py
+++ b/dartwrf/assim_synth_obs.py
@@ -282,21 +282,10 @@ def set_obserr_assimilate_in_obsseqout(oso, outfile="./obs_seq.out"):
# modify observation error of each kind sequentially
where_oso_iskind = oso.df.kind == kind
- if obscfg["error_assimilate"] == False:
- # osf_prior = obsseq.ObsSeq(cluster.dartrundir + "/obs_seq.final") # this file will be generated by `evaluate()`
-
- # where_osf_iskind = osf_prior.df.kind == kind
-
- # assim_err = get_parametrized_error(obscfg, osf_prior.df[where_osf_iskind])
- # oso.df.loc[where_oso_iskind, 'variance'] = assim_err**2
- #assert np.allclose(assim_err, oso.df['variance']**2) # check
+ try:
+ oso.df.loc[where_oso_iskind, 'variance'] = obscfg["error_assimilate"]**2
+ except:
pass
- else:
- # overwrite with user-defined values
- try:
- oso.df.loc[where_oso_iskind, 'variance'] = obscfg["error_assimilate"]**2
- except:
- pass
oso.to_dart(outfile)
--
GitLab