setup

analysis_only.py

@@ -18,6 +18,8 @@ prior_valid_time = dt.datetime(2008,7,30,12,30)
 assim_time = prior_valid_time
 
 
+cluster.setup()
+
 os.system(
     cluster.python+' '+cluster.scripts_rundir+'/assim_synth_obs.py '
                 +assim_time.strftime('%Y-%m-%d_%H:%M ')
@@ -28,4 +30,4 @@ os.system(
 
 
 
-# id_sat = create_satimages(time, depends_on=id)
\ No newline at end of file
+# id_sat = create_satimages(time, depends_on=id)

config/cfg.py

 from dartwrf import utils
 
 exp = utils.ExperimentConfiguration()
-exp.expname = "exp_v1.22_P2_rr_WV62_obs10_loc20_oe2" #"exp_v1.22_P2_rr_VIS+WV73_obs10_loc20_inf5"
+exp.expname = "exp_v1.22_P2_rr_VIS_obs10_loc20_inf5"
 exp.model_dx = 2000
 exp.n_ens = 40
 
 exp.filter_kind = 1
-exp.prior_inflation = 0
-exp.post_inflation = 4
+exp.prior_inflation = 5
+exp.post_inflation = 0
 exp.sec = True
 exp.reject_smallFGD = False
 exp.cov_loc_vert_km_horiz_km = (3, 20)
@@ -93,7 +93,7 @@ psfc = dict(plotname='SYNOP Pressure', plotunits='[Pa]',
             error_generate=50., error_assimilate=100.,
             cov_loc_radius_km=32)
 
-exp.observations = [wv62]
+exp.observations = [vis]
 exp.update_vars = ['U', 'V', 'W', 'THM', 'PH', 'MU', 'QVAPOR', 'QCLOUD', 'QICE', 'PSFC']
 #exp.update_vars = ['U', 'V', 'W', 'T', 'PH', 'MU', 'QVAPOR', 'PSFC']
 

cycled_exp.py

@@ -200,6 +200,7 @@ def verify_fast(depends_on=None):
 ################################
 if __name__ == "__main__":
     print('starting osse')
+    cluster.setup()
 
     timedelta_integrate = dt.timedelta(minutes=15)
     timedelta_btw_assim = dt.timedelta(minutes=15)

generate_free.py

@@ -16,6 +16,7 @@ from cycled_exp import *
 ################################
 print('starting osse')
 
+cluster.setup()
 id = None