.

free_forecast.py

@@ -8,7 +8,7 @@ import pandas as pd
 from dartwrf.workflows import WorkFlows
 
 
-w = WorkFlows(exp_config='nature.py', server_config='jet.py')
+w = WorkFlows(exp_config='exp_template.py', server_config='jet.py')
 id = None
 
 
@@ -26,7 +26,7 @@ if False: # generate_nature
     # id = w.create_satimages(begin, depends_on=id)
 
 
-if True:  # to continue a nature
+if False:  # to continue a nature
 
     start = dt.datetime(2008, 7, 30, 12)
     id = w.prepare_WRFrundir(start)  # create initial conditions
@@ -80,61 +80,62 @@ if True:  # to continue a nature
 
 
 
-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)
+if True:   # do a free run (all inits)
+    begin = dt.datetime(2008, 7, 30, 12,30)
+    #id = w.prepare_WRFrundir(begin)  # create initial conditions
+    # id = w.run_ideal(depends_on=id)
 
     #id = wrfinput_insert_wbubble(perturb=True, depends_on=id)
 
-    # restarts = pd.date_range(start=dt.datetime(2008, 7, 30, 10),
-    #                           end=dt.datetime(2008, 7, 30, 12),
-    #                           freq=dt.timedelta(minutes=60))
-    restarts = [dt.datetime(2008, 7, 30, 12, 30)]
+    restarts = pd.date_range(start=dt.datetime(2008, 7, 30, 13),
+                              end=dt.datetime(2008, 7, 30, 14,30),
+                              freq=dt.timedelta(minutes=30))
+    # restarts = [dt.datetime(2008, 7, 30, 12, 30)]
     
-    input_is_restart = False
+    input_is_restart = True
     time = begin
-    last_init = dt.datetime(2008, 7, 30, 9)  # dummy value
+    last_init = dt.datetime(2008, 7, 30, 12) 
+
     for i, next_restart in enumerate(restarts):
         print('run_WRF from', time, 'to', next_restart)
+
+        prior_path_exp = '/jetfs/scratch/lkugler/data/sim_archive/exp_v1.19_P2_noDA+1/' #w.cluster.archivedir
+        prior_init_time = last_init
+        prior_valid_time = time
+        id = w.prepare_IC_from_prior(prior_path_exp, prior_init_time, prior_valid_time, depends_on=id)
+
         id = w.run_ENS(begin=time, end=next_restart, 
                      input_is_restart=input_is_restart,
                      output_restart_interval=(next_restart-time).total_seconds()/60,
+                     first_minute=True,
                      #output_restart_interval=720,
                      depends_on=id)
 
         last_init = time
         time = next_restart
         input_is_restart = True
-        w.create_satimages(last_init, depends_on=id)
-
-        prior_path_exp = w.cluster.archivedir
-        prior_init_time = last_init
-        prior_valid_time = time
-        id = w.prepare_IC_from_prior(prior_path_exp, prior_init_time, prior_valid_time, depends_on=id)
+        id_sat = w.create_satimages(last_init, depends_on=id)
 
     # free run, no restart files anymore
-    end = dt.datetime(2008, 7, 30, 18)
-    print('run WRF from', time, 'until', end)
-    id = w.run_ENS(begin=time, end=end,
-             input_is_restart=input_is_restart,
-             #output_restart_interval=(next_restart-time).total_seconds()/60,
-             output_restart_interval=9999,
-             depends_on=id)
-    
+    # end = dt.datetime(2008, 7, 30, 18)
+    # print('run WRF from', time, 'until', end)
+    # id = w.run_ENS(begin=time, end=end,
+    #          input_is_restart=input_is_restart,
+    #          #output_restart_interval=(next_restart-time).total_seconds()/60,
+    #          output_restart_interval=9999,
+    #          depends_on=id)
     
     w.verify_wrf(depends_on=id)
-    id = w.create_satimages(time, depends_on=id)
-    w.verify_sat(depends_on=id)
+    #id = w.create_satimages(time, depends_on=id)
+    w.verify_sat(depends_on=id_sat)
 
 if False:  # to continue a free run
     start = dt.datetime(2008, 7, 30, 7)
-    end = dt.datetime(2008, 7, 30, 10)
 
     id = w.prepare_WRFrundir(start)  # create initial conditions
 
-    prior_path_exp = '/gpfs/data/fs71386/lkugler/sim_archive/exp_v1.19_P5_noDA' # w.cluster.archivedir
-    prior_init_time = dt.datetime(2008, 7, 30, 11)
+    prior_path_exp = '/jetfs/scratch/lkugler/data/sim_archive/exp_v1.19_P2_noDA' # w.cluster.archivedir
+    prior_init_time = dt.datetime(2008, 7, 30, 7)
     prior_valid_time = start
 
     id = w.prepare_IC_from_prior(prior_path_exp, prior_init_time, prior_valid_time, depends_on=id)