From 3fbd660214d4bcfd0941854dfa3526941ed606db Mon Sep 17 00:00:00 2001
From: lkugler <lukas.kugler@gmail.com>
Date: Wed, 24 May 2023 15:45:03 +0200
Subject: [PATCH] special nature run
---
free_forecast.py | 68 ++++++++++++++++++++++++++++++++++--------------
1 file changed, 48 insertions(+), 20 deletions(-)
diff --git a/free_forecast.py b/free_forecast.py
index 8cee19b..8106192 100755
--- a/free_forecast.py
+++ b/free_forecast.py
@@ -8,7 +8,7 @@ import pandas as pd
from dartwrf.workflows import WorkFlows
-w = WorkFlows(exp_config='cfg.py', server_config='srvx1.py')
+w = WorkFlows(exp_config='nature.py', server_config='jet.py')
id = None
@@ -26,33 +26,61 @@ if False: # generate_nature
# id = w.create_satimages(begin, depends_on=id)
-if False: # to continue a nature after spinup
-
- start = dt.datetime(2008, 7, 30, 7)
- end = dt.datetime(2008, 7, 30, 18)
+if True: # to continue a nature
+ start = dt.datetime(2008, 7, 30, 12)
id = w.prepare_WRFrundir(start) # create initial conditions
- id = w.run_ideal(depends_on=id)
+ # id = w.run_ideal(depends_on=id)
+
+ prior_path_exp = '/jetfs/scratch/lkugler/data/sim_archive//exp_v1.18_P1_nature+1/' # w.cluster.archivedir
- prior_path_exp = '/gpfs/data/fs71386/lkugler/sim_archive/exp_v1.19_P5_nat2' # w.cluster.archivedir
- prior_init_time = dt.datetime(2008, 7, 30, 7)
- prior_valid_time = dt.datetime(2008, 7, 30, 12)
+ time = start
+ prior_init_time = dt.datetime(2008, 7, 30, 6)
+ end = start + dt.timedelta(minutes=30)
+ restarts = pd.date_range(start=dt.datetime(2008, 7, 30, 12,30),
+ end=dt.datetime(2008, 7, 30, 14),
+ freq=dt.timedelta(minutes=30))
+
+ for i, next_restart in enumerate(restarts):
+
+ prior_valid_time = time
+ id = w.prepare_IC_from_prior(prior_path_exp, prior_init_time, prior_valid_time,
+ depends_on=id)
+
+ # integration time
+ end = next_restart
+
+ id = w.run_ENS(begin=start, end=end,
+ input_is_restart=True,
+ first_minute=True, # to get a +1 minute forecast after each restart
+ output_restart_interval=(end-start).total_seconds()/60,
+ #output_restart_interval=9999,
+ depends_on=id)
+
+ w.create_satimages(start, depends_on=id)
+
+ prior_init_time = time # this iteration's start
+ time = end # this iteration's end = next iteration's start
+
+ # after last restart
+ prior_valid_time = time
id = w.prepare_IC_from_prior(prior_path_exp, prior_init_time, prior_valid_time,
- new_start_time=start, # <---------- to start again after spinup
depends_on=id)
-
- id = w.run_ENS(begin=start, end=end,
- input_is_restart=False,
- output_restart_interval=(end-start).total_seconds()/60,
- #output_restart_interval=9999,
- depends_on=id)
+
+ end = time + dt.timedelta(minutes=5)
+ id = w.run_ENS(begin=time, end=end,
+ input_is_restart=True,
+ first_minute=True, # to get a +1 minute forecast after each restart
+ output_restart_interval=9999,
+ depends_on=id)
- id = w.create_satimages(start, depends_on=id)
+ w.create_satimages(start, depends_on=id)
-if True: # do a free run (all inits)
- begin = dt.datetime(2008, 7, 30, 7)
+
+if False: # do a free run (all inits)
+ begin = dt.datetime(2008, 7, 30, 12)
id = w.prepare_WRFrundir(begin) # create initial conditions
id = w.run_ideal(depends_on=id)
@@ -131,7 +159,7 @@ if False: # to continue a free run after spinup
prior_valid_time = dt.datetime(2008, 7, 30, 13,30)
id = w.w.prepare_IC_from_prior(prior_path_exp, prior_init_time, prior_valid_time,
- # new_start_time=start, # <---------- to overwrite start time
+ # new_start_time=start, # <---------- to overwrite start time / leads to a call to `create_updated_wrfinput_from_wrfout()`
depends_on=id)
#frequency_restart = (end-start).total_seconds()/60
--
GitLab