From 99dd6107e85e3c676668b7310d60b58e2508956b Mon Sep 17 00:00:00 2001
From: lkugler <lukas.kugler@gmail.com>
Date: Fri, 21 Apr 2023 08:42:03 +0200
Subject: [PATCH] tutorial
---
analysis_only.py | 4 ++--
config/cfg.py | 30 ++++++++++--------------------
config/srvx1.py | 2 +-
3 files changed, 13 insertions(+), 23 deletions(-)
diff --git a/analysis_only.py b/analysis_only.py
index 3076265..899241d 100755
--- a/analysis_only.py
+++ b/analysis_only.py
@@ -6,9 +6,9 @@ import datetime as dt
from dartwrf.workflows import WorkFlows
-prior_path_exp = '/mnt/jetfs/scratch/lkugler/data/sim_archive/exp_v1.19_P3_wbub7_noDA'
+prior_path_exp = '/users/students/lehre/advDA_s2023/data/sample_ensemble/'
prior_init_time = dt.datetime(2008,7,30,12)
-prior_valid_time = dt.datetime(2008,7,30,12,30)
+prior_valid_time = dt.datetime(2008,7,30,13)
assim_time = prior_valid_time
w = WorkFlows(exp_config='cfg.py', server_config='srvx1.py')
diff --git a/config/cfg.py b/config/cfg.py
index 46ea5e7..8a833a6 100755
--- a/config/cfg.py
+++ b/config/cfg.py
@@ -1,7 +1,7 @@
from dartwrf import utils
exp = utils.Experiment()
-exp.expname = "test_newcode"
+exp.expname = "obs-cold_localization-narrow"
exp.model_dx = 2000
exp.n_ens = 10
@@ -9,30 +9,20 @@ exp.filter_kind = 1
exp.prior_inflation = 0
exp.post_inflation = 4
exp.sec = True
-exp.reject_smallFGD = False
-exp.cov_loc_vert_km_horiz_km = (3, 20)
+exp.cov_loc_vert_km_horiz_km = (4, 40)
exp.superob_km = False # False or int (spatial averaging of observations)
exp.adjust_obs_impact = False
exp.use_existing_obsseq = False # False or pathname (use precomputed obs_seq.out files)
-#exp.use_existing_obsseq = '/jetfs/home/lkugler/data/sim_archive/NoImpactFactors/obs_seq_out/2008-07-30_%H:%M_obs_seq.out'
-#exp.use_existing_obsseq = '/jetfs/home/lkugler/data/sim_archive/exp_v1.21_P3_wbub7_VIS_obs10_loc20/obs_seq_out/2008-07-30_%H:%M_obs_seq.out'
-#exp.use_existing_obsseq = '/gpfs/data/fs71386/lkugler/sim_archive/exp_v1.21_P3_wbub7_REFL2D_obs10_loc20_oe5/obs_seq_out/2008-07-30_%H:%M_obs_seq.out'
-#exp.use_existing_obsseq = '/gpfs/data/fs71386/lkugler/sim_archive/exp_v1.21_P2_rr_VIS_obs20_loc4/obs_seq_out/2008-07-30_%H:%M_obs_seq.out'
+exp.use_existing_obsseq = '/users/students/lehre/advDA_s2023/dartwrf_tutorial/very_cold_observation.out'
-#exp.nature = '/mnt/jetfs/scratch/lkugler/data/sim_archive/exp_v1.19_P3_wbub7_nat/2008-07-30_12:00/1'
-exp.nature = '/mnt/jetfs/scratch/lkugler/data/sim_archive/exp_v1.18_P1_nature/2008-07-30_06:00/1'
+# path to the nature run, where we take observations from
+exp.nature = '/users/students/lehre/advDA_s2023/data/sample_nature/'
exp.input_profile = '/mnt/jetfs/home/lkugler/data/initial_profiles/wrf/ens/2022-03-31/raso.fc.<iens>.wrfprof'
-#exp.input_profile = '/gpfs/data/fs71386/lkugler/initial_profiles/wrf/ens/2022-03-31/raso.nat.<iens>.wrfprof'
-#exp.input_profile = '/gpfs/data/fs71386/lkugler/initial_profiles/wrf/ens/2022-05-18/raso.fc.<iens>.wrfprof'
-
-# localize vertically, if it has a vertical position
-# needs a horizontal scale too, to calculate the vertical normalization
-# since you can not specify different vertical localizations for diff. variables
-
-# n_obs= 22500: 2km, 5776: 4km, 121: 30km, 256:16x16 (20km); 961: 10km resoltn # radar: n_obs for each observation height level
+# n_obs can be 22500: 2km, 5776: 4km, 121: 30km, 256:16x16 (20km); 961: 10km resoltn
+# if radar: then n_obs is for each observation height level
vis = dict(plotname='VIS 0.6µm', plotunits='[1]',
kind='MSG_4_SEVIRI_BDRF', sat_channel=1,
@@ -82,15 +72,15 @@ q = dict(plotname='Specific humidity', plotunits='[kg/kg]',
t2m = dict(plotname='SYNOP Temperature', plotunits='[K]',
kind='SYNOP_TEMPERATURE', n_obs=1,
- error_generate=0.1, error_assimilate=1.,
- cov_loc_radius_km=20)
+ error_generate=0.1, error_assimilate=0.1,
+ cov_loc_radius_km=40)
psfc = dict(plotname='SYNOP Pressure', plotunits='[Pa]',
kind='SYNOP_SURFACE_PRESSURE', n_obs=1,
error_generate=50., error_assimilate=100.,
cov_loc_radius_km=32)
-exp.observations = [t]
+exp.observations = [t2m]
exp.update_vars = ['U', 'V', 'W', 'THM', 'PH', 'MU', 'QVAPOR', 'QCLOUD', 'QICE', 'PSFC']
#exp.update_vars = ['U', 'V', 'W', 'T', 'PH', 'MU', 'QVAPOR', 'PSFC']
diff --git a/config/srvx1.py b/config/srvx1.py
index 6e1cf0a..afbc74b 100755
--- a/config/srvx1.py
+++ b/config/srvx1.py
@@ -63,7 +63,7 @@ cluster.archive_base = utils.userhome+'/data/sim_archive/'
cluster.srcdir = '/users/staff/lkugler/AdvDA23/DART/WRF-4.3/run'
cluster.dart_srcdir = '/users/students/lehre/advDA_s2023/DART/models/wrf/work'
cluster.rttov_srcdir = '/users/students/lehre/advDA_s2023/RTTOV13/rtcoef_rttov13/'
-cluster.scriptsdir = utils.userhome+'/DART-WRF/dartwrf/'
+cluster.scriptsdir = utils.userhome+'/AdvDA23/DART-WRF/dartwrf/'
cluster.geo_em = '/users/students/lehre/advDA_s2023/data/geo_em.d01.nc'
# templates/run scripts
--
GitLab