.

dartwrf/assim_synth_obs.py

@@ -5,10 +5,10 @@ import numpy as np
 
 from config.cfg import exp, cluster
 from dartwrf.utils import symlink, copy, sed_inplace, append_file, mkdir, try_remove, print, shell
+from dartwrf.obs import error_models as err
 import dartwrf.create_obsseq as osq
 from dartwrf import wrfout_add_geo
 from dartwrf import obsseq
-from dartwrf.obs import error_models as err
 
 earth_radius_km = 6370
 

dartwrf/create_obsseq.py

@@ -237,8 +237,8 @@ def create_obs_seq_in(time_dt, list_obscfg,
 
         assert len(coords) == n_obs, (len(coords), n_obs)  # check if functions did what they supposed to
         for lat, lon in coords:
-            assert lat < 90 & lat > -90
-            assert lon < 180 & lon > -180
+            assert (lat < 90) & (lat > -90)
+            assert (lon < 180) & (lon > -180)
 
         kind = obscfg['kind']
         print('obstype', kind)

dartwrf/evaluate_posterior.py

@@ -7,7 +7,6 @@ from config.cfg import exp, cluster
 from dartwrf import assim_synth_obs as aso
 
 
-
 if __name__ == "__main__":
 
      assim_time = dt.datetime.strptime(sys.argv[1], "%Y-%m-%d_%H:%M")

dartwrf/evaluate_prior.py

@@ -6,6 +6,7 @@ import numpy as np
 from config.cfg import exp, cluster
 from dartwrf.utils import symlink, copy, sed_inplace, append_file, mkdir, try_remove, print, shell
 from dartwrf import assim_synth_obs as aso
+
 from osselyze.utils import get_prior_config
 
 if __name__ == "__main__":

dartwrf/obsseq_to_netcdf.py

 import os, sys, glob, warnings
+
 from config.cfg import exp, cluster
 import run_obs_diag as rod
 

dartwrf/prepare_namelist.py

@@ -11,6 +11,7 @@ Options:
 import os, sys, shutil, warnings
 import datetime as dt
 from docopt import docopt
+
 from config.cfg import exp, cluster
 from utils import sed_inplace, copy, symlink, mkdir
 

dartwrf/run_obs_diag.py

 import os, sys, shutil, glob
+
 from config.cfg import exp, cluster
 from utils import symlink, copy, sed_inplace, append_file, shell
 
 rundir_program = '/home/fs71386/lkugler/data/run_DART/'
 
-
 def prepare(obserr_iszero='.true.'):
     copy(cluster.scriptsdir+'/../templates/input.eval.nml',
             rundir_program+'/input.nml')

dartwrf/wrfout_add_geo.py

@@ -9,6 +9,7 @@ Example call:
 """
 import os, sys
 import netCDF4 as nc
+
 from config.cfg import exp, cluster
 
 fields_old = ["XLAT_M",   "XLONG_M",      "CLAT",