From 164914ad2a7921ab23e84d95a68b6731737e962d Mon Sep 17 00:00:00 2001
From: lkugler <lukas.kugler@gmail.com>
Date: Tue, 14 Dec 2021 02:48:30 +0100
Subject: [PATCH] update
---
scheduler.py | 8 ++++----
scripts/create_wbubble_wrfinput.py | 4 ++--
2 files changed, 6 insertions(+), 6 deletions(-)
diff --git a/scheduler.py b/scheduler.py
index 8d91ecc..4381ec7 100755
--- a/scheduler.py
+++ b/scheduler.py
@@ -134,8 +134,8 @@ def run_ENS(begin, end, depends_on=None, first_minute=True,
str(hist_interval),
'--radt='+str(radt),
'--restart='+restart_flag,]
- if restart_path:
- args.append('--rst_inname='+restart_path)
+ # if restart_path:
+ # args.append('--rst_inname='+restart_path)
if output_restart_interval:
args.append('--restart_interval='+str(int(float(output_restart_interval))))
@@ -212,7 +212,7 @@ def gen_obsseq(depends_on=None):
def verify(depends_on=None):
- s = my_Slurm("verify", cfg_update={"time": "240", "mail-type": "FAIL,END",
+ s = my_Slurm("verify-"+exp.expname, cfg_update={"time": "240", "mail-type": "FAIL,END",
"ntasks": "96", "ntasks-per-node": "96", "ntasks-per-core": "2"})
s.run(cluster.python_enstools+' /home/fs71386/lkugler/osse_analysis/analyze_fc.py '+exp.expname+' has_node',
depends_on=[depends_on])
@@ -266,7 +266,7 @@ if __name__ == "__main__":
# How long shall we integrate?
timedelta_integrate = timedelta_btw_assim
if time == dt.datetime(2008, 7, 30, 14): #this_forecast_init.minute in [0,]: # longer forecast every full hour
- timedelta_integrate = dt.timedelta(hours=3)
+ timedelta_integrate = dt.timedelta(hours=1)
# 3) Run WRF ensemble
id = run_ENS(begin=time, # start integration from here
diff --git a/scripts/create_wbubble_wrfinput.py b/scripts/create_wbubble_wrfinput.py
index a0ac07b..f452b34 100644
--- a/scripts/create_wbubble_wrfinput.py
+++ b/scripts/create_wbubble_wrfinput.py
@@ -7,7 +7,7 @@ import netCDF4 as nc
dx_km = 2
cr = 15 # km horizontal relaxation distance
-cz = 3000 # meters vertical relaxation distance
+cz = 2000 # meters vertical relaxation distance
perturbations = False
if len(sys.argv) > 1:
@@ -44,7 +44,7 @@ for iens in range(1, exp.n_ens+1):
xx, yy = np.meshgrid(dx, dy)
dr = np.sqrt(xx**2 + yy**2)[np.newaxis, :, :]
- pert = 3*np.exp(-(dr/cr)**2)*np.exp(-(z/cz)**2)
+ pert = 4*np.exp(-(dr/cr)**2)*np.exp(-(z/cz)**2)
ds.variables['T'][0,...] += pert
ds.variables['THM'][0,...] += pert
--
GitLab