From 9cbaaa6d701443e9ca31a73e2863a5357e621830 Mon Sep 17 00:00:00 2001
From: lkugler <lukas.kugler@gmail.com>
Date: Mon, 22 May 2023 10:35:18 +0200
Subject: [PATCH] .
---
config/cfg.py | 22 +++++++++++-----------
dartwrf/assim_synth_obs.py | 2 +-
templates/input.nml | 2 +-
3 files changed, 13 insertions(+), 13 deletions(-)
diff --git a/config/cfg.py b/config/cfg.py
index 59b40c4..c4b4a29 100755
--- a/config/cfg.py
+++ b/config/cfg.py
@@ -1,16 +1,16 @@
from dartwrf import utils
exp = utils.Experiment()
-exp.expname = "obs-T_inflation-1.1"
+exp.expname = "advDA_Raso-T_inf2-4"
exp.model_dx = 2000
-exp.n_ens = 20
+exp.n_ens = 40
exp.filter_kind = 1
-exp.prior_inflation = 3
-exp.inf_initial = 1.1
+exp.prior_inflation = 2
+exp.inf_initial = 4
exp.post_inflation = 0
exp.sec = True
-exp.cov_loc_vert_km_horiz_km = (4, 40)
+exp.cov_loc_vert_km_horiz_km = (2, 100)
exp.superob_km = False # False or int (spatial averaging of observations)
exp.adjust_obs_impact = False
@@ -41,7 +41,7 @@ wv62 = dict(plotname='Brightness temperature WV 6.2µm', plotunits='[K]',
wv73 = dict(plotname='Brightness temperature WV 7.3µm', plotunits='[K]',
kind='MSG_4_SEVIRI_TB', sat_channel=6,
n_obs=961, obs_locations='square_array_evenly_on_grid',
- error_generate=1., error_assimilate=3.,
+ error_generate=1., error_assimilate=1.,
cov_loc_radius_km=20)
ir108 = dict(plotname='Brightness temperature IR 10.8µm', plotunits='[K]',
@@ -59,16 +59,16 @@ radar = dict(plotname='Radar reflectivity', plotunits='[dBz]',
t = dict(plotname='Temperature', plotunits='[K]',
kind='RADIOSONDE_TEMPERATURE',
- n_obs=961, obs_locations='square_array_evenly_on_grid',
- #n_obs=1, obs_locations=[(45., 0.)],
+ #n_obs=961, obs_locations='square_array_evenly_on_grid',
+ n_obs=1, obs_locations=[(45., 0.)],
error_generate=0.2, error_assimilate=0.2,
- heights=[1000,], #range(1000, 17001, 2000),
- cov_loc_radius_km=40)
+ heights=range(700, 17701, 500),
+ cov_loc_radius_km=100)
q = dict(plotname='Specific humidity', plotunits='[kg/kg]',
kind='RADIOSONDE_SPECIFIC_HUMIDITY', n_obs=1,
error_generate=0., error_assimilate=5*1e-5,
- heights=[1000], #range(1000, 17001, 2000),
+ heights=[1500], #range(1000, 17001, 2000),
cov_loc_radius_km=0.1)
t2m = dict(plotname='SYNOP Temperature', plotunits='[K]',
diff --git a/dartwrf/assim_synth_obs.py b/dartwrf/assim_synth_obs.py
index d1cef07..3396f82 100755
--- a/dartwrf/assim_synth_obs.py
+++ b/dartwrf/assim_synth_obs.py
@@ -476,7 +476,7 @@ def prepare_inflation_2(time, prior_init_time):
if os.path.isfile(f_prior):
copy(f_prior, f_new)
print(f_prior, 'copied to', f_new)
- else:
+ else: # no prior inflation file at the first assimilation
warnings.warn(f_prior + ' does not exist. Using default file instead.')
copy(f_default, f_new)
diff --git a/templates/input.nml b/templates/input.nml
index cfca979..2c987da 100644
--- a/templates/input.nml
+++ b/templates/input.nml
@@ -63,7 +63,7 @@
inf_initial = <inf_initial>, 0.90,
inf_lower_bound = 1.0, 1.0,
inf_upper_bound = 1000000.0, 1000000.0,
- inf_damping = 0.9, 1.0,
+ inf_damping = 1.0, 1.0,
inf_sd_initial = 0.6, 0.0,
inf_sd_lower_bound = 0.6, 0.0,
inf_sd_max_change = 1.05, 1.05,
--
GitLab